автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Исследование эффективности применения международной модели ионосферы IRI-2001 для прогнозирования характеристик ВЧ радиосвязи

ШЛЮПКИН Александр Сергеевич

ИССЛЕДОВАНИЕ ЭФФЕКТИВНОСТИ ПРИМЕНЕНИЯ МЕЖДУНАРОДНОЙ МОДЕЛИ ИОНОСФЕРЫ 1М-2001 ДЛЯ ПРОГНОЗИРОВАНИЯ ХАРАКТЕРИСТИК ВЧ РАДИОСВЯЗИ

05.13.18 - Математическое моделирование, численные методы и комплексы программ

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

г. Ростов-на-Дону 2007 г.

003053257

Работа выполнена на кафедре радиофизики физического факультета Ростовского государственного университета

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

с. н. с. МАЛЬЦЕВА Ольга Алексеевна

Официальные оппоненты: доктор технических наук,

профессор ЖАК Сергей Вениаминович кандидат физико-математических наук, с. н. с. БРОНИН Андрей Геннадьевич

Ведущая организация: Институт прикладной геофизики имени

академика Федорова Е.К. (г. Москва)

Защита состоится « 1 » ^^Ч7^ 2007 г. в м часов на заседании диссертационного совета К.212.208.04 по физико-математическим и техническим наукам в Южном федеральном университете по адресу: 344090, г. Ростов-на-Дону, пр. Стачки 200/1, корпус 2, ЮГИНФО, к. 206.

С диссертацией можно ознакомиться в научной библиотеке ЮФУ по адресу: г. Ростов-на-Дону, ул. Пушкинская, 148.

Автореферат разослан « ^9» января 2007 г.

Ученый секретарь диссертационного Совета, кандидат физико-математических наук, доцент

Муратова Г.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы.

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

Задача прогнозирования решается путем использования математических моделей среды распространения ВЧ радиоволн - ионосферы. На основании таких моделей рассчитываются частотные и амплитудные характеристики радиотрасс. Достигнутые в последние годы успехи в ионосферных исследованиях и, как следствие, создание новых корректных моделей ионосферы, а также доступность применения ЭВМ со значительно возросшими вычислительными возможностями позволяют сегодня на качественно новой основе решать задачу прогнозирования характеристик ДКМ каналов.

На сегодняшний день в мире используется большое количество моделей ионосферы. Особое место среди них занимает модель IRI (International Reference of Ionosphere), которая, по мнению группы авторитетных экспертов, является наиболее разработанной и динамично развивающейся в своем классе. Последняя версия - IRI-2001 является глобальной медианной моделью ионосферы (т.е. позволяет строить долгосрочные прогнозы в любой точке земного шара). В ней предусмотрена возможность учета возмущенного состояния ионосферы, а также коррекции по данным текущей диагностики.

В то же время, созданная геофизиками, модель IRI не прошла еще достаточной проверки на трассах ВЧ связи. Так, в научной литературе мало данных о погрешностях прогнозирования иа базе IRI основной характеристики радиотрасс - максимальной применимой частоты (МПЧ). Причем это можно сказать не только о возмущенной, но и о спокойной ионосфере. Практически не исследованы адаптационные возможности модели, и, как следствие, мало информации о возможности применения модели для оперативного прогноза. Представляет также интерес оценка возможности использования для коррекции модели различных средств диагностики.

Из сказанного следует, что задача исследования точностных характеристик модели IRI, долгосрочного и оперативного прогноза условий распространения декаметровых волн ВЧ трасс на ее основе является актуальной, имеющей важное научно-прикладное значение.

Цели и задачи диссертаиии.

1. Исследовать эффективность применения математической модели "Международная справочная модель ионосферы IRI-2001" для долгосроч-

ного и оперативного прогнозирования условий распространения ДКМВ.

2. Осуществить выбор методов коррекции модели по данным текущей диагностики, позволяющих сократить погрешность прогнозирования характеристик ВЧ радиосвязи.

3. Разработать программную реализацию предложенных методов и оценить их эффективность.

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

- Провести сопоставительный анализ математических моделей ионосферы последних лет и определить место модели ERI-2001 среди них.

- Разработать программную реализацию модели, позволяющую производить прогнозирование параметров ионосферы при организации ВЧ радиосвязи с использованием программного модуля IRI-2001, доступного через сеть Internet.

- Осуществить поиск баз экспериментальных данных на радиотрассах для последующей статистически достоверной проверки прогнозов на основе IRI-2001.

- Разработать программную реализацию методов коррекции модели IRI-2001 по данным вертикального зондирования и зондирования со спутников.

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

- Оценить эффективность применения модели IRI-2001 для оперативного прогнозирования.

- Представить обоснование возможности использования спутниковой радионавигационной системы GPS для коррекции модели IRI-2001.

Научная новизна результатов исследования.

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

1. Количественные статистически обоснованные оценки погрешностей для односкачковых радиотрасс европейского региона для долгосрочного и оперативного прогнозов максимальной применимой частоты, выполненные на основе математической модели IRI

2. Количественные оценки эффективности коррекции модели по данным ВЗ (foF2) и измерениям ТЕС применительно к задачам оперативного прогнозирования МПЧ.

3. Метод корректировки модели IRI по данным ТЕС, учитывающий форму верхней части профиля ионосферы.

4. Алгоритмы и программные реализации корректировки модели IRI по данным ВЗ и ТЕС.

5. Алгоритм и программа получения ТЕС по результатам приема сигналов GPS.

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

Теоретические и практические результаты диссертационной работы использованы в НИР: «Стандерс-РГУ», «Крюшон-РГУ», «МАП-Б», выполненных по гособоронзаказам.

Полученные в ходе выполнения работы результаты, алгоритмы и программные продукты могут быть использованы при:

прогнозировании максимальных применимых и оптимальных рабочих частот (МПЧ и ОРЧ) на линиях ВЧ связи;

коррекции характеристик распространения ДКМВ по данным текущей диагностики;

создании оперативных систем диагностики ионосферы с использованием GPS;

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

расчетах ионосферных погрешностей абсолютных и относительных навигационных определений высотных (ионосферных) абонентов СРНС; оценке точностных характеристик систем ВЧ пеленгаторов.

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

Результаты работы были представлены и обсуждены на следующих конференциях: X и XI региональных конференциях по распространению радиоволн (С - Петербург, 2004, 2005), X всероссийской научной конференции студентов физиков (Москва, 2004), международной конференции по проблемам геокосмоса (С - Петербург, 2004), международном симпозиуме "Solar Extreme Events of 2003"( Москва, 2004), EGU-2005 (Вена, 2005), международной научной конференции "Излучение и рассеяние электромагнитных волн" (Таганрог, 2005), XI международном симпозиуме по ионосферным явлениям IES-2005 (Александрия, США, 2005), X международной конференции Ionospheric Radio Systems and Techniques (IRST) (Лондон, 2006), I международной конференции European Conference on Antennas & Propagation (Eu-CAP2006) (Ницца, 2006).

Публикации.

По теме диссертации опубликовано 20 печатных работ, из них 3 - в российских журналах, рекомендованных ВАК.

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

Диссертация состоит из введения, пяти глав материала и заключения с общим объемом в 145 листов, включая список цитируемой литературы и приложения.

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

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

Первая глава носит обзорный характер. Здесь анализируются существующие модели прогнозирования параметров ионосферы. Глава состоит из шести разделов. В первом разделе приводятся общие требования, предъявляемые к моделям ионосферы с позиций ВЧ радиосвязи. Они следующие:

1. Глобальность - позволяет прогнозировать параметры ионосферы в любой точке земного шара - отсутствуют широтно-долготные ограничения.

2. Адаптивность - возможность коррекции модели по результатам диагностики одного (или нескольких) из параметров ионосферы.

3. Возможность учета возмущенного состояния ионосферы.

4. Простота и высокая скорость расчетов.

5. Доступность - наличие модели в свободном доступе, справочных материалов и технической под держки.

Три следующих раздела посвящены краткому описанию теоретических, гибридных и эмпирических моделей соответственно. Показано, что применимость теоретических моделей для оперативного прогнозирования ВЧ каналов в условиях ограниченности ресурсов времени весьма проблематична. Примером такой модели может служить глобальная теоретическая модель ионосферы (Global Theoretical Ionospheric Model - GTIM). В качестве примера гибридных (полуэмпирических) математических моделей представлена отечественная полуэмпирическая модель ионосферы (ПЭМИ), признанная одной из лучших в мире. Однако упрощения, введенные по сравнению с теоретическими моделями, не позволяют намного снизить требования, предъявляемые к вычислительным ресурсам. Использование эмпирических математических моделей позволяет значительно упростить вычисления и тем самым повысить скорость работы прикладных программ расчета ВЧ каналов. Требования к вычислительным ресурсам при этом сводятся к минимуму. Однако эти модели носят медианный характер, т.е. их возможности ограничены долгосрочным прогнозированием. Для решения текущих задач оперативного прогноза требуется коррекция модели. Именно этот факт приводит к необходимости исследования возможности адаптации модели по результатам текущей диагностики состояния ионосферы.

Широко известны следующие эмпирические модели: IRI (International Reference of Ionosphere - международная справочная модель ионосферы); отечественная модель СМИ (Справочная Модель Ионосферы, НПО "Тайфун" г. Обнинск); европейская модель NeQuick.

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

В пятом разделе содержится краткая характеристика последней версии модели IRI (IRI-2001). Отличительные особенности модели, определив-

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

В последнем (шестом) разделе главы содержатся выводы.

Вторая глава посвящена детальному описанию модели IRI-2001. Показаны широкие возможности модели, позволяющие разработчикам использовать ее для решения различных задач.

В первом разделе приводятся общие положения модели, а также описываются входные и выходные параметры. Пакет программ распространяется NSSDC (National Space Science Data Centre) в виде FORTRAN подпрограмм и файлов коэффициентов CCIR(MKKP) и URSI, файлов документации. Доступ к программам организован через FIP-сервер.

Модель IRI является многопараметрической, включающей определение нескольких классов параметров, таких как: концентрация электронов и ионов и их температура; концентрация нейтральных компонент и их температура; характеристики магнитного поля и многие другие дополнительные величины.

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

Выходными параметрами являются следующие высотные профили: электронной плотности N(h) и Ne/NmF2, электронной температуры Те, ионной температуры Ti, нейтральной температуры Тп и процентного содержание ионов в диапазоне от 60 до 2000 км, а также ТЕС.

Второй раздел содержит основные положения алгоритма работы IRI. Распределение концентраций с высотой в IRI состоит из 2-х частей: модели нижней ионосферы для высот от начала ионосферы ho до высоты максимума hmF2 слоя F2; модели верхней части - для области выше hmF2.

Расчет исходных основных ионосферных параметров foF2 и M3000F2 в модели IRI может быть реализован двумя вариантами, выбор которых задан соответствующей входной опцией (1-URSI, 0-CCIR).

CCIR - это вычислительная программа прогноза ITU-R, построенная с использованием метода Лежандра-Фурье, в котором среднемесячные foF2 определяются коэффициентами разложения для любых значений широты, долготы и времени: n(<p,Ä,t) = a0(<p,Ä,t)+'£jaj((p,Z)-cos jt+bj{g>,X)-sin jf], где: (1)

jt=0

aj и bj - коэффициенты Фурье; Gk(cp, X) - сферические функции Лежандра; ф, X - географические координаты; п=76 для foF2; п=49 для M(3000)F2; m=6 для foF2; ш=4 для M(3000)F2.

Глобальные функции Gk((p, X) соответствуют каждому коэффициенту Фурье. Численные карты параметров foF2 и M(3000)F2 получаются с помощью 2-х наборов коэффициентов (URSI/CCIR) для низкой и высокой солнечной активности, для каждого месяца. Для вычисления foF2 используется индекс 1(312 солнечной активности.

Для наземных трасс важна ионосфера только ниже максимума области F. Для этих высот в модели IRI определяются следующие параметры вертикального распределения ионизации:

- максимальные (критические) частоты слоев: foF2, foFl, foE, foD;

- высоты максимумов слоев: hmF2, hmFl, hmE, hmD;

- полутолщины: Во (аналог утл для слоя F2), Cl (ym для слоя F1);

- параметры долины: высота hvt вышележащего слоя, на которой электронная концентрация равна NmE, высота Ь„ь и электронная концентрация Nvb в точке минимума долины.

Поскольку наибольшие ошибки моделей и прогнозов ионосферных параметров связаны с возмущенными условиями, в модели IRI введен учет возмущенного состояния ионосферы (ST-модель), механизм которого описан в третьем разделе. Учет возмущенных условий проводится как аддитивная поправка к спокойным значениям foF2. Модель дает корректирующий множитель CF - шторм-фактор, позволяющий перейти от foF2Mea. - критической частоты слоя F2 эмпирической медианной модели — к критической частоте на период возмущения: foF2Kopp. = CF-foF2MM.

В четвертом разделе рассматривается возможность вариации параметров IRI, позволяющая использовать модель в режиме адаптации по данным диагностики в реальном масштабе времени. Адаптивными параметрами могут быть индексы солнечной активности, критические частоты и высоты слоев. А наиболее перспективным средством текущей диагностики - использование ТЕС, восстановленного по результатам мониторинга ионосферы спутниками GPS.

В пятом разделе содержатся выводы второй главы. Несмотря на медианный характер, ограничивающий применение IRI для нужд радиосвязи, отмечаются широкие адаптационные возможности IR1 по данным текущей диагностики одного (или нескольких) из параметров, позволяющие компенсировать указанный недостаток.

Третья глава посвящена исследованию точности применения IRI-2001 для долгосрочного прогнозирования трасс ВЧ радиосвязи. Производится расчет количественных значений погрешностей прогнозирования для различных сезонов и уровней солнечной и геомагнитной активности.

В первом разделе главы описана методика проведения исследования и дается краткое представление о разработанном программном обеспечении, позволяющем производить необходимые для анализа вычисления. Основной

прием исследования - сопоставление суточных вариаций экспериментальных и модельных значений параметров foF2 и МПЧ. Экспериментальные foF2 были получены на 24 станциях ВЗ, большая часть из которых представлена на рисунке 2. Источниками получения этих данных были Internet-сайт SPIDR (Space Physics Interactive Data Resource), а также публикации в специальной литературе. Общий временной интервал подвергшихся анализу данных составил 7 лет. Точность прогнозирования МПЧ оценивалась по данным 6 европейских трасс НЗ, список которых и периоды измерений приведены в таблице 1.

О-станции НЗ; А.- станции ВЗ; станции ТЕС. Рис.2. Карта размещения использованных в работе станций.

Таблица 1. Список трасс наклонного зондирования.

№ Tpacca H3 Протяженность трассы, км Время наблюдения Индекс RZu

1 Chelveston - Linkoping 1216 1-31 марта 1993 67

2 Rome - Linkoping 1850 1-31 марта 1993 67

1-28 июня 1993 57

3 Inskip - Moscow 2551 1-30 апреля 2002 111

4 Cyprus - Moscow 2293

5 Inskip - Rostov 3042 11-15 апреля 2004 35

19-20 августа 2003 61

6 Inskip - Rome 1735 2-31 декабря 2003 43

Экспериментальные данные по трассам НЗ также были взяты из литературных источников, либо предоставлены лабораториями, проводившими

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

Во втором разделе главы приводятся результаты оценки точности долгосрочного прогнозирования с использованием различных моделей. В таблице 2 дается оценка точности прогноза ЛэР2 для различных условий, полученная ранее отечественными методами. Одним из критериев этой оценки является среднеквадратичное отклонение (СКО) о модельных значений от экспериментальных. В диссертации показано, что долгосрочный прогноз по Ш1 имеет более высокую точность, чем прогнозы с использованием других известных моделей.

Таблица 2. Оценка точности прогноза foF2.

широты уровень солн.акт. день ночь восход-заход

а(%) <т, МГц а, МГц а(%) <7, МГц

высокие макс(1979) 13-14 1,1-1,2 28,2-30,7 1,7-2,0 - -

мин(1987) 14,2-16,8 0,9-1,1 29,4-33,2 1,2-1,4 - -

средние мин-макс 10,9-12,4 1,2 17,4-20,0 0,9-1,0 15-18 1,3-1,6

низкие макс 27,9-39,9 4,6-6,4 60,0-76,9 4,7-5,2 28,7-9,5 5,5-5,6

Третий раздел содержит оценку долгосрочного прогноза йзР2 по модели ГО.1-2001, полученную в результате сравнения рассчитанных значений с экспериментальными данными для упомянутых выше станций ВЗ. Основной задачей была оценка влияния на прогноз шторм-фактора (БТ-модель). Из представленных в работе периодов наиболее типичные результаты получены в марте 1993 г., когда наряду со спокойными условиями наблюдались дни положительных и отрицательных возмущений. Геомагнитная обстановка в виде Бв^индекса для периода с 6 по 16 марта приведена на рисунке 3.

:: \f v^

.140 ---и--------

DST

Рис.3. Геомагнитная обстановка с 6 по 16 марта 1993 года.

Оценка точности IRI и эффективности шторм-фактора осуществлялась по экспериментальным значениям foE2, полученным на 15-ти средне-широтных станциях вертикального зондирования. Использование данных значительного числа станций ВЗ позволило сократить влияние на оценку эффективности шторм-фактора случайных погрешностей инструментального определения foF2. Эффективность С определялась по формуле:

С(%) = 100х G°s ~(2)

°off

Здесь стоЯ - среднеквадратичные отклонения (СКО) значений foF2, рассчитанных по IRI с опцией "off' (т.е. без шторм-факгора), от экспериментально наблюдаемых; аоп - СКО для foF2, найденные с опцией "on" (т.е. с учетом

6Маг 7Маг 8Маг ЗМаг 10 Маг 11 Маг 12Маг13Ма г 14 Маг 15Маг16Маг

\rW

шторм - фактора). В таблице 3 содержатся итоговые результаты исследования по экспериментальным данным за 1993 г, показывающие, что для дневных часов средние значения С *■ 39%, а для ночного времени С =23%. Общая статистика исследования показала, что при отрицательном возмущении применение шторм-фактора позволило сократить среднеквадратичную погрешность прогнозирования й)Р2 в 1,5+2,5 раза. В ночные часы эффективность шторм-фактора несколько ниже.

Таблица 3. Итоговые результаты исследования за 1993 год.

Условия День Ночь

Оогь МГц от,МГц С, % о0(Т ,МГц Соп ,МГц С, %

Спокойные 1,11 1,11 0 0,94 0,94 0

Отрицательная фаза 1,91 1,17 38,7 1,13 0,87 23,0

Среднее значения по всем возмущениям 1,26 1,07 15,1 1,09 1,01 7,3

Отметим, что применение шторм-фактора на положительной стадии возмущения постоянно сопровождается ухудшением соответствия между прогнозируемыми и наблюдаемыми значениями Л)Р2. Этот факт приводит к необходимости контроля БзЬиндекса при прогнозировании ионосферных параметров по ПИ-2001.

Четвертый раздел посвящен оценке точности ПУ и эффективности применения шторм-фактора на трассах НЗ (таблица 1) при прогнозировании МПЧ. В марте 1993 г. исследование проводилось на трассе №1. Наглядные результаты представлены на рисунке 4 для спокойной ионосферы (6.03.1993), отрицательного (9.03.1993) и положительного (16.03.1993) возмущений.

16.03.93

Сплошной линией показаны экспериментальные значения МНЧ, пунктирной - значения МПЧ, рассчитанные по модели Ш1-2001 без учета шторм

фактора, штрихпунктирной - МПЧ, найденные с учетом шторм-фактора (кривая, маркированная кружками - У8, описывается в четвертой главе). Использовано среднее европейское время (СЕТ).

Количественные оценки степени совпадения суточных зависимостей МПЧ и МНЧ, усредненные за март, приведены в таблице 4. Здесь о - среднеквадратичное отклонение прогнозируемых и экспериментальных значений МНЧ и МПЧ в МГц, о(%) - относительное среднеквадратичное отклонение. Результаты представлены отдельно для дня и ночи. Как и ранее, получено ухудшение соответствия при учете шторм-фактора для случая положительного возмущения и улучшение в случае отрицательного возмущения.

Таблица 4. Количественные оценки степени совпадения суточных зависимостей МПЧ и МНЧ в марте 1993 г. _____

O.MHz Ooff <3st Ovs doff, % Ost, % Ovs Cst, % Civs, %

День Спокойные условия 2,54 2,46 1,1 14,5 14,0 6,3 3,1 56,7

Отрицательные возмущения 5,43 2,96 2,47 57,2 31,2 26,0 45,5 54,5

Положительные возмущения 4,18 5,47 0,65 22,8 29,8 3,5 -30,9 84,4

Ночь Спокойные условия 0,92 0,56 0,34 15,2 9,3 5,6 39,1 63,0

Отрицательные возмущения 2,14 2,31 0,57 34,0 36,7 9,0 -7,9 73,4

Положительные возмущения 0,9 0,71 0,62 15,2 12,0 10,5 21,1 31,1

Для трассы №6 проведен анализ результатов обработки месячных данных НЗ - суточных ходов МНЧ обыкновенной компоненты, полученных в 15-минутных сеансах, с помощью chip sounder при односкачковом распространении от F2 слоя за 30 дней (с 2 по 31 декабря 2003 года). Экспериментальные данные получены QinetiQ (UK) и the National Institute of Geophysics and Vulcanology (INGV) в Риме и любезно предоставлены нам через Internet. Общий массив данных представлен на рисунке 5 (ломаные серые кривые), где также приведен суточный ход МНЧ, усредненный за 30 дней измерений в каждом из 96 временных разрезов, через 15 минут (толстая черная кривая). Здесь же представлен суточный ход МПЧ, рассчитанный по модели IRI-2001 (тонкая черная кривая), среднеквадратичная погрешность оказалась равной 9.7% относительно усреднения. Следует отметить, что декабрь 2003г. хотя и был относительно спокойным месяцем в геомагнитном отношении, однако, наблюдались слабые, в основном положительные, возмущения, поэтому шторм-фактор не только не дал существенного улучшения при описании экспериментальных МНЧ, но, наоборот, ухудшил ситуацию.

Декабрь 2003,

Рис.5. МПЧ и МНЧ для трассы №6 в декабре 2003 г.

Разброс экспериментальных данных, как видно из рисунка 5, достаточно велик. Для количественной оценки этого разброса были проведены расчеты как абсолютных — в МГц, так и относительных — в процентах, величин среднеквадратичного отклонения. Если рассматривать среднесуточную величину СКО, то она составляет 12.8%, наиболее изменчивым в процентном отношении является восходно-заходный период суток, когда разброс экспериментальных данных относительно среднего значения за месяц может превышать 20%, тогда как в полуденные часы он составляет всего 8-10%.

В пятом разделе содержатся выводы главы. Основные результаты сводятся к следующему. Результаты исследования прогноза 1ЪЕ2:

1. В спокойных условиях среднее значение а ~ 1 МГц.

2. Во время отрицательных возмущений получено о ~ 2 МГц. Процент улучшения при использовании шторм-фактора составил - 30%.

Результаты исследования прогноза МПЧ:

1. Значение относительной ошибки в спокойных условиях менее 15%.

2. Во время отрицательного возмущения относительная погрешность для дня более 50 %, для ночи — более 35%. Для положительной фазы соответственно более 23 и 15%. Использование шторм-фактора позволило снизить погрешность определения МПЧ на ~ 45 % днем, хотя практически не повлияло на погрешность в ночное время. В момент положительных возмущений применение шторм-фактора ухудшает прогноз.

Из этого следует, что для решения задач практической ВЧ связи, где требуется высокая устойчивость и гибкость, а ресурсы жестко ограничены, использование ГО1 без коррекции неприемлемо. Необходимо искать решение в использовании адаптационных свойств модели для сокращения погрешности.

Четвертая глава представляет исследование эффективности коррекции IRI-2001, т.е. эффективности оперативного прогнозирования ионосферных каналов связи. Корректировка модели при вычислении МПЧ трасс рассматривалась по двум возможным направлениям: 1) по экспериментальным данным foF2 станций ВЗ и 2) по экспериментальным значениям ТЕС станций TEC-GPS. При этом подразумевается, что станции должны быть расположены вблизи траектории распространения сигнала.

В первом разделе рассматривается метод корректировки модели по данным ВЗ. Его суть заключается в использовании данных foF2 для ввода в модель IRI-2001 в качестве входного параметра. Это позволяет изменить модельный профиль ионизации и остальные параметры в соответствии с экспериментальным значением foF2. На следующем этапе происходит расчет тра-екгорных характеристик и МПЧ трассы по скорректированной IRI-2001. Процедура определения МПЧ трассы требует знания пространственного распределения плотности ионизации вдоль траектории распространения сигнала. В общем случае, при корректировке модельного распределения ионизации вдоль радиотрассы опорные станции ВЗ могут находиться как внутри, так и вне зоны распространения волны. Тогда для получения параметров модели вдоль трассы возникает необходимость в построении пространственных карт этих параметров на момент связи по нескольким произвольно расположенным точкам с известными экспериментальными данными. Многолетние исследования привели к появлению множества достаточно точных для практического применения методов решения этой задачи, основанных на пространственной корреляции параметров ионосферы. В данной работе для восстановления пространственного распределения применялся более современный метод кригинга. Описание алгоритма, применительно к ионосферным задачам, представлено в диссертации.

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

Можно выделить следующие методы корректировки модели IRI: 1) модификация коэффициентов модели в соответствии с глобальными ионосферными картами GIM (GPS ТЕС); 2) варьирование эффективного значения параметра, характеризующего солнечную активность; 3) использование эффективной ширины ионосферы; 4) использование коэффициентов, корректирующих верхнюю часть ЩЬ)-профиля.

В методе 1 корректировке подвергаются коэффициенты разложения параметров йзГО и М(3000)Е2 в ряды Фурье в глобальном масштабе. Для корректировки модели вводится дополнительный множитель К( (р, А, г):

где символы low и high соответствуют низкой и высокой солнечной активности.

Скорректированные коэффициенты использовались для расчета ТЕС и обеспечили увеличение точности расчета ТЕС по модели IRI на 32.5%. Эта цифра соответствует спокойным условиям, для возмущенных условий соответствие намного хуже. Эффективность использования этого метода для корректировки foF2 неизвестна.

Метод 2 практически идентичен первому. Фактически, он является его прототипом. Он может использоваться в локальном и региональном вариантах и заключается в том, что подбирается такое эффективное значение ионосферного индекса RZ12 или IGi2 или само значение foF2, при котором рассчитанное значение ТЕС близко к экспериментальному. Это значение и является скорректированным.

Метод 3 связан с эффективной шириной ионосферы. Величину ТЕС можно представить в виде произведения Nm на эквивалентную ширину слоя

т:

Или, эквивалентно, - через ГоР2: ТЕС = 1.24-10"(^2)2 г, где (оР2 - в МГц, т -в км. Для определения т нужно пользоваться специальной климатологической моделью. Поскольку климатологическая модель х не обладает такой полнотой, как модели йэЮ и других параметров, пришлось использовать величины т, рассчитываемые из модели Ш1. Тогда процедуру корректировки можно свести к следующему. Для заданного момента времени по модели П11 определяются значения ТЕС, №пР2(Л>Р2) и т по формуле (4). Используя это значение т, по экспериментальному значению ТЕС определяется новое (скорректированное) значение №пР2, соответственно, ¡"оБ2.

В методе 4 значение т, определяемое в рамках модели ГО1 и являющееся медианным значением, корректируется по текущему значению ТЕС. В данной работе для такой коррекции предлагается использовать коэффициент К=т(э)/т(Ш.1), связывающий наблюдаемые значения х(э) с т(Ш1), который может быть представлен в виде множителей №п(ПИ)/Ыт(э) и ТЕС(э)./ТЕС(Ш1). Предполагая, что основной причиной отклонения реального Т от модельного является форма верхнего слоя ионосферы (Ь>ЬшР2), получаем выражение для коэффициента К=1-а(с-1)/с, где параметр а - вклад ТЕС верхнего М(Ъ)-профиля в полное значение ТЕС. Параметр с - коэффициент перехода от одной формы профиля к другой. Предлагаемый коэффициент обеспечивает корректировку по текущему значению ТЕС по формуле

-j -к{(рЛл)

(3)

T=TEC/NraF2

(4)

fc>F2(cK0p) = foF2(IKI)x —x ТЕС(э) (5)

^ 1]К TEC{IRI) K '

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

Числовое значение, характеризующее эффективность каждого из методов корректировки по отношению к долгосрочному прогнозу, оценивалось коэффициентом улучшения CI:

o(off)-o(cKoP)xm% 6)

<T(off)

Эффективность корректировки ПИ по данным ВЗ исследована в третьем разделе. Приведенная на рисунке 4 кривая, маркированная кружками, наглядно показывает улучшение прогноза МГГЧ как в спокойный период, так и в периоды возмущений. Таблица 4 содержит числовые значения точности коррекции по ВЗ oVs (как по абсолютной величине, так и в процентах), а также коэффициент улучшения CIVs- Результаты исследования по остальным трассам подтверждают факт улучшения прогноза МПЧ вне зависимости от сезона и характера возмущений.

Поскольку корректировка IRI по ТЕС-данным на первом этапе подразумевает восстановление foF2, в четвертом разделе рассматривается эффективность коррекции IRI предложенными методами. С помощью специальной программы были реализованы методы 2,3,4 (метод 1 не удалось реализовать из-за ограниченности ресурсов). Для тестирования методов корректировки данные по foF2 были взяты из SPIDR, а по ТЕС - из Internet-базы RAL (Rutherford Appleton Laboratory) по 25 станциям TEC-GPS.

Для иллюстрации результатов корректировки модели IRI по данным ТЕС описанными методами приводится пример расчетов по станции Rome для сентября 2003 г. Этот период характеризуются сильным возмущением межпланетного магнитного поля, повлекшим значительное изменение (уменьшение) NmF2 без изменения геомагнитных индексов. На рисунке 6 непрерывной толстой линией показана кривая экспериментальных значений, пунктир - расчет по исходной (без корректировки) модели IRI. Кривая, маркированная квадратами, показывает результаты корректировки методом №2, а кружками - методом №4.

Общая статистика показывает, что если отклонение модельных значений foF2 от реальных не превышает точности измерений (0.2-0.5 Мгц), что часто соответствует спокойным условиям, то методы №2 и №3 могут приводить к ухудшению соответствия, а метод №4 может оставить отклонение на исходном уровне или незначительно его улучшить. Во время возмущенных условий все методы улучшают соответствие (методы №2 и №3 - в 1.5-2 раза, метод №4 - в 2-3 раза). Таким образом, описанные методы позволяют использовать ТЕС-данные для дальнейшего прогноза МПЧ.

17 Sep 18 Sep 21 Sep 22 Sep

Рис.6. Результаты прогнозирования в возмущенных условиях 2003 г.

В пятом разделе исследуется эффективность коррекции IRI при прогнозировании МПЧ по результатам ТЕС-мониторинга. В исследовании использован метод №4 ТЕС-коррсктировки. В качестве примера, демонстрирующего эффективность корректировки, приводятся таблицы для описанных выше трасс. В таблице 5 содержится сравнение точностных характеристик долгосрочного прогноза (IRIoff), коррекции по ВЗ (IRIvs), коррекции по ТЕС (IRItec). в также коррекции модели СМИ (результаты взяты из литературных источников) по ВЗ (CMHys) на трассе №3 12.04.2002 г. Из таблицы видно, что, хотя корректировка по ВЗ дает наилучшие показатели, ТЕС-корректировка также позволяет получать улучшение по сравнению с долгосрочным прогнозом по IRI и оперативным прогнозом по СМИ.

Таблица 5. Таблица результатов корректировки 12.04.2002 г.

°мпч (МГц) °мпч (за сутки) °мпч (за день) стмпч (за ночь)

1 скачок 2 скачка 1 скачок 2 скачка 1 скачок 2 скачка

IRIofr 2,46 1,45 2,33 1,33 2,58 1,56

IRIvs 1,34 0,77 1,01 0,55 1,61 0,94

ШТЕС 1,46 1,08 1,12 0,76 1,74 1,33

CMHvs 2,04 1,24 2,34 1,52 1,69 0,87

8мпч (%) 8мпч (за сутки) 5мпч (за день) бщгч (за ночь)

1 скачок 2 скачка 1 скачок 2 скачка 1 скачок 2 скачка

IRIoff 11,75 10,81 10,27 9,34 13,43 12,40

IRIvs 6,40 5,74 4,45 3,86 8,38 7,47

ШтЕС 6,98 8,05 4,94 5,34 9,06 10,57

CMHvs 9,75 9,25 10,31 10,67 8,80 6,91

В таблице 6 приведены результаты ТЕС-коррекции для трассы КаЗ в спокойные и возмущенные периоды апреля 2002 г., а в таблице 7 - для трассы №5 в относительно спокойных условиях апреля 2004 г. Результаты по остальным трассам аналогичны.

Таблица 6. Результаты корректировки IRI для трассы №3 в апреле 2002 г.

Спокойные условия Возмущенные условия

Щ>(Т IRIvs ШтЕС IRIoff IRIsi IRIvs IRItec

0, MHz 2,55 1,47 2,12 6,03 3,71 3,58 3,24

а,% 14,5 8,4 12,1 34,3 21,1 20,4 18,4

С1,% 42,3 16,9 38,5 40,6 46,3

Таблица 7. Оценка эффективности

Шестой раздел содержит выводы главы. Результаты сравнения МПЧ, полученных с использованием 1Ш-2001 для европейских трасс, с экспериментальными МНЧ, показывают:

1) данные ВЗ обеспечивают стабильное сокращение ошибки прогноза на 25-40%, свидетельствуя о том, что текущее значение Л)Р2 является наиболее эффективным для контроля состояния ВЧ каналов;

2) ТЕС-корректировка модели Ш1 приводит к снижению погрешности прогноза значений Л>Р2 и МПЧ в среднем на 30%.

коррекции IRI-2001 для трассы №5.

Дата IRI IRIvs Штес

о, MHz 1,49 1,57 1,17

11.04.04 о, % 8,38 8,83 6,61

С1,% -5,4 21,5

с, MHz 2,18 0,98 1,31

12.04.04 а, % 13,65 6,14 8,2

С1,% 55,0 39,9

a, MHz 1,68 1,6 1,63

13.04.04 а,% 9,4 8,95 9,11

CI, % 4,8 3,0

а, MHz 2,09 1,18 1,58

14.04.04 а, % 11,26 6,36 8,5

С1,% 43,5 24,4

о, MHz 2,13 1,13 1,35

15.04.04 о, % 12,17 6,46 7,71

ci,% 46,9 36,6

О, MHz 1,91 1,29 1,41

среднее С, % 10,93 7,38 8,07

С1,% 32,5 26,2

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

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

Каждый спутник передает GPS-приемнику навигационное сообщение, содержащее параметры, необходимые для расчетов. Однако в исходном виде навигационное сообщение сложно использовать, поэтому на практике его преобразуют к строго регламентированному текстовому формату и сохраняют в виде нескольких специализированных файлов. Во втором разделе описывается наиболее распространенный из таких форматов - RJNEX (Receiver Independent Exchange Format — Формат передачи данных, не зависящий от приемника). В рамках диссертации была разработана компьютерная программа, извлекающая из RINEX-файлов необходимую для расчетов информацию.

Третий раздел посвящен алгоритму расчета ТЕС по данным, извлекаемым из RINEX-файлов, и дальнейшему использованию ТЕС для корректировки IRI. Общая схема алгоритма представлена на рисунке 7.

Рис.7. Схема использования GPS для оперативного прогноза.

Геометрия прохождения сигналов GPS изображена на рисунке 8. Ионосфера рассматривается в виде тонкой узкой оболочки без горизонтальных (широтно-долготных) вариаций на высоте hm = 350 км - метод тонкого слоя. Траектории строятся в приближении геометрической оптики. Измеряемыми величинами являются расстояния (дальности) или задержки при распространении сигнала между приемником и спутником. Расчетные уравнения сводятся к:

r.-W- * 1 ■ 1

/,2 + Л3

уь{А)+Ъ{р) , где (7)

Pi, Рг - псевдодальности на частотах fi и /2; ь{а) и ¿(d) - неизвестные межчастотные запаздывания приемника и спутника.

Рис.8. Геометрия прохождения сигналов GPS.

Величину b(A)+b{D) можно считать постоянной в рамках определенного интервала времени и оценивать раз в сутки при помощи Калибровки, требующей наличия RINEX-файлов за несколько часов. Информация об ионосфере заключена в измеряемой величине I полного электронного содержания ионосферы, которая является интегральной электронной концентрацией вдоль траектории луча между спутником и приемником - наклонный ТЕС, который связан с вертикальным ТЕС в подионосферных точках (В и С) функцией наклона.

В четвертом разделе производится проверка пригодности ТЕС, полученного по результатам GPS мониторинга, для корректировки модели IRI при проведении оперативного прогнозирования ионосферных трасс.

Пятый раздел - сводка результатов пятой главы. Сравнение рассчитанного ТЕС с доступными данными, использованными в качестве контрольного варианта, позволяет сделать вывод о возможности использования алгоритма для корректировки IRI-2001 при оперативном прогнозировании ионосферных трасс.

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

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

1. Оценки эффективности применения математической модели "Международная справочная модель ионосферы IRI-2001" для долгосрочного и оперативного прогнозирования условий распространения декаметровых волн для спокойного и возмущенного состоянии ионосферы.

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

3. Методы коррекции модели IRI по данным ВЗ и значениям полного электронного содержания (ТЕС) для оперативного прогнозирования состояния ВЧ каналов связи и количественная оценка их эффективности. Показано, что коррекция приводит к существенному сокращению погрешностей, прежде всего, для возмущенных периодов.

4. Способ получения ТЕС по результатам GPS мониторинга.

Личный вклад соискателя. Все основные результаты работы получены при активном участии автора. Автор принимал непосредственное участие в разработке 4-го метода ТЕС-корректировки. Программные реализации рассмотренных методов и алгоритмов, необходимые расчеты и их анализ выполнены самостоятельно.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ В СЛЕДУЮЩИХ РАБОТАХ

1. Барабашов Б.Г., Крашенинников И.В., Егоров И.Б., Мальцева О.А., Родионова В.Т., Черкашин Ю.Н., Шлюпкин А.С. Новые результаты наклонного зондирования и их использования для КВ радиосвязи //Труды международной научной конференции "Излучение и рассеяние электромагнитных волн" (ИРЭМВ-2005), Таганрог, 20-25 июня 2005. - С.383-385.

2. Барабашов Б.Г., Крашенинников И.В., Егоров И.Б., Мальцева О.А., Родионова В.Т., Черкашин Ю.Н., Шлюпкин А.С. Оценка эффективности использования модели IRI для планирования ВЧ радиосвязи //Сборник докладов XXI всероссийской научной конференции "Распространение радиоволн", Йошкар-Ола, 25-27 мая 2005. - Т. 2. - С.19-23.

3. Мальцева О.А., Полтавский О.С., Шлюпкин А.С. Использование GPS-данных в каналах КВ-связи //Электросвязь. - 2006. - №11.- С.47-49.

4. Мальцева О.А., Полтавский О.С., Шлюпкин А.С. Эффективность модели ионосферы IRI при определении условий распространения в ВЧ канале //Электромагнитные волны и электронные системы. - 2006. - №9. - С.21-24.

5. Мальцева О.А., Родионова В.Т., Шлюпкин А.С. Использование полного электронного содержания для текущей диагностики состояния ионосферного канала //Геомагнетизм и аэрономия. - 2005. - №4. - С.480-486.

6. Мальцева О.А., Родионова В.Т., Шлюпкин А.С. Об использовании параметра ПЭС для текущей диагностики состояния ионосферного канала //Тезисы X региональной конференции по распространению радиоволн, С-Петербург, 2004. - С.54-57.

7. Мальцева О.А., Шлюпкин А.С. Остаточная погрешность модели IRI и новый метод определения И(Ь)-профилей ионосферы //Труды XI региональной конференции по распространению радиоволн, С.- Петербург, 2005. -С.46-48.

8. Шлюпкин А.С. Использование навигационных спутников для радиосвязи //Тезисы 4-й международной научно-практической конференции "Телекоммуникационные технологии на транспорте России", ТелеКомТранс 2006, Сочи, 19-21 апреля 2006.- С.36-37.

9. Шлюпкин А.С., Мальцева О.А., Полтавский О.С. Использование информации навигационных спутников в каналах КВ радиосвязи //Сборник тезисов ВНКСФ-10, Москва, 1-7 апреля 2004. - С.989-991.

10. Barabashov В.G., Maltseva О.A., Rodionova V.T., Shlupkin A.S. Evaluation of HF Channel Characteristics on The Basis of The IRI Model //Proceedings of the llth International Ionospheric Effects Symposium (IES2005), Alexandria, USA, May 3-5,2005. - A127. - P.l-4.

11. Barabashov B.G., Maltseva O.A., Rodionova V.T., Shlupkin A.S. Evaluation of the IRI Model Efficiency for Operational Forecast of HF Propagation Conditions //The X EET International Conférence on IRST 2006, Jul. 2006, London: Proceedings Sériés IC517. - P.253-257.

12. Maltseva O.A., Barabashov B.G., Rodionova V.T., Shlyupkin A.S. One of possibilities to use TEC-GPS data for HF links //Proceedings of The first European Conference on Antennas and Propagation (EuCAP 2006), Nov. 6-10, 2006, Nice, France.-P02.12.-P.309.

13. Maltseva O.A., Barabashov B.G., Rodionova V.T., Shlyupkin A.S. One of possibilities to use TEC-GPS data for HF links //Special Publication, EuCAP 2006, Nov. 6-10, 2006, Nice, France. - SP-626 on CD. - P.l-5.

14. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Evaluation of real time TEC and STORM-TIME corrections in the IRI model //EGU2005, Apr. 25-29, 2005, Vienna, Austria: Geophysical Research abstracts. - V.7. - EGU05-A-03017. -P.l-3.

15. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Evaluation of real-time TEC and STORM-TIME correction in the IRI //Proceedings of the XXVIII General Assembly URSI, Oct. 23-29, 2005, New Daly, India. - 601b.2(0147) - P.244-249.

16. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Joint using TEC and IRI for studying ionospheric responses to geomagnetic disturbances //Program and Abstract Book International Symposium on Solar Extreme Events of 2003, Jul. 12-14, 2004, MSU, Moscow. - P.51.

17. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Space Weather for HF-communication links // Book of abstracts, 5th International Conference on Problems of Geocosmos, May 24-28,2004, St. Petersburg. - P.256.

18. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Space Weather for HF-communication links //Proceedings of the 5th International Conference "Problems of Geocosmos", 2004, Saint-Petersburg. - P.355-358.

19. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Using Total Electron Content of the Ionosphere for Forecasting Its Critical Frequency //International Conference on Problems of Geocosmos, May 23-27, 2006, St.-Petersburg: Book of abstracts. -P.160-161.

20. Maltseva O.A., Shlupkin A.S. The model IRI residual error and the new method of N(h) profile determination //EGU2006, Apr. 2-6,2006, Vienna, Austria: Geophysical Research abstracts. - V.8. - N 01180. - P.l-3.

Подписано в печать 26.О-/. О* I

Формат 60x84 '/и. Бумага офсетная. Печать офсетная. ;

ОбъемУ,Офиз. печ. л. Тираж/сО экз. Заказ Л«2- 55". ,

НПО РГПУ: !

344082, г. Ростов-на-Дону, ул. Большая Садовая, 33 1

Введение.

Глава 1. Обзор ионосферных моделей

1.1. Требования, предъявляемые к моделям ионосферы

1.2. Теоретические модели

1.3. Гибридные модели

1.4. Эмпирические модели

1.5. Краткая характеристика модели IRI

1.6. Выводы

Глава 2. Международная справочная модель ионосферы IRI

2.1. Общие положения 18 Входные параметры модели 19 Выходные параметры модели

2.2. Алгоритм модели

2.3. Учет возмущений (ST-модель)

2.4. Учет вариаций параметров

2.5. Выводы

Глава 3. Исследование точности применения IRI-2001 для долгосрочного прогнозирования foF2 и МПЧ

3.1. Методика исследования

3.2. Оценка точности прогноза

3.3. Результаты долгосрочного прогнозирования foF

3.4. Результаты долгосрочного прогнозирования МПЧ

3.5. Выводы

Глава 4. Исследование эффективности коррекции IRI

4.1. Метод корректировки по данным ВЗ

Восстановление пространственного распределения foF2 методом ^ кригинга

4.2. Методы корректировки по ТЕС 64 Точность определения ТЕС 65 Метод №1 69 Метод №2 71 Метод №3 72 Метод №

4.3. Оценка эффективности коррекции IRI-2001 по данным ВЗ

4.4. Оценка эффективности оперативного прогноза foF2 по ТЕС

4.5. Оценка эффективности оперативного прогноза МПЧ по ТЕС

4.6. Выводы

Глава 5. Использование GPS для расчета ионосферных трасс

5.1. Общие сведения о GPS

5.2. Описание формата RINEX

5.3. Алгоритм расчета ТЕС по данным GPS 113 Системы координат 114 Вычисление координат по эфемеридной информации спутников 115 Алгоритм обработки измерений для получения ТЕС 119 Калибровка

5.4. Опытная проверка возможности адаптации IRI-2001 по данным . ? . GPS

5.5. Выводы 127 Заключение 128 Список использованной литературы 130 Приложение 1. Формат данных RINEX 138 Приложение 2. Примеры RINEX-файлов

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

Задача прогнозирования решается путем использования математических моделей среды распространения ВЧ радиоволн - ионосферы. На основании таких моделей рассчитываются частотные и амплитудные характеристики радиотрасс. Достигнутые в последние годы успехи в ионосферных исследованиях и, как следствие, создание новых корректных моделей ионосферы, а также доступность применения ЭВМ со значительно возросшими вычислительными возможностями позволяют сегодня на качественно новой основе решать задачу прогнозирования характеристик ДКМ каналов.

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

На сегодняшний день в мире используется большое количество моделей ионосферы. Особое место среди них занимает модель IRI (International Reference of Ionosphere), которая, по мнению группы авторитетных экспертов [11], является наиболее разработанной и динамично развивающейся в своем классе. Последняя версия - IRI-2001 является глобальной медианной моделью ионосферы (т.е. позволяет строить долгосрочные прогнозы в любой точке земного шара). В ней предусмотрена возможность учета возмущенного состояния ионосферы, а также коррекции по данным текущей диагностики.

В то же время, созданная геофизиками, модель IRI не прошла еще достаточной проверки на трассах ВЧ связи. Так, в научной литературе мало данных о погрешностях прогнозирования на базе IRI основной характеристики радиотрасс - максимальной применимой частоты (МПЧ). Причем это можно сказать не только о возмущенной, но и о спокойной ионосфере. Практически не исследованы адаптационные возможности модели, и, как следствие, мало информации о возможности применения модели для оперативного прогноза. Представляет также интерес оценка возможности использования для коррекции модели различных средств диагностики.

Из сказанного следует, что задача исследования точностных характеристик модели IRI, долгосрочного и оперативного прогноза условий распространения декаметровых волн ВЧ трасс на ее основе является актуальной, имеющей важное научно-прикладное значение. Цели и задачи диссертации.

1. Исследовать эффективность применения математической модели "Международная справочная модель ионосферы IRJ-2001" для долгосрочного и оперативного прогнозирования условий распространения ДКМВ.

2. Осуществить выбор методов коррекции модели по данным текущей диагностики, позволяющих сократить погрешность прогнозирования характеристик ВЧ радиосвязи.

3. Разработать программную реализацию предложенных методов и оценить их эффективность.

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

- Провести сопоставительный анализ математических моделей ионосферы последних лет и определить место модели IR1-2001 среди них.

- Разработать программную реализацию модели, позволяющую производить прогнозирование параметров ионосферы при организации ВЧ радиосвязи с использованием программного модуля IRI-2001, доступного через сеть Internet.

- Осуществить поиск баз экспериментальных данных на радиотрассах для последующей статистически достоверной проверки прогнозов на основе IRI-2001.

- Разработать программную реализацию методов коррекции модели IRI-2001 по данным вертикального зондирования и зондирования со спутников, исследовать их эффективность.

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

- Оценить эффективность применения модели IRI-2001 для оперативного прогнозирования.

- Представить обоснование возможности использования спутниковой радионавигационной системы GPS для коррекции модели IRI-2001. Эти задачи решались на основе моделирования суточных зависимостей таких характеристик, как критическая частота ионосферы (foF2), полученная по данным вертикального зондирования (ВЗ), полное электронное содержание ионосферы в столбе (ПЭС или более употребительная в литературе аббревиатура ТЕС от английского Total Electron Content), рассчитываемое на основе навигационных сообщений глобальной навигационной системы (GPS), максимальная применимая частота (МПЧ), измеряемая методом наклонного зондирования (НЗ).

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

1. Количественные статистически обоснованные оценки погрешностей для од-носкачковых радиотрасс европейского региона для долгосрочного и оперативного прогнозов максимальной применимой частоты, выполненные на основе математической модели IRI.

2. Количественные оценки эффективности коррекции модели по данным ВЗ (foF2) и измерениям ТЕС применительно к задачам оперативного прогнозирования МПЧ.

3. Метод корректировки модели IRI по данным ТЕС, учитывающий форму верхней части профиля ионосферы.

4. Алгоритмы и программные реализации корректировки модели IRI по данным ВЗ и ТЕС.

5. Алгоритм и программа получения ТЕС по результатам приема сигналов GPS.

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

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

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

Результаты работы были представлены и обсуждены на следующих конференциях: X и XI региональных конференциях по распространению радиоволн (С

Петербург, 2004, 2005), => X всероссийской научной конференции студентов физиков (Москва, 2004), международной конференции по проблемам геокосмоса (С - Петербург, 2004), международном симпозиуме "Solar Extreme Events of 2003"( Москва, 2004),

EGU-2005 (Вена, 2005), => международной научной конференции "Излучение и рассеяние электромагнитных волн" (Таганрог, 2005), => XI международном симпозиуме по ионосферным явлениям IES-2005

Александрия, США, 2005), => X международной конференции Ionospheric Radio Systems and Techniques

IRST) (Лондон, 2006), => I международной конференции European Conference on Antennas & Propagation (EuCAP2006) (Ницца, 2006).

Пубчикации

По теме диссертации опубликовано 20 работ, из них 3 в российских журналах, рекомендованных ВАК.

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

1. Оценки эффективности применения математической модели "Международная справочная модель ионосферы IRI-2001" для долгосрочного и операгивного прогнозирования условий распространения декаметровых волн для спокойного и возмущенного состоянии ионосферы.

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

3. Методы коррекции модели IRI по данным ВЗ и значениям полного электронного содержания (ТЕС) для оперативного прогнозирования состояния ВЧ каналов связи и количественная оценка их эффективности. Показано, что коррекция приводит к существенному сокращению погрешностей, прежде всего, для возмущенных периодов.

4. Способ получения ТЕС по результатам GPS мониторинга.

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

Работа состоит из введения, пяти глав, заключения, списка литературы и двух приложений; изложена на 145 листах; содержит 35 рисунков и 39 таблиц.

Нумерация рисунков, таблиц и формул определяется номером главы и раздела.

Например, запись Рис.3.2.1. означает, что это первый рисунок второго раздела третьей главы. Список цитируемой литературы включает 134 наименования.

5.5. Выводы.

Спутники GPS системы, изначально не предусматривающей измерение параметра ТЕС, передают потребителю все необходимые данные для ТЕС-мониторинга в глобальном масштабе. В данной главе предложен успешный алгоритм извлечения информации навигационных спутников и вычисления ТЕС на ее основе, позволяющий выполнять коррекцию модели IRI-2001 при оперативном прогнозировании параметров ВЧ радиотрасс. Для каждого из пунктов реализованного алгоритма выполнена программная реализация.

ЗАКЛЮЧЕНИЕ

В диссертации исследована эффективность применения математической модели "Международная справочная модель ионосферы IRI-2001" для долгосрочного и оперативного прогнозирования условий распространения ДКМВ, предложены алгоритмы коррекции и разработана их программная реализация. Из полученных результатов, имеющих, прежде всего, прикладное значение, наибольшего внимания заслуживают следующие:

1. Выполнено исследование эффективности IRI-2001 для долгосрочного и оперативного прогнозирования условий распространения ДКМВ для спокойного и возмущенного состоянии ионосферы, различных сезонных и суточных вариаций.

2. По данным станций ВЗ и НЗ европейского региона, охватывающим большой временной интервал, получены количественные значения погрешности долгосрочного прогноза. Показано, что наибольшее значение ошибка имеет в период возмущенного состояния ионосферы, не позволяя использовать IRI-2001 в исходном виде для оперативного прогноза.

3. Предложены методы коррекции модели IRI-2001 по данным ВЗ и значениям ТЕС для оперативного прогнозирования состояния ВЧ каналов связи. Количественно оценена их эффективность. Показано, что коррекция приводит к существенному сокращению погрешностей, прежде всего, для возмущенных периодов.

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

5. Для всех описанных алгоритмов разработаны программные реализации, позволяющие использовать разработанную методику в условиях ограниченности временных и вычислительных ресурсов при осуществлении планирования ВЧ радиосвязи.

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

Автор считает своим приятным долгом выразить благодарность научному руководителю д.ф.-м. н. Мальцевой О.А. и куратору д.т.н., профессору Барабашову Б.Г. за помощь и постоянное внимание при выполнении работы, а также сотрудникам РГУ Анишину М.М. и Родионовой В.Т. за помощь в проведении исследований.

1. Advances of Space Research, 2001,- V.27. -№1. P. 1-180.

2. Angling M.J., Cannon P.S. Assimilation of radiooccultation measurements into background ionospheric models //Radio Sci. 2004. - V.39. - doi: 10.1029/2002RS002819.

3. Angling M.J., Khattatov B. Comparation study of two assimilative models of the ionosphere //Radio Sci. 2006. - V.41. - doi: 10.1029/2005RS003372.

4. Araujo-Pradere E.A., Fuller-Rowell T.J. STORM: An empirical storm-time ionospheric correction model. 2. Validation //Radio Sci. 2002. - V.37. - doi: 10.1029/2002RS002620.

5. Araujo-Pradere E.A., Fuller-Rowell T.J., Codrescu M.V. STORM: An empirical storm-time ionospheric correction model. 1.Model description //Radio Sci. 2002. -V.37.-doi: 10.1029/2001RS002620.

6. Belehaki A., Jakowski N. Studies of the ionospheric ionization over Athens using ground ionosonde and GPS-derived TEC values //Acta Geod. Geophys. Hung., 2002. -V.37.-P. 163-170.

7. Belehaki A., Tsagouri 1. Investigation of the relative bottomside/topside contribution to the total electron content estimates //Annals of Geophysics. 2002. - V.45. -№1. - P.73-86.

8. Bilitza D. International Reference Ionosphere Status 1995/1996. //Adv. Space Res. - 1997. - V.20. - №9. - P. 1751-1754.

9. Bilitza D. International Reference Ionosphere. //Radio Sci. 2001. - V.36. - №2.

10. Bilitza D. International Reference Ionosphere2000: Examples of improvements and new features //Adv. Space Res. 2003. - V.31. - №3. - P.757-767.

11. Bilitza D. The International Reference Ionosphere 1990 //NSSDC/WDC-R2S Report, National Space Science Center, Greenbelt, MD, 1990. P.90-92.

12. Bilitza D.(ed.) International Reference Ionosphere 1990 //NSSDC90-22, Green-belt, Maryland, 1990.

13. Bishop G.J., Daniel R.E.Jr., Rao S., Mazzella A.J.Jr., Bullett T.W. An application of parameterized real-time ionospheric specification model to regional ionospheric specification //Radio Sci. 2001. - V.36. - P. 1255-1263.

14. Breed A.M., Goodwin G.L. Ionospheric slab thickness and total electron content determined in Australia //Fremante AIP Congress paper, Sep/Oct 1998, Fremante. -P.1-17.

15. Broms M., Zolesi В. Variations in the ionosphere seen by oblique soundings //Proc. of COST238/PRIME Workshop, Sept. 1994, El Arenosillo, Spain. 1995. -P.45-55.

16. Brown L.D., Daniell R.E., Fox M.W., Klobuchar J.A., Doherty P.H. Evaluation of six ionospheric models as predictors of total electronic content. //Radio Sci. 1991. -V.26. -№4. - P. 1007-1015.

17. Cander L.R., Hickford J., Tsagouri I., Belehaki A. Real-time dynamic system for monitoring ionospheric propagation conditions over Europe //Electr. Lett. 2004. -V.40. - №4.

18. Chasovitin Yu.K., Gulyaeva T.L., Deminov M.G., Ivanova S.E. Russian standard model of ionosphere (SMI). //Proc. of the 2nd COST251 Workshop. Side. Turkey, March 1998. RAL. UK.

19. Ciraolo L. Evaluation of GPS L2-L1 biases and related daily TEC profiles. //Workshop on Modelling the Ionosphere for GPS applications, Sep. 29-30, 1993, Neustrelitz.

20. Ciraolo L., Spalla P. Comparison of ionospheric total electron content from the Navy Navigation Satellite System and GPS //Radio Sci. 1997. - V.32. - №3. -P.1071-1080.

21. Conkright R.O., Davies K., Musman S. Comparisons of ionospheric total electron contents made at Boulder, Colorado, using the Global Positioning System //Radio Sci. 1997. - V.32. - №4. - P. 1491 -1497.

22. Daniell R.E., Brown L.D., Anderson D.N., Fox M.W.,Doherty P.H., Decker D.T. Sojka J.J., Schunk R.W. Parametrized ionospheric model: a global ionospheric pa-rametrization based on first principles models //Radio Sci. 1995. - V.30. - P. 14991510.

23. Davies K., Hartmann G.K. Studying the ionosphere with the Global Positioning System //Radio Sci. 1997. - V.32. - №4. - P. 1695-1704 .

24. Dick M.I., Gulyaeva T.L. Testing of the accuracy of existing long-term ionospheric maps of foF2 //3rd COST251 Workshop Proc., El Arenosillo, Spain, COST251TD(99)003. 1999. - P.325-335.

25. Ducharme E.D., Petrie L.E., Eyfrig R. A method for predicting the F1-layer critical frequency. //Radio Sci. 1971. - V.6.

26. Essex E.A., Webb P.A., Horvath I., McKinnon C., Shilo N., Tate B. Monitoring the ionosphere/plasmasphere with low Earth orbit satellites: The Australian microsa-tellite FEDSAT. //ANARE Notes. 1999. - P. 1-6.

27. Ezquer R.G., Ortiz de Adler N., Heredia T. Predicted and measured total electron content at both peaks of the equatorial anomaly //Radio Sci. 1994. - V.29. - P.831-832.

28. Fact Sheet. U.S. Global Positioning System Policy, The White House, Office of Science and Technology Policy, National Security Council, March 29, 1996.

29. Feltens J. The International GPS service (IGS) ionosphere Working Group presented by P.H.Doherty, 2002.

30. Fuller-Rowell Т., Araujo-Pradere E., Cogrescu M. //Adv. Space Res. 2000. -V.25 - P.138-148.

31. Gaussiran T.L., Coker C., Coco D. FoF2 correlation with GPS derived TEC //Paper presented at National Radio Science Meeting, Union Radio Sci. Int., Boulder, Colo., Jan. 9-13, 1996.

32. Gulyaeva T.L. Progress in ionospheric informatics based on electron-density profile analysis of ionograms. //Adv. Space Res. 1987. - V.7. - №6.

33. Gulyaeva T.L., Huang X., Reinisch B.W. Ionosphere-plasmassphere model software for ISO. //Acta Geod. Geoph. Hung. 2002. - V.37. - P. 143-152.

34. Gulyaeva T.L., Kishcha P.V., Makarova N.R. AE-Storms and sub-storms occurrence: 1957-1990.

35. Hanbaba H.(ed.) Improved Quality of service in Ionospheric Telecommunication Systems Planning and Operation //COST action 251 Final Report, Space Research Center, Warsaw. 1999. - 302 p.

36. Hedin A.E. MSIS-86 Thermospheric Model. //J. Geophys. Res. 1987. - v.92. P .4649-4662.

37. Huang X., Reinisch B.W. Vertical electron content from ionograms in real time //Radio Sci. 2001. - V.36. - P.335-342.

38. Huang X., Reinisch B.W. Vertical electron density profiles from digisonde ionograms. The average representative profile-ARP. //Proc. of COST238/PRIME Workshop, Sep. 1994, El Arenosillo, Spain. 1995. - P.73-86.

39. Imel D.A. Evaluation of TOPEX/Poseidon dual-frequency ionosphere correction //J.Geophys. Res. 1994. - V.99. - NCI2. - P.24895-24906.

40. ION Newsletter. Spring 2000. - V. 10. - № 1.

41. ITU-R Reference Ionospheric Characteristics //Recommendation ITU-R P. 1239. 1997.-P.318-330.

42. Jaitinger., Zhang M.L.,Radicella S.M. An improved bottomside for the ionospheric electron density model NeQuick //Annals of Geophysics. 2005. - V.48. -№3. - P.525-533.

43. Jakowski N., I leise S., Wehrenpfennig A., Schluter S. TEC monitoring by GPS-A possible contribution to space weather monitoring //Phys. Chem. Earth C. 2001. -V.26. - №8. - P.609-613.

44. King J.W., Slater A.J. Errors in predicted values of foF2 and hmF2 compared with the ob-served day-to-day variability. //Telecomm. J. 1973. - V.40. - p.766-770.

45. Komjathy A., Langley R.B. Improvement of a global ionospheric model to provide ionospheric range error corrections for single-frequency GPS users. //Presented at the ION 52nd Annual Meeting, Jun. 19-21, 1996, Cambridge, MA. Юр.

46. Kouris S.S., Muggleton L.M. Diurnal variation in the E-layer ionization. //J. Atm. Terr. Phys. 1973.-V.35.

47. Liu R.Y., Smith P.A., King J.W. A new solar index which leads to improved foF2 predictions using the CCIR Attlas. //Telecom. J. 1983. - V.50. - №8. - P.408-414.

48. Lunt N., Kersley L., Bailey G.J. The influense of the protonosphere on GPS observations: model simulations. //Radio Sci. 1999. - V.34. -№3. - P.725-732.

49. Lunt N., Kersley L., Bishop G.J., Mazzella A.J. Bailey G.J. The effect of the protonosphere on the estimation of GPS total electron content: Validation using model simulations. //Radio Sci. 1999. - V.34. - №5. - P. 1261-1271.

50. Lunt N., Kersley L., Bishop G.J., Mazzella A.J. Bailey G.J. The protonospheric contribution to GPS total electron content: Two-station measurements. //Radio Sci. -1999. V.34 - №5. - P. 1281 -1285.

51. Lunt N., Kersley L., Bishop G.J., Mazzella A.J. The contribution of the protonosphere to GPS total electron content: Experimental measurements. //Radio Sci. -1999. V.34. - №5. - P. 1273-1280.

52. Maltseva O.A., Barabashov B.G., Rodionova V.T., Shlyupkin A.S. One of possibilities to use TEC-GPS data for HF links //Proceedings of The first European Conference on Antennas and Propagation (EuCAP 2006), Nov. 6-10, 2006, Nice, France. -P02.12. P.309.

53. Maltseva O.A., Barabashov B.G., Rodionova V.T., Shlyupkin A.S. One of possibilities to use TEC-GPS data for HF links //Special Publication, EuCAP 2006, Nov. 610,2006, Nice, France. SP-626 on CD. - P. 1-5.

54. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Evaluation of real time TEC and STORM-TIME corrections in the IRI model //EGU2005, Apr. 25-29, 2005, Vienna, Austria: Geophysical Research abstracts. V.7. - EGU05-A-03017. - P. 1-3.

55. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Evaluation of real-time TEC and STORM-TIME correction in the IRI //Proceedings of the XXVIII General Assembly URSI, Oct. 23-29, 2005, New Daly, India. 601 b.2(0147) - P.244-249.

56. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Space Weather for HF-communication links //Book of abstracts, 5th International Conference on Problems of Geocosmos, May 24-28,2004, St. Petersburg. P.256.

57. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Space Weather for HF-communication links //Proceedings of the 5th International Conference "Problems of Geocosmos", 2004, Saint-Petersburg. P.355-358.

58. Maltseva O.A., Poltavsky O.S., Shlupkin A.S. Using Total Electron Content of the Ionosphere for Forecasting Its Critical Frequency //International Conference on Problems of Geocosmos, May 23-27, 2006, St.-Petersburg: Book of abstracts. -P.160-161.

59. Maltseva O.A., Shlupkin A.S. The model IRI residual error and the new method of N(h) profile determination //EGU2006, Apr. 2-6, 2006, Vienna, Austria: Geophysical Research abstracts. V.8. - N 01180. - P. 1-3.

60. Miro G., de la Morena В., Marin D., Herraiz M., Farelo A.F., Gulyaeva T.L. Round-the-world propagation of the ionosphere disturbances at low middle latitudes. //3rd COST251 Workshop Proc, El Arenosillo, Spain. 1999. - CC)ST251TD(99)003.- P.169-174.

61. Miro G., Jakowski N., de la Morena B.A. Equivalent slab thickness of the ionosphere in middle latitudes based on TEC/foF2 observations over El Arenosillo. //3rd COST251 Workshop Proc. El Arenosillo , Spain. COST251 TD(99)003. - P.87-92.

62. Oliver M.A., Webster R. Kriging: A method of interpolation for geographical information systems. //Int. J. Geographical Information Systems. 1990. - V.4. - №3. -P.313-332.

63. Radicella S.M., Zhang M.L. The improved DGR analytical model of electron density height profile and total electron content in the ionosphere. //Annali di Geofisica. 1995. - V.38. -№1. - P.35-41.

64. Rishbeth H., Mendillo M. Patterns of F2-layer variability. //J.Atmos. Solar Terr. Phys. 2001. - V.63. -№15. - P.1661-1680.

65. Rush C.M., Miller D., Gibbs J. The relative daily variability of foF2 and hmF2 and their implications for HF radio propagation. //Radio Sci. 1974. - V.9. - №8/9. -P.749-756.

66. Rush C.M., Po Kempner M., Anderson D.N., Perry J., Stewart F.G.,Reasoner R. Maps of foF2 derived from observations and theoritical data. //Radio Sci. 1984. -V. 19.-№4.-P. 1083-1097.

67. Samardjiev Т., Bradley P.A., Cander Lj.R. and Dick M.I. Ionospheric mapping by computer contouring techniques. //Electronics Letters. 1993. - V.29. - №20. -P. 1794-1795.

68. Sardon E., Zarraoa N. Estimation of total electron content using GPS data: How stable are the differential satellite and receiver instrumental biases //Radio Sci. 1997.- V.32.-№5.-P.1899-1910.

69. Schunk R W, Sojka J J and Bowline M D Theoretical study of the electron temperature in the highlatitude ionosphere for solar maximum and winter conditions Hi Geophys Res. 1986. - V.91. - P. 12041.

70. Secan J.A., Wilkinson P.J. Statistical studies of an effective sunspot number. //Radio Sci. 1997. - V.37. - №4. - p. 1717-1724.

71. Stanislawska I., Bradley P.A., Juchnikowski G. Spatial correlation assessment of ionospheric parameters for limited-area mapping //GA URSI, Maastricht, Aug. 17-24, 2002.-№816.-4p.

72. Stanislawska I., Tulunay Y., Gulyaeva T.L. Transportable ionosonde in PRIME project. //Proc. of COST238/PRIME Workshop, Sep. 1994, El Arenosillo, Spain. -1995. P.91-100.

73. Stankov S.M. Evaluation of analytical ionospheric models used in electron density profile reconstruction //Acta Geod. Geoph. Hung. 2002. - V.37. - №2-3. -P.385-401.

74. Stankov S.M., Jakowski N., Heise S. Muntarov P., Kutiev I., Warnant R. A new method for reconstruction of the vertical electron density distribution in the upper ionosphere and plasmasphere //J.Geophys. Res. 2003. - V.108. - №A5.

75. Stewart F.G., Leftin M. Relationship between Ottawa 10.7cm solar radio noise flux and Zurich sunspot number //Telecomm. J. 1972. - V.39. - P. 159-169.

76. Sugiura M., Chapman S. The average morphology of geomagnetic storms with sudden commencement. //Abandl. Akad. Wiss. Gottingen Math. Phys. 1980. - Kl. -p.4.

77. Vice President Gore Announces, New Global Positioning System Modernization Initiative, The White House, Office of Vice President, January 25, 1999.

78. Warnant R. Reliability of the TEC computed using GPS measurements the problem of hardware biases //Acta Geod. Geoph. Hung. - 1997. - V.32. - №3-4. -P.451-459.

79. Webb P.A., Essex E.A. Electron density measurements of the plasmasphere -experimental observations and modeling studies. //Proc. Of WAPS'2000, La Trobe Univ., 2000. -P.98-103.

80. Wilkinson P., Wu J., Du J., Wang Y.-J. Real-time total electron content estimated using the international reference ionosphere //Adv. Space Res. 2001. - V.27. -P. 123-126.

81. Wilkinson P.J. Predictability of ionospheric variations for quiet and disturbed conditions. //J. Atm. Terr. Phys. 1995. - V.57. - P. 1469-1481.

82. Zarraoa N., Sardon E. Test of GPS for permanent ionospheric TEC monitoring at high latitudes //Ann. Geophysicae. 1996. -№14. - P. 11-19.

83. Ануфриева T.A., Шапиро Б.С. Геометрические параметры слоя F2 ионосферы. // М.: Наука, 1976. 90с.

84. Барабашов Б.Г., Анишин М.М. Прогнозирование характеристик декамет-ровых каналов связи //Труды НИИР. 2002. - С.99-101.

85. Барабашов Б.Г., Мальцева О.А. Возможности описания ионосферы в реальном времени и ее прогнозирования. //Изв. ВУЗов Сев.-Кавк. регион, естеств. науки. 2003. - №2. - С.33-38.

86. Барабашов Б.Г., Мальцева О.А. Ионосферное обеспечение однопозицион-ных пеленгаторов-дальномеров диапазона декаметровых волн //Труды НИИР, 2003. С.120-126.

87. Барабашов Б.Г., Мальцева О.А., Родионова В.Т.Контроль состояния ВЧ радиоканала в возмущенных условиях на основе модели IRI //Труды НИИР, 2005. С.140-145.

88. Васильев К.Н. и др. //Сб. "Ионосфера и солн.- земн. Связи", Алма-Ата. -Изд-во Наука. Каз. ССР. - 1977. - №3.

89. Голыгин В.А., Грозов В.П., Лесюта О.С., Сажин В.И. Коррекция критических частот в модели ионосферы по данным о полном электронном содержании. //БШФФ, 2003. -с.58-60.

90. ГОСТ 35645.140-86-ГОСТ 25645.142-86. Ионосфера Земли. Модели распределения электронной концентрации. 1986.

91. Иванов В.А., Рябова Н.В., Урядов В.П., Шумаев В.В. Прогнозирование и экстраполяция параметров КВ-радиоканала по данным наклонного зондирования ионосферы. //Радиотехника. 1997. - №7. - С.28-30.

92. Иванов В.Ф., Мысливцев Т.О., Ткачев Е.А., Троицкий Б.В. Мониторинг параметров ионосферы по данным навигационных измерений в СРНС ГЛО-HACC/GPS. //XXIY симпозиум по радиолокационному зондированию природных сред, 2006, С.Петербург. 8с.

93. Интерфейсный контрольный документ GPS, ICD-200C-002, 25.09.97, www.navcen.uscg.mil/geninfo/gpsdocuments/icd200/icd200c.pdf

94. Ионосфера Земли. Модель глобального распределения концентрации, температуры и эффективной частоты соударений электронов. //М., Изд-во стандартов, 1990.-812 с.

95. Котяшкин С.И. Определение ионосферной задержки сигналов в одночас-тотной аппаратуре потребителей спутниковой системы навигации NAVSTAR. //Зарубежная радиоэлектроника. 1989. - №5. - С.85-95.

96. Крашенинников И.В., Егоров И.Б., Коломийцев О.П., Черкашин Ю.Н. //Геомагнетизм и аэрономия. 2004. - Т.44. - №2. - С.221-226.

97. Мальцева О.А., Полтавский О.С., Шлюпкин А.С. Использование GPS-данных в каналах КВ-связи //Электросвязь. 2006. - №11.- С.47-49.

98. Мальцева О.А., Полтавский О.С., Шлюпкин А.С. Эффективность модели ионосферы IR1 при определении условий распространения в ВЧ канале //Электромагнитные волны и электронные системы. 2006. - №9. - С.21-24.

99. Мальцева О.А., Родионова В.Т., Шлюпкин А.С. Использование полного электронного содержания для текущей диагностики состояния ионосферного канала //Геомагнетизм и аэрономия. 2005. - №4. - С.480-486.

100. Мальцева О.А., Родионова В.Т., Шлюпкин А.С. Об использовании параметра ПЭС для текущей диагностики состояния ионосферного канала //Тезисы X региональной конференции по распространению радиоволн, С-Петербург, 2004.- С.54-57.

101. Мальцева О.А., Шлюпкин А.С. Остаточная погрешность модели IRI и новый метод определения Ы(Ь)-профилей ионосферы //Труды XI региональной конференции по распространению радиоволн, С.- Петербург, 2005. -С.46-48.

102. Математическая обработка результатов эксперимента. Справочное руководство. М., Наука. 1971.

103. Михайлов А.В., Булденкова С.Д., Михайлов В.В., Терехин Ю.Л. Сопоставление индексов солнечной активности в целях моделирования медианных значений foF2 //Геомагнетизм и аэрономия. 1990. -Т.30. -№1. -С.113-120.

104. Нишида А. Геомагнитный диагноз магнитосферы. //М.: Мир, 1980. 304с.

105. Ортиков М.Ю., Троицкий Б.В. Применение ионосферного индекса солнечной активности для расчета параметров ионосферы при обеспечении КВ-радиосвязи //Геомагнетизм и аэрономия. 2000. - Т.40. - №4. - С. 135-136.

106. Островский Г.И. Экспертная оценка точности определения параметров максимума Р2-области ионосферы. // Геомагнетизм и аэрономия. 1992. - Т.32. -№6. - С.78-83.

107. Поляков В.М., Суходольская В.Е., Ивельская М.К. и др. Полуэмпирическая модель ионосферы для широкою диапазона геофизических условий. М.: МЦД-Б, 1986.- 136 с.

108. Соловьев Ю.А. Спутниковая навигация и ее приложения. /ЭКО-ТРЕ11ДЗ, Москва, 2003.

109. Часовитин Ю.К. // Ионосф. исслед. М.: МГК, 1988.-№44. - С.6-13.

110. Часовитин Ю.К., Иванова С.Е., Максимова В.В., Сыкилинда Т.Н. Совре-мен-ные эмпирические модели ионосферы и пути их усовершенствования. //Ионосферные исследования. М., 1987 - №42. - С.38-51.

111. Шлюпкин А.С. Использование навигационных спутников для радиосвязи //Тезисы 4-й международной научно-практической конференции "Телекоммуникационные технологии на транспорте России", ТелеКомТранс 2006, Сочи, 19-21 апреля 2006. С.36-37.

112. Шлюпкин А.С., Мальцева О.А., Полтавский О.С. Использование информации навигационных спутников в каналах KB радиосвязи //Сборник тезисов ВНКСФ-10, Москва, 1-7 апреля 2004. С.989-991.

113. Автору данной диссертации принадлежат следующие публикации: 6, 7, 71 -79,105,106,119-123,133,134.