Бесплатный автореферат и диссертация по наукам о земле на тему
Отражение радиосигналов от слоистых структур в атмосфере и ионосфере
ВАК РФ 25.00.29, Физика атмосферы и гидросферы

Автореферат диссертации по теме "Отражение радиосигналов от слоистых структур в атмосфере и ионосфере"

МОСКОВСКИЙ ОРДЕНА ЛЕНИНА, ОРДЕНА ОКТЯБРЬСКОЙ РЕВОЛЮЦИИ, ОРДЕНА ТРУДОВОГО КРАСНОГО ЗНАМЕНИ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ им МВ ЛОМОНОСОВА

ФИЗИЧЕСКИЙ ФАКУЛЬТЕТ

На правах рукописи УДК 550 388 2+537 874 2

Голикова Елена Владимировна

Отражение радиосигналов от слоистых структур в атмосфере и ионосфере

Специальность 25 00 29 — физика атмосферы и гидросферы

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

ООЗ166621

Москва 2008

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

Научный руководитель

доктор физико-математических наук, профессор Куницын Вячеслав Евгеньевич

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

академик Маслов Виктор Павлович,

доктор физико-математических наук Гаврилов Борис Георгиевич

Ведущая организация

Институт земного магнетизма, ионосферы и распространения радиоволн им Н В Пушко-ва Российской Академии Наук

Защита состоится «17» апреля 2008 г в^О^часов на заседании диссертационного совета Д 501 001 63 при Московском государственном университете имени М В Ломоносова по адресу 119991, Москва, Ленинские горы, дом 1, стр 2, Физический факультет, аудитория СФ/1

С диссертацией можно ознакомиться в библиотеке физического факультета МГУ им М В Ломоносова

Автореферат разослан « марта 2008 г

Ученый секретарь

Диссертационного совета Д 501 001 63 кандидат физико-математических наук

В Б Смирнов

Общая характеристика работы

Актуальность темы Задача распространения радиоволн в слоистой ионосфере имеет большое значение для расчета радиотрасс и для методов дистанционного зондирования ионосферной плазмы Основным средством решения данной задачи долгое время были асимптотические методы, основанные на приближении геометрической оптики [1, 2, 3], что позволяло прогнозировать отражение радиоволн для плавно меняющихся профилей электронной концентрации Однако такой подход к решению данной задачи не применим при наличии тонких слоев, больших градиентов электронной концентрации, недостаточен для интерпретации многих экспериментальных данных, связанных с волновыми явлениями, и т д

С ростом современных компьютерных технологий стало возможным использовать эффективные численные методы, позволяющие решать задачи моделирования отражения радиосигналов от слоев, содержащих шум, от тонких спорадических слоев, слоев с разрывами функции плотности электронной концентрации и ее производных Первые результаты по расчету поля электромагнитной волны получены в работах [4, 5, 6] на основе численного решения краевых задач Однако в первых работах были проведены расчеты лишь для тонких (несколько длин волн) ионосферных слоев Численные результаты по отражению радиосигналов и расчету ионограмм для произвольных слоев были впервые получены в работах [7, 8, 9] Важной проблемой при этом является неустойчивость вычислений

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

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

ты и соотносить их с экспериментальными данными Для расчета произвольных слоев были предложены другие подходы к решению задачи отражения радиоволн, в частности, переход к новым переменным [7, 8], использование точных решений в случае представления потенциала в виде набора дельта-функций или набора прямоугольников [10, 11] Данные подходы позволили переформулировать задачу и добиться ее устойчивости В работах [12, 13] получены аналитические выражения для обобщенных коэффициентов отражения и прозрачности произвольной сферически-многослойной изотропной ионосферы и проведено моделирование сигнала, отраженного от многослойной ионосферы

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

Цели работы Сравнение численных методов расчета коэффициента отражения радиоволн на примере толстых немонотонных слоев Выбор эффективного численного метода решения прямой задачи расчета коэффициента отражения радиоволн от слоистых структур в атмосфере и ионосфере Исследование распространения радиоволн в ионосфере с учетом неоднородной тропосферы на модельных профилях показателя преломления Расчет коэффициента отражения радиоволн от нижней атмосферы при наклонном распространении волн, близком к горизонтальному Численное моделирование отражения радиоимпульсов от различных ионосферных профилей, в том числе с учетом случайного поля неоднородностей

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

Научная новизна Проведен анализ и сопоставление методов расчета коэффициента отражения на примере толстых немонотонных слоев

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

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

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

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

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

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

при наклонном распространении радиоволн, близком к горизонтальному (для углов порядка 89°) Увеличение угла падения радиоволны приводит к возрастанию амплитуды и частотного масштаба осцил-ляций Наблюдаемые волновые интерференционные эффекты соответствуют отражению радиоволны от неоднородных структур тропосферы

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

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

Апробация работы Результаты диссертации доложены на научных конференциях, симпозиуме

• VII Международная конференция студентов и аспирантов по фундаментальным наукам «Ломоносов-2000» Москва, 2000

• Ы1 Научная сессия, посвященная Дню Радио Москва, 2001

• XX Всероссийская конференция по распространению радиоволн Нижний Новгород, 2002

• Международный симпозиум «Атмосферная радиация» (МСАР-2006) Санкт-Петербург, 2006

• XVIII сессия Российского акустического общества Таганрог, 2006

Публикации Основные результаты диссертации изложены в 7 опубликованных работах, список которых приведен в конце автореферата

Структура и объем диссертации Диссертация состоит из введения, четырех глав, заключения и списка литературы, содержащего 104 наименования Объем диссертации составляет 111 страниц, включая 47 рисунков и таблицу

Содержание диссертации

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

В первой главе рассматриваются задачи распространения волн в слоистых средах и методы их решения Электрическое поле плоской гармонической волны описывается волновым уравнением, имеющим вид одномерного уравнения Гельмгольца

сР-Е

+ = 0, (1)

7 Ш

где к =--волновое число, ш - циклическая частота волны, с - скорость

с

света, е - диэлектрическая проницаемость В случае плоскослоистой изотропной плазмы задача распространения электромагнитной волны сводится к решению одномерного стационарного уравнения вида

j2 тр

a^[+(k*-q(x,u))E = 0 (2)

, Л e2N(x)

а(х,ш) =-q ~т~-г-г, (3)

' Едтс2 (1 + г^эфф/ш) v '

где N(x) - электронная концентрация, ео - диэлектрическая постоянная, ей m - заряд и масса электрона, г/эфф - эффективная частота соударений При отсутствии поглощения потенциал (3) становится действительным и не зависит от частоты Показатель преломления атмосферы п = у/е мало превышает единицу [14]

с потенциалом

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

Е(х, к) = егкх + Я{к)е-Лх при х -оо, Е(х, к) = Т(к)егкх при х -)• оо,

где Д(&) - комплексный коэффициент отражения содержит в себе как амплитуду, так и фазу отраженной волны, Т(к) - комплексный коэффициент прохождения Решением является нахождение коэффициентов отражения и прохождения, когда свойства среды известны

Сделан обзор моделей, для которых можно найти точное решение одномерного волнового уравнения (2) Рассмотрены существующие на сегодняшний день различные подходы к решению задачи отражения радиоволн от произвольного ионосферного профиля с учетом волновых явлений Приведен обзор работ по исследованию закономерности отражения радиоволн от различных ионосферных профилей

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

новой переменной У(х, к) = (¿Е (ж, к) /¿х) /Е (х, к)

^ + + (5)

с граничным условием позади слоя (х = а) У(а, к) — гк Аналогично, переход к усеченному потенциалу сводит волновое уравнение (2) к дифференциальному уравнению первого порядка вида

<Ш{Х, к) = _2гШ{х| к) _ Щ (1 + д ( к))2 (6)

ах 2 г к

с граничным условием Я(а, к) = О

Другим способом решения прямой задачи является использование точных решений, которые существуют в случае представления потенциала в виде набора дельта-функций д(х) = ^п (х — хп) или набора прямоугольников Тогда задача расчета коэффициента отражения сводится к вычислениям по рекуррентной формуле с граничным условием

Rn+i — 0 Расчет коэффициента отражения методом прямоугольников проводится по рекуррентной формуле

R _ (Рп - Рп+1) + Rn+ie2iPn+lAXn (рп + рп+1) п (Рп + Pn+i) + Rn+ (рп - pn+l) 1 '

где Ахп = xn+i -xn,pn= л/к2 - q~

Однако, при вычислении коэффициента отражения в области малых значений рп возникают трудности В этом случае предлагается коэффициент отражения рассматривать как Rn — — 1 + еп, где еп - малая добавка, равная

= 2рп (1 + Rn+ie2iPn+lAXn)

£п ~ (Рп+Рп+1) + iW2!p"+lA*" {Рп ~ Pn+l)

Исследование точности вычислений рассматриваемых методов для моделирования отражения радиоволн от ионосферы при вертикальном зондировании проводились путем численного эксперимента на примере полных ионосферных профилей (например слой F, величина потенциала которого составляет 1011 — 1012) За величину потенциала принята безразмерная величина, равная произведению характерного значения потенциала Q на квадрат его характерной ширины а , т е Qa2

Алгоритмы численных решений дифференциальных уравнений первого порядка (уравнения (5) и (6)) основаны на методе Рунге-Кутты 2-го и 4-го порядка точности Точность и эффективность численного расчета проверялась повторным расчетом коэффициента отражения с меньшим шагом сетки Отметим, что величина шага сетки бралась существенно меньше длины волны

На рис 1 представлена зависимость модуля погрешности расчета коэффициента отражения |ДД| от величины шага сетки Л/Л (здесь и далее Л - длина волны) для гдадкого унимодального (косинусного) слоя д(ж) = Q (1 — cos (2тгх/а)) /2 Порядок точности метода соответствует углу наклона кривой \AR(h/X)\ в логарифмических координатах

Приведенные результаты демонстрируют преимущество метода прямоугольников для данной задачи перед другими методами 2-го порядка точности Сравнение численных методов показало, что методы, основанные на уравнениях (5) и (6), имеют близкие значения погрешности коэффициента отражения

5

10"!

■4

ых

Рис. 1: Зависимость модуля погрешности расчета коэффициента отражения от величины шага для косинусного потенциала при (¡а = 1012, ка = 10е: 1, 2 - метод Рунге-Кутты 2-го и 4-го порядков для усеченного потенциала; 3,4- метод Рунге-Кутты 2-го и 4-го порядков для логарифмической производной; 5 - метод дельта-потенциалов; 6 - метод прямоугольников.

Оптимальный шаг сетки, необходимый для достижения погрешности коэффициента отражения не более 1% составляет /г/А ~ 10~2 для метода Рунге-Кутты второго порядка и для метода дельта-потенциалов, /г/А ~ 3 • 1СГ2 для Рунге-Кутты четвертого порядка и /г/А ~ Ю-1 для метода прямоугольников.

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

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

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

ления Показано, что заметная разница проявляется на частотах ниже 50 кГц На частотах порядка 10 МГц разница составляет не более 0,001% Учет тропосферы приводит к изменению фазы коэффициента отражения на величину порядка 7г

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

Ь, км Ь км

Рис 2 Примеры экспериментальных высотных профилей показателя преломления воздуха для радиоволн в летний период (14, 15]

Рис 3 Модуль коэффициента отражения при падении радиоволны на слой тропосферы под углом в = 85° и в = 89°

Изучены основные особенности наклонного распространения радиоволн в приземном слое тропосферы на примере высотных профилей показателя преломления, полученных по данным метеонаблюдений (рис 2) [14, 15]

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

По реальным профилям показателя преломления тропосферы проведены вычисления модуля коэффициента отражения радиоволн для наклонного зондирования На рис 3 представлены результаты расчета модуля коэффициента отражения |Д(/)| при падении радиоволны на слой тропосферы под углом в = 85° и в = 89° Увеличение угла падения радиоволны в приводит к возрастанию частотного масштаба осцилля-ций и амплитуды колебаний Величина модуля коэффициента отражения для в = 85° имеет порядок 10~3, а для в = 89° на порядок выше, 10~2 (рис 3) Чем более изрезанная форма профиля показателя преломления, тем сложнее вид имеет частотная зависимость коэффициента отражения

350 N. 1-1-ед

30,0 £ МГц

Рис. 4: Модельный высотный профиль показателя преломления воздуха в условиях температурных инверсий, модуль коэффициента отражения при разных углах падения радиоволны на слой тропосферы.

На рис. 4 приведены результаты расчета модуля коэффициента отражения |Л(/)| при разных углах падения радиоволны на модельный высотный профиль показателя преломления воздуха в условиях температурных инверсий. Пик коэффициента отражения соответствует случаю, когда длина волны совпадает с характерным размером неоднородностей. Для угла в — 88° максимум модуля коэффициента отражения лежит на частоте /накл ~ И-5 МГц и соответствующая ей длина волны А ~ 750 м совпадает с протяженностью неоднородностей Ь ¡=а 800 м (рис.4).

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

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

В работе промоделирован профиль электронной концентрации со

спорадическим слоем, имеющим форму параболы (рис. 5а). Особенностью параболического слоя является наличие разрыва первой производной электронной концентрации на его границах. Для спорадического слоя толщиной 500 м получена ионограмма вертикального зондирования (рис. 56). Ионограмма представляет собой двумерное распределение действующей высоты отражения К сигнала от несущей частоты /. Амплитуда отраженного импульса прямо пропорциональна интенсивности черного цвета на ионограмме.

Первый след на действующей высоте 100 км соответствует сигналу, отраженному от спорадического Ев-слоя. Второй и третий по высоте следы представляют сигналы, отраженные от Е- и ^-слоев. Поскольку в данной задачи поглощение не учитывается, на ионограмме наблюдаются следы от многократного переотражения сигнала в области Еа — Е долины (рис. 5 6). Уменьшение толщины спорадического слоя приводит к увеличению диапазона полупрозрачности. Большой диапазон полупрозрачности тонкого слоя может являться как следствием разрыва первой производной электронной концентрации, так и следствием большого градиента электронной концентрации. Таким образом, расчеты наглядно объясняют зависимость явлений, наблюдаемых при радиозондировании ионосферы, от известных свойств спорадических слоев.

Ь', км

Рис. 5: Модельный профиль электронной концентрации ионосферы с тонким спорадическим слоем в виде параболы (а) и ионограммы вертикального зондирования для профиля со спорадическим слоем шириной 500 м (6).

В работе рассмотрены потенциалы со случайным полем неоднородности, которые соответствуют представлению о характерных особенностях возмущений в ионосфере Для задания «шума» удобно выбрать следующий амплитудный пространственный спектр неоднородностей

' 0 при 0 < кп < ко,

Ь\3 /Ь

С(кп) =

2ккп) + Ъ \ 2ккп) ПРИ к°-кп-

Ь \~а

■^К I при кт < кп,

где а и Ъ - некоторые коэффициенты, Ь - интервал высот, на кото, 2тт

ром задан потенциал, кп — —=— - волновое число, соответствующее

Ь

п-й гармонике в спектре возмущений, константа а = 3 имеет смысл асимптотики энергетического спектра Спектр исключает низкие частоты, соответствующие масштабам самих ионосферных слоев, в интервале ко < кп < кт аппроксимируется полиномом третьей степени, причем функция С(&п) и ее первые производные в точках ко, кт непрерывны Созданное таким образом возмущение прибавляется к основному профилю потенциала

Профили электронной концентрации с добавлением «шума» различной амплитуды при масштабе неоднородностей 125 м и асимптотики энергетического спектра а = 3 представлены на рис 6 Величина «шума» составляет 0 5 — 7% от величины потенциала Модельные ионограммы на рис 7 наглядно демонстрируют зависимость структуры отраженных сигналов от величины неоднородностей

"Учет мелкомасштабных возмущений в ионосфере приводит к сложной структуре ионограмм При достаточно большой амплитуде неоднородностей (рис 6 б) наблюдается диффузное отражение (рис 7 б), которое принято связывать с горизонтальной неоднородностью ионосферных слоев Отметим, что диффузный характер отражения был получен нами в рамках модели плоскослоистой ионосферы, исключительно за счет вариаций вертикального профиля ионосферы Таким образом, оказывается, что диффузный характер сигналов вертикального зондирования можно объяснить наличием вертикальной неоднородности ионосферы

Рис. 6: Модельные профили электронной концентрации с добавлением «шума» различной амплитуды. Величина «шума» составляет 0.5% (а) и 7% (б) от величины потенциала.

Ь', км 1200-,

2.4 2,6 £ МГц

2.0 2,5 3,0

3.5 4.0

Г, МГц

Рис. 7: Ионограммы вертикального зондирования с характерным масштабом колебаний неоднородностей 125 м. Величина «шума» составляет 0.5% (а) и 7% (б) от величины потенциала.

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

Основные результаты и выводы

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

2 Промоделировано распространение радиоволн в ионосфере с учетом неоднородной тропосферы на модельных профилях показателя преломления Рассчитана зависимость производной фазы коэффициента отражения по волновому числу при учете влияния тропосферы Показано, что разница между величинами производной фазы коэффициента отражения с учетом и без учета тропосферы заметна на частотах ниже 50 кГц и составляет не более 0,001% на частотах порядка 10 МГц Учет тропосферы приводит к изменению фазы коэффициента отражения на величину порядка ж и модуля коэффициента отражения радиоволн на величину порядка Ю-4 Исследовано влияние неоднородностей показателя преломления воздуха на отражение радиоволн при вертикальном распространении волн в атмосфере Основной вклад в изменении коэффициента отражения вносят профили, содержащие большие градиенты показателя преломления на границах и разрывы функции показателя преломления

3 На примере высотных профилей показателя преломления, полученных по данным метеонаблюдений, продемонстрировано влияние приземного слоя тропосферы при наклонном распространении радиоволн Увеличение угла падения радиоволны приводит к возрастанию амплитуды и частотного масштаба осцилляций Величина модуля коэффициента отражения |Д(/)| при в = 89° составляет порядка 10~2 Наблюдаемые волновые интерференционные эффекты соответствуют многократному отражению радиоволны от неоднородных структур тропосферы При очень больших углах, в > 89°, становятся различимы интерференционные эффекты, связанные с отражением радиоволны от всего слоя тропосферы

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

5 Получен диффузный характер следов на ионограммах при достаточно большом по амплитуде поле неоднородностей Несмотря на то, что рассмотренная задача является одномерной по своей постановке (рассматривается плоскослоистая ионосфера), обнаруженная картина близка к той, которая наблюдается на реальных ионограммах и которая обычно объясняется поперечной неоднородностью отражающих слоев Оказывается, что диффузное отражение можно объяснить также и слоистой (вертикальной) структурой неоднородностей

Список публикаций

[1] Голикова Е В , Куницын В Е, Нестеров И А , Матвеев А С Моделирование отражения радиосигналов от слоистой атмосферы и ионосферы // Радиотехника и электроника 2005 Т50 №7 С 794-804

[2] Голикова Е В , Куницын В Е, Нестеров И А Моделирование ио-нограмм вертикального зондирования для немонотонных ионосферных слоев//Геомагнетизм и аэрономия 2008 Т48 №1 С 109-115

[3] Голикова Е В (Бочкова ЕВ) Сравнение точности методов расчета коэффициента отражения электромагнитных волн от слоистой плазмы // VII Международная конференция студентов и аспирантов по фундаментальным наукам "Ломоносов-2000" Секция "Физика" 2000 С 113-115

[4] Голикова ЕВ (Бочкова ЕВ), Нестеров И А Сравнение точности методов расчета комплексных коэффициентов отражения от слоистых сред // Ы1 Научная сессия, посвященная Дню Радио Тезисы докладов М 2001 Т2 С 283-285

[5] Голикова ЕВ (Бочкова ЕВ), Куницын ВЕ, Матвеев А С, Нестеров И А Моделирование отражения электромагнитных волн от слоистой ионосферы и нижней атмосферы // Труды XX всероссийской конференции по распространению радиоволн Нижний Новгород 2002

[6] Нестеров И А , Голикова Е В , Куницын В Е Распространение радиоволн в слоистой атмосфере и ионосфере // Международный симпозиум "Атмосферная радиация"(МСАР-2006) Сборник тезисов Санкт-Петербург 2006 С 121

[7] Голикова ЕВ К вопросу о влиянии слоистой атмосферы на отражение волн от атмосферных слоев // XVIII сессия Российского акустического общества Сборник трудов 2006 Т2 С 179-183

Список цитируемой литературы

[1] Гинзбург В Л Распространение электромагнитных волн в плазме M Наука, 1967

[2] Маслов В П Федорюк M В Квазиклассическое приближение для уравнений квантовой механики M Наука, 1976

[3] Маслов В П Асимптотические методы и теория возмущений M Наука, 1988

[4] Chessell CI Results of numerical calculations of reflection and transmission coefficients for thm highly ionized layers and their application to sporadic-i? reflections // Journal of Atmospheric and Terrestrial Physics 1971 V 33 P 1803-1822

[5] Miller К L , Smith L G Reflection of radio waves by sporadic-E layers // Journal of Atmospheric and Terrestrial Physics 1977 V 39 P 899911

[6] Nygren T A simple method for obtaining reflection and transmission coefficients and fields for electromagnetic wave m a horizontally stratified ionosphere // Planet Space Science 1981 V 29 №5 P 521528

[7] Куницын В E, Смородинов В А , Усачев А Б Отражение радиоволн от произвольного ионосферного слоя // Радиотехника и электроника 1989 ТЗЗ №2 С 233-240

[8] Куницын В Е, Усачев А Б Отражение радиоволн от немонотонных ионосферных слоев//Изв Вузов Радиофизика 1990 ТЗЗ №3 С 267-274

[9] Куницын В Е, Усачев А Б Расчет ионограмм вертикального зондирования слоистой ионосферы // Геомагнетизм и аэрономия 1993 ТЗЗ №1 С 145-148

[10] Бочкова ЕВ , Нестеров И А Сравнение точности методов расчета комплексных коэффициентов отражения от слоистых сред // Ы1 Научная сессия, посвященная Дню Радио Тезисы докладов М 2001 Т2 С 283-285

[11] Голикова Е В , Куницын В Е , Нестеров И А , Матвеев А С Моделирование отражения радиосигналов от слоистой атмосферы и ионосферы // Радиотехника и электроника 2005 Т 50 №7 С 794-804

[12] Михайлов С Я Отражение радиоволн от двухслойной изотропной ионосферы в областях прикритических частот // Радиотехника и электроника 2002 Т47 №10 С 1169-1179

[13] Михайлов С Я Моделирование сигнала, отраженного от произвольной сферически-слоистой изотропной ионосферы // Солнечно-земная физика 2005 Вып 7 С 100-109

[14] Казаков Л Я, Ломакин А Н Неоднородности коэффициента преломления воздуха в тропосфере М Наука, 1976

[15] Паршуков В А Влияние рефракции миллиметровых и сантиметровых волн в атмосфере на угол места источника излучения Дис канд ф -м наук М , 2004

Подписано к печати /V (73 08 Тираж 400 Заказ _

Отпечатано в отделе оперативной печати физического факультета МГУ

Содержание диссертации, кандидата физико-математических наук, Голикова, Елена Владимировна

Введение

1 Обзор задач распространения волн в слоистых средах

1.1 Распространение электромагнитных волн в слоистой плазме

1.2 Точные решения волнового уравнения

1.3 Известные результаты расчетов отражения радиоволн от ионосферы.

2 Методы решения прямой задачи отражения радиоволн и исследование области их применимости

2.1 Постановка прямой задачи отражения радиоволн.

2.2 Численные методы решения задачи отражения радиоволн от ионосферы.

2.3 Исследование погрешности расчета коэффициента отражения радиоволн от ионосферных слоев при вертикальном зондировании.

Введение Диссертация по наукам о земле, на тему "Отражение радиосигналов от слоистых структур в атмосфере и ионосфере"

Задачи распространения электромагнитных волн в слоистых средах имеют широкое применение в различных областях физики: в радиозондировании ионосферы и атмосферы, акустическом зондировании атмосферы, в сейсмологии, в акустике твердых тел и в подводной акустике. Для широкого класса задач электрическое поле плоской гармонической волны описывается волновым уравнением, имеющим вид одномерного уравнения Гельмгольца d2E + к2е(х,и;)Е = 0 , (1) где к =--волновое число, и> - циклическая частота волны, с - скорость света, е - дис электрическая проницаемость. В случае плазмы диэлектрическая проницаемость определяется выражением [1]: e2N(x) х,ш) — I--т-:-г,

0ти{ш + гиэфф) где N(x) - электронная концентрация, £о ~~ диэлектрическая постоянная, е и т - заряд и масса электрона, ^эфф - эффективная частота соударений. Для атмосферы показатель преломления п = \Je мало превышает единицу [2].

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

С ростом современных компьютерных технологий стало возможным использовать эффективные численные методы, позволяющие решать задачи моделирования отражения радиосигналов от слоев, содержащих шум, от тонких спорадических слоев, слоев с разрывами функции плотности электронной концентрации и ее производных. Первые результаты по расчету поля электромагнитной волны получены в работах [3]-[6] на основе численного решения краевых задач. Однако в первых работах были проведены расчеты лишь для тонких (несколько длин волн) ионосферных слоев. Численные результаты по отражению радиосигналов и расчету ионограмм для произвольных слоев были впервые получены в работах [7]-[11]. Важной проблемой при этом является неустойчивость вычислений, которая возникает в случае туннельного прохождения зондирующих волн. Для расчета произвольных слоев были предложены другие подходы к решению задачи отражения радиоволн, в частности, переход к новым переменным [8, 10], использование точных решений в случае представления потенциала в виде набора дельта-функций или набора прямоугольников [12]. Данные подходы позволили переформулировать задачу и добиться ее устойчивости. В работах [13, 14, 15] получены аналитические выражения для обобщенных коэффициентов отражения и прозрачности произвольной сферически-многослойной изотропной ионосферы и проведено моделирование сигнала, отраженного от многослойной ионосферы. Таким образом, предложенные методы дают возможность моделировать различные волновые эффекты и соотносить их с экспериментальными данными.

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

Цели диссертационной работы. Сравнение численных методов расчета коэффициента отражения радиоволн на примере толстых немонотонных слоев. Выбор эффективного численного метода решения прямой задачи расчета коэффициента отражения радиоволн от слоистых структур в атмосфере и ионосфере. Исследование распространения радиоволн в ионосфере с учетом неоднородной тропосферы на модельных профилях показателя преломления. Расчет коэффициента отражения радиоволн от нижней атмосферы при наклонном распространении волн, близком к горизонтальному. Численное моделирование отражения радиоимпульсов от различных ионосферных профилей, в том числе с учетом случайного поля неоднородностей.

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

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

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

Во второй главе рассмотрены приближенные методы решения волнового уравнения в случае плоскослоистой изотропной плазмы. Исследуется область применимости рассматриваемых методов для моделирования отражения радиоволн от ионосферы при вертикальном зондировании. Предложен оптимальный численный метод для данной задачи и на его основе проводятся дальнейшие расчеты.

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

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

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

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

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

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

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

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

По материалам диссертации опубликовано 7 научных работ [12], [16] - [21]. Результаты диссертации доложены на научных конференциях, симпозиуме:

• VII Международная конференция студентов и аспирантов по фундаментальным наукам «Ломоносов-2000». Москва, 2000.

• LII Научная сессия, посвященная Дню Радио. Москва, 2001.

• XX Всероссийская конференция по распространению радиоволн. Нижний Новгород, 2002.

• Международный симпозиум «Атмосферная радиация» (МСАР-2006). Санкт-Петербург, 2006.

• XVIII сессия Российского акустического общества. Таганрог, 2006.

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

Автор выражает благодарность профессору В. Е. Куницыну и старшему преподавателю И. А. Нестерову за помощь в работе и полезные обсуждения.

Заключение Диссертация по теме "Физика атмосферы и гидросферы", Голикова, Елена Владимировна

Основные результаты работы состоят в следующем:

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

2. Промоделировано распространение радиоволн в ионосфере с учетом неоднородной тропосферы на модельных профилях показателя преломления. Рассчитана зависимость производной фазы коэффициента отражения по волновому числу при учете влияния тропосферы. Показано, что разница между величинами производной фазы коэффициента отражения с учетом и без учета тропосферы заметна на частотах ниже 50 кГц и составляет не более 0,001% на частотах порядка 10 МГц. Учет тропосферы приводит к изменению фазы коэффициента отражения на величину порядка ж и модуля коэффициента отражения радиоволн на величину порядка Ю-4. Исследовано влияние неоднородностей показателя преломления воздуха на отражение радиоволн при вертикальном распространении волн в атмосфере. Основной вклад в изменении коэффициента отражения вносят профили, содержащие большие градиенты показателя преломления на границах и разрывы функции показателя преломления.

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

Величина модуля коэффициента отражения |Д(/)| при в = 89° составляет порядка Ю-2. Наблюдаемые волновые интерференционные эффекты соответствуют многократному отражению радиоволны от неоднородных структур тропосферы. При очень больших углах, в > 89°, становятся различимы интерференционные эффекты, связанные с отражением радиоволны от всего слоя тропосферы.

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

5. Получен диффузный характер следов на ионограммах при достаточно большом по амплитуде поле неоднородностей. Несмотря на то, что рассмотренная задача является одномерной по своей постановке (рассматривается плоскослоистая ионосфера), обнаруженная картина близка к той, которая наблюдается на реальных ионограммах и которая обычно объясняется поперечной неоднородностью отражающих слоев. Оказывается, что диффузное отражение можно объяснить также и слоистой (вертикальной) структурой неоднородностей.

Заключение

Библиография Диссертация по наукам о земле, кандидата физико-математических наук, Голикова, Елена Владимировна, Москва

1. Гинзбург В.Л. Распространение электромагнитных волн в плазме. М.: Наука, 1967.

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

3. Chessell C.I. Results of numerical calculations of reflection and transmission coefficients for thin highly ionized layers and their application to sporadic-i? reflections // Journal of Atmospheric and Terrestrial Physics. 1971. V.33. P. 1803-1822.

4. Chessell C.I., Thomas J.A., Bourne I.A. Experimental observations of the amplitudes of E.s and F-region reflections and their comparison with the thin-layer model for Es // Journal of Atmospheric and Terrestrial Physics. 1973. V.35. P. 545-561.

5. Miller K.L., Smith L.G. Reflection of radio waves by sporadic-^ layers // Journal of Atmospheric and Terrestrial Physics. 1977. V.39. P. 899-911.

6. Nygren T. A simple method for obtaining reflection and transmission coefficients and fields for electromagnetic wave in a horizontally stratified ionosphere // Planet. Space Science. 1981. V.29. №5. P. 521-528.

7. Куницын В.E., Смородинов В.А., Усачев А.Б. Коэффициент отражения радиоволн от немонотонного ионосферного слоя// Дифракция и распространение радиоволн. М.: МФТИ, 1985. С. 87-92.

8. Куницын В.Е., Смородинов В.А., Усачев А.Б. Амплитудно и фазочастотные харк-теристики отраженных от ионосферы волн // Радиотехника. 1987. №9. С. 61-63.

9. Куницын В.Е., Смородинов В.А., Усачев А.В. Отражение радиоволн от произвольного ионосферного слоя // Радиотехника и электроника. 1989. Т.ЗЗ. №2. С. 233-240.

10. Куницын В.Е., Усачев А.Б. Отражение радиоволн от немонотонных ионосферных слоев // Изв. Вузов. Радиофизика. 1990. Т.ЗЗ. №3. С. 267-274.

11. Куницын В.Е., Усачев А.Б. Расчет ионограмм вертикального зондирования слоистой ионосферы // Геомагнетизм и аэрономия. 1993. Т.ЗЗ. №1. С. 145-148.

12. Голикова Е.В., Куницын В.Е., Нестеров И.А., Матвеев А.С. Моделирование отражения радиосигналов от слоистой атмосферы и ионосферы // Радиотехника и электроника. 2005. Т.50. №7. С. 794-804.

13. Михайлов С. Я. Отражение радиоволн от двухслойной изотропной ионосферы в областях прикритических частот // Радиотехника и электроника. 2002. Т.47. №10. С. 1169-1179.

14. Михайлов С.Я., Ратовский К.Г., Медведев А.В. Моделирование сигнала, отраженного от двухслойной ионосферы в областях прикритических частот // Радиотехника и электроника. 2002. Т.47. №10. С. 1180-1186.

15. Михайлов С. Я. Моделирование сигнала, отраженного от произвольной сферически-слоистой изотропной ионосферы // Солнечно-земная физика. 2005. Вып. 7. С. 100-109.

16. Бочкова Е.В. Сравнение точности методов расчёта коэффициента отражения электромагнитных волн от слоистой плазмы // VII Международная конференция студентов и аспирантов по фундаментальным наукам "Ломоносов-2000". Секция "Физика". 2000. С. 113-115.

17. Бочкова Е.В., Нестеров И.А. Сравнение точности методов расчёта комплексных коэффициентов отражения от слоистых сред // LII Научная сессия, посвященная Дню Радио. Тезисы докладов. М. 2001. Т.2. С. 283-285.

18. Бочкова Е.В., Куницын В.Е., Матвеев А.С., Нестеров И.А. Моделирование отражения электромагнитных волн от слоистой ионосферы и нижней атмосферы // Труды XX всероссийской конференции по распространению радиоволн. Нижний Новгород. 2002.

19. Нестеров И.А., Голикова Е.В., Куницын В.Е. Распространение радиоволн в слоистой атмосфере и ионосфере // Международный симпозиум "Атмосферная радиация"(МСАР-2006). Сборник тезисов. Санкт-Петербург. 2006. С. 121.

20. Голикова Е.В. К вопросу о влиянии слоистой атмосферы на отражение волн от атмосферных слоев // XVIII сессия Российского акустического общества. Сборник трудов. 2006. Т.2. С. 179-183.

21. Голикова Е.В., Куницын В.Е., Нестеров И.А. Моделирование ионограмм вертикального зондирования для немонотонных ионосферных слоев // Геомагнетизм и аэрономия. 2008. Т.48. №. С. 109-115.

22. Дэвис К. Радиоволны в ионосфере. М.: Мир, 1973.

23. Алъперт Я. Л. Распространение электромагнитных волн и ионосфера. М.: Наука, 1972.

24. Fosterling К. Uber die Ausbreitung electromagnetischer Wellen in einem magnetisieren Medium bei senkrechter Inzidenz // Hochfrequenz-technik und Elelktroakustik. 1942. V. 59. Januar, p. 11-22.

25. Вабиков В.В. Метод фазовых функций в квантовой механике. М.: Наука, 1988.

26. Захаръев В.Н., Сузъко А.А. Потенциалы и квантовое рассеяние: Прямая и обратная задачи. М.: Энергоатомиздат, 1985.

27. Захаръев В.Н. Уроки квантовой интуиции. Дубна, ОИЯИ, 1996

28. Ландау Л.Д., Лившиц Е.М. Квантовая механика. М.: Наука, 1974.

29. Бреховских Л.М. Волны в слоистых средах. М.: Наука, 1973.

30. Руководство URSI по интерпретации и обработке ионограмм / Перевод с анг. под ред. Н.В. Медниковой. М.: Наука, 1977.

31. Гусев В.Д., Жидовленко И.Ю., Приходъко Л.И. Отражение и рассеяние радиоволн в ионосферном спорадическом слое Es // Радиотехника. 1986. №6. С. 71-73.

32. Гусев В.Д., Приходъко Л.И. Отражение радиоволн от неоднородных ионосферных слоев с поглощением // Распространение и дифракция электромагнитных волн в неоднородных средах. Тезисы докладов. М. 1994. С. 34-35

33. Усачёв А. Б. Отражение радиосигналов от слоистой ионосферы с учётом волновых явлений : Дис. . канд. ф.-м. наук. М., 1988.

34. Справочник по специальным функциям/ Под ред. М. Абрамовица и И. Стиган, М.: Наука, 1979.

35. Epstein P. Reflection of waves in an inhomogeneous absorbing medium // Proc. Nat. Acad. Science. USA. 1930. V.16. p. 627.

36. Анютин А.П. Об особенностях искажения AM и ЧМ сигналов, отражённых от слоя Эшитейна // Изв. вузов. Радиофизика. 1979. Т. 22. №6. С. 703-710.

37. Бреховских JI.M., Годин О.А. Акустика слоистых сред. М.: Наука, 1989.

38. Намазов С.А., Орлов Ю.И., Федоров Н.Н. Искажения радиоимпульсов при отражении от области максимума ионосферного слоя // Радиотехника и электроника. 1984. Т.29. №4. С. 609-619.

39. Орлов Ю.И., Федоров Н.Н. О границах применимости асимптотического описания полей радиоимпульсов вблизи критической частоты // Изв. Вузов. Радиофизика. 1984. Т.27. №. С. 1130-1135.

40. Керблай Т. С., Носова Г.Н. Применение аналитических моделей слоя Es при интерпретации ионограмм // Ионосферные модели: Сб. научн. трудов. М.:Наука, 1975. С. 169-175.

41. Гершман Б.Н., Игнатьев Ю. А., Каменецкая Г.Х. Механизмы образования ионосферного спорадического слоя Es на различных широтах. М.:Наука, 1976.

42. Фелсеп Л., Маркувиц Н. Излучение и рассеяние волн. Т.2. М.: Мир, 1978.

43. Jenkins В., Bailey G.J., Abdu М.А., Batista I.S., Balan N. Observations and model calculations of an additional layer in the topside ionosphere above Fortaleza, Brazil // Ann. Geophysicae. 1997. V.15. P. 753-759.

44. Zang D. Y. A new method of calculating the transmission and reflection coefficients and fields in a magnetized plasma layer // Radio Science. 1991. V.26. №6. P. 1415-1418.

45. Арсенин В.Я., Горячев В.А., Загонов В.П. и др. О методе расчёта коэффициентов отражения радиоволн СВ и KB диапазонов от ионосферы // Институт прикладной математики им. М.В. Келдыша. Препринт №134. 1978.

46. В.П. Крайнов, Л.П. Пресняков Фазовые функции потенциального рассеяния в оптике // Успехи физических наук. 1993. Т.163. №7. С. 85-92.

47. Ginocchio J.N. A class of exactly solvable potentials. One-dimensional Shrodinger equation // Annals of Physics. 1984. V.152. P. 203-219.

48. Маслов В.П., Федорюк М.В. Квазиклассическое приближение для уравнений квантовой механики. М.: Наука, 1976.

49. Маслов В.П. Асимптотические методы и теория возмущений. М.: Наука, 1988.

50. Крюковский А. С., Лукин Д. С. Краевые и угловые катастрофы в равномерной геометрической теории дифракции. М.: МФТИ, 1999.

51. G. Rumpker, J-M Kendall A Maslov-propagator seismogram for weakly anisotropic media // Geophys. J. Int. 2002. V.150. P.23-36.

52. X.Huang, J-M. Kendall, C.J. Thomson, G. F. West A comparison of the Maslov integral seismogram and the finite difference method // Geophys. J. Int. 1998. V.132. P.584-594.

53. J. Caddick, J-M. Kendall, D. G. Raymer Maslov shear-waveforms in highly anisotropic shales and implications for shear-wave splitting analyses // Revue de l'lnstitut Frangais du Petrole. 1998. V.53. №5. P.629-642.

54. G. Rumpker, A. Tommasi, J-M. Kendall Numerical simulations of depth-dependent anisotropy and frequency-dependent wave propagation effects // Journal of Geophysical Research. 1999. V.104. №B10. P.23,141-23,153.

55. Куличков C.H., Авилов К.В., Попов О.Е., Отрезов А.И., Буш Г.А., Барышников А.К. О некоторых результатах моделирования дальнего распространения инфразвука в атмосфере // Изв. РАН. Физика атмосферы и океана. 2004. Т.40. №2. С. 232-246.

56. Распространение звука во флуктуирующем океане / Под ред. С. Флатте. М.: Мир, 1982.

57. Нестеров И. А. Задачи зондирования слоистых сред и восстановление неунимодальных профилей : Дис. . канд. ф.-м. наук. М., 1998.

58. Михайлов С.Я. Коэффициент отражения радиоволн в областях прикритических частот многослойной ионосферы // Исследования по геомагнетизму аэрономии и физике солнца. Изд-во СО РАН, 1998. Вып. 109. Часть 1. С. 196-202.

59. Миркотан С.Ф., Бушулкин Н.И., Захаров В.И. Амплитудный и фазовый коэффициенты отражения и прохождения в области критической частоты тонкого ионосферного слоя.// В кн. Ионосферные исследования. Казань. 1997. №50. С. 130-135.

60. Гусев В.Д., Михайлова Е.Г, Приходъко Л.И. Расчет полей радиоволн и коэффициентов отражения при радиозондировании ионосферы // Радиотехника. 1999. №11. С. 22-25.

61. Булдырев B.C., Славянов С.Ю. Равномерные асимптотические разложения для решения уравнения типа Шредингера с двумя точками перехода / / Вестник ЛГУ: Физика, химия. 1968. Вып. 4. №22. С. 70-84.

62. Михайлов С.Я. Метод расчета пространственно-частотного распределения характеристик КВ-сигнала в трехмерно неоднородной ионосфере, основанный на волно-водпом подходе: Дис. . капд. ф.-м. наук. Иркутск, 1993.

63. Калиткин Н.Н. Численные методы. М.: Наука, 1978.

64. Паршуков В.А. Влияние рефракции миллиметровых и сантиметровых волн в атмосфере на угол места источника излучения : Дис. . канд. ф.-м. наук. М., 2004.

65. Паршуков В.А, Семенов А.А. Влияние облачности на рефракцию радиоволн // Изв. Вузов. Радиофизика. 1982. Т.25. №7. С. 844-846.

66. Шабелъников А.В. Рефракция электромагнитных волн в атмосфере Земли прина-личии слоистых неоднородностей // Радиотехника и электроника. 1968. Т.13. №12. С. 2115-2121.

67. Шабельников А.В. Рефракция миллиметровых волн в трехмерно-неоднородной атмосфере на приземных и космическихтрассах под малыми углами места. В кн.: Электромагнитные волны в атмосфере и космическом пространстве. М.: Наука, 1986. С. 25-32.

68. Van Vleck I.H. The absorption of microwaves by oxygen // Phys.Rev. 1947. V.71. №7. P. 413.

69. Brown R.R. On the poleward expansion of ionspheric absorption regions triggerd by sudden commencements of geomagnetic storms //J. Geophys. Res. 1978. V.83. ЖАЗ. P. 1169.

70. Franke S.I., Liu C.H. Modelling of equatorial multifreguency scintillation // Radio Science. 1985. V.20. №3. P. 403.

71. Le Vine D.M. Propagation delay in the atmosphere // Radio Science. 1972. V.7. №6. P. 625.

72. Askne I., Nordius H. Estimation of tropospheric delay for microwaves from surface wather data // Radio Science. 1987. V.22. №3. P. 379.

73. Стоцкий А.А. О флуктационных характеристиках тропосферы Земли // Радиофизика. 1973. Т.16. №5. С. 806.

74. Bennett I.A. The calculation of Doppler shifts due to a changing ionosphere // Journal of Atmospheric and Terrestrial Physics. 1967. V.29. №7. P. 887.

75. Михайлов С. Я. Многозначность восстановления профилей плазменной частоты по заданной ВЧХ и их различимость для наклонного распространения КВ-радиоволн в изотропной ионосфере // Изв. Вузов. Радиофизика. 2000. Т.63. №10. С. 855-872.

76. Куницын В.Е., Нестеров И.А. О задаче реконструкции профиля электронной концентрации слоистой плазмы // Фундаментальная и прикладная математика. 1999. Т.5. №2. С. 503-526.

77. Андрианов В.А. Распространение радиоволн в пограничном слое атмосферы // Итоги науки и техники. Сер. Радиотехника. М.:ВИНИТИ, 1994. Т.44. С. 3-80.

78. Яковлев О.И. Космическая радиофизика. М.: Научная книга, 1998.

79. Бин Б.Р., Даттон Е.Д. Радиометеорология / Перевод с анг. под ред. А.А. Семенова. JL: Гидрометеоиздат, 1971.

80. Пахомов JI.A., Пинус Н.Э., Шметер С.М. Аэрологические исследования изменчивости коэффициента преломления атмосферы для ультракоротких радиоволн. М.: Гидрометеоиздат, 1960.

81. Хргиан А.Х. Физика атмосферы. JI.: Гидрометеоиздат, 1978. Т.2.

82. Воронцов П.А. Аэрологические исследования пограничного слоя атмосферы. Гидрометеоиздат, 1960.

83. Вялъцева Э.Е. Изменчивость коэффициента преломления атмосферы для УКВ в пограничном слое // Метеорология и гидрология. 1972. №2. С. 8-14.

84. Типовые характеристики нижнего 300 метрового слоя атмосферы по измерениям на высотной мачте / Под ред. H.JI. Бызовой. М.: Гидрометеоиздат, 1982.

85. Метеорологический режим нижнего 300 метрового слоя атмосферы / Под ред. H.JI. Бызовой. М.: Гидрометеоиздат, 1984.

86. Крашенинников И.А., Лянной Б.Е., Снеговой А.А. Оперативный алгоритм пер-псчета ионограммы наклонного зондирования ионосферы в эффективную ионо-грамму вертикального зондирования //Дифракционнные эффекты коротких радиоволн. М.: ИЗМИРАН, 1982. С. 117-122.

87. Сажип В.И. Моделирование на ЭВМ распространения радиоволн в регулярной ионосфере. Иркутск: ИГУ, 1993.

88. Basler R.P., Scott T.D. Ionospheric structure from odliquebackscatter soundings // Radio Science. 1973. V.8. №5. P. 425-429.

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

90. Андреев М.Ю., Лукичева Т.Н., Мингалев B.C. Модельное исследование влияния главного ионосферного провала на наклонное распространение коротких радиоволн // Геомагнетизм и аэрономия. 2006. Т.46. №1. С. 99-105.

91. Грозов В.П., Куркин В.И., Ким А.Г., Матюшонок С.М., Подделъский И.Н., Розанов С.В., Смирнов В.Ф. Некоторые особенности поведения зимней среднеширот-ной ионосферы по данным НЗ и ВЗ // Солнечно-земная физика. 2004. Вып. 5. С. 120-123.

92. Ратовский К.Г. Исследование характеристик искажений отраженнного сигнала для исследования неоднородностей структуры ионосферы // Геомагнетизм и аэрономия. 2000. Т.40. №5. С. 83-88.

93. Bilitza D. International Reference Ionosphere 2000 // Radio Science. 2001. V.36. №2. P. 261-275.

94. Котович Г.В., Ким А.Г., Михайлов С.Я., Грозов В.П., Михайлов Я.С. Определение критической частоты f0F2 в средней точке трассы по данным наклонного зондирования на основе метода Смита // Геомагнетизм и аэрономия. 2006. Т.46. Ш. С. 547-551.

95. Вертоградов Г.Г., Вертоградов В.Г., Урядов В.П. Мониторинг волновых возмущений методом наклонного зондирования ионосферы // Изв. Вузов. Радиофизика. 2006. Т.49. №12. С. 1015-1029.

96. Рабинер Л., Гоулд Б. Теория и применение цифровой обработки сигналов. М.: Мир, 1978.

97. Сизиков B.C. Математические методы обработки результатов измерений. СПб.: Политехника, 2001.

98. Морозов В.А., Гребенников А.И. Методы решения некорректно поставленных задач. М.: Изд-во МГУ, 1992.

99. Куницын В.Е., Терещенко Е.Д. Томография ионосферы. М.: Наука, 1991.