Космические исследования, 2022, T. 60, № 6, стр. 471-478

Коэффициент поглощения дециметровых радиоволн (~19 см) в ионосфере Земли по результатам решения обратной задачи в радиозатменных спутниковых исследованиях во время магнитной бури в июне 2015 г.

В. Н. Губенко 1*, В. Е. Андреев 1, И. А. Кириллович 1, Т. В. Губенко 1, А. А. Павельев 1, Д. В. Губенко 1

1 Институт радиотехники и электроники им. В.А. Котельникова РАН
Фрязино, Россия

* E-mail: vngubenko@gmail.com

Поступила в редакцию 28.02.2022
После доработки 25.04.2022
Принята к публикации 04.05.2022

Полный текст (PDF)

Аннотация

По результатам анализа радиозатменных данных FORMOSAT-3/COSMIC обнаружено поглощение дециметровых (ДМ) радиоволн (длина волны ~19 см) в нижней высокоширотной ионосфере Земли. Предложен метод восстановления вертикальных профилей коэффициента поглощения путем решения обратной задачи радиозондирования в нижней ионосфере Земли. Этот метод является общим и может быть использован для различных диапазонов радиоволн и других сигналов глобальных навигационных спутниковых систем (ГНСС). Надежно идентифицированы слои поглощения в радиозатменных сеансах измерений, которые обусловлены мощными всплесками рентгеновского излучения и сильными изменениями геомагнитных условий во время бури. Найдено, что на высотах от ~90 до ~100 км величина коэффициента поглощения ДМ-радиоволн достигала значений (5.7 ± 1.4) ⋅ 10–3 дБ/км. Практическая значимость изучения эффектов поглощения радиоволн в D- и Е-областях ионосферы связана с обеспечением бесперебойной работы систем космической радиосвязи и навигации.

1. ВВЕДЕНИЕ

Ранее нами были проведены обработка и анализ ~100 сеансов радиозатменных спутниковых измерений FORMOSAT-3/COSMIC, выполненных во время геомагнитной бури 22–23.VI.2015 в высокоширотной ионосфере Земли [1, 2]. Было показано, что ионосферные возмущения характеристик радиоволн обусловлены как геомагнитными условиями, так и активностью мощных рентгеновских вспышек во время измерений. Для радиозатменных сеансов измерений, проведенных после начала (16.30 UT 22.VI.2015) мощных всплесков рентгеновского излучения и во время (18.36 UT 22.VI–02.00 UT 23.VI) главной фазы магнитной бури, были характерны повышенные уровни электронной плотности в D- и E-областях ионосферы планеты. Нами также был проведен поиск поглощения дециметровых (ДМ) радиоволн (длина волны сигнала ~19 см) на несущей GPS-частоте f1 = 1545.42 МГц. По результатам анализа измерений эйконала и интенсивности ДМ-сигнала спутников FORMOSAT-3/COSMIC, впервые было обнаружено поглощение радиоволн ДМ-диапазона в D- и Е-областях высокоширотной ионосферы Земли [1, 2]. Было найдено, что абсолютная величина интегрального поглощения на трассах радиозатменного зондирования составляет ~3 дБ в интервале высот ~50–90 км, а в отдельных случаях она достигает ~10 дБ на уровнях от ~90 до ~95 км.

В работах [3, 4] представлены выражения для определения рефракционного ослабления (X) мощности радиоволн:

(1)
$\begin{gathered} {\text{1}}~--~X\left( t \right) = m{{a}_{{{\psi }}}} = md{{~}^{{\text{2}}}}{{\psi /}}d{{t}^{{\text{2}}}}, \\ m = {{r}_{{{\psi }}}}{\text{/}}{{\left( {d{{p}_{0}}{\text{/}}dt} \right)}^{{\text{2}}}},~\,\,\,\,{{r}_{{{\psi }}}}~ = {{L}_{L}}{{L}_{{G~}}}{\text{/}}{{L}_{0}}, \\ \end{gathered} $
где aψ – ускорение эйконала (фазового пути ψ), р0 – прицельный параметр радиолуча, LL и LG – расстояния от приемника (LEO) до передатчика (GPS). Существует связь между рефракционным ослаблением X и нормированной мощностью сигнала (P), измеренной на низкоорбитальном спутнике FORMOSAT-3/COSMIC. Для нахождения безразмерной величины P, мощность PL принятого на низкоорбитальном спутнике сигнала нормировалась на среднее значение мощности (P0) радиоволн на высотах более 300 км, то есть P = PL/P0.

Интегральное поглощение (Г) радиозатменного ДМ-сигнала можно определить с помощью следующего соотношения:

(2)
$\begin{gathered} {\text{Г}}\,({\text{отн}}{\text{. ед}}{\text{.}}){\text{ = }}P\,({\text{отн}}{\text{. ед}}{\text{.}}){\text{/}}X\,({\text{отн}}{\text{. ед}}{\text{.}}), \\ {\text{Г}}\,({\text{дБ}}) = P\,({\text{дБ}}) - X\,({\text{дБ}}). \\ \end{gathered} $

Практическая значимость изучения эффектов поглощения радиоволн в D- и Е-областях ионосферы Земли связана с обеспечением бесперебойной работы систем космической радиосвязи и навигации. Возможность определения поглощения в атмосфере и ионосфере Земли ранее обсуждалось в работах [59]. Поглощение является функцией от метеопараметров и частоты и его можно использовать для формулировки обратной задачи восстановления метеопараметров. Измерения поглощения в трех частотных каналах (при достаточно высокой точности измерения амплитуды сигнала) позволяет восстанавливать профили давления, температуры и влажности без априорной информации [10].

Целью этой работы является восстановление высотных профилей коэффициента поглощения Z(h) дециметрового сигнала на высотах ниже ~100 км во время магнитной бури 2015 г. путем решения обратной задачи о поглощении радиоволн.

2. МЕТОД РЕШЕНИЯ ОБРАТНОЙ ЗАДАЧИ О ПОГЛОЩЕНИИ РАДИОВОЛН

Исходные экспериментальные данные о нормированной мощности P(h) радиоволн были сначала отфильтрованы методом скользящего среднего по 50 точкам. Таким образом, были получены усредненные вертикальные профили нормированной мощности Θ(h) дециметрового сигнала. Здесь h = r –R высота анализируемого слоя над поверхностью Земли, R = 6371 км – средний радиус планеты, r расстояние от центра планеты до ионосферного слоя. При вертикальной скорости опускания луча в ионосфере Земли ~2 км/с, это соответствует усреднению данных по вертикальному интервалу в ~2 км. При таком усреднении информация о мелкомасштабной структуре ионосферы практически не теряется, поскольку вертикальное разрешение при анализе методом геометрической оптики (ГО) равно 2(λLL)1/2 = 1.5 км, где λ – длина волны сигнала (~19 см), LL – расстояние от низкоорбитального спутника LEO до атмосферного лимба планеты (~3000 км). В геометрической оптике считается, что два физических луча (лучи вместе с их френелевскими объемами) являются различимыми, если они не пересекаются друг с другом. Анализ профилей рефракционного ослабления X(h) показал, что регулярное (осредненное) рефракционное ослабление 〈X〉 мощности радиоволн в сеансах измерений на высотах >50 км практически отсутствует. Поэтому мы полагаем, что наблюдаемое ослабление средней мощности 〈P(h)〉 сигнала может быть связано с поглощением радиоволн в нижней ионосфере Земли во время магнитной бури [1, 2, 11].

Рассмотрим более подробно, каким образом можно сделать обратное преобразование уравнения (3) для того, чтобы получить вертикальный профиль коэффициента поглощения Z(h):

(3)
$\Theta (h) = \exp \left[ { - \int\limits_s {ZdL} } \right]{\text{.}}$

Интегрирование в выражении (3) производится по траектории s зондирующего ионосферу радиолуча. Описываемая здесь процедура основана на том предположении, что ионосфера Земли является сферически симметричной. Геометрия задачи показана на рис. 1, где n лучей проходит через ионосферу планеты, состоящую из n сферических слоев. Рис. 1 иллюстрирует и объясняет способ решения обратной задачи о поглощении радиоволн, в котором не используется преобразование Абеля. Каждый слой ионосферы имеет постоянный коэффициент поглощения Z и постоянную толщину Δr = Δh = 2 км. Например, если анализируется интервал высот от 110 до 50 км, то число ионосферных слоев равно n = 30. Здесь удобно заменить уравнение (3) системой из n линейных уравнений вида:

(4)
$\begin{gathered} \Theta ({{h}_{i}}) = \exp \left[ { - \int\limits_s {ZdL} } \right] = \exp \left[ { - \sum\limits_{j = 1}^i {{{L}_{{ij}}}{{Z}_{j}}} } \right]{\text{,}} \\ {\text{где }}i{\text{ = 1}} \ldots n, \\ \end{gathered} $
Рис. 1.

Разрез ионосферы Земли, состоящей из n сферических слоев.

где Θ(hi) является величиной Θ для i-го луча, а Zj – коэффициент поглощения мощности радиоволн в j-ом слое. Матричный элемент Lij означает длину той части i-го луча, которая заключена внутри j-го слоя, как показано на рис. 1.

Выразим Θ(hi) в децибелах и развернем указанные выше выражения (4) в виде матричного уравнения:

(5)
$\begin{gathered} {{\Theta }_{{dB}}}({{h}_{i}}) = {\text{10}}{\kern 1pt} {\text{lg}}\left( {\Theta ({{h}_{i}}{\text{)}}} \right) = \left[ {\begin{array}{*{20}{c}} {{{\Theta }_{{dB}}}({{h}_{1}})} \\ {{{\Theta }_{{dB}}}({{h}_{2}})} \\ \vdots \\ {{{\Theta }_{{dB}}}({{h}_{i}})} \\ \vdots \\ {{{\Theta }_{{dB}}}({{h}_{n}})} \end{array}} \right] = - 10{\kern 1pt} \lg (e) \times \\ \times \,\,\left[ {\begin{array}{*{20}{c}} {{{L}_{{11}}}}&{}&{}&{}&{}&{} \\ {{{L}_{{21}}}}&{{{L}_{{22}}}}&{}&{}&{}&{} \\ \ldots & \ldots & \ldots &{}&{}&{} \\ {{{L}_{{i1}}}}&{{{L}_{{i2}}}}& \ldots &{{{L}_{{ii}}}}&{}&{} \\ \ldots & \ldots & \ldots & \ldots & \ldots &{} \\ {{{L}_{{n1}}}}&{{{L}_{{n2}}}}& \ldots &{{{L}_{{ni}}}}& \ldots &{{{L}_{{nn}}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{{Z}_{1}}} \\ {{{Z}_{2}}} \\ \vdots \\ {{{Z}_{i}}} \\ \vdots \\ {{{Z}_{n}}} \end{array}} \right]{\text{,}} \\ \end{gathered} $
где e = 2.718… – основание натуральных логарифмов; 10lg(e) = 4.343…; j = 1…i; i = 1…n; Zi – коэффициент поглощения радиоволн в слое i.

Двумерная матрица Lij связана с радиусом слоев ионосферы Земли как

(6)
${{L}_{{ij}}} = 2\left[ {\sqrt {r_{{j - 1}}^{2} - r_{i}^{2}} - \sqrt {r_{j}^{2} - r_{i}^{2}} } \right] = 2\left[ {\sqrt {{{{(R + {{h}_{{j - 1}}})}}^{2}} - {{{(R + {{h}_{i}})}}^{2}}} - \sqrt {{{{(R + {{h}_{j}})}}^{2}} - {{{(R + {{h}_{i}})}}^{2}}} } \right],$
для i > j и как
(7)
${{L}_{{ii}}} = 2\sqrt {r_{{i - 1}}^{2} - r_{i}^{2}} = 2\sqrt {{{{(R + {{h}_{{i - 1}}})}}^{2}} - {{{(R + {{h}_{i}})}}^{2}}} {\text{,}}$
для i = j.

Решая матричное уравнение (5) относительно Z, получаем:

(8)
$\begin{gathered} {{Z}_{1}} = - \frac{{{{\Theta }_{{dB}}}({{h}_{1}})}}{{4.343{{L}_{{11}}}}}, \\ {{Z}_{2}} = {{\left( { - \frac{{{{\Theta }_{{dB}}}({{h}_{2}})}}{{4.343}} - {{L}_{{21}}}{{Z}_{1}}} \right)} \mathord{\left/ {\vphantom {{\left( { - \frac{{{{\Theta }_{{dB}}}({{h}_{2}})}}{{4.343}} - {{L}_{{21}}}{{Z}_{1}}} \right)} {{{L}_{{22}}}}}} \right. \kern-0em} {{{L}_{{22}}}}}, \\ .{\text{ }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{.}} \\ {{Z}_{i}} = {{\left( { - \frac{{{{\Theta }_{{dB}}}({{h}_{i}})}}{{4.343}} - \sum\limits_{j = 1}^{i - 1} {{{L}_{{ij}}}{{Z}_{z}}} } \right)} \mathord{\left/ {\vphantom {{\left( { - \frac{{{{\Theta }_{{dB}}}({{h}_{i}})}}{{4.343}} - \sum\limits_{j = 1}^{i - 1} {{{L}_{{ij}}}{{Z}_{z}}} } \right)} {{{L}_{{ii}}}}}} \right. \kern-0em} {{{L}_{{ii}}}}}, \\ .{\text{ }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{. }}{\text{.}} \\ \end{gathered} $

Решение обратной задачи дает нам рекуррентные соотношения, которые явным образом связывают величины Z(hi) с ΘdB(hi).

Таким образом, нами определялись вертикальные профили Z(h) коэффициента поглощения радиоволн с масштабом высоты ~2 км и оценивались погрешности их восстановления δZ для каждого анализируемого радиозатменного сеанса измерений. Оценкой погрешности δZ для анализируемого профиля Z(h) служила абсолютная величина его максимальной невязки при решении обратной задачи о поглощении радиоволн. Проверка соответствия восстановленных профилей коэффициента поглощения исходным экспериментальным данным проводилась с помощью решения прямой задачи радиозондирования для слоeв с шириной ~2 км. В этом случае, вычисляемая мощность потока радиоволн PDP(hi) определялась по формуле:

(9)
${{P}_{{DP}}}{\text{(}}{{h}_{i}}{\text{)}} = - \left[ {\sum\limits_{j = 1}^{i - 1} {{{L}_{{ij}}}{{Z}_{j}} + {{L}_{{ii}}}{{Z}_{i}}} } \right] \times 4.343.$

3. АНАЛИЗ ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ И ЗАКЛЮЧЕНИЕ

С помощью решения обратной задачи радиозондирования были найдены вертикальные профили коэффициента поглощения ДМ-радиоволн и определены погрешности их восстановления в нижней ионосфере Земли. Анализ результатов показывает, что геомагнитные условия главной фазы бури и мощные всплески потоков рентгеновского излучения во время измерений сильно влияли на распространение радиоволн в E- и D-слоях ионосферы Земли. Проверка полученных результатов путем решения прямой задачи радиозондирования показала хорошее соответствие между расчетными и исходными экспериментальными данными:

(10)
${{P}_{{DP}}}{\text{(}}{{h}_{i}}{\text{)}} \approx {{\Theta }_{{dB}}}({{h}_{i}}).$

На рис. 2 показаны вертикальные профили нормированной мощности – P(h) (сплошная тонкая с изломами), нормированной мощности после фильтрации – ΘdB(h) (сплошная толстая), нормированной мощности из решения прямой задачи радиозондирования – PDP(h) (сплошная тонкая) и коэффициента поглощения ДМ-радиоволн Z(h), полученные в сеансах радиозатменных измерений 17.22 UT (панели (а) и (б)) и 19.01 UT (панели (в) и (г)) 22 июня 2015 г. в ионосфере Земли. Для каждого анализируемого сеанса указаны широта и долгота зондируемого района, дата и время проведения эксперимента, номера космических аппаратов. Справа указаны погрешности δZ восстановления профилей коэффициента поглощения Z(h). Измерения, относящиеся к профилям на панелях (а) и (б) рис. 2, были проведены после начала мощных всплесков рентгеновских лучей (~16.30 UT 22.VI.2015) и до прихода (18.36 UT 22.VI.2015) основного КВМ (коронального выброса массы) в магнитосферу планеты. Особенностью зондирования этого района с координатами 77.6° N и 66.1° E является то, что средняя мощность сигнала падает до уровня ~0.5 (–3 дБ) на высоте ~90 км и с уменьшением высоты остается на том же уровне (рис. 2а). Локальные максимумы коэффициента поглощения 2.5 ± 0.7 и 2.9 ± ± 0.7 дБ/км наблюдаются на высотах 92.9 и 88.6 км, соответственно (рис. 2б). Эти измерения могут быть подвержены влиянию всплесков рентгеновского излучения, а воздействие на них геомагнитных условий основной фазы шторма отсутствует. Измерения, относящиеся к радиозатменным профилям на панелях (в) и (г) рис. 2, выполненные в 19.01 UT 22.VI.2015, демонстрируют сильные ионосферные возмущения, обусловленные влиянием, как основной фазы магнитной бури, так и мощных потоков рентгеновского излучения. Здесь локальный максимум коэффициента поглощения (3.3 ± 1.3 дБ/км) расположен на высоте 63.8 км (рис. 2г). В табл. 1 представлены высотные интервалы поглощения ДМ-радиоволн в ионосфере Земли, высоты hmax максимумов коэффициента поглощения и значения максимумов Zmax для анализируемых в этой работе сеансов измерений.

Рис. 2.

Вертикальные профили нормированной мощности – P(h) (сплошная тонкая с изломами), нормированной мощности после фильтрации – ΘdB(h) (сплошная толстая), нормированной мощности из решения прямой задачи радиозондирования – PDP(h) (сплошная тонкая) и коэффициента поглощения ДМ-радиоволн Z(h), полученные в сеансах радиозатменных измерений 17.22 UT (панели (а) и (б)) и 19.01 UT (панели (в) и (г)) 22.VI.2015 в ионосфере Земли. Справа указаны погрешности δZ восстановления профилей коэффициента поглощения Z(h).

Таблица 1.  

Высотные интервалы поглощения ДМ-радиоволн в ионосфере Земли, высоты hmax максимумов коэффициента поглощения и значения максимумов Zmax

Дата измерений Время измерений Номер спутника LEO Номер спутника GPS Широта и долгота района измерений
(h ~ 60 км)
Высотный интервал поглощения, км Высота максимума коэф. поглощения hmax, км Значение максимума коэф. поглощения (Zmax ± δZ) · 10–3, дБ/км
22.VI.2015 17.22 UT 6 22 77.6° N
66.1° E
96–72 92.9
88.6
81.4
74.2
2.5 ± 0.7
2.9 ± 0.7
1.4 ± 0.7
1.6 ± 0.7
22.VI.2015 19.01 UT 6 16 67.5° N
40.5° E
70–62 67.8
63.8
1.6 ± 1.3
3.3 ± 1.3
22.VI.2015 19.41 UT 1 30 78.1° N
95.0° E
100–91 95.8
91.8
2.8 ± 1.4
5.7 ± 1.4
22.VI.2015 23.00 UT 5 4 76.4° N
137.5° E
99–43 88.7
84.2
45.4
1.0 ± 0.9
1.0 ± 0.9
1.2 ± 0.9
23.VI.2015 05.06 UT 6 22 70.7° N
118.0° W
104–74 101.7
80.8
3.3 ±1.7
2.0 ± 1.7
23.VI.2015 07.28 UT 1 1 69.2° N
123.3° W
95–79 91.8 2.2 ± 1.1

Наиболее надeжно поглощение сигналов ДМ-диапазона было обнаружено в радиозатменном сеансе измерений, представленном на рис. 3а и 3б. Эти измерения проведены в 19.41 UT 22.VI.2015 и они демонстрируют наиболее сильные возмущения, обусловленные как мощными всплесками рентгеновского излучения, так и изменением геомагнитных условий во время главной фазы шторма. Как видно из представленных данных, при радиозондировании данного района (78.1° N; 95.0° E) полярной шапки Земли (луч опускается сверху вниз), сначала мощность дециметровых радиоволн понижается до уровня ~0.1 (–10 дБ) на высотах от 101.5 до 90.3 км, затем увеличивается до ~0.5 (–3 дБ), и далее продолжает оставаться на том же уровне при уменьшении высоты (рис. 3а). Локальный максимум коэффициента поглощения в этом сеансе достигает наибольшего значения (для всех анализируемых сеансов) 5.7 ± 1.4 дБ/км и расположен на высоте 91.8 км (рис. 3б). Поглощение в радиозатменном сеансе измерений, представленном на рис. 3 (панели (в) и (г)), выражено менее ярко по сравнению с сеансами, рассмотренными выше.

Рис. 3.

Профили нормированной мощности – P(h), нормированной мощности после фильтрации – ΘdB(h), нормированной мощности из решения прямой задачи радиозондирования – PDP(h) и коэффициента поглощения ДМ-радиоволн Z(h), полученные в сеансах радиозатменных измерений 19.41 UT (панели (а) и (б)) и 23.00 UT (панели (в) и (г)) 22.VI.2015 в ионосфере планеты. Здесь обозначения такие же, как и на рис. 2.

На рис. 4 представлены данные о поглощении радиоволн в ионосфере Земли, полученные в сеансах радиозатменных измерений 05.06 UT (панели (а) и (б)) и 07.28 UT (панели (в) и (г)) 23.VI.2015 в ионосфере Земли. Указанные времена измерений соответствуют концу главной фазы бури (~03–05 UT 23.VI.2015) [1]. Как видно из представленных на рис. 4 данных, поглощение радиоволн наблюдается довольно четко в этих сеансах. Локальные максимумы коэффициента поглощения 3.3 ± 1.7 (рис. 4б) и 2.2 ± 1.1 дБ/км (рис. 4г) наблюдаются на высотах 101.7 и 91.8 км, соответственно.

Рис. 4.

Профили нормированной мощности – P(h), нормированной мощности после фильтрации – ΘdB(h), нормированной мощности из решения прямой задачи радиозондирования – PDP(h) и коэффициента поглощения ДМ-радиоволн Z(h), полученные в сеансах радиозатменных измерений 05.06 UT (панели (а) и (б)) и 07.28 UT (панели (в) и (г)) 23.VI.2015 в ионосфере Земли. Обозначения такие же, как на рис. 2 и 3.

Таким образом, нами предложен метод восстановления вертикальных профилей коэффициента поглощения с помощью решения обратной задачи радиозондирования в нижней ионосфере Земли. Этот метод является общим и может быть использован для различных диапазонов радиоволн и других ГНСС-сигналов. Надежно идентифицированы слои поглощения в радиозатменных сеансах измерений во время магнитной бури 22–23.VI.2015. Найдено, что на высотах от ~90 до ~100 км величина коэффициента поглощения ДМ-радиоволн достигала значений (5.7 ± 1.4) · 10–3 дБ/км. Практическая значимость изучения эффектов поглощения радиоволн в D- и Е-областях ионосферы Земли связана с обеспечением бесперебойной работы систем космической радиосвязи и навигации.

Работа выполнена в рамках государственного задания.

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

  1. Губенко В.Н., Андреев В.Е., Кириллович И.А. и др. Радиозатменные исследования возмущений в ионосфере земли во время магнитной бури 22–23.VI.2015 // Геомагнетизм и аэрономия. 2021. Т. 61. № 6. С. 713–722. https://doi.org/10.31857/S0016794021060067

  2. Губенко В.Н., Андреев В.Е., Кириллович И.А. и др. Вариации параметров радиоволн в высокоширотной ионосфере Земли на трассах спутник–спутник во время геомагнитной бури 22–23.VI.2015 // Космич. исслед. 2021. Т. 59. № 3. С. 191–195. https://doi.org/10.31857/S0023420621030055

  3. Gubenko V.N., Pavelyev A.G., Kirillovich I.A., Liou Y.-A. Case study of inclined sporadic E layers in the Earth’s ionosphere observed by CHAMP/GPS radio occultations: Coupling between the tilted plasma layers and internal waves // Advances in Space Research. 2018. V. 61. № 7. P. 1702–1716. https://doi.org/10.1016/j.asr.2017.10.001

  4. Губенко В.Н., Кириллович И.А. Модуляция спорадических Е-слоев мелкомасштабными атмосферными волнами в высокоширотной ионосфере Земли // Солнечно-земная физика. Т. 5. № 3. С. 116–129. https://doi.org/10.12737/szf-53201912

  5. Liou Y.-A., Pavelyev A.G., Liu S.-F. et al. FORMOSAT-3/COSMIC GPS radio occultation mission: preliminary results // IEEE T. Geosci. Remote. 2007. V. 45. P. 3813–3826.

  6. Pavelyev A.G., Liou Y.A., Wickert J. et al. New Applications and Advances of the GPS Radio Occultation Technology as Recovered by Analysis of the FORMOSAT-3/COSMIC and CHAMP Data-Base / Eds. Steiner A., Pirscher B., Foelsche U., Kirchengast G. Springer-Verlag, Berlin. 2009. P. 165–178. https://doi.org/10.1007/978-3-642-00321_9

  7. Pavelyev A.G., Liou Y.-A., Zhang K. et al. Identification and localization of layers in the ionosphere using the eikonal and amplitude of radio occultation signals // Atmos. Meas. Tech. 2012. V. 5. P. 1–16. https://doi.org/10.5194/amt-5-1-2012

  8. Pavelyev A.G., Liou Y.-A., Matyugov S.S. et al. Application of the locality principle to radio occultation studies of the Earth’s atmosphere and ionosphere // Atmos. Meas. Tech. 2015. V. 8. P. 2885–2899. https://doi.org/10.5194/amt-8-2885-2015

  9. Wickert J., Pavelyev A.G., Liou Y.A. et al. Amplitude scintillations in GPS signals as a possible indicator of ionospheric structures // Geophys. Res. Lett. 2004. V. 31. L24801. https://doi.org/10.1029/2004GL020607

  10. Горбунов М.Е. Физические и математические принципы спутникового радиозатменного зондирования атмосферы Земли. М.: ГЕОС, 2019.

  11. Andreev V.E., Gubenko V.N., Pavelyev A.A. et al. Absorption of decimetre radio waves in the Earth’s high-latitude ionosphere during a geomagnetic storm in June 2015 // J. Physics: Conference Series. 2020. V. 1632. P. 012008. https://doi.org/10.1088/1742-6596/1632/1/012008