КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ

Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.

 

ОБЩИЕ СВЕДЕНИЯ


Номер 17-77-20009

НазваниеФормирование, диагностика и мониторинг плазменных структур различных масштабов в Арктическом и Антарктическом регионах ионосферы Земли, их влияние на качество радиосвязи, радиолокации и работу навигационных спутниковых систем

РуководительКлименко Максим Владимирович, Доктор физико-математических наук

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

Период выполнения при поддержке РНФ 07.2017 - 06.2020 

Конкурс№24 - Конкурс 2017 года по мероприятию «Проведение исследований научными группами под руководством молодых ученых» Президентской программы исследовательских проектов, реализуемых ведущими учеными, в том числе молодыми учеными.

Область знания, основной код классификатора 07 - Науки о Земле, 07-608 - Ионосфера

Ключевые словаполярная ионосфера, модель ионосферы, распространение радиоволн, радиолокация, радиосвязь, навигация, плазменные неоднородности и неустойчивости, космическая погода Арктики и Антарктики, ГЛОНАСС, навигационные спутниковые системы

Код ГРНТИ37.15.29


СтатусУспешно завершен


 

ИНФОРМАЦИЯ ИЗ ЗАЯВКИ


Аннотация
Изучение динамики полярной структуры ионосферы является актуальным разделом геофизики, имеющим большое прикладное значение, особенно для России. В отличие от экваториальной и среднеширотной ионосферы, полярная область характеризуется сравнительно низкими значениями электронной концентрации, являясь при этом сильно структурированной средой из-за наличия большого количества неоднородностей различных масштабов, от метров до сотен километров. Неоднородности возникают вследствие процессов магнитосферно-ионосферного взаимодействия, связанных с авроральной активностью, в частности, с высыпаниями энергичных частиц, а также процессов переноса и разогрева плазмы за счет магнитосферных электрических полей и продольных токов. При таком механизме формирования ионосферы высоких широт, неоднородности потоков частиц передаются из магнитосферы в ионосферу, что создает большие помехи для радиосвязи в ионосфере Земли за счет интерференции радиоволн на плазменных неоднородностях. Кроме того, понимание особенностей полярной ионосферы необходимо для корректного использования спутниковых навигационных систем в высоких широтах. Основные ошибки определения местоположения и времени с помощью систем GPS (Global Positioning System) и ГЛОНАСС (Глобaльная навигационная спутниковая система) связаны с прохождением сигналов через ионосферу. Существующие же в настоящее время эмпирические модели ионосферы, чаще всего используемые в качестве моделей среды для задач распространения радиоволн, недостаточно точно описывают параметры среды, в том числе в периоды геомагнитных возмущений. Это относится и к наиболее часто используемой международной справочной модели IRI, особенно в высокоширотной области. Увеличение темпов развития программ по освоению Арктики и наметившийся рост в использовании высоких широт для перевозки пассажиров воздушным транспортом делают задачу обеспечения надежности радиосвязи одной из наиболее важных. Таким образом, создание глобальной эмпирической модели, включающей высокоширотную ионосферу, является особенно актуальным и приоритетным направлением в исследовании ионосферы Земли. Основные цели данного проекта, позволяющие получить результаты мирового уровня, следующие: 1. Построить глобальную эмпирическую модель параметров максимума F области ионосферы и полного электронного содержания более точную, чем существующие мировые и российские аналоги; 2. Используя улучшенную модель ионосферы и набор моделей расчета лучевых траекторий и радиотрасс КВ-диапазона, которые будут доработаны в ходе выполнения данного проекта, создать на их основе комплекс программ прогноза применимых частот, который может быть использован для уточнения параметров работы передающих станций, систем радиолокации и противовоздушной обороны; 3. Провести на основе созданного комплекса программ фундаментальные исследования в области физики ионосферы (по глобальной структуре и долготному эффекту в распределении электронной концентрации, влиянию солнечной и геомагнитной активности в различные сезоны в выделенных широтных и долготных секторах, и т.д.) и распространении радиоволн (исследование функционала оптической длины пути, влияние магнитного поля на распространение КВ радиоволн, и т.д.); 4. Развить и апробировать методику адаптации параметров модели высокоширотной ионосферы по данным наклонного полного электронного содержания (ПЭС), предложенную и используемую на практике для среднеширотной ионосферы. Также в ходе выполнения проекта планируется осуществить детальное исследование формирования и динамики плазменных неоднородностей различных масштабов в полярной и приполярной ионосфере Земли, поскольку плазменные неоднородности могут порождать новые и изменять существующие ионосферно-магнитосферные связи и оказывать существенное влияние на качество радиосвязи, радиолокации и работу навигационных спутниковых систем. В рамках проекта планируется исследовать влияние высокоширотных неоднородностей различных масштабов на качество работы радиолокационных систем. Научная новизна данного проекта состоит: (1) в создании более точной, чем имеющиеся в настоящее время, глобальной эмпирической модели максимума слоя F2 на основе которой могут быть получены новые знания о широтной и долготной структуре полярной ионосферы; (2) в разработке новых методов расчетов радиотрасс в трехмерно-неоднородной ионосфере; (3) в отказе от упрощений, связанных с предположением об однородности, бесконечности и стационарности плазмы, при изучении неоднородностей в ионосферной плазме.

Ожидаемые результаты
Одним из основных результатов данного проекта станет создание более точной (особенно в высокоширотном регионе), чем имеющиеся в настоящее время, глобальной эмпирической модели максимума слоя F2, на основе которой и модели IRI будет создан комплекс программ прогноза оптимально применимых частот в высокоширотном Арктическом и Антарктическом регионе. На основе созданной модели могут быть получены новые знания о широтной и долготной структуре полярной ионосферы. Впервые в рамках одного проекта будут реализованы и применены две методики расчета радиотрасс в выбранной модели ионосферы: решение уравнения эйконала методом характеристик с использованием метода стрельбы и прямой вариационный метод. Создание комплекса программ, включающих модели среды и распространения радиоволн позволяет проводить оперативную оценку параметров сигнала для осуществления эффективной КВ, УКВ и спутниковой связи в различных регионах Земли. Также будет разработан пакет программ и сервисов, позволяющих использовать модель ионосферы и программы расчета радиотрасс в интерактивном режиме. Постановка задачи разрабатываемого метода расчета радиотрасс на основе вариационного принципа универсальна и в перспективе применима к широкому спектру актуальных научных задач (распространение сейсмических волн в земной коре, ультразвуковая томография и т.д.). Поэтому решение проблемы нахождения всех видов решений данным способом обуславливает важное фундаментальное значение решаемой задачи. Результаты проекта могут быть использованы для решения практических задач, в том числе, радиосвязи, радиолокации и спутниковой навигации в Арктике и Антарктике. Знание характеристик плазменных неоднородностей позволит учесть в аппаратуре амплитудные и фазовые искажения, возникающие при распространении сигнала в пространственно неоднородной и нестационарной среде. В результате выполнения проекта будет развит метод диагностики высокоширотной ионосферы в Арктическом регионе, который будет заключаться в адаптации параметров эмпирической модели ионосферы по данным наклонного ПЭС (полное электронное содержание) с учетом азимутального разделения коррекции по наблюдаемым спутникам. До настоящего времени подобные исследования проводились только в среднеширотном регионе и без учета азимутального разделения.


 

ОТЧЁТНЫЕ МАТЕРИАЛЫ


Аннотация результатов, полученных в 2017 году
В ходе выполнения проекта показано, что международная справочная модель ионосферы IRI, наиболее часто используемая для решения научных и прикладных задач радиосвязи, навигации и радиолокации, некорректно описывает основные морфологические особенности распределения ионосферных параметров в области высоких и субавроральных широт, то есть в Арктическом и Антарктическом регионах. Нами создан полный банк ионосферных спутниковых и наземных данных. На языке программирования Python была написана программа, позволяющая группировать, сортировать и анализировать данные ионосферных наблюдений. Реализовано автоматическое пополнение массива данных новыми данными. Результатом работы программы является матрица, содержащая географические координаты, средние арифметические и медианные значения ионосферных параметров, количество вертикальных профилей электронной концентрации и количество отбракованных данных по каждому из выбранных критериев. Различные спутниковые данные были тщательно проанализированы, очищены от некорректных значений, состыкованы и теперь могут использоваться для построения глобальной эмпирической модели параметров F2 области ионосферы. На создаваемую нами модель накладываются требования воспроизведения ею пространственной структуры и крупномасштабной динамики высокоширотной ионосферы. Проведены теоретические исследования механизмов формирования высокоширотных крупномасштабных неоднородностей на основе результатов расчетов Глобальной Самосогласованной Модели Термосферы, Ионосферы, Протоносферы (ГСМ ТИП), которые позволили лучше понять особенности высокоширотной ионосферы. Проведена количественная оценка эффективности методики адаптации справочных моделей ионосферы по данным наклонного полного электронного содержания для высокоширотного региона путем сравнения модельных и экспериментальных результатов. Результаты апробации показали, что для марта уменьшение ошибки определения ключевого ионосферного параметра - критической частоты F2 слоя, foF2, при коррекции моделей наблюдается как в ночное, так и в дневное время. В остальные месяцы после коррекции модели среднеквадратическая ошибка определения foF2 увеличивается по сравнению с результатами расчетов моделей без коррекции. Отмечается улучшение модельных результатов при коррекции в освещенное время (за исключением декабря месяца). Выяснение причин улучшений и ухудшений результатов модельного описания ионосферы после коррекции позволит определить вектор дальнейших действий для решения прикладной задачи коррекции модели среды в режиме реального времени. Решение данной задачи важно для корректной работы радиолокационных и навигационных радиофизических систем. В ходе выполнения проекта параллельно развиваются две модели распространения радиоволн, которые в конечном итоге будут дополнять друг друга при решении различных научных и прикладных задач. Одна из моделей распространения радиоволн, построенная на основе решения уравнения эйконала, за отчетный период была дополнена программной реализацией трехмерного алгоритма “пристрелки”. Это позволило нам исследовать различные характеристики радиотрасс в трехмерной слабо неоднородной анизотропной среде как в спокойных условиях, так и в периоды геомагнитных бурь. Это позволило изучить влияние магнитного поля Земли на формирование лучевых траекторий и характеристики радиотрасс как в спокойных условиях, так и во время геомагнитной бури. В ходе проведения вычислительных экспериментов получено, что учет магнитного поля Земли при расчете радиотрасс приводит: 1) к сильной азимутальной расходимости лучей; 2) к значительным изменениям максимально-применимых частот для E и F областей ионосферы (до 1-2 МГц). Оба этих результата говорят о важности учета анизотропии среды в задачах расчета радиотрасс. Другим перспективным методом расчета радиотрасс между двумя объектами, положение которых зафиксировано, является использование прямого вариационного метода. В ходе исследований в рамках данного проекта была доказана возможность раздельного поиска двух типов радиолучей (верхних и нижних лучей) для различных рабочих частот на основе вариационного метода. Реализован комбинированный подход глобальной оптимизации, основанный на принципе Ферма и последовательном поиске верхних и нижних лучей, который обладает высокой степенью сходимости и оптимальным временем проведения расчетов, а также возможностью определения полного набора решений. Важными преимуществами разработанного метода являются: 1) решение проблемы задания начального приближения, задаваемого в виде линии, проведенной на поверхности Земли, соединяющей пункты передачи и приема для последовательного нахождения всех односкачковых трасс; 2) осуществление процедуры глобального поиска всех радиолучей между пунктами передачи и приема; 3) точное попадание всех найденных радиолучей из точки передачи в точку приема. На основе двух описанных выше моделей распространения радиоволн создано программное обеспечение для построения и визуализации ионограмм наклонного зондирования на языке программирования C++ и среды Gnuplot. Предусмотрена возможность одновременного построения модельных результатов и данных наблюдений ионограмм наклонного зондирования, что позволяет осуществлять сопоставление и анализ полученных результатов. В рамках выполнения проекта произведена обработка и анализ различных спутниковых и радарных измерений параметров плазмы на различных высотах ионосферы и плазмосферы, а также характеристик аврорального овала, полярного ветра и языка ионизации, оказывающих значительное влияние на взаимосвязь процессов между ионосферой и магнитосферой и на прикладные аспекты, связанные с распространением радиоволн в высокоширотном регионе, как в спокойных условиях, так и особенно в периоды геомагнитных бурь. Впервые GPS-наблюдения на борту метеорологической миссии MetOp были использованы для выявления неоднородностей плазмы на высотах выше 835 км. Были проведены оценки проявления ионосферных неоднородностей различных масштабов в различных высотных интервалах, в том числе на высотах плазмосферы. Наши результаты показывают, что при сильных геомагнитных бурях на высотах внешней ионосферы могут наблюдаться значительные неоднородности плотности плазмы, эффекты которых могут проявляться и выше. Быстрые изменения фазовой задержки сигнала во внешней ионосфере могут значительно ухудшать эффективность работы GPS на низкоорбитальных спутниках, когда бортовые GPS приемники могут временно испытывать полную/частичную потерю сигналов, проходящих через значительные градиенты плотности плазмы. Мульти-инструментальный анализ наземных и спутниковых данных выявил, что ионосферные неоднородности, обнаруженные на высотах внешней ионосферы совпадают с положением крупномасштабных ионосферных неоднородностей. На основании результатов расчетов модели ГСМ ТИП произведена интерпретация формирования крупномасштабных ионосферных неоднородностей в высоких широтах на различных фазах геомагнитной бури. Сформулирована теория, а также были созданы и протестированы вычислительные алгоритмы для решения задачи о влиянии плазменных неоднородностей на формирования широкополосной электростатической турбулентности (шума) в полярных областях Земли. Проведено численное моделирование для изучения роли различных ионосферных источников энергии в генерации широкополосной низкочастотной турбулентности. Неоднородные электрические поля и плазменные неоднородности могут приводить к возникновению регулярно наблюдаемого широкополосного шума в верхней ионосфере, а также к нагреву ионосферных ионов и их последующему оттоку в магнитосферу. Для задания неоднородных профилей электрического поля, концентрации плазмы и продольных скоростей заряженных частиц в численных расчетах использовались спутниковые данные и впервые в качестве исходных данных использовались ракетные измерения, полученные при изучении ионосферы высоких широт и обратных течений в каспе. Показано, что электростатическая широкополосная турбулентность может быть идентифицирована как разновидность электростатических ионно-циклотронных или ионно-акустических волн, возбуждаемых неоднородным распределением плотности энергии волны.

 

Публикации

1. Захаренкова И., Черняк Ю. Underutilized spaceborne GPS observations for space weather monitoring Space Weather, Space Weather, 16. https://doi.org/10.1002/2017SW001756 (год публикации - 2018) https://doi.org/10.1002/2017SW001756

2. Ильясов А. А., Чернышов А. А., Могилевский М. М., Головчанская И. В., Козелов Б. В. Неустойчивость, вызванная неоднородным распределением плотности энергии, как возможный источник электростатического широкополосного шума Химическая физика, том 37, № 5 (год публикации - 2018)

3. Клименко М.В., Клименко В.В., Дэспирак И.В., Захаренкова И.Е., Козелов Б.В., Черняков С.М., Андреева Е.С.,Терещенко Е.Д., Веснин А.М., Коренькова Н.А., Гомонов А.Д., Васильев Е.Б., Ратовский К.Г. Disturbances of the thermosphere-ionosphere-plasmasphere system and auroral electrojet at 30°E longitude during the St. Patrick's Day geomagnetic storm on 17–23 March 2015 Journal of Atmospheric and Solar-Terrestrial Physics, - (год публикации - 2018) https://doi.org/10.1016/j.jastp.2017.12.017

4. Котова Д.С., Оводенко В.Б., Ясюкевич Ю.В., Клименко М.В., Мыльникова А.А., Козловский А.Е. Коррекция модели NeQuick в высокоширотном регионе с использованием данных наклонного полного электронного содержания Physics of Auroral Phenomena, Т. 40, № 1 (41), С. 120-124 (год публикации - 2017)

5. Носиков И.А., Клименко М.В., Бессараб П.Ф., Жбанков Г.А. Применение прямого вариационного метода поиска верхних и нижних лучей в задаче расчета КВ-радиотрасс в ионосфере Сборник трудов БШФФ-2017 и XV Конференции молодых ученых, Иркутск. 2017, С. 129 – 131 (год публикации - 2017)

6. Чернышов A.A., Спичер A., Ильясов A.A., Милош В., Клаусен Л., Саито Ю., Джин Я., Моен Й. Studies of small-scale plasma inhomogeneities in the cusp ionosphere using sounding rocket data Physics of Plasmas, Volume 25, Issue 4 (год публикации - 2018) https://doi.org/10.1063/1.5026281

7. Чирик Н.В., Клименко М.В., Клименко В.В., Карпачев А.Т., Ратовский К.Г., Коренькова Н.А. Принципы обработки и отбора данных радиозатменных наблюдений для исследования F2-слоя ионосферы Химическая физика, Т. 36. № 12, С. 66–74. (год публикации - 2017) https://doi.org/10.7868/S0207401X17120081

8. Чугунин Д.В., Клименко М.В., Клименко В.В. Характеристики потоков полярного ветра на высотах ~20000 км Химическая физика, том 37, № 5 (год публикации - 2018)

9. - Ученые активно ищут способ нейтрализовать магнитные бури газета «Страна Калининград», - (год публикации - )


Аннотация результатов, полученных в 2018 году
На основе банка данных прямых спутниковых измерений, внешнего зондирования ионосферы и радиозатменных наблюдений, собранных в 2017 году, созданы две эмпирические модели параметров максимума электронной концентрации в F2 области ионосферы для всех моментов местного времени с часовым разрешением и всех уровней солнечной активности: 1) глобальная модель для всех сезонов (по месяцам); 2) Модель Главного Ионосферного Провала (МГИП) с разбиением на четыре месяца зимнего солнцестояния в северном и южном полушариях. На основе результатов эмпирических моделей выделены долготные вариации положения главного ионосферного провала. Показано, что амплитуда долготных вариаций положения ГИП в дневное время больше чем ночью. На основе статистической обработки данных четко выделены ГИП и высокоширотный провал при высокой солнечной активности в утренние часы. Модельные исследования показали, что долготные вариации положения ГИП в дневное время определяются долготными изменениями зенитного угла Солнца и долготного распределения плотности атомарного кислорода. Формирование долготных изменений в ночном ГИП связано с долготным эффектом ионизации высокоэнергичных высыпающихся частиц, состава нейтральной атмосферы и высокоширотного электрического поля. На основе модельных расчетов, спутниковых и наземных данных выявлена задержка вариаций дневных значений foF2 и ПЭС относительно изменений F10.7. Эта задержка для высоких широт больше для foF2 по сравнению с ПЭС. Согласно результатам моделирования, изменение отношения O / N2 представляется основным механизмом формирования этой задержки. Проведено исследование причин ухудшения модельного описания foF2 после коррекции эмпирических моделей ионосферы IRI-plas и NeQuick на основе данных наклонного ПЭС при улучшении модельного описания вертикального ПЭС в некотором регионе. Некорректное воспроизведение высотной структуры ионосферы является одной из основных причин неэффективности работы используемого метода коррекции ионосферы для описания foF2. Предварительный сравнительный анализ модельных и экспериментальных карт ПЭС показал, что выбранные модели ионосферы наилучшим образом воспроизводят широтно-долготную структуру ПЭС в марте и хуже всего в солнцестояние. Это и является причиной выявленной нами сезонной зависимости эффективности коррекции foF2 по данным наклонного ПЭС. В рамках изучения источников, свойств и физических механизмов широкополосной турбулентности в полярных областях ионосферы показано, что плазменная неустойчивость, вызванная неоднородным распределением плотности энергии может успешно объяснить результаты ракетных измерений. Показана принципиальная возможность применимости метода локальной перемежаемости для изучения нерегулярных турбулентных процессов и нерегулярных структур электрического поля в полярных областях ионосферы. В нашем исследовании наблюдается качественное согласие между высокими значениями инкремента плазменной неустойчивости, вызванной неоднородным распределением плотности энергии и областей с высоким значением параметра перемежаемости в полярных областях. Это говорит о том, что развитие данного рода плазменных неустойчивостей связано с турбулентными негауссовыми (перемежаемыми) процессами. Кроме того, показано, что развитие данной неустойчивости существенно зависит не только от наличия неоднородностей в электрическом поле, но также сильно зависит от продольных скоростей заряженных частиц в высокоширотной области ионосферы. Проанализированы случаи, когда на спутнике Интербол-2 наблюдалось пересечение полярной границы аврорального овала (Альфвеновская аврора) и когда были одновременные изображения авроральных свечений на спутнике POLAR. На основе выбранных случаев были получены зависимости скорости поперечного нагрева ионов от скорости движения полярной границы аврорального овала к полюсу. Показано, что чем выше скорость границы аврорального овала к полюсу, тем наблюдается более интенсивный поперечный нагрев окружающей фоновой плазмы Выявлено пространственное и временное совпадение между появлением интенсивной радиоавроры по данным УКВ радара в районе Красноярска и увеличением индекса ROT по данным приемника ГНСС сигналов на авроральной станции Норильск во время геомагнитной бури 27-28 мая 2017 года. Это указывают на наличие связи между неоднородностями с размерами ~10-100 км и неоднородности в Е слое с пространственными размерами ~35 см. В областях формирования ионосферных неоднородностей наблюдаются сильные электрические поля и увеличенные дрейфы ионов. Показано, что время возникновения ионосферных неоднородностей согласуется с пространственно-временными изменениями продольных токов второй зоны. На основании анализа результатов модели ГСМ ТИП и данных наблюдений показано, что во время геомагнитной бури St. Patrick 2015 года в приполюсных областях, где скорости дрейфа плазмы максимальны, язык ионизации значительно подавляется за счет повышенной плотности молекулярного азота. На высотах плазмосферы область, подобная языку ионизации, занимает всю полярную шапку. Это можно объяснить меньшим влиянием рекомбинационных процессов на распределение электронной плотности с ростом высоты. Нами проанализированы и выделены основные особенности и области применения разрабатываемых методов расчета радиотрасс. Вариационный метод расчета лучевых траекторий радиоволн имеет следующие преимущества и особенности: а) Положения передатчика и приемника зафиксированы в соответствии с граничными условиями, что позволяет избежать процедуры пристрелки и сопутствующие ей проблемы сходимости, возникающие особенно остро при нахождении верхних и волноводных решений граничной задачи в трехмерно-неоднородной среде; б) Имеются перспективы решения проблемы нахождения всех решений при многолучевости на основе созданного нами алгоритма глобальной оптимизации. Для определения стационарных решений любого типа (минимумы, седловые точки и максимумы функционала оптической длины радиолуча) нами разработан универсальный метод оптимизации функционала оптической длины пути, названный методом обобщенной силы. На основе чередования различных вариантов метода обобщенной силы разработана, реализована и апробирована методика глобальной оптимизации для определения полного числа решений граничной задачи. С помощью этого алгоритма мы избавляемся от необходимости задания шага трассировки лучей от которого зависит число найденных решений с помощью трассировки и пристрелки. Применение решения уравнения эйконала с начальными условиями является оптимальным при нахождения поля волны от антенны передатчика с заданной диаграммой направленности. Нами было показано, что учет анизотропии среды приводит к существенной расходимости лучевых траекторий, и, как следствие, к усложнению процедуры трехмерной пристрелки. Показано, что длина оптического пути верхнего луча в случае распространения в анизотропной среде в целом больше, чем в изотропном случае. Затухание сигнала в изотропной среде больше, чем при учете анизотропии. С помощью сравнения результатов расчетов лучевых траекторий, полученных вариационным методом и методом на основе решения уравнения эйконала, установлена взаимосвязь между типом стационарных точек функционала и каустиками. Радиолучи, не образующие каустических особенностей, являются минимумами функционала оптической длины пути (верхние лучи). Лучи, касающиеся каустику 1 раз, соответствуют седловой точке 1 рода (нижние лучи). Аналогично, лучи, образующие каустики и касающиеся ее N раз, соответствуют седловой точке N-го рода (пример: лучи в межслоевых каналах).

 

Публикации

1. Chirik N.V., Klimenko M.V., Karpachev A.T., Klimenko V.V., Pustovalova L.V. Radio Occultation Data for High-Latitudinal F2-Layer Empirical Modeling Proceedings of VI International conference “Atmosphere, Ionosphere, Safety”, - (год публикации - 2018)

2. Карпачев А.Т., Клименко М.В., Клименко В.В. Longitudinal variations of the ionospheric trough position Advances in Space Research, V. 63 (2), P. 950-966 (год публикации - 2019) https://doi.org/10.1016/j.asr.2018.09.038

3. Карпачев А.Т., Чирик Н.В., Клименко В.В., Клименко М.В. Эмпирические модели электронной концентрации в максимуме высокоширотного F2 слоя ионосферы на основе спутниковых измерений для различных уровней солнечной активности Сборник конференции “Астрономия-2018” Том 2. Солнечно-земная физика – современное состояние и перспективы, C. 98-101 (год публикации - 2018) https://doi.org/10.31361/eaas.2018-2.025

4. Клименко М., Чирик Н., Котова Д., Носиков И., Жбанков Г., Ратовский К., Карпачев А., Пустовалова Л., Бессараб Ф., Бессараб П., Клименко В., Захаренкова И., Благовещенский Д. Development of improved ionospheric empirical model and software for HF ray tracing 2018 2nd URSI Atlantic Radio Science Meeting (AT-RASC), S-G02-09 (год публикации - 2018) https://doi.org/10.23919/URSI-AT-RASC.2018.8471348

5. Котова Д., Оводенко В., Ясюкевич Ю., Мыльникова А., Клименко М. Ground-Based GNSS Data for the Ionosphere Model Correction at high-latitudes 2018 2nd URSI Atlantic Radio Science Meeting (AT-RASC), S-G02-12 (год публикации - 2018) https://doi.org/10.23919/URSI-AT-RASC.2018.8471472

6. Котова Д.С. Алгоритм трехмерной «пристрелки» на основе численной модели распространения коротких радиоволн в ионосфере Proceedings of the 41st Annual Seminar “Physics of Auroral Phenomena”, C. 129-130 (год публикации - 2018) https://doi.org/10.25702/KSC.2588-0039.2018.41.129-130

7. Котова Д.С., Клименко М.В., Клименко В.В., Благовещенский Д.В., Захаров В.Е. Анализ высокоширотных ионограмм наклонного зондирования в период геомагнитной бури 17 марта 2015 г. Proceedings of the 41st Annual Seminar “Physics of Auroral Phenomena”, C. 131-134 (год публикации - 2018) https://doi.org/10.25702/KSC.2588-0039.2018.41.131-134

8. Котова Д.С., Клименко М.В., Клименко В.В., Благовещенский Д.В., Захаров В.Е., Веснин А.М. Model Simulation of Oblique Sounding Ionograms between Lovozero and Gorkovskaya during St. Patrick's Day 2015 Geomagnetic Storm Proceedings of VI international conference Atmosphere, Ionosphere and Safety, Part 2, P. 44-49 (год публикации - 2018)

9. Котова Д.С., Оводенко В.Б. Метод коррекции параметров эмпирических моделей ионосферы по данным наклонного полного электронного содержания Сборник конференции “Астрономия-2018” Том 2. Солнечно-земная физика – современное состояние и перспективы, С. 114-117 (год публикации - 2018) https://doi.org/10.31361/eaas.2018-2.029

10. Котова Д.С., Оводенко В.Б., Ясюкевич Ю.В., Клименко М.В., Мыльникова А.А., Козловский А.Е., Гусаков А.А. Коррекция эмпирических ионосферных моделей IRI-Plas и NeQuick в высокоширотном регионе с использованием данных отдаленных приемников сигналов глобальных навигационных спутниковых систем Химическая физика, Т. 37, № 7, с. 87–92 (год публикации - 2018) https://doi.org/10.1134/S0207401X18070129

11. Носиков И.А., Клименко М.В., Бессараб П.Ф. Identification of Low and High Ionospheric Rays by a Direct Variational Method Radio Science Conference (URSI AT-RASC), 2018 2nd URSI Atlantic. Gran Canaria, Spain, S-B03-P (год публикации - 2018) https://doi.org/10.23919/URSI-AT-RASC.2018.8471460

12. Носиков И.А., Клименко М.В., Бессараб П.Ф., Жбанков Г.А. Features of Radio Ray Identification by a Direct Variational Method Proceedings of VI International conference “Atmosphere, Ionosphere, Safety”, Part 2, P. 49–53 (год публикации - 2018)

13. Оводенко В., Захаренкова И., Клименко М., Тютин И., Успенский М., Котова Д., Ратовский К., Чирик Н., Клименко В., Рахматулин Р., Пашинин А., Дмитриев А., Суворова А. Ionospheric Irregularities over Norilsk during the 27-28 May 2017 Geomagnetic Storm 2018 2nd URSI Atlantic Radio Science Meeting (AT-RASC), G05-5 (год публикации - 2018) https://doi.org/10.23919/URSI-AT-RASC.2018.8471415

14. Сомина Е.Р., Носиков И.А., Клименко М.В., Бессараб П.Ф., Жбанков Г.А., Иванова В.А., Подлесный А.В. Oblique Sounding Ionogram Simulation for HF Ray Traces over Siberian and Far Eastern regions Proceedings of VI International conference “Atmosphere, Ionosphere, Safety”, Part 2, P. 133–137 (год публикации - 2018)

15. Чернышов А., Спичер А., Ильясов А., Милош В. Small-scale plasma inhomogeneities in the Earth's ionosphere of the Arctic region Proceedings of VI International conference “Atmosphere, Ionosphere, Safety”, Part 2, P. 62-67. (год публикации - 2018)

16. Чернышов А., Спичер А., Ильясов А., Милош В., Могилевский М. Role of plasma inhomogeneities in the generation of broadband waves in the polar ionosphere 2nd URSI AT-RASC, Gran Canaria, H02-5 (год публикации - 2018) https://doi.org/10.23919/URSI-AT-RASC.2018.8471639

17. Черняк Ю.В., Захаренкова И.Е. Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations Advances in Space Research, Volume 63, Issue 6, Pages 1845-1859 (год публикации - 2018) https://doi.org/10.1016/j.asr.2018.10.036

18. Чирик Н.В., Клименко М.В., Карпачев А.Т., Ратовский К.Г., Клименко В.В., Лещенко В.С., Коренькова Н.А. Оптимальный индекс солнечной активности связанный с F10.7, для эмпирической модели F2-слоя ионосферы Химическая физика, T. 37, № 7, C. 93–96 (год публикации - 2018) https://doi.org/10.1134/S0207401X18070051

19. Чугунин Д.В., Могилевский М.М. Ion heating on the poleward moving of the auroral oval polar boundary Proceedings of VI International conference “Atmosphere, Ionosphere, Safety”, Part 2, P. 92-96. (год публикации - 2018)

20. - Ученые активно ищут способ нейтрализовать магнитные бури газета «Страна Калининград», - (год публикации - )


Аннотация результатов, полученных в 2019 году
На основе тестирования созданных в ходе выполнения проекта эмпирических моделей высокоширотной ионосферы показано, что созданные модели воспроизводят основную структурную особенность суб-авроральной ионосферы (главный ионосферный провал) и его долготную вариацию схожим образом с канадской моделью E-CHAIM в отличие от справочной модели IRI. Не очень хорошее воспроизведение структуры главного ионосферного провала моделью IRI приводит к тому, что широтно-долготное распределение полного электронного содержания хуже всего воспроизводится этой моделью в зимнее солнцестояние. Для этого периода согласно результатам наших статистических исследований коррекция эмпирических моделей по данным ГНСС приемников не улучшает, а ухудшает воспроизведение foF2 и полного электронного содержания (ПЭС) в регионе, тем самым вводя ограничения на использование данного алгоритма коррекции IRI-подобных моделей при решении практических задач радиолокации и КВ-радиосвязи в суб-авроральном и высокоширотном регионе. Была выявлена сезонная зависимость эффективности коррекции моделей IRI-plas и NeQuick по данным приемников ГНСС сигналов в средних и субавроральных широтах. Наилучшая эффективность коррекции foF2 и полного электронного содержания получена в период весеннего равноденствия, а наихудшая - в зимнее солнцестояние. Наиболее вероятным объяснением сезонной зависимости эффективности коррекции IRI-подобных моделей является сезонная зависимость точности описания моделью наблюдаемого пространственного распределения ПЭС. Информация по результатам этих исследований была представлена на новостном ресурсе indicator.ru: https://indicator.ru/earth-science/modeli-ionosfery-priemnikov-sputnikov-27-02-2020.htm. Во время геомагнитных возмущений методика коррекции модели по данным ГНСС приемников не ухудшает свою эффективность, а в некоторых случаях эффективнее корректирует модель, чем для спокойных условий. Показано, что в целом коррекция модели улучшает описание ПЭС для всех месяцев за исключением декабря. Улучшение модельного описания ПЭС за счет коррекции по выбранной методике наилучшим образом проявляется в окрестности ГНСС приемника, но в целом имеет достаточно широкий радиус действия. На основе кросс-корреляционного анализа данных канадского радара некогерентного рассеяния, результатов теоретической модели ГСМ ТИП и индекса солнечной активности F10.7 выявлен высотный профиль зависимости задержки вариаций электронной концентрации относительно вариаций индекса солнечной активности. Вертикальная структура этой задержки содержит максимум на высотах F2 слоя ионосферы, резко уменьшается при уменьшении высоты и довольно слабо уменьшается с ростом высоты. Исследования показали, что эта задержка больше для foF2 по сравнению с ПЭС и для высоких широт по сравнению со средними широтами. С помощью результатов модели ГСМ ТИП выявлено, что основным возможным механизмом формирования указанной выше задержки является изменение отношения O/N2. Данное исследование указывает на необходимость продолжения исследований солнечно-циклических, сезонных, долготных и суточных изменений задержек вариаций электронной концентрации относительно вариаций солнечной активности для дальнейшего включения данного факта в эмпирические модели ионосферы. За отчетный период были осуществлены следующие улучшения моделей распространения КВ-радиоволн в ионосфере: 1) в модели, основанной на решении системы бихарактеристических уравнений, осуществлен учет расходимости радиолучей при расчете энергии, приходящего радиосигнала; 2) в модели расчета траекторных характеристик радиотрасс на основе вариационного метода доработан и программно реализован алгоритм систематического поиска множества лучевых траекторий с заданными граничными условиями, получивший название глобальной оптимизации (об этом сообщалось на новостных ресурсах indicator.ru, РИА новости и др.). Обе численные модели расчета распространения радиоволн дополнены блоками расчета волновых фронтов и ионограмм наклонного зондирования по найденным семействам лучевых траекторий. Для расширения возможностей вариационного метода получены аналитико-численные выражения для функционала фазового пути в анизотропной среде, ее первой и второй производных. Реализован предварительный вариант оптимизационной схемы расчета “верхних” лучей (лучи Педерсена) с учетом геометрии луча и ориентации волнового фронта в анизотропной ионосфере, требующий дальнейшего развития и тестирования. Численный анализ собственных значений матрицы Гессе волноводных и многоскачковых лучевых траекторий показал, что эти решения соответствуют седловым точкам порядка N функционала фазового пути. Порядок седловых точек N равен равен числу каустик, которых касается луч. Этот результат дает возможность на основании метода обобщенной силы разработать и реализовать метода поиска волноводных и многоскачковых лучей на основе вариационного метода. Исследования влияния учета анизотропии при расчете распространения радиоволн позволили объяснить на ионограммах наклонного зондирования более пологий след от лучей Педерсена, которые располагаются ближе к следам от нижних лучей, по сравнению с результатами при учете анизотропии. Учет магнитного поля приводит к сильной азимутальной расходимости лучей, что затрудняет их нахождение в условиях трехмерно неоднородной анизотропной среды. Основываясь на спутниковых данных ИНТЕРБОЛ-2 низкочастотного электрического поля и функции распределения ионов от 1 эВ до 20 кэВ было показано, что поперечный нагрев ионов на полярной границе аврорального овала во время суббури происходит быстрее, чем ионно-циклотронный нагрев. Основываясь на изображениях аврорального овала в ультрафиолетовом спектре показано, что существенную роль в поперечном нагреве ионов играют индуцированные электрические поля, вызванные быстрым движением/расширением аврорального овала к полюсу. Нагретые таким образом ионы вследствие первого адиабатического инварианта ускоряются вдоль магнитного поля и составляют существенную часть оттока ионосферных ионов в магнитосферу. Показано, что сбои в навигационных сигналах GPS связаны с высыпанием заряженых частиц в высокоширотной ионосфере. Подтверждены результаты других работ, что частота L2 имеет большую вероятность сбоев, чем L1, как в спокойных условиях, так и во время геомагнитной активности в полярных областях. В результате исследования впервые продемонстрировано, что вероятность больших скачков ПЭС в несколько раз выше, чем вероятность сбоя по фазе на частотах L1 и L2. Максимум возникновения сбоев приходится на фазу восстановления геомагнитной суббури. Рост Bz компоненты межпланетного поля на фазе восстановления суббури, является косвенным индикатором роста турбулизации параметров плазмы. Проведенный в работе корреляционный анализ взаимосвязей сбоев параметров L1, L2 и скачков ПЭС показал их значительную корреляцию с динамикой развития магнитной суббури и развитием оптических явлений, связанных с авроральными высыпаниями. Анализ развития ионосферных неоднородностей, полученных по данным УКВ и КВ радара и ROT по данным цепочки ГНСС приемников во время геомагнитной бури 27-28 мая 2017 года показал: (а) одновременное начало появления радиоавроры в УКВ диапазоне и интенсивных вариаций ROT для ГНСС приемника вблизи области зоны действия УКВ радара; (б) одновременное начало появления радиоавроры в КВ диапазоне и интенсивных вариаций ROT для ГНСС приемников вблизи области зоны действия КВ радара; (в) синхронное снижение интенсивности радиоавроры по данным УКВ и КВ радара и ROT в период ослабления интенсивности продольных токов в полуденные часы MLT при переходе от утренней к вечерней структуре продольных токов первой и второй зоны. Временной интервал усиления продольных токов первой и второй зоны и их смещения в сторону экватора совпадает с усилением интенсивности радиоавроры в УКВ диапазоне и увеличением флуктуаций ROT. Это говорит о том, что данные ГНСС приемников и магнитометров могут использоваться как сопутствующие данные при решении задачи радиолокации. Показана принципиальная возможность проводить измерения концентрации плазмы и ее неоднородностей в ионосфере Земли, в том числе и на полярных широтах, на сверхмалых космических аппаратах (наноспутники, кубсаты) с применением радиофизических методов, позволяющие по свойствам принимаемого излучения определять характеристики среды, через которую проходило излучение. Оценки расстояния между двумя спутниками показали, что возможно проводить измерения на расстояниях от 100 м до 100 км. Такие подходы, предложенные участниками проекта, могут быть использованы для исследования ионосферных неоднородностей не только естественного происхождения, но и искусственной ионосферной турбулентности.

 

Публикации

1. Клименко М.В., Захаренкова И.Е., Клименко В.В., Лукьянова Р.Ю., Черняк Ю.В. Simulation and observations of the polar tongue of ionization at different heights during the 2015 St. Patrick's Day storm Space Weather, V. 17, 1073–1089 (год публикации - 2019) https://doi.org/10.1029/2018SW002143

2. Клименко М.В., Карпачев А.Т., Ратовский К.Г.,, Котова Д.С., Оводенко В.Б., Ясюкевич Ю.В., Ясюкевич А.С., Жбанков Г.А., Теменс Д., Клименко В.В. Ionosphere as a Medium of Radio Wave Propagation in Different Applied Tasks IEEE 2019 Russian Open Conference on Radio Wave Propagation (RWP), P. 117-120 (год публикации - 2019) https://doi.org/10.1109/RWP.2019.8810245

3. Клименко М.В., Клименко В.В., Бессараб Ф.С., и др. Последние достижения и проблемы в моделировании системы термосфера–ионосфера ТРУДЫ МЕЖДУНАРОДНОЙ БАЙКАЛЬСКОЙ МОЛОДЕЖНОЙ НАУЧНОЙ ШКОЛЫ ПО ФУНДАМЕНТАЛЬНОЙ ФИЗИКЕ и Конференции молодых ученых «Взаимодействие полей и излучения с веществом», Лекция. БШФФ-2019, С. 11-16 (год публикации - 2019)

4. Котова Д.С., Захаренкова И.Е., Клименко М.В., Оводенко В.Б., Тютин И.В., Чугунин Д.В., Чернышов А.А., Ратовский К.Г., Чирик Н. В., Успенский М. В., Клименко В. В., Рахматулин Р. А., Пашин А. Ю., Дмитриев А. В., Суворова А. В. Formation of Ionospheric Irregularities in the East Siberian Region during the Geomagnetic Storm of May 27–28, 2017 Russian Journal of Physical Chemistry B, Vol. 14, No. 2, pp. 377–389 (год публикации - 2020) https://doi.org/10.1134/S1990793120020232

5. Котова Д.С., Оводенко В.Б., Ясюкевич Ю.В., Клименко М.В., Ратовский К.Г., Мыльникова А.А., Андреева Е.С., Козловский А.Е., Коренькова Н.А., Нестеров И.А., Туманова Ю.А. Efficiency of updating the ionospheric models using total electron content at mid- and sub-auroral latitudes GPS Solutions, V. 24, N. 25 (год публикации - 2020) https://doi.org/10.1007/s10291-019-0936-x

6. Носиков И.А., Клименко М.В., Бессараб П.Ф. Глобальная оптимизация как способ исключения проблемы задания начальных условий при модельном расчете радиотрасс Сборник трудов XXVI Всероссийской открытой научной конференции «Распространение радиоволн» (РРВ-26), Том II, с. 426–429 (год публикации - 2019)

7. Носиков И.А., Клименко М.В., Жбанков Г.А., Подлесный А., Иванова В., Бессараб П.Ф. Generalized Force Approach to Point-to-Point Ionospheric Ray Tracing and Systematic Identification of High and Low Rays IEEE Transactions on Antennas and Propagation, - (год публикации - 2019) https://doi.org/10.1109/TAP.2019.2938817

8. Носиков И.А., Толченников А.А., Клименко М.В., Доброхотов С.Ю. Вариационный метод расчета лучевых траекторий и фронтов волн цунами, порожденных локализованным источником Журнал вычислительной математики и математической физики, Т. 60, №. 8. (год публикации - 2020)

9. Чернышов А.А., Милош В.Й., Жин Й., Захаров В.И. Relationship between TEC jumps and auroral substorm in the high-latitude ionosphere Scientific Reports, 10, 6363 (год публикации - 2020) https://doi.org/10.1038/s41598-020-63422-9

10. - Ученые придумали новый способ расчета распространения радиоволн Indicator.ru, 14 ФЕВРАЛЯ В 16:17 ТЕХНИЧЕСКИЕ НАУКИ (год публикации - )

11. - Российские ученые сделали открытие в области радио- и сейсмоволн РИА Новости, 09:18 04.02.2020 (год публикации - )

12. - Физики научились точно предсказывать распространение коротких радиоволн сайт Балтийского федерального университета им. И. Канта, 4 Февраля 2020 (год публикации - )

13. - Калининградский учёный рассказал об открытии в области радио- и сейсмоволн новостной ресурс klops.ru, 04.02.2020 09:06 (год публикации - )

14. - Ученые придумали новый способ расчета распространения радиоволн РНФ, Новости Президентской программы исследовательских проектов (год публикации - )

15. - Модели ионосферы Земли скорректировали по данным приемников сигналов спутников навигации Журнал: “Все о Космосе”, 3:24 28/02/2020 (год публикации - )

16. - Модели ионосферы Земли скорректировали по сигналам спутников навигации indicator.ru, 27-02-2020 (год публикации - )

17. - Модели ионосферы Земли скорректировали по сигналам спутников навигации Вестник Отделения Наук о Земле РАН, 28 февраля 2020 г. (год публикации - )


Возможность практического использования результатов
Распределение различных ионосферных параметров, и, в особенности, электронной концентрации, оказывают влияние на условия прохождения радиосигналов, что важно для работы систем связи, телекоммуникационных систем, загоризонтной и надгоризонтной радиолокации и глобальных навигационных систем (включая отечественную спутниковую систему ГЛОНАСС). Существующие эмпирические модели ионосферы, чаще всего используемые в качестве моделей среды для задач распространения радиоволн, недостаточно точно описывают параметры среды, в том числе в периоды геомагнитных возмущений. Это относится и к наиболее часто используемой международной справочной модели IRI, особенно в высокоширотной области. В ходе выполнения проекта созданы глобальные эмпирические модели, включающие высокоширотную ионосферу. Увеличение темпов развития программ по освоению Арктики и наметившийся рост в использовании высоких широт воздушно-космическими силовыми структурами и для перевозки пассажиров воздушным транспортом делают задачу обеспечения надежности радиосвязи одной из наиболее важных. Одним из перспективных подходов к расчету КВ-радиотрасс является прямой вариационный метод, позволяющий проводить расчеты между приёмной и передающей станциями, положение которых зафиксировано. Кроме того, стоит отметить, что разработанный вариационный метод расчета радиотрасс универсален и может применяться в любых задачах, где справедливо приближение геометрической оптики, например для расчета распространения волн цунами и сейсмических волн. Для обеспечения надежной и непрерывной радиосвязи нужно решать задачу уточнения частотного диапазона приемно-передающих станций радиосвязи в конкретных гелио-геофизических условиях в конкретных регионах. Решению таких задач было уделено особое внимание в ходе выполнения данного проекта. Результаты выполненной работы будут интересны для разных служб и ведомств России. Об этом свидетельствует недавний проект, выполняемый ОАО НПК НИИДАР для обеспечения надежной работы систем радиолокации и использующий результаты исследований эффективности корректировки модели ионосферы, полученные в ходе выполнения данного проекта.