Бесплатный автореферат и диссертация по наукам о земле на тему
Крупномасштабные явления в солнечном ветре и магнитосферно-ионосферная активность
ВАК РФ 25.00.29, Физика атмосферы и гидросферы

Автореферат диссертации по теме "Крупномасштабные явления в солнечном ветре и магнитосферно-ионосферная активность"

На правах рукописи

д/

003452834

Ревунов Сергей Евгеньевич

КРУПНОМАСШТАБНЫЕ ЯВЛЕНИЯ В СОЛНЕЧНОМ ВЕТРЕ И МАГНИТОСФЕРНО-ИОНОСФЕРНАЯ АКТИВНОСТЬ

25.00.29 - физика атмосферы и гидросферы

АВТОРЕФЕРАТ

диссертации на соискание учёной степени кандидата физико-математических наук

1

.. J

Нижний Новгород, 2008 г.

003452834

Работа выполнена в Нижегородском государственном педагогическом университете.

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

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

доктор физико-математических наук, профессор Деминов Марат Гарунович кандидат физико-математических наук Ягодкина Оксана Ивановна

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

Объединенный институт физики Земли РАН, (ОИФЗ), г. Москва

Защита диссертации состоится « 2 » декабря 2008 г. в 16 часов 00 мин, на заседании диссертационного совета Д 002.237.01 в Учреждении Российской Академии наук Институте земного магнетизма, ионосферы и распространения радиоволн им. Н.В. Пушкова РАН по адресу: 142190, г. Троицк, Московской области (проезд автобусом №398 от станции метро «Тёплый стан» до остановки «ИЗМИРАН»),

С диссертацией можно ознакомиться в библиотеке ИЗМИР АН. Автореферат разослан «31» октября 2008 а

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

диссертационного совета Д 002.237.01 доктор физико-математических наук: г ^_^ /Михайлов Ю.М./

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

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

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

Актуальность проблемы

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

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

Еще одно направление космических исследований получило развитие по мере накопления спутниковой информации о солнечной активности и состоянии межпланетной

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

»

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

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

модельного описания энергетического бюджета магнитосферы и пространственно-временного распределения магнитного поля в ней.

В настоящее время одновременно с развитием радиосвязи с использованием космических аппаратов продолжается использование коротковолновой (КВ) радиосвязи через среднеширотную ионосферу. Для повышения надежности такой радиосвязи необходимо совершенствование оперативного прогнозирования ионосферных параметров на время от 30 мин до нескольких часов. Оно может быть осуществлено на основе прогнозирования критической частоты ионосферного слоя ¥2 на различные интервалы времени с помощью ИНС с учетом геофизических параметров.

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

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

Выполнено исследование связи основных параметров ионосферного коротковолнового радиоканала - максимально наблюдаемой частоты (МНЧ) ионосферного распространения КВ радиосигнала и критической частоты №2 среднеширотной ионосферы с геоэффективными параметрами межпланетной среды.

Цели и задачи работы '

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

1. Исследование возможности выполнения классификации вариаций индекса геомагнитной активности Бз! в периоды магнитосферных возмущений в зависимости от вызывающего их солнечного источника на основе сопоставления

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

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

3. ' Анализ связи максимально наблюдаемой и критической частоты ионосферного

распространения коротковолнового радиосигнала с параметрами солнечного ветра и межпланетного магнитного поля для повышения надежности ионосферной KB радиосвязи.

»

Научная новизна

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

На основе одножидкостного магнитогидродинамического подхода с применением технологии искусственных нейронных сетей разработана методика автоматического разделения скачков в регистрируемых параметрах космической плазмы и магнитного поля на классы, отвечающие известным типам разрывов. Для этого была разработана классификационная ИНС вида «слой Кохонена» и выполнена классификация скачков параметров зарегистрированных на КА WIND в 1996-1999 гт. по двум разным алгоритмам. Для разрывов установленных классов определены ориентации плоскостей их поверхностей в трехмерном пространстве на основе одномерных наблюдений. Также проведено специальное исследование эволюционности найденных ударных волн.

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

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

Разработан алгоритм прогнозирования на интервал 0,5-3 часа максимально наблюдаемой и критической частоты на основе технологии искусственных нейронных сетей. Методика прогнозирования применена к данным наклонного зондирования на трассе 1пБ1ар (Англия) - Ростов-на-Дону, которые получены для различных геофизических условий. Определены характерные времена ионосферной реакции на изменения параметров солнечного ветра и межпланетного магнитного поля, определяющих последовательность развития магнитосферно-ионосферных возмущений, приводящих к изменениям КВ радиосвязи.

Научная и практическая ценность

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

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

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

На основе анализа данных наклонного зондирования с использованием линейно-частотной модуляции (ЛЧМ), полученных на трассе Inskip (Англия) - Ростов-на-Дону, установлена связь максимально наблюдаемой и критической частоты с изменениями параметров солнечного ветра и межпланетного магнитного поля, определяющих последовательность развития магнитосферно-ионосферных возмущений, что приводит к изменению условий KB радиосвязи. Практическая ценность выполненной работы заключается в применении ее результатов для повышения надежности ионосферной KB радиосвязи.

Степень достоверности полученных результатов

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

»

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

Экспериментальные данные

В диссертационном исследовании использовались следующие экспериментальные данные:

1. Параметры солнечного ветра и ММП с разным усреднением, полученные из базы спутниковой системы OMNI [http://www.ngdc.noaa.gov]

2. Данные по индексам глобальной геомагнитной активности [http://www.ngdc.noaa.gov; http://cdaweb.gsfc.nasa.gov/istp_public;

. http://swdcwww.kugi.kyoto-u.ac.jp/wdc/Sec3.html]

3. Данные по состоянию ионосферного слоя F2, предоставленные НИРФИ

4. Данные о состоянии межпланетной среды, полученные с космических аппаратов (КА) Wind и Асе [http://www.srl.caltech.edu/ACE/ASC/rtsw.html; http://cdaweb.gsfc.nasa.gov/istp_public]

5. Данные о состоянии межпланетной среды, полученные из базы SPIDR [http://clustl.wdcb.ru/spidr/index.jsp]

6. Данные по разрывам параметров межпланетной среды, предоставленные ИКИ РАН.

На защиту выносятся следующие положения:

1.' Классификация Dst-индекса на основе параметров солнечного ветра и межпланетного магнитного поля, относящихся к конкретным явлениям солнечной активности: вспышки, корональные дыры, волокна, корональные выбросы вещества и их сочетания; классификация разрывов в межпланетной плазме на

основе метода нейронных сетей и определение трехмерных ориентаций их плоскостей на основе одномерных спутниковых наблюдений.

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

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

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

Результаты диссертационной работы докладывались на следующих семинарах и

конференциях:

1. IX Международная научно-техническая конференция «Радиолокация, навигация и связь», 22-24 апреля 2003 г., Воронеж

2. XXIII General Assembly, Sapporo, June 30-Му 11,2003

3. VIII Нижегородская сессия молодых ученых. Естественнонаучные дисциплины, 2025 апреля 2003 г., Дзержинск

4. Десятая Всероссийская научная конференция студентов-физиков и молодых ученых, 1-7 апреля 2004 г., Москва

5. European Geosciences Union, General Assembly, Nice, France, 25-30 April 2004

6. Ill Международная конференция «Солнечно-земные связи и электромагнитные предвестники землетрясений», 16-21 августа 2004, Институт космофизических исследований и распространения радиоволн (ИКИР) ДВО РАН (с. Паратунка, Камчатская обл.)

7. Международная конференция «Солнечно-земная физика», Иркутск, 20-25 сентября 2004 г.

8. Семинар «Физика авроральных явлений», 1-4 марта 2005 г., Полярный геофизический институт

9. 21 Всероссийская конференция по распространению радиоволн, Йошкар-Ола, 2427 мая 2005 г.

10. 10th IAGA Scientific Assembly, Toulouse-France, July 18-29, 2005

11. 9th Asian-Pacific Regional International Astronomical Union (IAU) Meeting (APRIM-2005), Bali, Indonesia, 26-29 July 2005

12. Всероссийская конференция «Экспериментальные и теоретические исследования основ прогнозирования гелиогеофизической активности» (10-15 октября 2005 г., г. Троицк)

13. XXVIIIURSI General Assembly New Delhi, India, 23 to 29 October 2005

14. 29 Annual Seminar «Physics of auroral phenomena», Apatity, 27 February-3 March 2006, PGIKSC RAS

15. XI Нижегородская сессия молодых ученых. Математические науки, 22-25 мая 2006 г., Нижний Новгород

16. 36th COSPAR SCIENTIFIC ASSEMBLY, Beijing, China, 16-23 July 2006

17. 11 Пулковская международная конференция по физике Солнца «Физическая природа солнечной активности и прогнозирование ее геофизических проявлений», ГАО РАН, Пулково, Санкт-Петербург, 2-7 июля 2007

18. European Geosciences Union, General Assembly 2007, Vienna, Austria, 15-20 April

2007

19. 31 Annual Seminar «Physics of auroral phenomena», Apatity, 26-29 February 2008, PGI KSC RAS

20. European Geosciences Union, General Assembly 2008, Vienna, Austria, 13-18 April

2008

21. Семинары НИРФИ, ИЗМИР АН и НГПУ Публикации

По теме диссертации автором опубликовано 6 статей в рецензируемых журналах и 27 работ в сборниках и трудах конференций.

Личный вклад соискателя

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

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Диссертация состоит из введения, четырех глав и заключения, содержит 123 страницы, 44 рисунка, 2 схемы и 5 таблиц. Библиография включает в себя 100 наименований.

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

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

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

В третьем параграфе проводится классификации геомагнитных возмущений, характеризующихся индексом геомагнитной активности Ов^ с привлечением данных о концентрации плазмы солнечного ветра, ее скорости и вертикальной составляющей вектора межпланетного магнитного поля. Концентрация и скорость плазмы характеризуют кинетическое давление солнечного ветра, а вертикальная Вг компонента межпланетного магнитного поля является наиболее геоэффективным параметром, от знака и амплитуды которой зависит поток энергии внутрь магнитосферы из межпланетной среды: при Вг<0 поток усиливается. Для проведения классификации используется метод ИНС, который позволяет выявить ноЕые и уточнить уже установленные закономерности связи параметров геомагнитной возмущенности с параметрами солнечной активности.

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

В четвертом параграфе экспериментально устанавливается число классов событий космической погоды наиболее адекватное для данной задачи классификации, а в пятом проводится отождествление полученных классов с вызвавшими их явлениями солнечной активности. Классификация выполнялась по набору параметров: Dst, N, V, Bz. Данный набор отвечает физической сути задачи, моделируя причинно-следственные связи. Этим демонстрируется возможность объединения параметров, характеризующих причины (источники солнечной активности и структура солнечного ветра) и следствие (геомагнитное возмущение, представляемое Dst(t) вариацией). Таким образом описывается взаимодействие подошедшего потока солнечной плазмы с магнитосферой.

На Рис. 1, в качестве примера, приведены два из восьми найденных классоб, выделенных с помощью разработанного метода классификации. Показаны относящиеся к этим классам конфигурации индекса Dst(t), концентрации плазмы солнечного ветра N(t), скорости солнечного ветра V(t) и вертикальной компоненты вектора ММП Bz(t).

17Л2тОШ> 17Я2М-23М 05/08/98-10:00 >0вЛ8/9«-10:00 10/1 СЛ8-00:00 > te/f0i»&-14:00 13/01/99-0? 00 >13/01 #9-18:00

Источник - потоки СВ

с выраженными ударными волнами

5 10 15 20 26 hours

5 10 16 20 hours

ь)

д 1V во 50 40 30 к 550 500 460 400 10 6 0 -5

\ \ 20 10 1 d- / V 360 300 -16 •20

5 10 15 20 25 hours

5 10 15 20 25 hours

5 10 15 20 25 hours

>8-20:00 16/04/99-04:00 > 17/04/99-03:00 12,09/99-01;00> 12Д9Я9-10:00

Источник - KBB s форме магнитного облака

5 10 16 20 26 hours

Рис. 1. Классы различных сильных магнитных бурь, вызванных разними источниками - ударными волнами (а) и КВВ в форме магнитных облаков (Ь).

Заметим, что разделение приведенных классов выполнено не только на основе амплитуды индекса Бб!, но и с учетом формы, т.е. динамики Бз^вариации: на Рис. 1а индекс имеет отчетливый провал, а на Рис. 1Ь провалу предшествует повышение. Участвующие в классификации параметры К, V и Вг так же имеют свои уникальные особенности внутри класса, на основе которых можно сделать вывод о солнечном источнике. Так, колебания

скорости с повышением, рост плотности и напряженности ММГГ на фоне невозмущенного солнечного ветра указывают на наличие ударной волны в потоке (см. Рис. 1а), а волновые возмущения в записях ММП совместно с ростом скорости потока и наличием пика плотности могут соответствовать магнитному облаку (см. Рис. 1Ь). Таким образом разработан комплексный подход к изучению геомагнитных возмущений как следствия солнечно-земных связей. Комплекс событий космической погоды (отклонения параметров солнечного ветра и ММП от фоновых значений и последующие изменения в индексе Эв!) является основой новой классификации геомагнитных возмущений содержащей их причину. Классифицируются не просто типы возмущений в солнечном ветре и амплитуда индекса Бэ!, а события в совокупности по схеме «причина+следствие».

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

Nfcm3) |B|(iiTI) 2¡P(Pa)x10'18

timo, sec

Рис. 2. Параметры события 04/04/96 типа ударная волка, полученные на основе измерений на КА (а); его расчетная плоскость разрыва в солнечно-эклиптической системе координат SE в условных единицах (б).

Седьмой параграф содержит результаты исследования по поиску оценки связи Кр индекса геомагнитной активности с электрическим полем солнечного ветра Есв и с проекцией этого поля на геомагнитный диполь. Показано, что для всех исследуемых лет, суммарно по сезонам года, максимум классического коэффициента корреляции Л(Кр,УВз), где V - скорость солнечного ветра, а Вб - южная компонента вектора ММП, наблюдается для интервалов осеннего и весеннего равноденствия и равен 0,52. Многофакторный коэффициент корреляции ДКр, Е-1СВ* Е) | св), где Ехсв — компонента вектора электрического поля солнечного ветра, перпендикулярная геомагнитному диполю, а Е| | св — компонента, параллельная ему, максимален для летнего периода (0,48) и минимален дня зимнего (0,38). Подтвержденный нами максимум сезонного изменения коэффициентов корреляции свидетельствует о влиянии ориентации электрического поля солнечного ветра относительно земного магнитного диполя на состояние геомагнитной активности. При этом результаты распределения уровня корреляции по сезонам наиболее резко проявляются при представлении Есв в системах координат, связанных с геомагнитным диполем - в солнечно-магнитосферной (БМ) и солнечно-магнитной (ЭМАО).

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

Рис. 3. Результат классификации на три класса спектров модуля геомагнитного поля за 26 августа 1998 года. Представлен вид со стороны северной полярной шапки. Точкой в центре отмечен геомагнитный полюс Земли, а) - овал, определяемый на основе его современной модели ОМЭР; Ь) - овал, определяемый на основе классификации спектров модуля поля.

Каждый сегмент на рисунке соответствует определенной станции. Цветом показаны три выделенных класса. Темные сегменты-станции соответствуют среднеширотной зоне (класс 1), серые - овалу полярных сияний (авроральная зона, класс 2) и светлые - полярной шапке (полярная зона, класс 3). Вычисляя коэффициент подобия, учитывающий площадь сравниваемых частей сегментов и их относительное расположение, было получено, что согласованность сетки-ответа нейронной сети с сеткой модельной авроральной зоны для разных дней составила 61-65%.

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

В пятом параграфе выполнен поиск физических причин разделения спектральных характеристик магнитных возмущений в полярных областях. Показано, что в спектрах станций, принадлежащих разным зонам, присутствуют характерные пики частот. Так, в авроральной зоне отмечен пик на частоте 7,5*10"3 Гц. Данная частота соответствует плазменным колебаниям - устойчивым пульсациям Рс4. Амплитуда таких пульсаций

обычно возрастает с приближением к зоне полярных сияний. В субавроральной зоне обнаружен пик на частотах 0,7-1,1*10"3 Гц. Этот низкочастотный пик согласно установившемуся мнению может быть вызван собственньми колебаниями замкнутых силовых линий магнитного поля Земли. Выполненный корреляционный анализ спектров геомагнитного поля на станциях, отнесенных ИНС к полярной шапке, и спектров ММП в солнечном ветре показал их высокую согласованность (корреляция до 0,96). Это дает основание полагать, что за вариации в полярной шапке ответственны колебания в солнечном ветре. Таким образом при работе нейронная сеть устанавливает спектральные особенности вариаций геомагнитного поля, измеренного меридиональной цепью наземных станций, и делает вывод о принадлежности станции к той или иной полярной области. Далее на этой основе нейросеть выполняет оценку границ полярной шапки и аврорального овала.

Четвертая глава посвящена исследованию среднеширотной ионосферы. Рассмотрена связь максимально наблюдаемой и критической частоты слоя F2 с параметрами околоземного космического пространства. В первом параграфе приводятся результаты базового корреляционного исследования связи частот с наиболее геоэффективными параметрами межпланетной среды. Во втором параграфе представлены результаты анализа временного шага между временем наблюдения параметров межпланетной среды и прогнозируемым временем начала использования этой частоты. Вычисление корреляции критической частоты foF2 и максимально наблюдаемой с параметрами межпланетной среды необходимо для обоснования использования при обучении ИНС дополнительных солнечно-земных параметров. При этом учитывается задержка на время переноса возмущения в солнечном ветре от КА до магнитосферы и на время реакции магнитосферы на изменение ПСВ. Так, в результате проведенного корреляционного анализа, была установлена оптимальная задержка, при которой достигаются максимальные коэффициенты корреляции, равная трем часам. Перенос возмущенных параметров солнечного ветра и ММП от патрульного КА до границ магнитосферы занимает в среднем 50-60 мин (по данным Space Science Center, Institute of Geophysics and Planetary Physics, University of California, Los Angeles, http://www-ssc.igpp.ucla.edu). Согласно установившимся представлениям оставшиеся два часа тратятся на разогрев нейтральной атмосферы вследствие джоулевых потерь под действием вновь пришедших и высыпавшихся частиц. Данная задержка связана с инерцией атмосферы. На высотах больше 250 км (высота слоя F2 - 300-400 км) скорость химической рекомбинации ионов сравнима со скоростью их диффузии в нейтральном газе и механизмы переноса (диффузия) начинают влиять на величину ионной и электронной

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

становятся ионы О*:

о*+о2-+о;+о

О'+Л^ ->МГ+ ЛГ

. Кроме того, протекают многие другие

реакции, однако, приведенные показывают, как могут образоваться ионы N0* и 0\,

являющиеся на ряду с О* наиболее многочисленными.

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

В четвертом параграфе продемонстрирована методика поиска дополнительных входных параметров, а так же комбинаций параметров для повышения качества прогнозирования. Рассмотрен вопрос о необходимости разделения технологии прогнозирования МНЧ для спокойных и возмущенных интервалов. Пятый параграф посвящен вопросу повышения качества восстановления профиля МНЧ при добавлении производных высокого порядка ряда МНЧ в число входных параметров при обучении и дальнейшей работе нейросети. Показано, что с помощью разветвленной ИНС можно проводить прогнозирование МНЧ с заблаговременностью 0,5-3 часов с точностью не ниже 85%. Диаграмма, представленная на Рис. 5, иллюстрирует динамику качества прогноза в каждой серии численных экспериментов для майского интервала данных. В этих экспериментах в качестве дополнительного параметра добавлялся индекс БвЪ или рентгеновское излучение, или Вг компонента вектора ММП, задержанная на 3 часа относительно частоты. Эти параметры давали приблизительно одинаковый вклад в изменение качества прогноза.

22-03-2002 23-03-2002 24-03-202 25-03-2002

часы. ЦТ &л)

Рис. 4. Результат прогноза на I час для интервала с 22 марта 00:00 по 25 марта 23:59. Входной массив -последовательность критических частот и её производная. Сплошной линией показана реальная последовательность, пунктиром - прогноз. РЕ=92%, 11=0.96.

■ ([3] Hf+di*df2+d0*add[31

0f+df[3] B1+df+df2+dG+ad<l[4)

я hdfcdQ (3] ■ f+d1+df2+dr3tadd[4'l *T+dtKli2+dG[3]

05 1 i .5 2 2i 3

Рис. 5. Динамика эффективности прогнозирования МНЧ (FKp) на майском интервале при последовательном добавлении производных высокого порядка и дополнительного параметра при обучении разветвленной сети с помощью механизма контрольной проверки.

Столбцы на рисунке 5 соответствуют последовательному добавлению производных (df, d£2, dfi) и дополнительного параметра (add) в обучающий массив. В квадратных скобках указан тип сети. Знак «*» означает применение дополнительного алгоритма коррекции к ответному сигналу ИНС.

В заключении представлены основные результаты и выводы проведенного исследования:

1. Выполнена классификация скачков параметров солнечного ветра по данным КА WIND за трехлетний период их измерений (1996-1999 гг.) используя одножидкостный МГД

-18-

подход с помощью специально созданной классификационной искусственной нейронной сети вида «слой Кохонена». Для выделенных классов скачков разработан и применен способ определения ориентации плоскостей их поверхностей в солнечном ветре по одномерным наблюдениям на КА. Выполнен анализ эволюционное™ найденных ударных волн, который продемонстрировал их потенциальную неустойчивость. 2. Выделены классы структур ОхЬвариаций, соответствующие магнитным бурям определенной интенсивности и конкретным типам возмущений параметров солнечного ветра и межпланетного магнитного поля, которые связаны с указанными солнечными плазменными потоками:

Класс №1 - слабые бури, являющиеся следствием взаимодействия с магнитосферой Земли коротирующих с Солнцем штоков (согс^й^-зй-еатз).

Класс №2 - слабые бури, являющиеся следствием взаимодействия с магнитосферой Земли медленных гелиосферных потоков с низкой концентрацией плазмы. Класс №3 - умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли коротирующих потоков, которые до этого взаимодействовали с НСЭ-потоком или гелиосферным токовым слоем.

Класс №4 - умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли быстрых, высокоинтенсивных соп^и^-потоков.

Класс №5 - умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли высокоскоростного потока, прошедшего мимо Земли и задевшего ее своим флангом.

Класс №6 - умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли сразу нескольких гелиосферных потоков.

Класс №7 - сильные бури, являющиеся следствием взаимодействия с магнитосферой Земли уединенных, не имеющих периодической повторяемости, потоков с выраженными ударными волнами.

Класс №8 - сильные бури, являющиеся следствием взаимодействия КВВ в форме магнитного облака с магнитосферой Земли.

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

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

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

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

1. Бархатов Н.А., Валов ВЛ., Макаров А.В., Понятов А.А., Ревунов С.Е., Сахаров С.Ю., Урядов В.П., Метод искусственных нейронных сетей и его применение к прогнозированию ионосферного KB радиоканала.' Труды IX Международной научно-технической конференции «Радиолокация, навигация и связь», 22-24 апреля 2003 г., Воронеж, Том III, С.1853-1864.

2. Barkhatov N.A., Ponyatov А.А., Uryadov V.P., Barkhatova O.A., Revunov S.E., Forecasting of ionospheric hf channel by the artificial neural networks method // Abstracts of IUGG XXIII General Assembly, Sapporo, Juhe 30-July 11,2003, GAII.03/02A/A13-005, C.A298

3. Бархатов H.A., Ревунов C.E., Прогнозирование критической частоты ионосферного слоя F2 методом искусственных нейронных сетей. VIII Нижегородская сессия молодых ученых. Естественнонаучные дисциплины, 20-25 апреля 2003 г., Дзержинск, Тезисы докладов, Н. Новгород, 2003, С.98-99

4. Barkhatov N.A., Levitin А.Е., Revunov S.E. Comprehensive classification of geomagnetic disturbances. Geophysical Research Abstracts, V.6, 04796, European Geosciences Union, General Assembly 2004, ST5„Nice, France, 25-30 April 2004

5. Ревунов C.E., Дудкина И.С. Повышение эффективности прогнозирования критической частоты ионосферного слоя F2. Сборник тезисов Десятой Всероссийской научной конференции студентов-физиков и молодых ученых: Тезисы докладов: В 2 т. - Москва, издательство АСФ России, 2004, Т.2. С. 915

6. Barkhatov N.A., Revunov S.E., Uryadov V.P. Improvement of ionospheric layer f2 critical frequency forecasting. Geophysical Research Abstracts, V.6, 04789, European Geosciences Union, General Assembly 2004, ST14, Nice, France, 25-30 April 2004

7. Бархатов Н.А., Левитин А.Е., Ревунов С.Е. Космическая классификация глобальных геомагнитных возмущений. Тезисы докладов III Международной конференции «Солнечно-земные связи и электромагнитные предвестники землетрясений», 16-21 августа 2004, Институт космофизических исследований и распространения радиоволн (ИКИР) ДВО РАН (с. Паратунка, Камчатская обл.), С.46

8. Бархатов Н.А., Левитин А.Е., Ревунов С.Е. Классификация глобальных геомагнитных возмущений на основе комплекса событий космической погоды. Тезисы докладов международной конференции «Солнечно-земная физика», Иркутск, 20-25 сентября 2004 г., С.46

9. Barkhatov N.A., Revunov S.E., Uryadov V.P. (2004), Forecasting of the critical frequency of the ionosphere F2 layer by the method of artificial neural networks. Int. J. Geomagn. Aeron., 5, GI2010, DOI: 10.1029/2004GI000065

10. Бархатов H.A., Ревунов C.E., Урядов В.П. Технология искусственных нейронных сетей для прогнозирования критической частоты ионосферного слоя F2. // Известия ВУЗов «Радиофизика», 2005, T.4S. С.1-15

11. Бархатов Н.А., Ревунов С.Е., Урядов В.П., Вертоградов Г.Г., Вертоградов В.Г., Валов В.А. Прогнозирование МНЧ с помощью искусственных нейронных сетей. Материалы 21 Всероссийской конференции по распространению радиоволн, Йошкар-Ола, 24-27 мая 2005 г.

12. Barkhatov N., Levitin A., Revunov S. Location of high latitude zones of geomagnetic activity by the data from the Greenland chain of magnetic stations. // Abstracts of 10-th IAGA Scientific Assembly, Toulose-France, July 18-29,2005, IAGA2005-A-01244

13. Barkhatov N.A., Revunov S.E., Ulibina R.I. Classification of discontinuities in parameters of space plasma using magnetohydrodynamic approach. Abstracts of 9th Asian-Pacific Regional International Astronomical Union (IAU) Meeting (APRIM_2005), Bali, Indonesia, 26-29 July 2005, P. 17

14. Barkhatov N.A., Gromova L.I., Dremuhina L.A, Levitin A.E., Revunov S.E., Ulibina R.I. The analysis of magnetospheric disturbance power according to its solar source. Abstracts of 9th Asian-Pacific Regional International Astronomical Union (IAU) Meeting (APRIM_2005), Bali, Indonesia, 26-29 July 2005, P.l 8

15. Бархатов H.A., Левитин A.E., Ревунов С.Е. Установление границ областей полярной шапки и овала полярных сияний по данным Гренландской цепочки магнитных станций. Тезисы Всероссийской конференции «Экспериментальные и теоретические исследования основ прогнозирования гелиогеофизической активности» (10-15 октября 2005 г., г. Троицк), С.7

16. Бархатов Н.А., Громова Л.И., Дремухина Л.А., Левитин А.Е., Ревунов С.Е., Улыбина Р.И. Анализ энергетики магнитосферного возмущения в зависимости от его солнечного источника. Тезисы Всероссийской конференции «Экспериментальные и теоретические исследования основ прогнозирования гелиогеофизической активности» (10-15 октября 2005 г., г. Троицк), С.44

17. Бархатов Н.А., Ревунов С.Е., Урядов В.П., Вертоградов Г.Г, Вертоградов В.Г. Учет проявлений солнечной активности при прогнозировании максимально наблюдаемой частоты ионосферного KB радиоканала. Тезисы Всероссийской конференции «Экспериментальные и теоретические исследования основ прогнозирования гелиогеофизической активности» (10-15 октября 2005 г., г. Троицк), С.59

18. Barkhatov N.A., Revunov S.E., Uryadov V.P., Vertogradov G.G., Vertogradov V.G. Forecasting of hf propagation with use of artificial neural networks. XXVIIIth URSI General Assembly New Delhi, India, 23 to 29 October 2005

19. Бархатов H.A., Левитин A.E., Ревунов C.E. Классификация глобальных геомагнитных возмущений на основе комплекса событий космической погоды. В сб. «Солнечно-земная физика», вып.8, Изд. СО РАН, Новосибирск, 2005, С.129-131

20. Ревунов С.Е., Улыбина Р.И. Нейросетевая классификация разрывов параметров космической плазмы, XI Нижегородская сессия молодых ученых, Математические науки, 22-25 мая 2006 г., Нижний Новгород, Тезисы докладов, 2006, С.36-37

21. Barkhatov N.A., Levitin A.E., Revunov S.E. Monitoring of polar cap spatial dynamics basing on data from Greenland chain of magnetic observatories. Abstracts of 29 Annual Seminar «Physics of auroral phenomena», Apatity, 27 February-3 March 2006, PGIKSC RAS, P.31

22. Barkhatov N.A., Revunov S.E. Classification of abrupt changes of the interplanetary medium parameters by a method of artificial neural networks. Abstracts of 29 Annual Seminar «Physics of auroral phenomena», Apatity, 27 February-3 March 2006, PGI KSC RAS, P.51

23. Бархатов H.A., Ревунов C.E., Вертоградов Г.Г, Вертоградов В.Г, Валов В.А., Урядов В.П. Прогнозирование максимально наблюдаемой частоты ионосферного КВ-радиоканала методом искусственных нейронных сетей. // Геомагнетизм и аэрономия, 2006, Т.46, N 1,С. 88-98

24. Barkhatov N.A., Dremukhina L.A., Gromova L.I., Levitin A.E., Revunov S.E., Ulibina R.I. «Magnetosphere power budget role in the task of classification of magnetospheric activity sources» 36th Cospar scientific assembly, Beijing, China, 16-23 July 2006. Session D2.I/E3.1, A-00873

25. Бархатов H.A., Громова Л.И., Дремухина Л.А., Левитин А.Е., Ревунов С.Е., Улыбина Р.И. Учет энергетического бюджета магнитосферы в задаче классификации источников магнитосферной активности. Изв.РАН (сер.физ), 2006, 70, N10, C.1535-I537

26. • Barkhatov N.A., Levitin А.Е., Revunov S.E. Complex Classification of Global Geomagnetic

Disturbances. // Kosmicheskie Issledovaniya, 2006, Vol. 44, No. 6, pp. 488-499

27.. Barkhatov N.A., Levitin A.E., Revunov S.E. Monitoring of polar cap spatial dynamics basing on data from Greenland chain of magnetic observatories. Proceedings Qf 29 Annual Seminar «Physics of auroral phenomena», Apatity, 2007, PGI KSC RAS, P.75

28. Barkhatov N.A., Revunov S.E. Classification of abrupt changes of the interplanetary medium parameters by a method of artificial neural networks. Proceedings of 29 Annual Seminar «Physics of auroral phenomena», Apatity, 2007, PGI KSC RAS, P. 175

29. Бархатов H.A., Левитин A.E., Ревунов C.E. Контроль пространственной динамики полярной шапки и полярного овала по данным гренландской сети магнитных станций. // Геомагнетизм и аэрономия, 2007, Т.47, N2, С. 230-235

30. Бархатов Н.А., Громова Л.И., Левитин А.Е., Ревунов С.Е. Анализ влияния электрического поля солнечного ветра на геомагнитную активность. Тезисы докладов 11 Пулковской международной конференции по физике Солнца: Физическая природа солнечной активности и прогнозирование ее геофизических проявлений. ГАО РАН, Пулково, Санкт-Петербург, 2-7 июля 2007, С. 17

31. Бархатов Н.А., Левитин А.Е., Ревунов С.Е., Смирнова А.С., Снегирев С.Д., Хвиюзова Т.А. Классификация событий космической погоды. Тезисы докладов 11 Пулковской международной конференции по физике Солнца: Физическая природа солнечной активности и прогнозирование ее геофизических проявлений. ГАО РАН, Пулково, Санкт-Петербург, 2-7 июля 2007, С. 19

32. Barkhatova О.М., Revunov S.E., Barkhatov N.A., Levitin A.E. Nonlinear connection of global and polar current systems on the main phase of a geomagnetic storm. Abstracts of 31 Annual Seminar «Physics of auroral phenomena», Apatity, 26-29 February 2008, PGI KSC RAS, P.12

33. Barkhatov N., Levitin A., Revunov S. Connection of geomagnetic activity with orientation of the solar wind electric field concerning the terrestrial magnetic moment. Geophysical Research Abstracts, Vol. 10, EGU2008-A-00000,2008 EGU General Assembly 2008

Оглавление диссертации Введение

ГЛАВА 1. Отражение динамнкн параметров межпланетной среды в мапштосферно-ионосферной активности

§1.1 Крупномасштабные явления солнечной активности §1.2 Физические причины геомагнитной активности

§1.3 Современные методы исследования природы магнитосферно-ионосферной активности

§1.4 Постановка решаемых задач

ГЛАВА 2. Разработка комплексного подхода к изучению геомагнитных возмущений как следствия солнечно-земных связей

§2.1 Геомагнитные возмущения как следствие крупномасштабных явлений солнечной активности

§2.2 Нейросетевой классификационный подход

§2.3 Задачи классификации геомагнитных возмущений на основе солнечно-земных связей §2.4 Принципы отбора глобальных геомагнитных вариаций и поиск адекватного числа их классов

§2.5 Определение космических источников геомагнитных возмущений

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

§2.7 Анализ влияния электрического поля солнечного ветра на геомагнитную активность

§2.8 Выводы

ГЛАВА 3. Особенности высокоширотных геомагнитных возмущений §3.1 Необходимость контроля пространственной динамики полярной шапки и аврорального овала

§3.2 Разработка методики классификации пространственного положения станций по результатам Фурье-анализа геомагнитных данных

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

§3.4 Обработка результатов классификации и оценка уровня согласованности

классификационного и традиционного модельного подходов

§3.5 Поиск физических причин разделения спектральных характеристик магнитных

возмущений в полярных областях

§3.6 Выводы

ГЛАВА 4. Исследование связи максимально наблюдаемой и критической частоты слоя Г2 среднеширотной ионосферы с геоэффективными параметрами межпланетной среды

§4.1 Корреляционное исследование связи критической частоты слоя Р2 и МНЧ с

ключевыми геоэффективными параметрами межпланетной среды

§4.2 Определение времени реакции критической частоты слоя Р2 ионосферы и МНЧ на

взаимодействие межпланетной среды с магнитосферой

§4.3 Поиск оптимальной нейросетевой архитектуры для прогноза ионосферных

параметров

§4.4 Определение дополнительных входных параметров нейросети для повышения

качества прогнозирования критической частоты слоя Р2 и МНЧ

§4.5 Восстановление критической частоты и МНЧ на разветвленной нейросети методом

добавления производных высокого порядка

§4.6 Выводы

Заключение

Литература

Благодарности

При получении результатов, представленных в диссертации, я опирался на знания и профессиональный опыт моего учителя и научного руководителя доктора физико-математических наук, профессора Николая Александровича Бархатова; сотрудников лаборатории геомагнитных вариаций ИЗМИРАН Анатолия Ефимовича Левитана, Людмилы Ивановны Громовой, Лидии Алексеевны Дремухиной; сотрудника ИКИ РАН Георгия Наумовича Застенкера; сотрудника НИРФИ Валерия Павловича Урядова. Большую помощь и поддержку мне оказали также сотрудники Научно-исследовательской лаборатории Физики Солнечно-Земных связей Нижегородского Государственного Педагогического Университета. Всем им выражаю свою признательность и благодарность.

© С.Е. Ревунов (Научно-исследовательская лаборатория Физики Солнечно-Земных связей НГПУ, http://spacelabnnov.l10mb.com)

Крупномасштабные явления в солнечном ветре и магнитосферно-ионосферная активность

Подписано к печати qkjt>2008 г. У сл. печ. л. 1,5. Заказ 4Z7 Тираж 100 экз.

Отдел полиграфии AHO «МУК НГПУ». 603950, Нижний Новгород, ГСП-37, ул. Ульянова, 1

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

ВВЕДЕНИЕ.

ГЛАВА 1. Отражение динамики параметров межпланетпой среды в магнитосферно-ионосферной активности

§1.1 Крупномасштабные явления солнечной активности.

§1.2 Физические причины геомагнитной активности.

§1.3 Современные методы исследования природы магнитосферно-ионосферной активности.

§ 1.4 Постановка решаемых задач.

ГЛАВА 2. Разработка комплексного подхода к изучению геомагнитных возмущений как следствия солнечно-земных связей

§2.1 Геомагнитные возмущения как следствие крупномасштабных явлений солнечной активности.

§2.2 Нейросетевой классификационный подход.

§2.3 Задачи классификации геомагнитных возмущений на основе солнечно-земных связей.

§2.4 Принципы отбора глобальных геомагнитных вариаций и поиск адекватного числа их классов.

§2.5 Определение космических источников геомагнитных возмущений.

§2.6 Классификация резких скачков в параметрах межпланетной среды.

§2.7 Анализ влияния электрического поля солнечного ветра на геомагнитную активность.

§2.8 Выводы.

ГЛАВА 3. Особенности высокоширотных геомагнитных возмущений

§3.1 Необходимость контроля пространственной динамики полярной шапки и аврорального овала.

§3.2 Разработка методики классификации пространственного положения станций по результатам Фурье-анализа геомагнитных данных.

§3.3 Разработка методики классификации пространственного положения станций по результатам вейвлет-анализа геомагнитных данных.

§3.4 Обработка результатов классификации и оценка уровня согласованности классификационного и традиционного модельного подходов.

§3.5 Поиск физических причин разделения спектральных характеристик магнитных возмущений в полярных областях.

§3.6 Выводы.

ГЛАВА 4. Исследование связи максимально наблюдаемой и критической частоты слоя F2 среднеширотной ионосферы с геоэффективными параметрами межпланетной среды

§4.1 Корреляционное исследование связи критической частоты слоя F2 и МНЧ с ключевыми геоэффективными параметрами межпланетной среды.

§4.2 Определение времени реакции критической частоты слоя F2 ионосферы и МНЧ на взаимодействие межпланетной среды с магнитосферой.

§4.3 Поиск оптимальной нейросетевой архитектуры для прогноза ионосферных параметров.

§4.4 Определение дополнительных входных параметров нейросети для повышения качества прогнозирования критической частоты слоя F2 и МНЧ.

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

§4.6 Выводы.

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

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

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

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

Методы, используемые в исследовании

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

Цель работы

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

Научная новизна

Применение ИНС для классификации Dst-вариаций геомагнитного поля помогло выявить характерные особенности в реакции геомагнитного поля на приход к Земле солнечных возмущений конкретной природы: вспышки, корональные дыры, активизировавшиеся волокна (протуберанцы), корональные выбросы вещества и их сочетания. Анализ полученных результатов позволил отождествить каждый выделенный класс с конкретным типом возмущений параметров солнечного ветра и межпланетного магнитного поля.

На основе одножидкостного магнитогидродинамического подхода с применением технологии искусственных нейронных сетей разработана методика автоматического разделения скачков в регистрируемых параметрах космической плазмы и магнитного поля на классы, отвечающие известным типам разрывов. Для этого была разработана классификационная ИНС вида «слой Кохонена» и выполнена классификация скачков параметров зарегистрированных на КА WIND в 1996-1999 гг. по двум разным алгоритмам. Для разрывов установленных классов определены ориентации плоскостей их поверхностей в трехмерном пространстве на основе одномерных наблюдений. Также проведено специальное исследование эволюционности найденных ударных волн.

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

Разработан алгоритм прогнозирования на интервал 0,5-3 часа максимально наблюдаемой и критической частоты ионосферного слоя F2 на основе технологии искусственных нейронных сетей. Методика прогнозирования применена к данным наклонного зондирования, полученным на трассе Inskip (Англия) — Ростов-на-Дону. Выполненные исследования охватили временной интервал, содержащий различные геофизические условия. Определены характерные времена ионосферной реакции на изменения ключевых параметров солнечного ветра и межпланетного магнитного поля, определяющих последовательность развития магнитосферно-ионосферных возмущений в свою очередь приводящих к изменениям коротковолновой (KB) радиосвязи.

Научная и практическая ценность

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

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

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

Экспериментальные данные

В диссертационном исследовании использовались следующие экспериментальные данные:

1. Параметры солнечного ветра и межпланетного магнитного поля, полученные из базы спутниковой системы OMNI [http://www.ngdc.noaa.gov]

2. Данные по индексам глобальной геомагнитной активности [http://www.ngdc.noaa.gov; http://cdaweb.gsfc.nasa.gov/istppublic; http://swdcwwvv.kugi.kyoto-u.ac.jp/wdc/Sec3.html]

3. Данные по состоянию ионосферного слоя F2, предоставленные НИРФИ

4. Данные о состоянии межпланетной среды, полученные с КА Wind и Асе [http://www.srl.caltech.edu/ACE/ASC/rtsw.html; http://cdaweb.gsfc.nasa.gov/istpj3ublic]

5. Данные о состоянии межпланетной среды, полученные из базы SPIDR [http://clvistl.\vdcb.ruyspidr/index.jsp]

6. Данные по разрывам параметров межпланетной среды, предоставленные ИКИ РАН.

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

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

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

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

Краткое содержание глав

Во введении сформулированы проблемы и изложены основные задачи, решаемые в диссертационном исследовании.

В первой главе диссертации представлен обзор крупномасштабных явлений солнечной активности, влияние динамики параметров межпланетной среды на магнитосферно-ионосферную активность. Рассмотрены геоэффективные проявления таких событий. Уделено внимание современным подходам к решению геофизических задач. В первом параграфе приведено описание наиболее активных и, как следствие, наиболее геоэффективных последствий активности Солнца. Представлены основные свойства таких активных событий, как вспышки, волокна, корональные дыры, корональные выбросы вещества (КВВ). Второй параграф посвящен описанию геомагнитной активности, существующим индексам, характеризующим глобальное геомагнитное состояние. Рассмотрены особенности развития геомагнитных возмущений среднеширотной и высокоширотной зон. Представлен физический механизм развития геомагнитной активности. В третьем параграфе приводятся различные современные подходы к решению геофизических задач. Среди них: построение трехмерных картин геоэффективных структур в межпланетном пространстве с привлечением «multi-spacecraft method», применение вейвлет-анализа - масштабируемого амплитудно-временного разложения геомагнитных данных, использование нелинейного корреляционного инструмента — искусственных нейронных сетей. В четвертом параграфе ставятся задачи по исследованию причин геомагнитной активности с привлечением современных методик решения.

Вторая глава посвящена исследованию возможности создания комплексной классификации глобальных геомагнитных возмущений с привлечением данных о параметрах околоземного космического пространства. Первый параграф посвящен описанию вариаций глобальных геомагнитных возмущений, описываемых индексом Dst. Второй параграф описывает инструмент классификации - приводятся архитектура нейросети, описан алгоритм обучения и работы. В третьем параграфе, при учете физики механизмов развития геомагнитных бурь и причинно-следственных связей, ставится задача классификации геомагнитных возмущений. В качестве инструмента классификации используется технология самообучающейся нейросети. В четвертом параграфе приводятся выкладки по отбору Dst-вариаций и сопутствующих данных о параметрах межпланетной среды, участвующих в классификации. Экспериментально устанавливается наиболее адекватное для данной задачи число классов событий. В пятом параграфе делается попытка отождествить полученные классы вариаций геомагнитного поля с вызвавшими их солнечными источниками. В шестом параграфе с привлечением отработанной технологии ИНС-классификации приводятся результаты исследования по разделению скачков на классы, отвечающие известным типам магнитогидродинамических разрывов. Предложен алгоритм установления ориентации поверхностей разрывов по одномерным наблюдениям скачков параметров солнечного ветра на космических аппаратах. Седьмой параграф содержит результаты исследования по поиску оценки связи Кр индекса геомагнитной активности с электрическим полем солнечного ветра Есв и с проекцией этого поля на геомагнитный диполь. Показано, что основным геоэффективным параметром, как это уже неоднократно отмечено многими исследователями, является южная (Bz<0) компонента вектора межпланетного магнитного поля. Именно эта компонента поля играет главную роль и в наличии связи Кр индекса с взаимной ориентацией Есв и магнитного момента Земли. Таким образом, демонстрируется возможность объединения всех параметров характеризующих как причину, т.е. структуру солнечного ветра, так и следствие, т.е. глобальное геомагнитное возмущение.

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

Заключительная, четвертая, глава посвящена исследованию среднеширотной ионосферы. Рассмотрена связь максимально наблюдаемой и критической частоты слоя F2 с параметрами околоземного космического пространства. В первом параграфе приводятся результаты базового корреляционного исследования связи частот с ключевыми геоэффективными параметрами межпланетной среды. Это исследование вызвано необходимостью ограничения числа солнечно-магнитосферно-ионосферных связей для выполнения успешного прогноза частоты методом ИНС. Установлено, что разные параметры по-разному отражаются в профиле частот в возмущенные и спокойные интервалы. Во втором параграфе представлены результаты поиска оптимального времени задержки для параметров ММП и ПСВ. Данное исследование проводилось на разветвленной ИНС и показало, что для нейросети является предпочтительной задержка в 2,5-3,5 часа. Обучение и тестирование при этом происходит более корректно и эффективно. Этот результат согласуется с предыдущим исследованием корреляционной связи параметров. Третий параграф главы посвящен поиску оптимальной архитектуры ИНС, на которой в дальнейшем проводиться основная часть численных экспериментов. Описаны особенности процесса обучения нейросети, разрабатывается методика прогнозирования частот. В четвертом параграфе продемонстрирована методика поиска дополнительных входных параметров, а так же комбинаций параметров для повышения качества прогнозирования. Выясняется вопрос о необходимости разделения технологии прогнозирования МНЧ для спокойных и возмущенных интервалов. Пятый параграф посвящен вопросу повышения качества восстановления профиля частот при добавлении их производных высокого порядка в число входных параметров при обучении и дальнейшей работе нейросети. Показано, что на разветвленной ИНС можно проводить прогнозирование МНЧ и критической частоты на интервалы 0,5-3 часа с точностью не ниже 85%.

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

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

§4.6 Выводы

Данная часть диссертационного исследования посвящена разработке алгоритма прогнозирования на интервал 0,5-3 часа максимально применимой частоты и критической частоты ионосферного слоя F2 на основе технологии искусственных нейронных сетей. Методика прогнозирования применена к данным наклонного JT4M зондирования. полученным на трассах Inskip (Англия) - Ростов-на-Дону и Inskip - Нижний Новгород. Выполненные исследования охватили временной интервал, содержащий различные геофизические условия. Определены характерные времена ионосферной реакции на изменения ключевых параметров солнечного ветра и межпланетного магнитного поля, определяющих последовательность развития магнитосферно-ионосферных возмущений в свою очередь приводящих к изменениям KB радиосвязи.

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

Разветвленные сети работают более стабильно — при перезапуске обучения такой сети качество ее обучения меняется менее заметно и, как следствие, качество прогноза от опыта к опыту практически не изменяется, а неудовлетворительные результаты появляются крайне редко. Обучение подобных мощных сетей требует страховки от переобучения — механизма контрольной проверки. Это несколько увеличивает время обучения, но дает высокое качество прогноза. Другим достоинством разветвленной сети является применимость к ее ответному сигналу методики коррекции. Максимальные результаты (до 99%) при прогнозировании частот на 0,5-1 час достигнуты на такой разветвленной ИНС. Та же сеть выполняет прогноз на 2 часа с точностью 91% и выше. Прогноз на 2 часа выполняется с 88%.

Качество прогноза на возмущенном интервале при обучении сети только по предыстории низкое, так как профиль частот МНЧ и foF2 сильно изрезан. Необходимо дополнить поступающую на вход сети информацию сопутствующими процессами, такими как сведения о текущей глобальной геомагнитной обстановке, концентрации солнечного ветра. Для возмущенного интервала введение в обучающую последовательность данных о частицах различных энергий способствовало улучшению восстановления профиля МНЧ и повышению точности трехчасового прогноза на ДРЕ^5%. Однако, этого не достаточно для разделения технологии прогнозирования на спокойных и возмущенных интервалах. Для прогнозирования на спокойных и возмущенных интервалах на 0,5-3 часа следует придерживаться единой схемы прогноза, так как обучающая последовательность обычно уже включает и спокойные и возмущенные интервалы. Обученная в таких условиях сеть имеет равноценные возможности для прогноза интервалов МНЧ и foF2 любой степени возмущенности.

Достоинством метода ИНС при восстановлении временных рядов является их применимость в целях прогноза в режиме «on line». Получая непрерывный поток необходимых данных (например, через Internet) сеть способна с определенной дискретностью создавать прогнозируемый ряд и при необходимости корректировать свой уровень обучения. Но, так как обучаемость сети в одних и тех же условиях может быть разной, то для достоверного прогноза следует пользоваться несколькими (3-5) параллельно работающими одинаковыми сетями (группой). Достоверный ответ тот, который с малой расходимостью получен на нескольких сетях из группы. Остальные сети, дающие другой, существенно отличный результат следует переобучить (с учетом новых данных) и снова включить в группу. Такая стратегия обеспечит максимально достоверный прогноз параметра в режиме реального времени.

ЗАКЛЮЧЕНИЕ

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

Перечислим важнейшие результаты, пути их получения и положения, выносимые на защиту:

1. Выполнена классификация скачков параметров солнечного ветра по данным КА WIND за трехлетний период их измерений (1996-1999 г.) используя одножидкостный МГД подход с помощью специально созданной классификационной искусственной нейронной сети вида «слой Кохонена». Для выделенных классов скачков разработан и применен способ определения ориентации плоскостей их поверхностей в солнечном ветре по одномерным наблюдениям на КА. Выполнен анализ эволюционности найденных ударных волн, который продемонстрировал их потенциальную неустойчивость.

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

Класс №1 — слабые бури, являющиеся следствием взаимодействия с магнитосферой Земли коротирующих с Солнцем потоков (corotating-streams).

Класс №2 — слабые бури, являющиеся следствием взаимодействия с магнитосферой Земли медленных гелиосферных потоков с низкой концентрацией плазмы.

Класс №3 - умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли коротирующих потоков, которые до этого взаимодействовали с HCS-потоком или гелиосферным токовым слоем.

Класс №4 — умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли быстрых, высокоинтенсивных corotating-потоков.

Класс №5 - умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли высокоскоростного потока, прошедшего мимо Земли и задевшего ее своим флангом.

Класс №6 - умеренные бури, являющиеся следствием взаимодействия с магнитосферой Земли сразу нескольких гелиосферных потоков.

Класс №7 - сильные бури, являющиеся следствием взаимодействия с магнитосферой Земли уединенных, не имеющих периодической повторяемости, потоков с выраженными ударными волнами.

Класс №8 - сильные бури, являющиеся следствием взаимодействия с магнитосферой Земли корональных солнечных выбросов в форме магнитного облака.

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

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

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

Сопоставление классификационных результатов на основе ИНС с аналитическими вычислениями положения южной границы овала на базе эмпирических моделей, основанных на корреляционной связи спутниковых данных с индексами геомагнитной активности AL, Dst [76], показывает, что они согласуются (совпадение 60% и выше). Метод определения динамики пространственного положения полярной шапки и овала полярных сияний по данным Гренландской цепочки магнитных станций можно применять к данным других подобных цепочек геомагнитных измерений. При этом возможна оценка границ полярной шапки и аврорального овала в реальном масштабе времени. Определены физические причины разделения спектральных характеристик магнитных возмущений в полярных областях

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

3. Реализовано прогнозирование на интервал 0,5-3 часа максимально наблюдаемой и критической частоты ионосферного слоя F2 на основе метода искусственных нейронных сетей с использованием параметров солнечного ветра и межпланетного магнитного поля. Определены характерные времена реакции ионосферы на изменения ключевых параметров солнечного ветра и межпланетного магнитного поля.

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

Разветвленные сети работают более стабильно - при перезапуске обучения такой сети качество ее обучения меняется менее заметно и, как следствие, качество прогноза от опыта к опыту практически не изменяется, а неудовлетворительные результаты появляются крайне редко. Обучение подобных мощных сетей требует страховки от переобучения — механизма контрольной проверки. Это несколько увеличивает время обучения, но дает высокое качество прогноза. Другим достоинством разветвленной сети является применимость к ее ответному сигналу методики коррекции. Максимальные результаты (до 99%) при прогнозировании частот на 0,5-1,5 часа достигнуты на такой разветвленной ИНС. Та же сеть выполняет прогноз на 2 часа с точностью 91% и выше. Прогноз на 3 часа выполняется с эффективностью 88%.

Установлено, что качество прогноза на возмущенном интервале при обучении сети только по предыстории низкое, так как профиль частот МНЧ и foF2 сильно изрезан. Необходимо дополнить поступающую на вход сети информацию сопутствующими процессами, такими как сведения о текущей глобальной геомагнитной обстановке, концентрации солнечного ветра. Для прогнозирования на спокойных и возмущенных интервалах на 0,5-3 часа следует придерживаться единой схемы прогноза, так как обучающая последовательность обычно уже включает и спокойные и возмущенные интервалы. Обученная в таких условиях сеть имеет равноценные возможности для прогноза интервалов МНЧ и f0F2 любой степени возмущенности.

Предложена групповая нейросетевая методика восстановления временных рядов в режиме «on line». Получая непрерывный поток необходимых данных через Internet сеть способна с определенной дискретностью создавать прогнозируемый ряд и при необходимости корректировать свой уровень обучения. Для достоверного прогноза следует пользоваться несколькими (3-5) параллельно работающими одинаковыми сетями (группой). Достоверный ответ тот, который с малой расходимостью получен на нескольких сетях из группы. Остальные сети, дающие другой, существенно отличный результат следует переобучить (с учетом новых данных) и снова включить в группу. Такая стратегия обеспечит максимально достоверный прогноз параметра в режиме реального времени.

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

Степень достоверности полученных результатов

Результаты, изложенные в диссертационной работе, согласуются с соответствующими данными экспериментальных наблюдений. Их сопоставление с аналитическими и численными исследованиями других авторов показали, что полученные результаты дополняют и уточняют эти исследования. Все результаты, представленные в диссертации опубликованы в рецензируемых научных журналах «Космические исследования», «Геомагнетизм и аэрономия», «Известия РАН — Физическая серия», «Известия ВУЗов —

Радиофизика», а также представлены на Российских и международных конференциях и научных семинарах НИРФИ и ИЗМИР АН.

Личный вклад соискателя

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

Благодарности

При получении результатов, представленных в диссертации, я опирался на знания и профессиональный опыт моего учителя и научного руководителя доктора физико-математических наук, профессора Николая Александровича Бархатова; сотрудников лаборатории геомагнитных вариаций ИЗМИР АН Анатолия Ефимовича Левитина, Людмилы Ивановны Громовой, Лидии Алексеевны Дремухиной; сотрудника ИКИ РАН Георгия Наумовича Застенкера; сотрудника НИИРФИ Валерия Павловича Урядова. Большую помощь и поддержку мне оказали также сотрудники Научно-исследовательской лаборатории Физики Солнечно-Земных связей Нижегородского Государственного Педагогического Университета. Всем нм выражаю свою признательность и благодарность.

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

1. Бархатов Н.А., Валов В.А., Макаров А.В., Понятов А.А., Ревунов С.Е., Сахаров С.Ю., Урядов В.П., Метод искусственных нейронных сетей и его применение к прогнозированию ионосферного KB радиоканала. Труды IX Международной научно-технической конференции «Радиолокация, навигация и связь», 22-24 апреля 2003 г., Воронеж, Том III, С. 1853-1864.

2. Barkhatov N.A., Ponyatov А.А., Uryadov V.P., Barkhatova O.A., Revunov S.E., Forecasting of ionospheric hf channel by the artificial neural networks method // Abstracts of IUGG XXIII General Assembly, Sapporo, Juhe 30-July 11, 2003, GAII.03/02A/AI3-005, C.A298

3. Бархатов H.A., Ревунов C.E., Прогнозирование критической частоты ионосферного слоя F2 методом искусственных нейронных сетей. VIII Нижегородская сессия молодых ученых. Естественнонаучные дисциплины, 20-25 апреля 2003 г., Дзержинск, Тезисы докладов, Н. Новгород, 2003, С.98-99

4. Barkhatov N.A., Levitin A.E., Revunov S.E. Comprehensive classification of geomagnetic disturbances. Geophysical Research Abstracts, V.6, 04796, European Geosciences Union, General Assembly 2004, ST5, Nice, France, 25-30 April 2004

5. Ревунов С.E., Дудкина И.С. Повышение эффективности прогнозирования критической частоты ионосферного слоя F2. Сборник тезисов Десятой Всероссийской научной конференции студентов-физиков и молодых ученых: Тезисы докладов: В 2 т. - Москва, издательство АСФ России, 2004, Т.2. С. 915

6. Barkhatov N.A., Revunov S.E., Uryadov V.P. Improvement of ionospheric layer f2 critical frequency forecasting. Geophysical Research Abstracts, V.6, 04789, European Geosciences Union, General Assembly 2004, ST 14, Nice, France, 25-30 April 2004

7. Бархатов H.A., Левитин A.E., Ревунов C.E. Космическая классификация глобальных геомагнитных возмущений. Тезисы докладов III Международной конференции «Солнечно-земные связи и электромагнитные предвестники землетрясений», 16-21 августа 2004, Институт космофизических исследований и распространения радиоволн (ИКИР) ДВО РАН (с. Паратунка, Камчатская обл.), С.46

8. Бархатов Н.А., Левитин А.Е., Ревунов С.Е. Классификация глобальных геомагнитных возмущений на основе комплекса событий космической погоды. Тезисы докладов международной конференции «Солнечно-земная физика», Иркутск, 20-25 сентября 2004 г., С.46

9. Barkhatov N.A., Revunov S.E., Uryadov V.P. (2004), Forecasting of the critical frequency of the ionosphere F2 layer by the method of artificial neural networks. Int. J. Geomagn. Aeron., 5, GI2010, D01:10.1029/2004GI000065

10. Бархатов H.A., Ревунов C.E., Урядов В.П. Технология искусственных нейронных сетей для прогнозирования критической частоты ионосферного слоя F2. // Известия ВУЗов «Радиофизика», 2005, Т.48. С.1-15

11. Бархатов Н.А., Ревунов С.Е., Урядов В.П., Вертоградов Г.Г., Вертоградов В.Г., Валов В.А. Прогнозирование МНЧ с помощью искусственных нейронных сетей. Материалы 21 Всероссийской конференции по распространению радиоволн, Йошкар-Ола, 24-27 мая 2005 г.

12. Barkhatov N., Levitin A., Revunov S. Location of high latitude zones of geomagnetic activity by the data from the Greenland chain of magnetic stations. // Abstracts of 10-th IAGA Scientific Assembly, Toulose-Franee, July 18-29, 2005, IAGA2005-A-01244

13. Barkhatov N.A., Revunov S.E., Ulibina R.I. Classification of discontinuities in parameters of space plasma using magnetohydrodynamic approach. Abstracts of 9th Asian-Pacific Regional International Astronomical Union (IAU) Meeting (APR1M2005), Bali, Indonesia, 26-29 July 2005, P. 17

14. Barkhatov N.A., Gromova L.I., Dremuhina L.A, Levitin A.E., Revunov S.E., Ulibina R.I. The analysis of magnetospheric disturbance power according to its solar source. Abstracts of 9th Asian-Pacific Regional International Astronomical Union (IAU) Meeting (APR1M2005), Bali, Indonesia, 26-29 July 2005. P. 18

15. Бархатов H.A., Левитин A.E., Ревунов C.E. Установление границ областей полярной шапки и овала полярных сияний по данным Гренландской цепочки магнитных станций. Тезисы Всероссийской конференции «Экспериментальные и теоретические исследования основ прогнозирования гелиогеофизической активности» (10-15 октября 2005 г., г. Троицк), С.7

16. Бархатов Н.А., Громова Л.И., Дремухина Л.А., Левитин А.Е., Ревунов С.Е., Улыбина Р.И. Анализ энергетики магнитосферного возмущения в зависимости от его солнечного источника. Тезисы Всероссийской конференции «Экспериментальные и теоретические исследования основ прогнозирования гелиогеофизической активности» (10-15 октября 2005 г., г. Троицк), С.44

17. Бархатов Н.А., Ревунов С.Е., Урядов В.П., Вертоградов Г.Г, Вертоградов В.Г. Учет проявлений солнечной активности при прогнозировании максимально наблюдаемой частоты ионосферного KB радиоканала. Тезисы Всероссийской конференции «Экспериментальные и теоретические исследования основ прогнозирования гелиогеофизической активности» (10-15 октября 2005 г., г. Троицк), С.59

18. Barkhatov N.A., Revunov S.E., Uryadov V.P., Vertogradov G.G., Vertogradov V.G. Forecasting of hf propagation with use of artificial neural networks. XXVIIIth URSI General Assembly New Delhi, India, 23 to 29 October 2005

19. Бархатов H.A., Левитин A.E., Ревунов C.E. Классификация глобальных геомагнитных возмущений на основе комплекса событии космической погоды. В сб. «Солнечно-земная физика», вып.8, Изд. СО РАН, Новосибирск, 2005, С. 129-131

20. Ревунов С.Е., Улыбина Р.И. Нейросетевая классификация разрывов параметров космической плазмы, XI Нижегородская сессия молодых ученых, Математические науки, 22-25 мая 2006 г., Нижний Новгород, Тезисы докладов, 2006, С.36-37

21. Barkhatov N.A., Levitin А.Е., Revunov S.E. Monitoring of polar cap spatial dynamics basing on data from Greenland chain of magnetic observatories. Abstracts of 29 Annual Seminar «Physics of auroral phenomena», Apatity, 27 February-3 March 2006, PGI KSC RAS, P.31

22. Barkhatov N.A., Revunov S.E. Classification of abrupt changes of the interplanetary medium parameters by a method of artificial neural networks. Abstracts of 29 Annual Seminar «Physics of auroral phenomena», Apatity, 27 February-3 March 2006, PGI KSC RAS, P.51

23. Бархатов H.A., Ревунов C.E., Вертоградов Г.Г, Вертоградов В.Г, Валов В.А., Урядов В.П. Прогнозирование максимально наблюдаемой частоты ионосферного КВ-радиоканала методом искусственных нейронных сетей. // Геомагнетизм и аэрономия, 2006, Т.46, N 1, С. 88-98

24. Barkhatov N.A., Dremukhina L.A., Gromova L.I., Levitin A.E., Revunov S.E., Ulibina R.I. «Magnetosphere power budget role in the task of classification of magnetospheric activity sources» 36th Cospar scientific assembly, Beijing, China, 16-23 July 2006. Session D2.1/E3.1, A-00873

25. Бархатов H.A., Громова Л.И., Дремухина Л.А., Левитин А.Е., Ревунов С.Е., Улыбина Р.И. Учет энергетического бюджета магнитосферы в задаче классификации источников магнитосферной активности. Изв.РАН (сер.физ), 2006, 70. N10, С. 1535-1537

26. Barkhatov N.A., Levitin А.Е., Revunov S.E. Complex Classification of Global Geomagnetic Disturbances. // Kosmicheskie Issledovaniya, 2006, Vol. 44, No. 6, pp. 488-499

27. Barkhatov N.A., Levitin A.E., Revunov S.E. Monitoring of polar cap spatial dynamics basing on data from Greenland chain of magnetic observatories. Proceedings of 29 Annual Seminar «Physics of auroral phenomena», Apatity, 2007, PGI KSC RAS, P.75

28. Barkhatov N.A., Revunov S.E. Classification of abrupt changes of the interplanetary medium parameters by a method of artificial neural networks. Proceedings of 29 Annual Seminar «Physics of auroral phenomena», Apatity, 2007, PGI KSC RAS, P. 175

29. Бархатов H.A., Левитин A.E., Ревунов C.E. Контроль пространственной динамики полярной шапки и полярного овала по данным гренландской сеги магнитных станций. // Геомагнетизм и аэрономия, 2007, Т.47, N2, С. 230-235

30. Бархатов Н.А., Громова Л.И., Левитин А.Е., Ревунов С.Е. Анализ влияния электрического поля солнечного ветра на геомагнитную активность. Тезисы докладов 11 Пулковской международной конференции по физике Солнца: Физическая природа солнечной активности и прогнозирование ее геофизических проявлений. ГАО РАН, Пулково, Санкт-Петербург, 2-7 июля 2007, С. 17

31. Бархатов Н.А., Левитин А.Е., Ревунов С.Е., Смирнова А.С., Снегирев С.Д., Хвиюзова Т.А. Классификация событий космической погоды. Тезисы докладов 11 Пулковской международной конференции по физике Солнца: Физическая природа солнечной активности и прогнозирование ее геофизических проявлений. ГАО РАН, Пулково, Санкт-Петербург, 2-7 июля 2007, С. 19

32. Barkhatova О.М., Revunov S.E., Barkhatov N.A., Levitin A.E. Nonlinear connection of global and polar current systems on the main phase of a geomagnetic storm. Abstracts of 31 Annual Seminar «Physics of auroral phenomena», Apatity, 26-29 February 2008, PGI KSC RAS, P. 12

33. Barkhatov N., Levitin A., Revunov S. Connection of geomagnetic activity with orientation of the solar wind electric field concerning the terrestrial magnetic moment. Geophysical Research Abstracts, Vol. 10, EGU2008-A-00000, 2008 EGU General Assembly 2008

Библиография Диссертация по наукам о земле, кандидата физико-математических наук, Ревунов, Сергей Евгеньевич, Нижний Новгород

1. Витинский Ю. И. Солнечная активность. -М.: Наука, 1969. - С. 192.

2. Харгривс Д.К. Верхняя атмосфера и солнечно-земные связи. Введение в физику околоземной космической среды. Л.: Гидрометеоиздат., 1982.

3. Хундхаузен А. Расширение короны и солнечный ветер. М., Изд. Мир, 1976

4. Смольков Г.Я., Фомичев В.В., Снегирев С.Д. Радиоизлучение Солнца и космическая погода // Солнечно-земная физика. 2002. - Вып. 2. - С. 31.

5. Гибсон Э. Спокойное Солнце / Пер. с англ. М.: Мир, 1977.

6. Акасофу С.И., Чепмен С. Солнечно-земная физика. 4.1,2 М.: Мир, 1974. - С. 444.

7. Витинский Ю.И., Копецкий М., Куклин Г.В. Статистика пятно-образовательной деятельности Солнца —М.: Наука, 1986.

8. Ишков В. Н. Солнечные геоэффективные явления: как и когда они воздействуют на околоземное космическое пространство // Солнечно-земная физика, Иркутск, ИСЗФ, 2002.-Вып2.-С.10.

9. Куликовский А.Г., Любимов Г.А. Магнитная гидродинамика. М.: Физматгиз, 1962, 246 с.

10. Баранов В.Б., Краснобаев К.В. Гидродинамическая теория космической плазмы. М.: Наука, 1977, 336 с.

11. Пудовкин М.И., Семенов B.C. Теория пересоединения и взаимодействие солнечного ветра с магнитосферой Земли. М. Наука, 1985.

12. Брюнелли Б.Е., Намгаладзе А.А. Физика ионосферы. М.Наука. 1988. 527с.

13. Пудовкин М.И., Распопов О.М., Клейменова Н.Т. Возмущения электромагнитного поля Земли. Л.: Изд-во ЛГУ, 1976. С. 247.

14. Sonnerup B.U.O, Н. Hasegawa Orientation and motion of two-dimensional structures in a space plasma, J. Geophys. Res., 2005, Vol. 110, N A06208, doi:10.1029/2004JA010853.

15. Milan S.E., Lester M., Cowley S.W.H., Oksavik K., Brittnacher M., Greenwald R.A., Sofko G., Villain J.-P. Variations in the polar cap area during two substorm cycles // Ann. Geophys. Vol. 21. P. 1121.2003

16. Blanchard G.T., Lyons L.R., Samson J.C., Rich F.J. Locating polar cap boundaiy from observations of 6300A auroral emission //J. Geoph. Res. Vol. 100. No. A5. P. 7855. 1995

17. Newell P.T., Liou Kan, Sotirelis Т., Meng Ching-I. Polar Ultraviolet Imager observations of global auroral power as a function of polar cap size and magnetotail stretching // J. Geoph. Res. Vol. 106. No. A4. P. 5895. 2001

18. Gerard, J.-C., Hubert В., Grard A., Meurant M., and Mende S. B. Solar wind control of auroral substorm onset locations observed with the IMAGE-FUV imagers, J. Geoph. Res., 2004, V. 109, A03208, doi: 10.1029/2003JA010129

19. Misiti M., Misiti Y., Oppenheim G., Poggi J.-M. Wavelet Toolbox for use with MATLAB. The Math Works Inc. 941 p. 2000

20. Дьяконов В.П., Круглов В.В. Математические пакеты расширения MATLAB. Специальный справочник. СПб.: Питер, 2001, 480с.

21. Дьяконов В.П. MATLAB 6. Учебный курс. СПб.: Питер. 593 с. 2001

22. Круглов В.В., Борисов В.В. Искусственные нейронные сети. Теория и практика. М.: Горячая линия-Телеком. 384 с. 2000

23. Медведев B.C., Потёмкин В.Г. Нейронные сети. MATLAB 6. М.: Диалог-МИФИ, 2002, 496с.

24. Бархатов Н.А., Королев А.В., Пономарев С.М., Сахаров С.Ю. Долгосрочное прогнозирование индексов солнечной активности методом искусственных нейронных сетей. Известия ВУЗов "Радиофизика", 2001, Т.44, № 9, С.806

25. Бархатов Н.А., Беллюстин Н.С., Левитин А.Е., Сахаров С.Ю. Сравнение эффективности предсказания индекса геомагнитной активности Dst искусственными нейронными сетями. // Изв. ВУЗов "Радиофизика", 2000, Т.43, № 5, С.385

26. Бархатов Н.А., Левитин А.Е., Сахаров С.Ю. Метод искусственных нейронных сетей как способ восстановления пробелов в записях отдельных магнитных обсерваторий по данным других станций. // Геомагнетизм и аэрономия, 2002, Т.42, С. 195

27. Munsami, V. Determination of the effects of substorms on the storm-time ring current using neural networks. J. Geophys. Res., 2000, Vol. 105, №. A12, P.27833

28. Gleisner H., Lundstedt H. A neural network-based local model for prediction of geomagnetic disturbances. J. Geophys. Res., 2001, Vol. 106, № A5, P.8425

29. Lundstedt H., Wintoft P. Prediction of geomagnetic storms from solar wind data with the use of a neural network. Ann. Geophys., 1994, Vol. 12, № 1, P. 19

30. Gleisner H., Lundstedt H., Wintoft P. Predicting geomagnetic storms from solar-wind data using time-delay neural networks. Ann Geophys., 1996, Vol. 14, № 7, P.679

31. Sutcliffe P.R. Substorm onset identification using neural networks and Pi2 pulsations. Ann Geophys., 1997, Vol. 15, № 10, P.1257

32. Sutcliffe P.R. The development of a regional geomagnetic daily variation model using neural networks. Ann Geophys., 2000, Vol. 18, № 1, P. 120

33. Conway, A.J.; Macpherson, K.P.; Blacklaw, G.; Brown, J.C. 1998 A neural network prediction of solar cycle 23. J. Geophys. Res. Vol. 103, № A12, P.29, 733

34. Fessant F., Bengio S., Collobert D. On the prediction of solar activity using different neural network models. Ann Geophysic., 1996, Vol. 14, № 1, P.20

35. Wintoft, P.; Lundstedt, H. A neural network study of the mapping from solar magnetic fields to the daily average solar wind velocity. J. Geophys. Res. 1999, Vol. 104, № A4, P.6729

36. Gonzalez W.D., Tsurutani B.T. Criteria of interplanetary parameters causing intense magnetic storms (Dst<-100 nT), Planetary and Space Science, 1987, Vol. 36, P. 1101

37. Echer E. and W.D. Gonzalez. Geoeflectiveness of interplanetary shocks, magnetic clouds, sector boundary crossings and their combined occurrence, Geophys. Res. Lett., 2004, 31, L09808, 10.1029/2003 GL019199

38. Watari S., M. Vandas and T. Watanabe, Formation of a strong southward IMF near the solar maximum of cycle 23, Annales Geophysicae, 2004, Vol. 22, P. 673

39. Иванов К.Г. Солнечные источники потоков межпланетной плазмы на орбите Земли. // Геомагнетизм и аэрономия, Т.36, № 2, С. 19, 1996.

40. Бархатов Н.А., Гольберг К.Ю., Зырянова М.С., Иванов К.Г. МГД-моделирование взаимодействия высокоскоростного потока и медленного солнечного ветра. // Геомагнетизм и аэрономия, 2004, Т.44, N 1, С. 28-36.

41. Gonzalez W.D., Joselyn J.A., Kamide Y., Kroehl H.W., Rostoker G., Tsurutani B.T., Vasyliunas V.M. What is geomagnetic storm? J. Geophys. Res., 1994, Vol. 99, № A4, P. 5771

42. Tsurutani B.T., Gonzalez W.D., Gonzalez Alicia, Tang Frances, Arballo J.K., Okada Masaki. Interplanetary origin of geomagnetic activity in the declining phase of the solar cycle. J. Geophys. Res., 1995, Vol. 100, № Al 1, P. 21717

43. Нишида А. Геомагнитный диагноз магнитосферы. M.: "Мир", 1980, 300 с.

44. Tu C.Y., Marsch Е. MHD structures, waves and turbulence in the solar wind: Observations and theories, Space Sci. Rev., 1995, V.73, P.l

45. Russell C.T., McPherron R.L. Semiannual variation of geomagnetic activity // J. Geophys. Res. 1973,V.78, N1, P.92-107.

46. Яновский Б.М. Земной магнетизм, 1978; Л.: Изд. Ленинградского университета, 591 С.

47. Lyatsky W., Tan A. Latitudinal effect in semiannual variation of geomagnetic activity // J.Geophys. Res., 2003, V.108.N5, 1204, doi:10.1029/2002JA009467.

48. Кузнецова T.B., Лаптухов А.И., Кузнецов В.Д. Учет геометрии взаимодействия солнечного ветра с магнитным полем Земли в геоэффективных параметрах и прогнозе геомагнитной активности // Астрономический вестник, 2006, Т. 40, №6, с. 561-568

49. Rice W.K. М., Zank J.P., Lee G. Particle acceleration and coronal mass ejection driven shocks: Shocks of arbitrary strength. J. Geophys. Res., 2003, Vol. 108, № A10, 1369, 10.1029/2002JA009756

50. Lee M.A., An analytical theory of morphology, flows. And shock compressions at corotating interaction regions in solar wind. J. Geophys. Res., 2000, Vol. 105, № A5, P. 10491

51. Бархатов H.A., Гольберг К.Ю., Зырянова M.C., Иванов К.Г. Локальное МГД моделирование взаимодействия высокоскоростного потока и медленного солнечного ветра. «Геомагнетизм и аэрономия», 2003, Т.43, № , С. 28

52. Физика магнитосферы. Под ред. Грингауза К.И., Жулина И.А., М.: Мир, 1972, 592 с.

53. Л.ДЛандау, Е.М.Лифшиц. Теоретическая физика: Гидродинамика. Том 6, М.: Наука, 1986, 736 с.

54. Ройхваргер З.Б. Эволюционность и ширина ударных волн, Письма в ЖЭТФ, 1974, Т. 19, Вып. 8, с. 496-499

55. Lindsay G.M., Russell С.Т., Luhmann J.G. Predictability of Dst index based upon solar wind conditions monitored inside 1 AU. J. Geophys. Res., 1999, Vol. 104, № A5, P. 10335

56. Temerin M., Li Xinlin. A new model for the prediction of Dst on the basis of the solar wind. J. Geophys. Res., 2002, Vol. 107, № A12, 1472, 10.1029/2001JA007532

57. Zhou Xiaoyan and Tsurutani Bruce T. Interplanetary shock triggering of nightside geomagnetic activity: Substorms, pseudobreakups, and quiescent events. J. Geophys. Res., 2001, Vol. 106, № A9, P. 18957

58. Emilia K., Huttunen J., Koskinen Hannu E.J., Schwenn Rainer. Variability of magnetospheric storms driven by different solar wind perturbations. J. Geophys. Res., 2002, Vol. 107, № A7, 1121, 10.1029/2001JA900171

59. Webb D.F., Cliver E.W., Crooker N.U., Cyr O.C. St., Thompson B.J. Relationship of halo coronal mass ejections, magnetic clouds, and magnetic storms. J. Geophys. Res., 2000, Vol. 105, №A4, P. 7491

60. Wu Chin-Chun and Lepping R.P. Effects of magnetic clouds on the occurrence of geomagnetic storms: The first 4 years of Wind. J. Geophys. Res., 2002, Vol. 107, № A10, 1314, 10.1029/2001JA000161

61. Nakai H., Kamide Y. Solar cycle variations in the storm-substorm relationship. J. Geophys. Res., 1999, Vol. 104, № A10, P. 22695

62. Vennerstroem S. Interplanetary sources of magnetic storms: A statistical study. J. Geophys. Res., 2001, Vol. 106, № A12, P. 29175

63. Lyatsky W. and Tan A. Solar wind disturbances responsible for geomagnetic storms. J. Geophys. Res., 2003, Vol. 108, № A3. 1134, 10.1029/2001 JA005057

64. Jurac S., Richardson J.D. The dependence of plasma and magnetic field correlations in the solar wind on geomagnetic activity. J. Geophys. Res., 2001, Vol. 106, № A12, P. 2919505

65. Wang Y.M., Ye P.Z., Wang S., Zhou G.P., Wang J.X. A statistical study on the geoeffectiveness of Earth-directed coronal mass ejections from March 1997 to December 2000. J. Geophys. Res., 2002, Vol. 107, № A11, 1340, 10.1029/2002JA009244

66. Leamon Robert J., Canfield Richard C., Pevtsov Alexei A. Properties of magnetic clouds and geomagnetic storms associated with eruption of coronal sigmoids. J. Geophys. Res., 2002, Vol. 107, №A9, 1234, 10.1029/2001JA000313

67. Stampler R., Lockwood V., Wild M.N., Clark T.D.G. Solar causes of the long-term increase in geomagnetic activity // J. Geophys. Res., 1999, V. 104, N12, P. 28325-28342.

68. Crooker N.U., Cliver E.W., Tsurutani B.T., The semiannual variation of great geomagnetic storms and the postshock Russell-McPherron effect preceding coronal mass ejecta // Geophys. Res. Lett., 1992, V. 19, P. 429.

69. Crooker N.U., Cliver E.W., A postmodern view of M-regions // J. Geophys. Res., 1994, V.99, P. 23383

70. Tsurutani B.T., Gonzalez W.D., Tang P., Lee Y.T., Great magnetic storms, // Geophys. Res. Lett., V. 19, N73, 1992.

71. Исаев С. И., Пудовкин М.И. Полярные сияния и процессы в магнитосфере Земли. Изд. Наука, Ленингр. отд., 1972, С. 244.

72. Воробьев В.Г, Ягодкина О.И. Влияние магнитной активности на глобальное распределение зон авроральных вторжений // Геомагнетизм и аэрономия. Т. 45. № 4. С.43 8-448. 2005

73. Иванов В.А., Куркин В.И., Носов В.Е., Урядов В.П., Шумаев В.В. ЛЧМ ионозонд и его применение в ионосферных исследованиях // Изв. ВУЗов "Радиофизика", 2003, Т.46, № 11, С.919.

74. Зевакина Р.А., Киселева М.В. О возможности количественного прогноза интенсивности отрицательных ионосферных возмущений. // Сб. Прогнозирование ионосферы и условий распространения радиоволн. М. Наука 1985, с.39-43.

75. Бархатов Н.А., Левитин А.Е., Рябкова Г.А. Искусственные нейронные сети для прогнозирования индексов геомагнитной активности по параметрам околоземного космического пространства. В сб. "Солнечно-земная физика", вып.2 (115), Иркутск,2002, С. 104

76. Baker D.N., Pulkkinen T.I., Hesse М., McPherron R.L. A quantitative assessment of energy storage and release in the Earth's magnetotail // J. Geophys. Res. Vol. 102. No. A4. P. 7159 (96JA03961). 1997

77. Burlaga, L.F., C. Wang, J.D. Richardson and N.F. Ness. Evolution of the multiscale statistical properties of corotating streams from 1 to 95 AU, J. Geophys. Res., 2003, Vol. 108, № A7, 1305, 10.1029/2003JA009841

78. Chin-Chun Wu and R.P. Lepping. Effect of solar wind velocity on magnetic cloud-associated magnetic stonn intensity J. Geophys. Res., 2002, Vol. 107, №. All, 1346, 10.1029/2002JA009396

79. Falconer D.A., R.L. Moore and G.A. Gary. A measure from line-of-sight magnetograms for prediction of coronal mass ejections, J. Geophys. Res., 2003. Vol. 108, № A10, 1380, 10.1029/2003JA010030

80. Gerard J.-C., Hubert В., Grard A., Meurant M., Mende S.B. Solar wind control of auroral substorm onset 13. P. 2165-2174. 1976

81. Huang, C.-S., G.D. Reeves, J.E. Borovsky, R.M. Skoug, Z.Y. Pu, and G. Le. Periodic magnetospheric substorms and their relationship with solar wind variations, J. Geophys. Res.,2003, Vol 108, №A6, 1255, 10.1029/2002JA009704

82. Kahler S.W. Hudson H.S. Origin and development of transient coronal holes, J. Geophys. Res., 2001, Vol. 100,№A12, P. 29239

83. Kan J.R., Lee L.C. Energy coupling function and solar wind-magnetosphere dynamo // Geophys. Res. Lett. Vol. 6. No 7. P. 577-580. 1979

84. Lam H.-L., Boteler D.H., Trichtchenko L. Case studies of space weather events from their launching on the Sun to their impacts on power systems on the Earth // Ann. Geophys. Vol. 20. P.1073-1079.2002

85. Luhmann J.G., У. Li, C.N. Arge, P.R. Gazis and R. Ulrich. Solar cycle changes in coronal holes and space weather cycles J. Geophys. Res., 2002, Vol. 107, № A8, 1154, 10.1029/2001JA00755092.