Бесплатный автореферат и диссертация по наукам о земле на тему

Ионосферные возмущения на различных фазах 23-го цикла солнечной активности по данным глобальной сети GPS

ВАК РФ 25.00.29, Физика атмосферы и гидросферы

Автореферат диссертации по теме "Ионосферные возмущения на различных фазах 23-го цикла солнечной активности по данным глобальной сети GPS"

Российская академия наук Дальневосточное отделение Институт космофизических исследований ираст

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

Живетьев Илья Валерьевич

Ионосферные возмущения на различных фазах 23-го цикла солнечной активности по данным глобальной сети GPS

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

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

Иркутск-2007

003069276

Диссертация выполнена в Институте космофизических исследований и распространения радиоволн Дальневосточного отделения Российской академии наук

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

Доктор физико-математических наук, профессор Афраймович Эдуард Леонтьевич

Официальные оппоненты

Доктор физико-математических наук, профессор Пономарев Евгений Александрович

Доктор физико-математических наук, профессор Попов Георгий Васильевич

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

Томский государственный университет

Защита состоится 22 мая 2007 г в 14 часов на заседании диссертационного совета Д 003 034 01 при Институте солнечно-земной физики СО РАН (664033, Иркутск, ул Лермонтова, 126, ИСЗФ СО РАН)

С диссертацией можно ознакомиться в библиотеке ИСЗФ СО РАН

Автореферат разослан "апреля 2007 г

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

Поляков В И

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

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

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

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

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

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

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

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

На протяжении многих лет огромные усилия исследователей предпринимаются для поисков научного обоснования построения системы предсказания ЗТ на основе мониторинга состояния геосфер Существует множество указаний на различные признаки готовящихся ЗТ в изменении локальных параметров нейтральной и заряженной компонент ионосферы, а также в акустическом и электромагнитом поле излучения эпицентров будущих ЗТ Считается, что на ускорение подготовки ЗТ могут влиять различные глобальные проявления динамики солнечной активности, обусловленные регулярными изменениями в течение солнечного цикла, а также вспышками, выбросами корональной массы, вызывающими магнитные и ионосферные бури [8, 9] Обсуждается гипотеза о том, что не только литосферные процессы могут быть источником ионосферных (атмосферных) возмущений, но и глобальные изменения в атмосфере сами влияют на характер протекания литосферных процессов, при этом атмосферные катаклизмы могут служить «спусковым» механизмом высвобождения сейсмической энергии [7] Эта гипотеза также требует надежных экспериментальных доказательств

В настоящее время широкое распространение получили исследования литосферно-ионосферных процессов на основе изучения возмущений полного электронного содержания (ПЭС) в ионосфере Земли, регистрируемых при зондировании ионосферы сигналами современных спутниковых навигационных систем - американской GPS и российской ГЛОНАСС

В некоторых работах [5, 12] отмечено отклонение от регулярного суточного хода ПЭС вблизи эпицентра ЗТ за несколько часов и дней до главного толчка Предпринимались неоднократные попытки обнаружить в вариациях ПЭС предсказанное теорией изменение спектра возмущений ионосферной плазмы, обусловленное генерацией в районе эпицентра акустико-гравитационных волн (АГВ) [10], которые регистрируются в форме перемещающихся ионосферных возмущений (ПИВ) Однако до сих пор не получено надежных свидетельств подготовки катастрофических землетрясений по вариациям ПЭС, как, впрочем, и по другим характеристикам сигналов радиофизических средств мониторинга сейсмоопасных районов

Для надежной идентификации признаков подготовки ЗТ необходимо хорошо знать глобальные и локальные характеристики «фоновых» возмущений ПЭС, обусловленных изменениями солнечной и геомагнитной активности Эта проблема составляет предмет настоящей диссертации

В диссертации изучение глобальных свойств ионосферы основано на концепции глобального электронного содержания (ГЭС), предложенной в 2005 г Э Л Афраймовичем Основное внимание в изучении локальных характеристик уделено Камчатскому региону Сейсмическая деятельность у Тихоокеанского

побережья Камчатки, Курильских островов и северо-восточной Японии достигает наивысшего на Земле уровня В районе Камчатского залива и Командорских островов расположена область стыка Курило-Камчатской и Алеутской островных дуг, которая относится к наиболее важным и интересным тектоническим узлам мира Сюда подходит и крупный подводный хребет Императорских гор В этом районе находится самое резкое пересечение глубоководных желобов в мире, а также один из магматических центров мира — Ключевская группа вулканов

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

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

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

1 Разработка программно-алгоритмического комплекса, предназначенного для определения глобального электронного содержания по данным международной сети GPS и его моделирования с использованием ионосферной модели IRI-2001

2 Формирование базы данных глобальной и региональной Камчатской сетей GPS за период 1996-2006 гг, необходимой для изучения динамики глобального электронного содержания и фоновых перемещающихся ионосферных возмущений

3 Исследование динамики глобального электронного содержания в течение 23-го цикла солнечной активности в сравнении с динамикой ультрафиолетового излучения и радиоизлучения Солнца

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

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

Научная новизна исследования:

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

анализе данных за 23-й цикл солнечной активности получены уточненные параметры связи динамики ионосферы, индекса солнечной активности F10 7 и потока УФ-излучения Солнца

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

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

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

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

Личный вклад автора

Основные результаты работы являются оригинальными и получены либо лично автором, либо при его непосредственном участии

Автору принадлежат

1 Разработка методов и программных средств определения глобального электронного содержания с помощью глобальных ионосферных карт и ионосферной модели IRI-2001

2 Создание и обработка баз данных глобальной и Камчатской сетей приемников GPS за 1996-2006 гг

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

4 Комплексное исследование пространственно-временных характеристик абсолютной и относительной амплитуды вариаций ПЭС в спокойных и возмущенных геомагнитных условиях в период с 1999 по 2005 г

5 Исследование ионосферного отклика на главный толчок Кроноцкого землетрясения 5 12 1997 г

Автор принимал непосредственное участие в разработке концепции глобального электронного содержания по данным измерений ПЭС на глобальной сети GPS и моделировании с использованием модели IRI-2001, в сравнении вариаций глобального электронного содержания с вариациями потока ультрафиолетового излучения Солнца и индекса солнечной активности F10 7

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

Основные результаты и выводы, приведенные в диссертации, докладывались на Байкальской научной молодежной школе по фундаментальной физике БШФФ-2004, БШФФ-2005, БШФФ-2006, Иркутск, XXI Всеросийской научной конференции «Распространение радиоволн», Йошкар-Ола, 2005, III международной конференции «Солнечно-Земные связи и электромагнитные предвестники землетрясений», с Паратунка, Камчатская обл, 2004, XXVIIIth General Assembly of International Union of Radio Science (URSI), New Delhi, 2005, VI международной конференции "Problems of Geocosmos", С-Петербург, май 2006, Международной конференции IGS-2006, Дармштадт, Германия, май 2006, Международной школе по космической погоде, Триест, Италия, май 2006, XIII международном симпозиуме «Оптика атмосферы и океана Атмосферная физика», Томск, ИОА СО РАН, июль 2006, 36-ой сессии COSPAR, Пекин, Китай, июль 2006, IRI 2006 Workshop, Buenos Aires, октябрь 2006, 3 европейской неделе по космической погоде, Брюссель, ноябрь 2006, 30-м Апатигском семинаре «Физика авроральных явлений», Апатиты, февраль 2007

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

1 Новый ионосферный параметр - глобальное электронное содержание, отличающийся от известных параметров боШлыпим соответствием отображения глобальных характеристик ионосферы в эксперименте и моделировании, метод определения и программный комплекс Уточненные параметры связи динамики ионосферы, индекса солнечной активности F10 7 и потока УФ-излучения Солнца во время 23-го солнечного цикла

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

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

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

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

Диссертация состоит из введения, трех глав, заключения и библиографического указателя, содержащего 183 ссылки Общий объем диссертации - 187 страниц, включая 6 таблиц и 53 рисунка

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

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

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

Во второй главе дано описание разработанного в диссертации метода вычисления глобального электронного содержания на основе использования экспериментальных данных и моделирования (раздел 2 1) Метод основан на использовании GIM-карт [14], вычисляемых по данным международной сети приемников GPS в нескольких научных центрах (JPL, США, CODE, Швейцария и др ) Для моделирования использовалась ионосферная модель IRI-2001 Приводится краткое описание программного комплекса GREC, который позволяет рассчитывать ГЭС на основе эмпирических и модельных входных данных для различных фаз цикла солнечной активности

Пространственный диапазон карт GIM в стандартном формате IONEX -±180° по долготе и ±90° по широте, размер элементарной ячейки GIM - 5° по долготе и 2 5° по широте, общее количество ячеек GIM равно 5184 Для каждого момента времени с двухчасовым временным разрешением из файлов IONEX известны значения вертикального ПЭС 1ц, где индексы i, j указывают на

координаты (широту и долготу) ячейки GIM Глобальное электронное содержание G(t) вычисляется путем суммирования по всем ячейкам GIM значений ПЭС, умноженных на площадь соответствующей ячейки GIM

G = TSI (i)

Значения ГЭС на основе IRI-2001 определялись в два этапа сначала рассчитывались значения ПЭС по всему земному шару в узлах широтно-долготной сетки с шагом 5° по долготе и 2 5° по широте с временным шагом 2 часа, затем полученные модельные значения ПЭС суммировались с использованием формулы (1) Итогом являлась временная зависимость смоделированных значений ГЭС M(t) Для оценки значений ГЭС в диссертации предложена единица измерения GECU = 1032 электронов

В разделах 2 2, 2 3 рассмотрена динамика ГЭС в течение 23-го цикла солнечной активности и сравнение эксперимента с результатами моделирования Временная зависимость ГЭС G(t) за 1996—2006 гг (рис 1а, черная кривая) демонстрирует существенную изменчивость — от 0 5 до 3 2 GECU Сравнение с вариациями числа солнечных пятен Rsn и с индексом солнечной активности F10 7 F(t), равным потоку солнечного радиоизлучения на длине волны 10 7 см в единицах сеп (10 22 Втм 2Гц-1), показало, что зависимость G(t) адекватно отражает влияние изменения солнечной активности на ионосферу Особенно четко это видно при сравнении зависимостей G(t), F(t) и Rsn(t), сглаженных с годовым временным окном (рис 16) Сглаженные ряды экспериментальных G(t) и модельных M(t) значений ГЭС (рис 16), демонстрируют хорошее соответствие Среднее значение относительного отклонения [G(t)-M(t)]/G(t), %, между сглаженными рядами G(t) и M(t) для данного периода не превышают 1 5%, среднеквадратичное отклонение (СКО) в пределах 5%

Регрессионная зависимость ГЭС от индекса F10 7 за период 1998-2006 гг демонстрирует линейную зависимость со значительным разбросом (СКО 50%) при относительно малых значениях F10 7 и тенденцию к насыщению при высоких значениях F10 7 (рис 26) Это согласуется с результатами, полученными в работе [13] при анализе длинных рядов суточных значений F10 7 и значений электронной концентрации NmF2 в максимуме Р2-слоя, измеренных на 20 ионосферных станциях Авторы [13] обнаружили эффект насыщения NmF2 при высоких значениях F10 7 В то же время значения ГЭС и потока УФ-излучения Солнца в диапазоне длин волн 0 1-50 нм за период 19962006 гг связаны линейной зависимостью с незначительным разбросом (СКО около 20%, рис 2а)

В диссертации установлено, что для ГЭС характерны сильные сезонные вариации с относительной амплитудой от 10% на подъеме и спаде солнечной активности до 30% в период максимума Это иллюстрируют представленные на рис 1в вариации относительной амплитуды dG(t)/G(t) (черная кривая), отфильтрованные в диапазоне периодов от 100 до 300 дней

G(t), GECU

35

G(t); M(t) GECU

25 —I -IRI2001

- — - F10 7

2 — -ГЭС

1 5 _ —--R$n S

- - * ""* jf

1 —

05

Rsn(t), F(t), сел

i- 200 ■ 160 ■ 120 80 40 0

dG/G; dM/M, %

40 —i

20 0 -20 -40

IRI-2001 ГЭС

R(t)

38 36 34

32 —|

3 28 26

IRI-2001 ГЭС

Г.

~1-1-1-1-1-1-1-1-1-1-1

1996 1997 1998 1999 2000 2001 2002 20032004 2005 2006 Год

Рис 1 а - вариации экспериментальных значений ГЭС (черная кривая), относительной амплитуды 27-дневных вариаций ГЭС (серая кривая), б -сглаженные с годовым временным окном вариации экспериментальных (черная кривая) и расчетных (серая кривая) значений ГЭС, индекса F10 7 и относительного числа солнечных пятен Rs„, в - сезонные вариации экспериментальных (черная кривая) и расчетных (серая кривая) значений ГЭС, г - отношение R(t) значений ГЭС дневной стороны Земли к значениям ГЭС ночной стороны (экспериментальные значения - черная кривая, расчетные -серая кривая)

Максимальные значения сезонных вариаций ГЭС приходятся на периоды равноденствия; это вполне согласуется с тем, что плотность нейтральной атмосферы и значения электронной концентрации максимальны в апреле и октябре fl]. Для рядов M(t) также характерны полугодовые колебания (рис, ¡в -серая кривая), однако сезонные вариации экспериментальных и модельных рядов ГЭС находятся не в фазе. Так, разница между моментами времени максимумов экспериментальных значении ГЭС и данных моделирования достигает 20^10 дней.

, GECU

1998-2006

M » И I 1 I 1 !

О 50 100150200250 F10.7, с.е.п.

0 0.2 0.4 0.6 0.8 EUV, фотонов cftf'tr' 1011

G, GECU

3.5 п

í 996-2006

Рис. 2. Регрессионные зависимости значений ГЭС от УФ (а) и значений ГЭС от индекса FIO,7 (б): аппроксимирующие зависимости показаны серыми кривыми.

Для понимания глобальных ионосферных процессов важны оценки отношения дневной и ночной ионизации. Анализ отношения R значений ГЭС дневной стороны Земли к ГЭС ночной стороны показал, что временная зависимость R(t) испытывает сильные сезонные вариации с максимумами в периоды солнцестояния. ] 1а рис ir изображены ряды экспериментальных Ro(t) и модельных RmU) данных, сглаженных с временным окном 81 день. В периоды роста и максимума солнечной активности диапазон изменения и среднее значение вариаций RM(t) заметно меньше (от 2.7 до 3,1 и 2.9 соответственно), чем соответствующие значения R(;(t) (от 2.8 до 3.5 и 3.2). Это означает, что в модели 1RI-2001 степень ионизации ночной стороны завышена. Кроме того, сезонные вариации экспериментальных и модельных рядов отношения R находятся не в фазе.

Одним из важных факторов влияния солнечного излучения на ионосферу является 27-дневная цикличность, обусловленная вращением Солнца. Однако 27-дневная модуляция ионосферных параметров, измеренных в одном пункте, маскируется множеством других факторов, что затрудняет получение достоверных количественных характеристик. ГЭС более отчетливо отражает эту особенность динамики ионосферы (раздел 2,2.4). На рис. За в показаны

27-дневные вариации dG(t), dF(t), dU(t), полученные с помощью фильтрации исходных рядов ГЭС G(t), индекса F10 7 F(t) и потока УФ излучения U(t) в диапазоне периодов 20-40 дней Для иллюстрации выбран период, который включает октябрь-ноябрь 2003 г, когда были зарегистрированы сильные гелиосферные возмущения [3]

Корреляционный анализ данных для периода 1999-2006 гг выявил высокую степень подобия 27-дневных вариаций рядов dG(t), dF(t) и dU(t) Коэффициент корреляции между рядами dG(t) и dU(t) достигает 0 94, между dG(t) и dF(t), а также dF(t) и dU(t) несколько ниже - 0 83 и 0 79 соответственно Установлено, что 27-дневные вариации ГЭС запаздывают в среднем на два дня по сравнению с соответствующими вариациями F10 7 и УФ (рис 36) Однако разброс между рядами ГЭС и F10 7 существенно больше разброса между ГЭС и УФ

На рис Зг (и За) представлены зависимости огибающей (касательной к максимумам) 27-дневных вариаций dG(t), dF(t) и dU(t), показанных на рис За-в В начале активного периода интенсивность 27-дневных вариаций F10 7, УФ и ГЭС начала возрастать практически синхронно, с разницей в несколько дней Спад интенсивности вариаций dG(t) и dU(t) также синхронен, однако наблюдается запаздывание около 30 дней относительно вариаций dF(t) Это запаздывание может быть связано с разницей между скоростью угасания яркости локальных источников на Солнце в УФ- и радиодиапазоне, а также с различными условиями «видимости» образований на лимбе Солнца, когда область генерации УФ-излучения находится уже за лимбом, а соответствующая область радиоизлучения еще не затемнена Солнцем [16]

На рис 1а серой кривой показана временная зависимость огибающей 27-дневных вариаций ГЭС G27(t), сглаженная с годовым временным окном Максимальное отклонение относительной амплитуды G27(t) меняется от 8% на подъеме и спаде солнечной активности до 2% в периоды максимума Такая закономерность обусловлена особенностями активных образований на поверхности Солнца, число которых возрастает с увеличением солнечной активности [6, 15], при этом относительная амплитуда 27-дневной модуляции потока ультрафиолетового излучения уменьшается Напротив, в периоды роста и особенно спада солнечной активности на поверхности Солнца образуется несколько интенсивных активных образований, которые могут существовать в течение нескольких оборотов Солнца [6, 15], вызывая глубокую модуляцию потока УФ-излучения

В разделе 2 2 5 приведены результаты исследования динамики регионального электронного содержания (РЭС) на протяжении 23-го цикла солнечной активности РЭС рассчитывалось для трех поясов экваториального и двух среднеширотных в Северном и Южном полушариях (±30° и от 30° до 65° геомагнитной широты соответственно)

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

два раза Отношение 11(4) РЭС дневной стороны Земли к ночной на экваторе в 1 5 раза ниже, чем дня ГЭС, что свидетельствует о более высокой степени ночной ионизации, особенно в период максимума солнечной активности

dG(t), GECU 03-1

dF(t), сеп

dG(t), GECU 03

dU(t), фот c*f2 с' 10"

0 12g.

dU(t), фот см"2 с' 10" 0 12

dF(t), сеп

г- 80

dll(t), фот cut с110", EdG(t), GECU ■ 03

EdF(t), сел 80

2003

2004

Номер дня Год

Рис 3 27-дневные вариации ГЭС и индекса F10 7 (а), ГЭС и УФ (б), УФ и F10 7 (в), огибающая 27-дневных вариаций ГЭС, индекса F10 7 и УФ (г)

Для средних широт характерен выраженный годовой цикл с максимальным значением R(t) около 8 0 для Южного полушария зимой и для Северного полушария летом Такие изменения могут быть связаны с естественным годовым циклом освещенности полушарий

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

В разделе 3 1 для оценки потенциальной возможности регистрации эффектов сейсмо-ионосферного взаимодействия было проведено исследование

ионосферного отклика на Кроноцкое ЗТ, которое произошло 5 декабря 1997 г в 1126 UT у Тихоокеанского побережья Камчатки Это самое сильное землетрясение в данном районе, зарегистрированное за время детальных сейсмологических наблюдений, его магнитуда составила М=7 9 [4] Тем не менее, впервые проведенный в диссертации анализ данных Камчатской региональной сети GPS показал, что амплитуда ионосферного отклика на главный толчок ЗТ в виде вариаций ПЭС с характерным периодом около 500 с не превышала 0 02 TECU (1016 ж2 эл), что близко к уровню фоновых вариаций ПЭС

Для оценки потенциала обнаружения возможных предвестников ЗТ проведено исследование основных характеристик относительной амплитуды вариаций ПЭС при различных уровнях геомагнитной и солнечной активности Исследование проводилось для 12 спокойных дней (индекс КР<3) и 40 возмущенных суток (Кр>3), в период с 1999 по 2005 г Для получения статистически значимых результатов были выбраны три широтных региона, представленных в сети Интернет наибольшим количеством GPS-станций высокоширотный (50-80° N, 200-300° Е, географические координаты, 59 станций), среднеширотный (20-50° N, 200-300° Е, 817 станций) и экваториальный (±20° N, 0-360° Е, 76 станций)

На первом этапе обработки данных значения «наклонного» ПЭС Is(t), вычисленные на основе представленных в сети Интернет RINEX-файлов, преобразуются в соответствующие «вертикальные» значения Iv(t) В результате фильтрации рядов Iy(t) в диапазонах периодов, соответствующих КМ ПИВ (2060 мин) и среднемасштабным (СМ) ПИВ (2-10 мин), получаются вариации абсолютного значения ПЭС 5I(t) Затем вычисляются среднее значение Imjl = <IV> и среднеквадратичное отклонение (СКО) вариаций ПЭС dlm>n = <(812)>1/2 для всех рядов Iv(t)mn и 5I(t)mn, полученных для всех m станций GPS и для всех п ИСЗ, видимых с данной станции в течение выбранного интервала времени, общее число рядов составило около 106 Относительная амплитуда вариаций ПЭС (dl/l)min определяется путем нормировки величины СКО на фоновое значение 1о, в качестве которого используется значение абсолютного вертикального ПЭС, полученное по GIM-картам

Для выбранной территории значения Imjn, dlm,n и dl/lm n усреднялись с целью получения средних значений ПЭС <1>, абсолютной <dl> и относительной амплитуд <dl/l> вариаций ПЭС в диапазоне периодов КМ ПИВ и СМ ПИВ Среднее число рядов вариаций ПЭС с хорошим качеством данных за сутки составляло не менее 18000

На рис 4 изображены вариации <dl/l> относительной амплитуды вариаций ПЭС в течение дня, когда наблюдалась магнитная буря 30 10 2003 г (слева, суммарное значение КР = 56), а также суммарные суточные зависимости усредненных значений <dl/l> для 12 спокойных дней (справа, КР < 3) Зависимости <dl/l> для СМ ПИВ показаны тонкими линиями, для КМ ПИВ -жирными При низком уровне геомагнитной возмущенности относительная амплитуда возмущений ночью существенно превышает дневные значения

на экваторе (рис 4г) и высоких широтах (рис 4е) в 3-5 раз, на средних широтах (рис 4д) - в 2 раза

Магнитная буря 30 10 2003 г

<с11/1>, %

60 п

50 :

40 -

30 :

20 : -Ъ

10 :

0 - 111

12 споконых суток

Высокие широты (50°-80вМ, 200°-300°Е) - <с!1/1>, % Э- 25Т

1 5 1

05 0

<с!1/1>, %

8 -|

6 -

Средние широты (20°-50Х 200°-300°Е) б.

<сН/|>, %

1 6 1

1 2 - П

08 -

04 -

- —

0 -I —|—I—Г"

и-

<сШ1>, % 4

Экватор (±20 N. 0 -360 Е) в.

11

■ 11111111111 0 4 8 1216 20 24 Время, !_Т

<с!1/!>, %

1 6 п

1 2 - \Л

08 -

04 -0 - мм

0 4 8 1216 20 24 Время, ЬТ

Рис 4 Вариации относительной амплитуды вариаций ПЭС во время магнитной бури 30 10 2003 г (а-в) и суммарные суточные зависимости относительной амплитуды для 12 спокойных дней (г-с) Шкала местного времени соответствует долготе 255° Е

В возмущенных условиях отмеченная закономерность суточного хода относительной амплитуды вариаций ПЭС нарушается, если главная фаза магнитной бури приходится на местное дневное время, как это произошло 30 10 2003 г (рис 4а-в) Таким образом, геомагнитный контроль амплитуды вариаций ПЭС на высоких и средних широтах оказывается более существенным, чем регулярные суточные изменения Этот вывод согласуется с полученными

ранее результатами [11] На экваторе амплитуда возмущений слабо зависит от уровня геомагнитной активности (рис 4в)

На рис 5 представлены регрессионные зависимости относительной амплитуды КМ вариаций ПЭС от значений индекса КР (а-в) и от индекса F10 7 (г-е), точками показаны усредненные в интервале 3 ч значения <dl/l>, линиями -аппроксимирующие прямые В среднем относительная амплитуда вариаций ПЭС пропорциональна значению геомагнитного индекса Кр, наиболее выражена эта зависимость на высоких широтах (коэффициент пропорциональности к=0 37), слабее на средних широтах (к=0 2) и меньше всего на экваторе (k<0 1) Следует заметить, что средние значения относительной амплитуды КМ- и СМ-возмущений получены при усреднении по большой территории и по всем «видимым» в данный интервал времени ИСЗ GPS Как следствие, среднее значение <dl/l> существенно ниже максимально возможной амплитуды возмущений

•а.

Высокие широты (50°-80°N, 200°-300°Е)

<dl/l>, %

6 -

5 4

3 -2 1

Кр

<dl/l> % г

10 - Г.

8 -6 - • • • • • •

4 -2 - • #

0 - —г I I ! I I I 1

Средние широты (20°-50°N, 200°-300°Е) <dl/l>, % р.

Экватор (±20°N, 0°-360°Е) <d!/l>, %

<dl/l> %

10 -

8 -

6 -4 - •

2 -

0 -

8 -6 -4 -2 -О -

100 200 300 F10 7, сел

100 200 300 F10 7, с ел

-1 I I —I-г

100 200 300 400 F10 7, сел

Рис 5 Регрессионные зависимости относительной амплитуды КМ-вариаций ПЭС от значений индекса Кр (а-в) и от индекса F10 7 (г-е)

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

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

Таким образом, фоновые ионосферные возмущения геомагнитного происхождения могут оказаться основным препятствием для обнаружения относительно слабых волновых возмущений, предположительно генерируемых на стадии подготовки ЗТ [10].

о, км Большая магнитная буря,

2000 issc UTj 29 октября 2003 г.

i

400 - ' ' f i—|'| ' I—1—i—1

6 6.г 6.4 6,6 6.8 7 7.2 7.4 7.5 7.В в Чэс.ит

Рис. 6. Пространственное распределение положения минимумов и максимумов крупномасштабных вариаций 11ЭС.

С другой стороны, КМ ПИВ большой амплитуды, сравнимой е суточными изменениями ПЭС, могут маскировать возможные отклонения от нормального суточного хода, с которыми связывают процессы подготовки ЗТ [12], Такие

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

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

На рис 6 для момента времени 07 00 UT (спустя —50 минут после внезапного начала магнитной бури 29 октября 2003 г в 06 11 UT) жирными кривыми черного и серого цвета показаны соответственно приблизительные положения минимумов и максимумов вариаций ПЭС Для примера тонкой черной кривой показано меридиональное распределение амплитуды КМ вариаций ПЭС в секторе С (Камчатская региональная сеть станций GPS), черными точками -положения максимума и минимумов вариаций Волновые векторы возмущений для различных меридиональных секторов А-Е, показанные стрелками, ориентированы в экваториальном направлении Кажущийся центр кольцевого возмущения находится в авроральной области Динамика КМ-возмущения вдоль направления распространения в секторе С представлена на верхней панели рис 6

При анализе данных для всех наборов GPS-станций впервые зарегистрировано КМ ПИВ в форме уединенной волны с фазовым фронтом «кольцевого» типа, образовавшееся в авроральной зоне после внезапного начала магнитной бури 29 10 2003 г Возмущение длительностью около 40 мин и с длиной волны до 3000 км распространялось в экваториальном направлении со скоростью до 1500 м/с в утреннем, дневном и вечернем секторах Северного полушария и 700-1000 м/с в ночном секторе

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

1 Впервые разработан метод и программы для вычисления по данным GPS и ионосферной модели IRI-2001 нового ионосферного параметра -глобального электронного содержания Установлено, что в течение 23-го цикла солнечной активности ГЭС менялось в пределах 0 5-3 2 GECU в соответствии с динамикой УФ-излучения Солнца и индекса F10 7 Для ГЭС характерны сильные (до 30%) сезонные вариации с максимумом в периоды равноденствия, а также сезонные вариации отношения ГЭС на дневной и ночной сторонах Земли с максимумом во время летнего и зимнего солнцестояния 27-дневные вариации ГЭС подобны соответствующим

изменениям FIO 7 и УФ (коэффициент корреляции больше 0 9), но запаздывают в среднем на два дня Относительная амплитуда 27-дневных вариаций уменьшается с 8% в периоды роста и спада солнечной активности до 2% в период максимума, что соответствует динамике активных образований на Солнце в течение 11-летнего цикла активности

2 Установлено, что в целом экспериментальные и модельные значения ГЭС демонстрируют хорошее соответствие, однако обнаружены и заметные расхождения Относительная разность между рядами рассчитанного по модели IRI-2001 и экспериментального ГЭС увеличиваются от 1 5% до 9 0 % по мере уменьшения временного окна сглаживания от одного года до 10 дней Фаза сезонных вариаций в модели различается (с отставанием или опережением) на 20—40 дней относительно экспериментальных значений В период максимума расчетное отношение ГЭС для дневной и ночной сторон Земли меняется в более узком интервале (от 2 7 до 3 1), чем в эксперименте (от 2 8 до 3 5)

3 Установлено, что среднее значение относительной амплитуды ПИВ увеличивается до 8% пропорционально значению планетарного индекса КР, коэффициент пропорциональности на высоких широтах к=0 37, на средних к=0 2, на экваторе k<0 1 При низком уровне геомагнитной активности относительная амплитуда возмущений ночью существенно превышает дневные значения на экваторе и на высоких широтах в 3-5 раз, на средних широтах в 2 раза Во время магнитных бурь геомагнитный контроль амплитуды вариаций ПЭС на высоких и средних широтах оказывается более существенным, чем регулярные суточные изменения На экваторе амплитуда возмущений практически не зависит от уровня геомагнитной активности

4 Впервые для определения фазовой скорости и формы волнового фронта крупномасштабных возмущений в различных долготных секторах Северного полушария использованы пять наборов GPS-станций, максимально распределенных по широте, включая Камчатскую региональную сеть Получено экспериментальное доказательство существования крупномасштабных возмущений с кольцевой формой волнового фронта Кажущийся центр кольцевого возмущения находится в авроральной области Такие возмущения появляются во время максимальных изменений магнитного поля Земли, имеют характерную длительность порядка 40 мин, относительную амплитуду до 40% и перемещаются в направлении, близком экваториальному со скоростью до 1500 м/с в утреннем, дневном и вечернем секторах Северного полушария и 700-1000 м/с в ночном секторе

Основные публикации автора по теме диссертации:

1 Афраймович Э Л, Живетьев И В , Бузевич А В Крупномасштабные перемещающиеся ионосферные возмущения во время магнитной бури 29-30 октября 2003 г по данным камчатской сети станций GPS // Геомагнетизм и аэрономия 2005 Т 45 № 1 С 123-126

2 Афраймович Э JI, Живетьев И В Ионосферный отклик главного толчка Кроноцкого землетрясения 5 12 1997 г по данным Камчатской региональной сети станций GPS // Вулканология и сейсмология 2005 № 4 С 56-62

3 Афраймович Э Л, Живетьев И В Динамика крупномасштабных перемещающихся ионосферных возмущений аврорального происхождения Реакция на внезапное начало магнитной бури 29 октября 2003 // Труды XXI Всероссийской научной конференции «Распространение радиоволн» 2005 Т 1 С 129-133

4 Афраймович Э Л, Живетьев И В , Косогоров Е А Относительная амплитуда вариаций полного электронного содержания по данным глобальной сети GPS // Труды VIII конференции молодых ученых «Астрофизика и физика околоземного космического пространства» Иркутск 2005 С 116-117

5 Афраймович Э Л , Астафьева Э И , Живетьев И В Солнечная активность и глобальное электронное содержание // Доклады Академии наук 2006 Т 409 № 3 С 399-402

6 Афраймович Э Л , Астафьева Э И , Ясюкевич Ю В , Живетьев И В Глобальное электронное содержание Метод определения и динамика в течение 23-го цикла солнечной активности // Труды IX конференции молодых ученых «Физические процессы в космосе и околоземной среде» БШФФ-2006 Изд-во ИСЗФ СО РАН Иркутск. 2006 С 109-111

7 Afraimovich Е L, Astafieva ЕI, Oinats А V , Yasukevich Y V , Zhivetiev I V Global electron content and solar activity comparison with IRI modeling results // Proceedmgs of IGS-2006, Darmstadt, 8-11 May http //nng esoc esa de/ws2006/Papers/p_Afraimovich_GEC_IRI_IGS_update pdfgz

8 Afraimovich E L , Astafyeva EI, Oinats A V , Yasukevich Yu V , Zhivetiev I V Global electron content as a new index of solar activity Comparison with IRI modelmg results //IRI News 2006 V 13 N 1 October 2006 A5

9 Voeykov S V , Afraimovich E L, Kosogorov E A, Perevalova N P , Zhivetiev IV Relative amplitude of medium-scale traveling ionospheric disturbances as deduced from global GPS network //

http //www cosis net/abstracts/COSPAR2006/00122/COSPAR2006-A-00122 pdf

ЦИТИРУЕМАЯ ЛИТЕРАТУРА

1 Брюнелли Б Е , Намгаладзе А А Физика ионосферы 1988 M Наука 499 с

2 Вернадский В И Научная мысль как планетное явление 1977 M Наука 191 с

3 Веселовский И С , Панасюк M И, Авдюшин СИ и др Солнечные и гелиосферные явления в октябре-ноябре 2003 г причины и следствия // Космические исследования 2004 Т 42 № 5 С 453-508

4 Гусев А А, Левина В И , Салтыков В А, Гордеев Е И Сильное Кроноцкое землетрясение 5 декабря 1997 года основные данные, сейсмичность очаговой зоны, механизм очага Макросейсмический эффект // Кроноцкое землетрясение на Камчатке 5 декабря 1997 года предвестники, особенности, последствия 1998 Петропавловск-Камчатский С 32-54

5 Захаренкова И Е Использование измерений сигналов системы GPS для обнаружения ионосферных предвестников землетрясений // Автореферат диссертации на соискание ученой степени кандидата физико-математических наук Калининград 2007 20 с

6 Мордвинов А В , Плюснина Л А Когерентные структуры в динамике крупномасштабного магнитного поля Солнца // Астрономический журнал 2001 Т 78 №8 С 753-760

7 Морозова Л И Тайфуны и сейсмичность // Доклады Академии Наук 2006 Т 410 №3 С 397^100

8 Соболев Г А, Шестопалов И П, Харин Е П Геоэффекгивные солнечные вспышки и сейсмическая активность Земли // Физика Земли 1998 № 7 С 85-90

9 Сытинский АД О связи землетрясений с солнечной активностью // Изв АН СССР Физика Земли 1989 №2 С 13-30

10 Afraimovich EL, Astafieva EI, Gokhberg MB, Lapshin VM, Permyakova V E , Steblov G M, Shalimov S L Variations of the total electron content m the ionosphere from GPS data recorded during the Hector Mme earthquake of October 16, 1999, California // Russian Journal of Earth Sciences 2004 V 6 N 5 P 339-354

11 Afraimovich E L , Kosogorov E A, Lesyuta O S , Yakovets A F , Ushakov 11 Geomagnetic control of the spectrum of traveling ionospheric disturbances based on data from a global GPS network // Ann Geophys 2001 V 19 N 7 P 723-731

12 Liu J Y, Chuo Y J, Shan S J, Tsai Y B , Chen YI, Pulinets S A, Yu S B Pre-earthquake ionospheric anomalies registered by continuous GPS TEC measurements // Ann Geophys 2004 V 22 P 1585-1593

13 Liu L , Wan W, Nmg B , Pirog O M , Kurkin VI Solar activity variations of the ionospheric peak electron density // J Geophys Res V 111 doi 10 1029/2006JA011598

14 Mannucci A J , Wilson B D , Yuan D N , Ho C H , Lindgwister U J , Runge T F A global mappmg technique for GPS-derived ionospheric TEC measurements//RadioSci 1998 V 33 N3 P 565-582

15 Mordvinov A V , Willson R C Effect of large-scale magnetic fields on total solar irradiance//Solar Phys 2003 V 215 P 5-16

16 Nusinov A A, Katyushina V V Lyman-alpha line intensity as a solar activity index in the far ultraviolet range//Solar Phys 1994 V 152 P 201-206

Отпечатано в издательском отделе ИСЗФ СО РАН Заказ № 80 от 11 апреля 2007 г Объем 22 с Тираж 150 экз

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

Список таблиц.

Список рисунков.

Список обозначений и сокращений.

Введение.

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

1.1. Ионосферные предвестники землетрясений.

1.1.1. Аномалии полного электронного содержания.

1.1.2. Изменения спектра акустико-гравитационных волн.

1.2. Вариации уровня солнечной активности и ионосферные параметры.

1.3. Геомагнитные возмущения.

1.4. Методы мониторинга ионосферных возмущений сигналами спутниковой системы GPS.

1.4.1. Общие сведения о спутниковой радионавигационной системе GPS. Глобальная сеть двухчастотных приемников GPS.

1.4.2. Определение полного электронного содержания по данным GPS.

1.4.3. Соответствие пространственно-временных характеристик вариаций ПЭС и параметров локальных неоднородностей электронной концентрации. Ракурсный эффект.

1.4.4. Определение пространственно-временных характеристик ионосферных возмущений. Метод SADM-GPS.

1.4.5. Глобальные карты абсолютного значения «вертикального» полного электронного содержания.

1.4.6. Спектр вариаций полного электронного содержания.

Глава 2. Глобальное электронное содержание.

2.1. Метод определения глобального электронного содержания, программный комплекс.

2.1.1. Общие сведения.

2.1.2. Метод определения глобального электронного содержания.

2.1.3. Ограничения анализа.

2.1.4. Программный комплекс.

2.1.5. Сравнение значений ГЭС, вычисленных по данным различных лабораторий.

2.2. Солнечная активность и ионосфера на протяжении 23-го солнечного цикла.

2.2.1. Цикл солнечной активности.

2.2.2. Полугодовые вариации.

2.2.3. Вариации «день-ночь».

2.2.4. 27-дневные вариации.

2.2.5. Электронное содержание в геомагнитных координатах.

2.3. Моделирование глобального электронного содержания.

2.3.1. Расчет ГЭС с помощью ионосферной модели IRI-2001.

2.3.2. Сравнение экспериментальных и модельных значений ГЭС.

2.4. Обсуждение результатов.

Глава 3. Ионосферные возмущения, обусловленные геомагнитной и сейсмической активностью.

3.1. Ионосферный отклик на землетрясения: Кроноцкое землетрясение

5 декабря 1997 г.

3.1.1. Общие сведения.

3.1.2. Исходные данные.

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

3.1.4. Результаты и обсуждение.

3.2. Относительная амплитуда возмущений полного электронного содержания.

3.2.1. Общие сведения об эксперименте и метод обработки.

3.2.2. Суточная зависимость амплитуды вариаций ПЭС.

3.2.3. Распределение вероятности относительной амплитуды вариаций ПЭС и зависимость от значения индексов Кр и F10.7.

3.2.4. Обсуждение.

3.3. Характеристики волнового фронта перемещающихся ионосферных возмущений аврорального происхождения.

3.3.1. Методы исследования и результаты.

3.3.2. Обсуждение результатов.

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

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

Между тем, начиная с классических работ академика Вернадского [29] в мировой литературе все большее число публикаций посвящено разработке предложенной им концепции изучения геофизических явлений, как единого процесса взаимодействия геосфер - литосферы, атмосферы, ионосферы, магнитосферы, межпланетной среды и Солнца, составляющих единую систему «Солнце-Земля». Одной из геосфер, наиболее чувствительных к процессам, сопровождающим катастрофические явления в твердой земле и атмосфере, является ионосфера.

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

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

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

На протяжении многих лет огромные усилия исследователей предпринимаются для поисков научного обоснования построения системы предсказания ЗТ на основе мониторинга состояния геосфер. Существует множество указаний на различные признаки готовящихся ЗТ в изменении локальных параметров нейтральной и заряженной компоненты ионосферы, а также в акустическом и электромагнитном поле излучения эпицентров будущих ЗТ. Считается, что на ускорение подготовки ЗТ могут влиять различные глобальные проявления динамики солнечной активности, обусловленные регулярными изменениями в течение солнечного цикла, а также вспышками, выбросами корональной массы, вызывающими магнитные и ионосферные бури [29, 64, 68]. Обсуждается гипотеза о том, что не только литосферные процессы могут быть источником ионосферных (атмосферных) возмущений, но и глобальные изменения в атмосфере сами влияют на характер протекания литосферных процессов, при этом атмосферные катаклизмы могут служить «спусковым» механизмом высвобождения сейсмической энергии [55]. Эта гипотеза также требует надежных экспериментальных доказательств.

В настоящее время широкое распространение получили исследования литосферно-ионосферных процессов на основе изучения возмущений полного электронного содержания (ПЭС) в ионосфере Земли, регистрируемых при зондировании ионосферы сигналами современных спутниковых навигационных систем - американской GPS и российской ГЛОНАСС.

В некоторых работах [39, 149] отмечено отклонение от регулярного суточного хода ПЭС вблизи эпицентра ЗТ за несколько часов и дней до главного толчка. Предпринимались неоднократные попытки обнаружить в вариациях ПЭС предсказанное теорией изменение спектра возмущений ионосферной плазмы, обусловленное генерацией в районе эпицентра акустико-гравитационных волн (АГВ) [75]. Однако до сих пор не получено надежных свидетельств подготовки катастрофических землетрясений по вариациям ПЭС, как, впрочем, и по другим характеристикам сигналов радиофизических средств мониторинга сейсмоопасных районов.

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

В диссертации изучение глобальных свойств ионосферы основано на концепции глобального электронного содержания (ГЭС), предложенной в 2005 г. Э.Л. Афраймовичем. Основное внимание в изучении локальных характеристик уделено Камчатскому региону. Сейсмическая деятельность у тихоокеанского побережья Камчатки, Курильских островов и северо-восточной Японии достигает наивысшего на Земле уровня. В районе Камчатского залива и Командорских островов расположена область стыка Курило-Камчатской и Алеутской островных дуг, которая относится к наиболее важным и интересным тектоническим узлам мира. Сюда подходит и крупный подводный хребет

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

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

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

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

1. Разработка программно-алгоритмического комплекса, предназначенного для определения глобального электронного содержания по данным международной сети GPS и для его моделирования с использованием ионосферной модели IRI-2001.

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

3. Исследование динамики глобального электронного содержания в течение 23-го цикла солнечной активности в сравнении с динамикой ультрафиолетового и радиоизлучения Солнца.

4. Исследование глобальных характеристик крупномасштабных перемещающихся ионосферных возмущений (КМ ПИВ), генерируемых в авроральной зоне во время мощных магнитных бурь.

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

Научная новизна исследования:

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

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

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

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

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

Личный вклад автора

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

Автору принадлежат:

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

2. Создание и обработка баз данных глобальной и камчатской сетей приемников GPS за 1996-2006 г.г.

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

4. Комплексное исследование пространственно-временных характеристик абсолютной и относительной амплитуды вариаций ПЭС в спокойных и возмущенных геомагнитных условиях в период с 1999 по 2005 г.г.

5. Исследование ионосферного отклика на главный толчок Кроноцкого землетрясения 5.12.1997 г.

Автор принимал непосредственное участие в разработке концепции глобального электронного содержания по данным измерений ПЭС на глобальной сети GPS и моделировании с использованием модели IRI-2001, в сравнении вариаций глобального электронного содержания с вариациями потока ультрафиолетового излучения Солнца и индекса солнечной активности F10.7.

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

Основные результаты и выводы, приведенные в диссертации, докладывались на Байкальской научной молодежной школе по фундаментальной физике БШФФ-2004, БШФФ-2005, БШФФ-2006, Иркутск; XXI Всеросийской научной конференции "Распространение радиоволн", Йошкар-Ола, 2005; III международной конференции "Солнечно-Земные связи и электромагнитные предвестники землетрясений", с. Паратунка, Камчатская обл., 2004; XXVIIIth General Assembly of International Union of Radio Science (URSI), New Delhi, 2005; VI международной конференции "Problems of Geocosmos", С-Петербург, май 2006; Международной конференции IGS-2006, Дармштадт, Германия, май 2006; Международной школе по космической погоде, Триест, Италия, май 2006; XIII международном симпозиуме «Оптика атмосферы и океана. Атмосферная физика», Томск, ИОА СО РАН, июль 2006; 36-ой сессии COSPAR, Пекин, Китай, июль 2006; IRI 2006 Workshop, Buenos Aires, октябрь 2006; 3 европейской неделе по космической погоде, Брюссель, ноябрь 2006; 30-м Апатитском семинаре «Физика авроральных явлений», Апатиты, февраль 2007.

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

1. Новый ионосферный параметр - глобальное электронное содержание, г отличающийся от известных параметров большим соответствием отображения глобальных характеристик ионосферы в эксперименте и моделировании; метод определения и программный комплекс. Уточненные параметры связи динамики ионосферы, индекса солнечной активности F10.7 и потока УФ излучения Солнца во время 23-го солнечного цикла.

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

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

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

Диссертация состоит из введения, трех глав, заключения и библиографического указателя, содержащего 183 ссылки. Общий объем диссертации - 187 страниц, включая 6 таблиц и 53 рисунка.

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

Заключение

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

1. Впервые разработан метод и программы для вычисления по данным GPS и ионосферной модели IRI-2001 нового ионосферного параметра - глобального электронного содержания. Установлено, что в течение 23-го цикла солнечной активности ГЭС менялось в пределах 0.5-3.2 GECU в соответствии с динамикой УФ излучения Солнца и индекса F10.7. Для ГЭС характерны сильные (до 30%) сезонные вариации с максимумом в периоды равноденствия, а также сезонные вариации отношения ГЭС на дневной и ночной сторонах Земли с максимумом во время летнего и зимнего солнцестояния. 27-дневные вариации ГЭС подобны соответствующим изменениям F10.7 и УФ (коэффициент корреляции больше 0.9), но запаздывают в среднем на 2 дня. Относительная амплитуда 27-дневных вариаций уменьшается с 8% в периоды роста и спада солнечной активности до 2% в период максимума, что соответствует динамике активных образований на Солнце в течение 11-летнего цикла активности.

2. Установлено, что в целом экспериментальные и модельные значения ГЭС демонстрируют хорошее соответствие, однако обнаружены и заметные расхождения. Относительная разность между рядами расчетного по модели IRI-2001 и экспериментального ГЭС увеличиваются от 1.5% до 9.0 % по мере уменьшения временного окна сглаживания от одного года до 10 дней. Фаза сезонных вариаций в модели различается (с отставанием или опережением) на 20-40 дней относительно экспериментальных значений. В период максимума расчетное отношение ГЭС для дневной и ночной сторон Земли меняется в более узком интервале (от 2.7 до 3.1), чем в эксперименте (от 2.8 до 3.5).

3. Установлено, что среднее значение относительной амплитуды ПИВ увеличивается до 8% пропорционально значению планетарного индекса Кр; коэффициент пропорциональности на высоких широтах к=0.37, на средних к=0.2, на экваторе к<0.1. При низком уровне геомагнитной активности относительная амплитуда возмущений ночью существенно превышает дневные значения: на экваторе и на высоких широтах в 3-5 раз, на средних широтах в 2 раза. Во время магнитных бурь геомагнитный контроль амплитуды вариаций ПЭС на высоких и средних широтах оказывается более существенным, чем регулярные суточные изменения. На экваторе амплитуда возмущений практически не зависит от уровня геомагнитной активности.

4. Впервые для определения фазовой скорости и формы волнового фронта крупномасштабных возмущений в различных долготных секторах северного полушария использованы пять наборов GPS-станций, максимально распределенных по широте, включая Камчатскую региональную сеть. Получено экспериментальное доказательство существования крупномасштабных возмущений с кольцевой формой волнового фронта. Кажущийся центр кольцевого возмущения находится в авроральной области. Такие возмущения появляются во время максимальных изменений магнитного поля Земли, имеют характерную длительность порядка 40 мин, относительную амплитуду до 40% и перемещаются в направлении, близком к экваториальному со скоростью до 1500 м/с в утреннем, дневном и вечернем секторах северного полушария и 700-1000 м/с в ночном секторе.

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

Автор искренне благодарит своего научного руководителя доктора физико-математических наук, профессора Афраймовича Эдуарда Леонтьевича.

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

Диссертация выполнена при частичной поддержке научно-исследовательского проекта СО РАН - ДВО РАН - НАН Украины "Развитие радиофизических методов диагностики околоземного космического пространства и исследование механизмов взаимодействия геосферных оболочек" N 3.24.

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

1. Аароне Д. Глобальная морфология ионосферных мерцаний // ТИИЭР. 1982. Т.70. №4. С.45-65.

2. Акасофу С.И., Чепмен С. Солнечно-земная физика. Т.1. 1975. М.: Мир. 512с.

3. Алимов О.А., Гохберг М.Б., Липеровская Е.В., и др. Эффект резких уменьшений плотности спорадического слоя Es ионосферы предвестник землетрясений// ДАН СССР. 1989. Т.305. №6. С.1335-1339.

4. Андрианов В.А., Мосин B.JL, Смирнов В.М. Применение радиосигналов спутниковой навигационной системы для зондирования ионосферы Земли //Радиотехника и электроника. 1996. Т.41. №9. С. 1029-1032.

5. Архив навигационных файлов центра SOP АС // ftp://lox.ucsd.edu/pub/nav.

6. Архив файлов IONEX // ftp://cddisa.gsfc.nasa.gov/pub/gps/products/ionex.

7. Архив файлов RINEX центра SOPAC // ftp://lox.ucsd.edu/pub/rinex.

8. Афраймович Э.Л., Абрамов В.А., Варшавский И.И., Величанский Б.Н., и др. Наблюдение ионосферных эффектов наземного промышленного взрыва методами радиозондирования // Изв. АН СССР. Физика Земли. 1985. №11. С.99-103.

9. Афраймович Э.Л., Астафьева Э.И., Живетьев И.В. Солнечная активность и глобальное электронное содержание // Доклады Академии Наук. 2006. Т.409. №3. С.399-402.

10. Афраймович Э.Л., Астафьева Э.И., Кирюшкин В.В. Ионосферное возмущение в ближней зоне эпицентра землетрясения на о. Хоккайдо 25 сентября 2003 г. //Изв. Вузов. Радиофизика. 2005. Т.48. №4. С.299-313.

11. Афраймович Э.Л., Воейков С.В. Экспериментальное доказательство существования уединенной внутренней гравитационной волны в атмосфБере земли во время мощной магнитной бури 30.10.2003 г. // Доклады Академии Наук. 2004. Т.399. №6. С.683-686.

12. Афраймович Э.Л., Вугмейстер Б.О., Гохберг М.Б., Альперович Л.С., и др. Акустическая волна взрыва // Изв. АН СССР. Физика Земли. 1985. №11. С.32-42.

13. Афраймович Э.Л., Живетьев И.В. Ионосферный отклик главного толчка Кроноцкого землетрясения 5.12.1997 г. по данным Камчатской региональной сети станций GPS // Вулканология и сейсмология. 2005. №4. С.56-62.

14. Афраймович Э.Л., Живетьев И.В., Бузевич А.В. Крупномасштабные перемещающиеся ионосферные возмущения во время магнитной бури 2930 октября 2003 г. по данным камчатской сети станций GPS // Геомагнетизм и аэрономия. 2005. Т.45. №1. С.123-126.

15. Афраймович Э.Л., Косогоров Е.А., Лесюта О.С., Ушаков И.И. Спектр перемещающихся ионосферных возмущений по данным глобальной сети GPS // Изв. ВУЗов. Радиофизика. 2001. T.XLIV. №10. С.828-839.

16. Афраймович Э.Л., Косогоров Е.А., Плотников А.В. Ударно-акустические волны, генерируемые при запусках ракет и землетрясениях // Космические исследования. 2002. Т.40. №3. С.383-393.

17. Афраймович Э.Л., Косогоров Е.А., Плотников А.В., Уралов A.M. Параметры ударно-акустических волн, генерируемых при землетрясениях // Физика Земли. 2001. №6. С. 16-28.

18. Афраймович Э.Л., Перевалова Н.П. Моделирование измерений полного электронного содержания на GPS радиоинтерферометре // Солнечно-земная физика. 2004. Т.4. С.71-78.

19. Ахмедов P.P., Куницын В.Е. Моделирование ионосферных возмущений, вызванных землетрясениями и взрывами // Геомагнетизм и аэрономия. 2004. Т.44. №1. С.105-112.

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

21. Бучаченко А.Л., Ораевский В.Н., Похотелов О.А., и др. Ионосферные предвестники землетрясений // Успехи физических наук. 1996. Т. 166. №9. С.1023-1029.

22. Вернадский В.И. Научная мысль как планетное явление. 1977. М.: Наука. 191 с.

23. Веселовский И.С., Панасюк М.И., Авдюшин С.И., и. д. Солнечные и гелиосферные явления в октябре-ноябре 2003 г.: причины и следствия // Космические исследования. 2004. Т.42. №5. С.453-508.

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

25. Витинский Ю.И., Оль А.И., Сазонов Б.И. Солнце и атмосфера Земли. 1976. Л.: Гидрометеоиздат. 353 с.

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

27. Гохберг М.Б., Моргунов В.А., Похотелов О.А. Сейсмоэлектромагнитные явления. 1988. М.: Наука. 174 с.

28. Гохберг М.Б., Шалимов С.Л. Литосферно-ионосферная связь и ее моделирование // Российский журнал наук о Земле. 2000. Т.2. №2.

29. Дружин Г.И. Опыт прогноза Камчатских землетрясений на основе наблюдений за электромагнитным ОНЧ излучением // Вулканология и сейсмология. 2002. №6. С.51-62.

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

31. Захаренкова И.Е. Использование измерений сигналов системы GPS для обнаружения ионосферных предвестников землетрясений // Автореферат диссертации на соискание ученой степени кандидата физико-математических наук. Калининград. 2007. С.20.

32. Захаренкова И.Е., Шагимуратов И.И., Лаговский А.Ф. Модификация ионосферы в период подготовки землетрясений по данным спутниковой системы GPS // Труды XXI Всероссийской научной конференции "Распространение радиоволн". 2005. Т.1. С.194-198.

33. Захаренкова И.Е., Шагимуратов И.И., Лаговский А.Ф., Кранковский А. Исследование ионосферных предвестников для Землетрясений класса м~5.0 // Электронный журнал «Исследовано в России». 2006. Т.39. С. 1047-1055.

34. Иванов-Холодный Г.С., Никольский Г.М. Солнце и ионосфера. 1969. М.: Наука. 456 с.

35. Исходный текст модели IRI-2001 // ftp://nssdcftp.gsfc.nasa.gov/models.

36. Калинин Ю.К., Сергеенко Н.П. Движущиеся уединенные макронеоднородности, возникающие в ионосфере за несколько часов до катастрофических землетрясений // Доклады Академии Наук. 2002. Т.387. №1. С. 105-107.

37. Карты ПЭС высокоширотных областей "NOR" // http://www.kn.nz.dlr.de/daily/tec-np.

38. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. 1977. М.: Наука. 832 с.

39. Кравцов Ю.А., Фейзулин З.И., Виноградов А.Г. Прохождение радиоволн через атмосферу земли. 1983. М.: Радио и связь. 224 с.

40. Кривицкий Б.Х., ed. Справочник по теоретическим основам радиоэлектроники. Т.2. 1977. М.: Энергия. 305-308.

41. Кринберг И.А., Тащилин А.В. Ионосфера и плазмосфера. 1984. М.: Наука. 189 с.

42. Левин В.Е., Гордеев Е.И., Бахтиаров В.Ф., Касахара М. Предварительные результаты GPS мониторинга на Камчатке и Командорских островах // Вулканология и сейсмология. 2002. №1. С.3-11.

43. Леонович Л.А., Афраймович Э.Л., Портнягина О.Ю. Скорости и направления перемещения крупномасштабных возмущений полного электронного содержания во время больших магнитных бурь // Геомагнетизм и аэрономия. 2004. Т.44. №2. С.166-173.

44. Липеровский В.А., Алимов О.А., Шалимов С.Л., и др. Исследование F-области ионосферы перед землетрясениями // Физика Земли. 1990. №12. С.77-85.

45. Липеровский В.А., Похотелов О.А., Шалимов С.Л. Ионосферные предвестники землетрясений. 1992. М.: Наука.

46. Мордвинов А.В., Плюснина Л.А. Когерентные структуры в динамике крупномасштабного магнитного поля Солнца // Астрономический журнал. 2001. Т.78. №8. С.753-760.

47. Морозова Л.И. Тайфуны и сейсмичность // Доклады Академии Наук. 2006. Т.410. №3. С.397-400.

48. Намгаладзе А.А., Захаров Л.П. Возмущение состава нейтральной атмосферы и термосферных ветров на F-область ионосферы // Исследование ионосферной динамики. 1979.М.: ИЗМИР АН. С.84-95.

49. Намгаладзе А.А., Латышев К.С. Влияние верхних граничных условий на моделируемые ионосферные параметры // Геомагнетизм и аэрономия. 1976. Т. 16. №1. С.43-49.

50. Нусинов А.А. Ионосфера как природный детектор для исследования долговременных изменений потоков солнечного геоэффективного излучения // Геомагнетизм и аэрономия. 2004. Т.44. №6. С.779-786.

51. Перцев Н.Н., Шалимов C.JT. Генерация атмосферных гравитационных волн в сейсмически активном регионе и их влияние на ионосферу // Геомагнетизм и Аэрономия. 1996. Т.36. №2. С.111-118.

52. Пулинец С.А., Легенька А.Д., Зеленова Т.И. Зависимость сейсмо-ионосферных вариаций в максимуме слоя от местного времени // Геомагнетизм и аэрономия. 1998. Т.38. №3.

53. Ришбет Г., Гарриот O.K. Введение в физику ионосферы. 1975. Л.: Гидрометеоиздат. С.304.

54. Смертин В.М., Намгаладзе А.А. Исследование зависимости характеристик внутренних гравитационных волн от параметров источника // Геомагнетизм и аэрономия. 1981. Т.21. №2. С.302-308.

55. Смертин В.М., Намгаладзе А.А. Ионосферные эффекты внутренних гравитационных волн, генерируемых авроральной электроструей // Ионосферное прогнозирование. 1982.М.: Наука. С.107-113.

56. Соболев Г.А., Шестопалов И.П., Харин Е.П. Геоэффективные солнечные вспышки и сейсмическая активность Земли // Физика Земли. 1998. №7. С.85-90.

57. Соловьев Ю.А. Системы спутниковой навигации. 2000. М.: Эко-Трендз. 267 с.

58. Сорокин В.М., Чмырев В.М. Электродинамическая модель ионосферных предвестников землетрясений и некоторых видов катастроф // Геомагнетизм и аэрономия. 2002. Т.42. №6. С.821 -830.

59. Сытинский А.Д. Связь сейсмичности Земли с солнечной активностью и атмосферными процессами. 1987. Л.: Гидрометеоиздат. 99 с.

60. Сытинский А.Д. О связи землетрясений с солнечной активностью // Изв. АН СССР. Физика Земли. 1989. №2. С.13-30.

61. Физика космоса: Маленькая энциклопедия. 1986. М.: Советская энциклопедия. 783 с.

62. Харгривс Д.К. Верхняя атмосфера и солнечно-земные связи 1982. Л.: Гидрометеоиздат. С.352

63. Хэррис Ф.Д. Использование окон при гармоническом анализе методом дискретного преобразования Фурье // ТИИЭР. 1978. Т.вв. №1. С.60-95.

64. Afraimovich E.L. Dynamics and anisotropy of traveling ionospheric disturbances as deduced from transionospheric sounding data. Statistical angle-of-arrival and doppler method (SADM) // Preprint ISTP. 1995. №5-95. P.54.

65. Afraimovich E.L. Statistical angle-of-arrival and doppler method (SADM) for determining characteristics of the dynamics of the transionospheric radio signal interference pattern //Acta Geod. Geoph. Hung. 1997. V.32. №3-4. P.461-468.

66. Afraimovich E.L., Astafieva E.L, Kirushkin V.V. Localization of the source of ionospheric disturbance generated during an earthquake // International Journal of Geomagnetism and Aeronomy V.6. doi: 10.1029/2004GI000092

67. Afraimovich E.L., Astafyeva E.I., Oinats A.V., Yasukevich Y.V., Zhivetiev I.V. Global electron content as a new index of solar activity. Comparison with IRI modeling results // IRI News. 2006. V.13. №1.

68. Afraimovich E.L., Boitman O.N., Zhovty E.I., Kalikhman A.D., Pirog T.G. Dynamics and anisotropy of traveling ionospheric disturbances as deduced from transionospheric sounding data // Radio Sci. 1999. V.34. №2. P.477-487.

69. Afraimovich E.L., Chernukhov V.V., Demyanov V.V. Updating the ionospheric delay model for single-frequency equipment of users of the GPS // Radio Sci. 2000. V.35. №1. P.257-262.

70. Afraimovich E.L., Karachentsev V.A. Testing of the transionospheric radiochannel using data from the global GPS network // Annals of Geophysics. 2003. V.46. №6. P. 1229-1246.

71. Afraimovich E.L., Kosogorov E.A., Lesyuta O.S., Yakovets A.F., Ushakov I.I. Geomagnetic control of the spectrum of traveling ionospheric disturbances based on data from a global GPS network // Ann. Geophys. 2001. V.19. №7. P.723-731.

72. Afraimovich E.L., Minko N.P., Fridman S.V. Spectral and dispersion characteristics of medium-scale travelling ionospheric disturbances as deduced from transionospheric sounding data // J. Atmos. Terr. Phys. 1994. V.56. №11. P.1431-1446.

73. Afraimovich E.L., Palamartchouk K.S., Perevalova N.P. GPS radio interferometry of travelling ionospheric disturbances // J. Atm. Terr. Phys. 1998. V.60. №12. P.1205-1223.

74. Afraimovich E.L., Palamartchouk K.S., Perevalova N.P. Statistical angle-of-arrival and doppler method for GPS interferometry of TIDs // Adv. Space Res. 2000. V.26. №6. P.1001-1004.

75. Afraimovich E.L., Perevalova N.P., Plotnikov A.V., Uralov A.M. The shock-acoustic waves generated by the earthquakes // Ann. Geophys. 2001. V.19. №4. P.395-409.

76. Beach T.L., Kelley M.C., Kintner P.M. Total electron content variations due to nonclassical traveling ionospheric disturbances: Theory and Global Positioning System observations // J. Geophys. Res. 1997. V.102. P.7279-7292.

77. Bertel L., Bertin F., Testud J. De la mesure du contenu electronique integre appliquee a l'observation des ondes de gravite de moyenne echelle // J. Atm. Terr. Phys. 1976. V.38. P.261-270.

78. Beynon W.J.G., Brown G.M. Region E and solar activity // J. Atm. Phys. 1959. V.15. P.168-174.

79. Bhuyan P.K., Borah R.R. TEC derived from GPS network in India and comparison with models, in GA URSI'2005. 2005. Delhi, India. http://www.ursi.org/Proceedings/ProcGA05/pdf/GFl a.6(0754).pdf.

80. Bilitza D. International Reference Ionosphere 1990 // NSSDC. Greenbelt, Maryland. 1990. P.90-22.

81. Bilitza D. International reference ionosphere // Radio Sci. 2001. V.36. №2. P.261-275.

82. Bilitza D. A correction for the IRI topside electron density model based on Alouette/ISIS topside sounder data // Adv. Space Res. 2004. V.33. №6. P.838-843.

83. Blanc E., Jacobson A.R. Observation of ionospheric disturbances following a 5-kt chemical explosion. 2. Prolonged anomalies and stratifications in the lower thermosphere after shock passage // Radio Sci. 1989. V.24. P.739-746.

84. Bowman G.G. A review of some recent work on midlatitude spread F occurrence as detected by ionosondes // J. Geomagn. Geolectr. 1990. V.42. P.109-138.

85. Bowman G.G. Nighttime mid-Latitude travelling ionospheric disturbances associated with mild spread-F conditions // J. Geomagn. Geolectr. 1991. V.43. №8. P.899-920.

86. Bowman G.G. Upper atmosphere neutral-particle density variations compared with spread-F occurrence rates at locations around the world // Ann. Geophys. 1992. V.10. P.676-682.

87. Brekke A., Doupnik J.D., Banks P.M. Incoherent scatter measurements of E-region conductivities and currents in the auroral zone // J. Geophys. Res. 1974. V.79. P.3773-3790.

88. Calais E., Minster B.J. GPS detection of ionospheric perturbations following the January 1994, Northridge earthquake // Geoph. Res. Lett. 1995. V.22. P. 1045-1048.

89. Calais E., Minster B.J., Hofton M.A., Hedlin M.A.H. Ionospheric signature of surface mine blasts from Global Positioning System measurements // Geophys. J. Int. 1998. V.132. P. 191-202.

90. Carlos R.C., Massey R.S. The Los Alamos beacon receiver Array // IEEE Transactions on Geosciences and Remote Sensing // 1994. V.32. P.954-958.

91. Center for Orbit Determination in Europe, University of Berne, Switzerland (CODG) // http://www.cx.unibe.ch.

92. Coisson P., Radicella S.M. The IRI topside parameters // Adv. Radio Sci. 2004. V.2. P.249-251.

93. Coisson P., Radicella S.M. Ionospheric topside models compared with experimental electron density profiles // Ann. Geophys. 2005. V.48. №3. P.497-503.

94. Coisson P., Radicella S.M., Leitinger R., Nava B. New option for IRI topside electron density profile using NeQuick model // IRI NewsLetter. 2005. №12. P.5-7.

95. Coisson P., Radicella S.M., Nava B. Comparisons of experimental topside electron concentration profiles with IRI and NeQuick models // Ann. Geophys. 2002. V.45. №1. P.l 11-116.

96. Daniell J.R.E., Brown L.D., Anderson D.N., Fox M.W., Doherty P.H., Decker D.T., Sojka J.J., Schunk R.W. Parameterized ionospheric model: A global ionospheric parameterization based on first principles models // Radio Sci. 1995. V.30. P.1499-1510.

97. Davies K. Recent progress in satellite radio beacon studies with particular emphasis on the ATS-6 radio beacon experiment // Space Science Review. 1980. V.25. №4. P.357-430.

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

99. Davis M.J. The integrated ionospheric response to internal atmospheric gravity waves // J. Atm. Terr. Phys. 1973. V.35. P.929-959.

100. Fitzgerald T.J. Observations of total electron content perturbations in GPS signals caused by a ground level explosion // J. Atmos. Terr. Phys. 1997. V.59. №7. P.829-834.

101. Gallagher D.L., Craven P.D., Comfort R.H. An empirical model of the Earth's plasmasphere //Adv. Space Res. 1988. V.8. P. 15-24.

102. Geodetic Survey Division of Natural Resources Canada (EMRG) // http://www.nrcan-rncan.gc.ca/.

103. GEOPACK-2005 // http://modelweb.gsfc.nasa.gov/magnetos/data-based/Geopack 2005.doc.

104. Georges T.M. HF Doppler studies of TID's // J. Atm. Terr. Phys. 1968. V.30. №5. P.735-746.

105. Georges T.M., Hooke W.H. Wave-induced fluctuations in ionospheric electron content: a model indicating some observational biases // J. Geophys. Res. 1970. V.75. P.6295-6308.

106. Giovanni G., Radicella S.M. An analytical model of the electron density profile in the ionosphere // Adv. Space Res. 1990. V.10. №11. P.27-30.

107. Grup Universitat Politecnica de Catalunya (UPCG) // http://www.upc.es.

108. Gurtner W. The Receiver Independent Exchange Format Version 2 // http://igscb.ipl.nasa.gov/igscb/data/format/rinex2.txt.

109. Hao J., Tang Т., Li D. Progress in the research on atmospheric electric field anomaly as an index for short-impending prediction of earthquakes // J. Earthquake. Pred. Res. 2000. V.8. №3. P.241-255.

110. Haykowicz L.A., Hunsucker R.D. A simultaneous observation of large-scale periodic TIDs in both hemispheres following an onset of auroral disturbances // Planet. Space Sci. 1987. V.35. №6. P.785-791.

111. Hegai V.V., Kim V.P., Liu J.Y. The ionospheric effect of atmospheric gravity waves excited prior to strong earthquake // Adv. Space Res. 2006. V.37. P.653-659.

112. Hegai V.V., Kim V.P., Nikiforova L.I. A possible generation mechanism of acoustic-gravity waves in the ionosphere before strong earthquakes // J. Earthquake Pred. Res. 1997. V.6. №4. P.584-589.

113. Hines C.O. Internal atmospheric gravity waves at ionospheric heights // Can. J. Phys. 1960. V.38. P.1441-1481.

114. Но C.M., Mannucci A.J., Lindqwister U.J., Pi X., Tsurutani B.T. Global ionosphere perturbations monitored by the worldwide GPS network // Geophys. Res. Lett. 1996. V.23. №22. P.3219-3222.

115. Носке К., Schlegel К. A review of atmospheric gravity waves and travelling ionospheric disturbances: 1982-1995 // Ann. Geophys. 1996. V.14. P.917-940.

116. Hofmann-Wellenhof В., Lichtenegger H., Collins J. Global Positioning System: Theory and Practice. 1992. New York: Springer-Verlag Wien. P.327.

117. Huang X., Reinisch B.W., Bilitza D., Benson R.F. Electron density profiles of the topside ionosphere // Ann. Geophys. 2002. V.45. №1. P.125-130.

118. Hunsucker R.D. Atmospheric gravity waves generated in the high-latitude ionosphere. A review // Review of Geophysics. 1982. V.20. №2. P.293-315.

119. IRI-2001 // ftp://nssdcftp.gsfc.nasa.gov/models/space/model/ionos/iri.html.

120. ITU-2004: Ionospheric propagation data and prediction methods required for the design of satellite services and systems // Recommendation ITU-R. 2004. P.531-7.

121. Jacobson A.R., Carlos R.C., Massey R.S., Wu G. Observations of traveling ionospheric disturbances with a satellite-beacon radio interferometer: Seasonal and local time behavior//J. Geophys. Res. 1995. V.100. P. 1653-1665.

122. Jakowski N., Fichtelmann В., Jungstand A. Solar activity control of Ionospheric and thermospheric processes // J. Atmos. Terr. Phys. 1991. V.53. P.l 125-1130.

123. Jakowski N., Heise S., Wehrenpfennig A., Schluter S., Reimer R. GPS/GLONASS-based TEC measurements as a contributor for space weather // J. Atmos. Sol.-Terr. Phys. 2002. V.64. №5-6. P.729-735.

124. Jakowski N., Sardon E., Schluter S. GPS-based TEC observations in comparison with IRI95 and the european TEC model NTCM2 // Adv. Space Res. 1998. V.22. P.803-806.

125. Jet Propulsion Laboratory of California Institute of Technology (JPLG) // http://www.jpl.nasa.gov.

126. Kelley M.C., Miller C.A. Electrodynamics of midlatitude spread F 3. Electrohydrodynamic waves? A new look at the role of electric fields in thermospheric wave dynamics // J. Geophys. Res. 1997. V.102. P.l 153911548.

127. Klobuchar J.A. Ionospheric time-delay algorithm for single-frequency GPS users // IEEE Transactions on Aerospace and Electronics System. 1986. V.23. №3. P.325-331.

128. Kondo G. The variation of the atmospheric field at the time of earthquake // Mem. Kakioka Magnet. Observatory. 1968. V.13. №1. P.l 1-13.

129. Kotake N., Otsuka Y., Tsugawa Т., Ogawa Т., Saito A. Climatological study of GPS total electron content variations caused by medium-scale travelingionospheric disturbances // J. Geophys. Res. V.lll. doi: 10.1029/2005JAO11418

130. Kumar P., Foufoula-Georgiou E. Wavelet Applications in Geophysics: A Review// Reviews of Geophysics. 1997. V.35. №4. P.385-412.

131. Lawrence R.S., Little C.G., Chivers J.A. The influence of the ionosphere upon radio wave propagation «Earth-space» // Proc. IEEE. 1964. V.52 №4. P.5-30.

132. Leitinger R., Zhang M., Radicella S.M. An improved bottomside for the ionospheric electron density model NeQuick // Ann. Geophys. 2005. V.48. №3. P.525-534.

133. Li Y.Q., Jacobson A.R., Carlos R.C., Massey R.S., Taranenko Y.N., Wu G. The blast wave of the Shuttle plume at ionospheric heights // Geophys. Res. Lett. 1994. V.21. P.2737-2740.

134. Liu J.Y., Chen Y.I., Chuo Y.J., Tsai H.F. Variations of ionospheric total electron content during the Chi-Chi earthquake // Geophys. Res. Lett. 2001. V.28. №7. P.1383-1386.

135. Liu J.Y., Chuo Y.J., Shan S.J., Tsai Y.B., Chen Y.I., Pulinets S.A., Yu S.B. Pre-earthquake ionospheric anomalies registered by continuous GPS TEC measurements//Ann. Geophys. 2004. V.22. P.1585-1593.

136. Liu L., Wan W., Ning В., Pirog O.M., Kurkin V.I. Solar activity variations of the ionospheric peak electron density // J. Geophys. Res. V.lll. doi: 10.1029/2006JAO 11598

137. Mannucci A.J., Wilson B.D., Yuan D.N., Но C.H., Lindgwister U.J., Runge T.F. A global mapping technique for GPS-derived ionospheric TEC measurements // Radio Sci. 1998. V.33. №3. P.565-582.

138. Mercier C. Observations of atmospheric gravity waves by radiointerferometry // J. Atm. Terr. Phys. 1986. V.48. №7. P.605-624.

139. Mercier C. Some characteristics of atmospheric gravity waves observed by radio-interferometry // Ann. Geophys. 1996. V.14. P.42-58.

140. Mercier С., Jacobson A.R. Observations of atmospheric gravity waves by radio interferometry: are results biased by the observational technique // Ann. Geophys. 1997. V.15. P.430-442.

141. Mordvinov A.V., Willson R.C. Effect of large-scale magnetic fields on total solar irradiance// Solar Phys. 2003. V.215. P.5-16.

142. Nusinov A.A., Katyushina V.V. Lyman-alpha line intensity as a solar activity index in the far ultraviolet range // Solar Phys. 1994. V.152. P.201-206.

143. Oliver W.L., Otsuka Y., Sato M., Takami Т., Fukao S. Climatology of F region gravity waves propagation over the middle and upper atmosphere radar //J. Geophys. Res. 1997. V.102. P. 14499-14512.

144. Ondoh T. Seismo-ionospheric effects, Atmospheric and ionospheric electromagnetic phenomena associated with earthquakes // ed. M. Hayakawa, Tokyo. 1999. P.789-803.

145. Oraevsky V.N., Ruzhin Y.Y., Shagimuratov I.I. Anomalies of ionospheric TEC above Turkey before two strong earthquakes at 1999 // Proceed. 15th Wroclaw EMC Symposium. 2000.

146. Oyama K.I., Noguchi K., Izawa M., Saito A., Otsuka Y., Tohyama F. Local time, annual, latitude, and seasonal variations of total electron content over Japan // Institute of space and astronautical science. ISAS Research Note 796. 2005.

147. Parrot M. Statistical study of ELF/VLF emissions recorded by a low-altitude satellite during seismic events // J. Geophys. Res. 1994. V.99. №12. P.23339-23347.

148. Perkins F. Spread F and ionospheric currents // J. Geophys. Res. 1973. V.78. P.218-226.

149. Pulinets S.A. Seismic activity as a source of the ionospheric variability // Adv. Space Res. 1998. V.22. №6. P.903-906.

150. Pulinets S.A., Boyarchuk K. Ionospheric precursors of earthquakes. 2004. Berlin, Heidelberg: Springer-Vergal.

151. Pulinets S.A., Legen'ka A.D., Gaivoronskaya T.V., Depuev V.K. Main phenomenological features of ionospheric precursors of strong earthquakes // J. Atmos. Sol.-Terr. Phys. 2003. V.65. №16-18. P.1337-1347.

152. Radicella S.M., Leitinger R. The evolution of the DGR approach to model electron density profiles // Adv. Space Res. 2001. V.27. №1. P.35-40.

153. Richmond D., Matsushita S. Thermospheric response to a magnetic substorm // J. Geophys. Res. 1975. V.80. №19. P.2839-2850.

154. Row R.V. Acoustic-gravity waves in the upper atmosphere due to a nuclear detonation and an earthquake //J. Geophys. Res. 1967. V.72. P. 1599-1610.

155. Schaer S., Beutler G., Rothacher M. Mapping and predicting the ionosphere. // Proceedings of the IGS AC Workshop, Darmstadt, Germany. February 9-11. 1998. P.307-320.

156. Schaer S., Gurtner W. IONEX: The IONosphere Map EXchange Format Version 1 // Proceedings of the IGS AC Workshop, Darmstadt, Germany, February 9-11. 1998.

157. SEM/ SOHO // http://www.usc.edu/dept/space science/semdatafolder/.

158. SIDC Solar Influences Data Analysis Center // http://sidc.oma.be/sunspot-data/.

159. Solar Radio Monitoring Programme // http://www.drao-ofr.hia-iha.nrc-cnrc.gc.ca.

160. Spoelstra T.A., Kelder H. Effects produced by the ionosphere on radio interferometry//Radio Sci. 1984. V.19. P.779-788.

161. Torrence C., Compo G.P. A Practical Guide to Wavelet Analysis // Amer. Meteor. Soc. 1998. V.19. P.61-78.

162. Tsugawa Т., Saito A., Otsuka Y. A statistical study of large-scale traveling ionospheric disturbances using the GPS network in Japan // J. Geophys. Res. V.109. doi: 10.1029/2003J AO 10302

163. Tsugawa Т., Shiokawa K., Otsuka Y., Ogawa Т., Saito A., Nishioka M. Geomagnetic conjugate observations of large-scale traveling ionospheric disturbances using GPS networks in Japan and Australia // J. Geophys. Res. V.lll. doi: 10.1029/2005JA011300

164. Voeykov S.V., Afraimovich E.L., Kosogorov E.A., Perevalova N.P., Zhivetiev I.V. Relative amplitude of medium-scale traveling ionospheric disturbances as deduced from global GPS network // 2006.

165. Williams P.J.S., Virdi T.S., Lewis R.V., Lester M., Rodger A.S., McCrea I.W., Freeman K.S.C. Worldwide atmospheric gravity wave study in the European sector 1985-1990//J. Atmos. Terr. Phys. 1993. V.55. P.683-696.

166. Wilson B.D., Mannucci A.J., Edwards C.D. Subdaily northern hemisphere maps using the IGS GPS network // Radio Sci. 1995. V.30. P.639-648.

167. Yeh K.C., Liu C.H. Radio wave scintillations in the ionosphere // Proceedings of the IEEE. 1982. V.70. №4. P.324.