УЛЬТРАНИЗКОЧАСТОТНЫЕ ВОЛНЫ В КОСМОСЕ И НА ЗЕМЛЕ Текст научной статьи по специальности «Физика»

УЛЬТРАНИЗКОЧАСТОТНЫЕ ВОЛНЫ В КОСМОСЕ И НА ЗЕМЛЕ

В.А. Пилипенко, заведующий лабораторией1, главный научный сотрудник2

Д.Д. Позднякова1, инженер

Н.В. Савельева1, младший научный сотрудник

1Институт физики Земли РАН

^Геофизический Центр РАН

(Россия, г. Москва)

DOI:10.24412/2500-1000-2024-9-3-163-205

Работа над обзором была поддержана грантом РНФ №24-77-10012 (ПДД).

Аннотация. Дан краткий обзор современных представлений о физической природе волновых вариаций геомагнитного поля в ультранизкочастотном (УНЧ) диапазоне - от долей мГц до первых Гц. Шумы и колебания УНЧ диапазона отражают турбулентный характер взаимодействия плазмы солнечного ветра с земной магнитосферой. Возбуждение некоторых типов УНЧ волн являются проявлением резонансного взаимодействия волн и частиц в околоземной плазме. Физика УНЧ волн опирается на теоретический аппарат физики плазмы и магнитной гидродинамики, а экспериментальной базой исследований являются данные наблюдений на спутниках и наземных станциях. Приведены морфологические характеристики основных типов геомагнитных пульсаций, а затем рассмотрены физические механизмы их генерации и возможная роль в протекании плазменных процессов в околоземной среде. Геомагнитные УНЧ волны оказываются не только одним из основных средств диагностики околоземного пространства, но и активным фактором воздействия как на околоземную среду, так и на наземные и спутниковые технологические системы. Изложение рассчитано на читателей, работающих в смежных областях (космические технологии, геофизика), которые хотели бы получить общее представление о физике УНЧ волн в околоземном пространстве.

Ключевые слова: магнитосфера Земли, УНЧ волны, солнечно-земные связи, радиационный пояс Земли, геоиндуцированные токи.

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

Все околоземное пространство заполнено ультранизкочастотными (УНЧ, и11ха-1о,^й"едиепсу - по западной терминологии) волнами. УНЧ волны занимают диапазон, покрывающий три порядка по частоте, от первых мГц до нескольких Герц. УНЧ волны проявляются в данных спутниковых и наземных магнитометров как геомагнитные пульсации, т. е. короткопе-риодные колебания магнитного поля. Верхняя частота диапазона геомагнитных пульсаций соответствует гирочастоте протонов в магнитосфере (~3-5 Гц), нижняя определяется основным периодом колебаний магнитосферной полости в целом (~доли мГц). УНЧ волны, пронизывающие все околоземное космическое пространство, доносят до земной поверхности информацию о свойствах окружающей Зем-

лю плазмы [Walker, 1993; Guglielmi and Pokhotelov, 1996].

УНЧ пульсации геомагнитного поля были, по существу, первыми электромагнитными волнами, зарегистрированными человечеством с помощью простейших магнитометров. Энергетически УНЧ вол-

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

Рис. 1. Общая структура магнитосферы Земли

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

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

г 4

\

Рис. 2. Спектральные плотности мощности различных типов природных излучений в околоземном пространстве.

Морфологические свойства и классификация геомагнитных пульсаций

Геомагнитные пульсации разделяются на два больших класса [Троицкая и Гулье-льми, 1969]: непрерывные (Pc - pulsations continuous) пульсации от Рс1 до Рс5, которые характеризуются квазисинусоидальной формой, и иррегулярные (Pi - pulsations irregular) пульсации от Pi1 до Pi3, имеющие вид отдельных всплесков или широкополосных вариаций. Как правило, непрерывные (Рс) пульсации характерны для дневного времени, а иррегулярные пульсации (Pi) - для ночного. Генерация пульсаций типа Pc определяется крупномасштабной динамикой околоземной плазмы и в большинстве случаев связана с взаимодействием солнечного ветра с магнитосферой. Иррегулярные пульсации типа Pi, как правило, связаны с развитием

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

Основные типы геомагнитных пульсаций исторически классифицируются согласно диапазону частот/периодов, в котором они наблюдаются [Пудовкин и др., 1976]. Как правило, каждый тип пульсаций отличается от другого не только диапазоном периодов, но и физической природой возбуждения волн. Установленные классификацией границы периодов для того или иного типа пульсаций носят условный характер, т. к. параметры околоземной среды крайне динамичны. Характерная амплитуда геомагнитных пульсаций резко возрастает с увеличением их периодов: так, если амплитуда короткопериодных пульсаций (Рс1-2, РИ) составляет сотые и десятые доли нТл, то амплитуда длинно-

периодных пульсаций (Рс5, Pi3) может достигать нескольких сотен нТл.

Геомагнитные пульсации Рс5

Пульсации Рс5 (Т ~ 150-600 с, f ~ 1,56,7 мГц) заполняют низкочастотную часть спектра УНЧ волн. В высоких широтах амплитуды Рс5 обычно составляют ~40-100 нТл. Пульсации Рс5 являются, в ос-

новном, высокоширотным явлением и типичны для восстановительной фазы суббури. На рисунке 3 показан пример регистрации Рс5 пульсаций в авроральных широтах. Колебания Рс4 диапазона (Т~40-150 с) как правило являются гармоникой Рс5 колебаний.

Рис. 3. Пример записи пульсаций Рс5 в авроральных широтах (Х-компонента)

Максимум амплитуды Рс5 наблюдается в авроральной зоне и смещается к более низким широтам при увеличении магнитной активности. В суточном ходе частоты появления пульсаций отмечается два максимума: основной - в утренние часы (0709 МЬТ), и второй, несколько меньший - в послеполуденные часы (15-17 МЬТ). Период наблюдаемых пульсаций возрастает с увеличением широты. Направление вращения вектора поляризации волн различно на широтах выше и ниже амплитудного максимума на данном меридиане.

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

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

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

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

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

го сектора (рис. 4). Амплитуда глобальных Рс5 пульсаций почти на порядок больше, чем амплитуда обычных Рс5 пульсаций и может достигать ~300-400 нТл.

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

Наблюдения на спутниках показали, что на больших расстояниях (>6Яе, где Яе -радиус Земли) пульсации диапазона Рс5 являются типичным явлением во внешних областях возмущенной магнитосферы. С наблюдательной точки зрения магнито-сферные УНЧ волны в частотном диапазоне Pc4-5 можно разделить по их характерной поляризации на основные типы: тороидальные, компрессионные, полои-дальные. В свою очередь поляризация МГД волн тесно связана с их пространственной структурой.

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

- Полоидальные пульсации являются азимутально-мелкомасштабными колебаниями с волновыми числами m>>1. Эта разновидность колебаний получила название storm-time Pc5 waves (рис. 5). Хотя амплитуда этих колебания велика (на геостационарной орбите достигает 1/3 от по-

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

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

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

дальными, т. е. имеют малый азимутальный масштаб.

Модельная пространственная структура тороидальных и полоидальных пульсаций показана на рисунке 6.

Рис. 6. Иллюстрация пространственной структуры тороидальных (малые азимутальные числа m) и полоидальных (большие m) колебаний

Геомагнитные пульсации Рс2-4

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

с ростом широты точки наблюдения. Так, амплитуда колебаний с периодом Т~20-30 с в средних широтах составляет десятые доли нТл, а в высоких - единицы и первые десятки нТл. В магнитосфере на спутниках геомагнитные пульсации диапазона РсЗ-4 также наблюдаются регулярно. Пример пульсаций Рс3, зарегистрированных на спутнике ММБ, дан на рисунке 7.

Рис. 7. Пример магнитограммы пульсаций Рс3, зарегистрированных трехкомпонентным

магнитометром спутника ММБ

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

синфазное в ^-компоненте и противофазное в О-компоненте поля.

С ростом магнитной активности период Рс3 колебаний, регистрируемых в данной точке, уменьшается и сдвигается в диапазон Рс2, а при уменьшении магнитной активности - в диапазон Рс4. В широтном распределении амплитуд РсЗ пульсаций четко выделяется два максимума: средне-

широтный (геомагнитные широты 55-60°) и высокоширотный (геомагнитные широты 70-74°) по обе стороны проекции плаз-мопаузы. Поляризация волн вблизи максимумов линейная, знак вращения вектора поляризации меняется на широте амплитудного максимума.

Характеристики Рс2-4 пульсаций тесно связаны с параеметрами турбулентной области перед фронтом ударной волны (форшок), где происходит возбуждение УНЧ волн (upstream waves) пучками отраженных протонов на фронте ударной волны. Частота возбуждаемых волн связана с величиной межпланетного магнитного поля (В) линейным соотношением: f [мГц]~6В [нТл]. Надежно установлен эффект модуляции интенсивности Рс3 пульсаций ориентацией межпланетного магнитного поля (ММП), а именно конусного угла □ □- угла между направлением ММП и линией Земля-Солнце. Интенсивность

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

Геомагнитные пульсации Рс1

Геомагнитные пульсации Рс1 (диапазон частот от долей Гц до нескольких Гц) наблюдаются как в виде неструктурированного излучения, так и в виде структурированной последовательности квазимонохроматических волновых пакетов. Структурированные Рс1 получили поэтическое название «жемчужины» (pearls) за их характерную модуляцию амплитуды в виде отдельных цугов длительностью 1 -4 мин., повторяющихся через 1-2 мин., и наблюдающихся непрерывными сериями длительностью от нескольких минут до нескольких часов. Колебания могут одновременно регистрироваться в большом диапазоне долгот (до 120°) и несколько сотен км по широте. Пример магнитограммы и спектрограммы Рс1 пульсаций «жемчужин» показан на рисунке 8.

Рис. 8. Магнитограммы (вверху) и динамический спектр (внизу) структурированных Рс1

пульсаций (жемчужин)

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

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

Подклассом Рс1 пульсаций являются колебания с убывающим периодом (КУП) или IPDP (irregular pulsations with diminishing periodj, которые на земной поверхности наиболее часто наблюдаются в послеполуденном и вечернем секторах в возмущенных условиях. На магнитограммах они видны как серии отдельных цугов, подобных колебаниям Pcl, но с монотонно увеличивающейся частотой от долей Гц до нескольких Гц (рис. 9, верхняя панель).

f, Гц 1.5

1.0

О,?

-

-

■J рч m f

j jj ,4 1,- £ m

m Щ Ш EHUtl -УЙ В ''flfffe' - L. <5 r z j

16

17

1»

19 IT

19 20 21

Рис. 9. Динамические спектры геомагнитных пульсаций КУП (вверху) и геомагнитных

пульсаций PiB и PiC (внизу).

Чем сильнее магнитная буря, тем вероятнее возбуждение IPDP и тем круче нарастает частота колебаний. Другим важным отличием IPDP пульсаций от «жемчужин» является их синхронное появление в сопряженных точках. Генерация КУП

наблюдается на фоне интенсивных магни-тосферных суббурь и усиления кольцевого тока.

Иррегулярные геомагнитные пульсации Р

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

Короткопериодные иррегулярные пульсации (РИ) с ^40 с разделяются на два типа: широкополосные импульсные (РiВ) и шумовые (РС) в конечной частотной полосе ^~5-10 с). В возмущенных условиях оба типа пульсаций могут регистрироваться одновременно. Пример динамического спектра этих пульсаций приведен на рисунке 9 (нижняя панель).

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

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

Геомагнитные пульсации П2-3

Во время геомагнитной суббури в ночном секторе магнитосферы и ионосферы

разыгрываются целый комплекс волновых процессов. Неотъемлемым элементом суббури являются Pi2 геомагнитные пульсации (T~50-150 с), представляющие собой цуг затухающих колебаний длительностью порядка нескольких периодов [Keiling and Takahashi, 2011]. На земной поверхности эти глобальные пульсации наблюдаются на всех широтах, от экватора до ночной стороны полярной шапки. На низких широтах Pi2 регистрируются на всех долготах, включая даже дневную сторону, в средних и высоких широтах они отмечаются наиболее часто в ночном секторе. Амплитуда Pi2 пульсаций возрастает с широтой, на низких широтах она составляет единицы и даже десятые доли нТл, а в авроральной зоне может достигать значений ~100 нТл.

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

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

2015 1015 (256), LT GEO 0.61

22:10 22:12 22:14 22:16 22:18 22:20

UT

Рис. 10. Пример одновременных Pi2 пульсаций на станциях наземной сети AMBER (Африка) и на низкоорбитальном спутнике SWARM (положение станций и проекция орбиты показаны на отдельной вставке). Верхние три панели - магнитограммы компонент Bt, By, Bp на спутнике, 4-я панель - географическая широта спутника Ф, следующие шесть панелей - магнитограммы наземных станций, нижняя панель - индекс авроральной активности АЕ

На фоне суббури часто наблюдаются длиннопериодные иррегулярные пульсации Pi3 (или Ps6) с ^10-20 минут. Эти пульсации являются тонкой структурой

суббури, и представляют собой квазипериодическую последовательность импульсов с амплитудами до нескольких сотен нТл (рис. 11).

06 - 07 APR 2000

06/04

Рис. 11. Пример суббури с одновременными Pi3 пульсациями (Х-компонента). Верхняя панель показывает вариации SYM-H индекса, характеризующего интенсивность магнитной бури. Красная черта указывает период времени с аномалиями в работе систем железнодорожной сигнализации.

Длиннопериодные волны в дневной пограничной области

Дневная высокоширотная область характеризуется увеличенной УНЧ активностью, включающей длиннопериодные квазипериодические вариации и Pc5 волны. Нерегулярные низкочастотные (T-доли мГц) пульсации наблюдаются почти каждый день, когда наземная обсерватория оказывается в пределах нескольких часов от местного полудня (рис. 12). Эти высокоширотные квазипериодические возмущения получили название IPCL (Irregular Pulsations at Cusp Latitudes). Ранние исследования дневной УНЧ-активности в высоких широтах давали надежду, что эти длиннопериодные нерегулярные вариации могут быть использованы в качестве простого индикатора положения дневного каспа и границы полярной шапки. Предположительно, граница замкнутых и открытых силовых линий должна проявлять-

ся в переходе от колебаний Pc5 к IPCL колебаниям с большим периодом.

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

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

Events), аномалиями горячего потока (Hot Flow Anomalies), перемежаемостью МГД

турбулентности магнитослоя.

Рис. 12. Пример IPCL пульсаций на высокоширотных станциях сети IMAGE. Местный

полдень соответствует 09 UT.

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

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

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

Физические механизмы УНЧ волн МГД колебания (пульсации Рс3-5) Хотя исследование плазменных и электромагнитных волновых явлений в околоземном пространстве опирается на математический аппарат и теоретические подходы, разработанные в физике плазмы и магнитной гидродинамике, нельзя сказать, что космическая геофизика занимается только приложением физических представлений, выработанных в общей физике плазмы. Ряд принципиально новых физических идей, получивших затем широкое развитие, взят физиками из космической геофизики: пересоединение силовых линий и аннигиляция магнитных полей, бес-столкновительные ударные волны [Solar Wind Sources of Magnetospheric Ultra-Low-

Frequency Waves, 1994]. Сюда же следует добавить и введенные Альвеном представления о магнитной гидродинамике и МГД волнах, отмеченные Нобелевской премией [Pilipenko, 1990]. Эти волны диапазона от первых единиц до сотни мГц регистрируются спутниковыми и наземными магнитометрами в виде Рс3-5 пульсаций геомагнитного поля [Гульельми и Троицкая, 1973]. Конверсия МГД волн в области аль-веновского резонанса, впервые предложенная для интерпретации структуры Рс3-5 волн в магнитосфере Земли, ныне широко используется в работах по управляемому термоядерному синтезу для радиочастотного нагрева плазмы.

МГД описывает крупномасштабную динамику плазмы как медленные движения проводящей жидкости и основана на самосогласованной системе уравнений гидродинамики, дополненных уравнениями Максвелла. Масштабы МГД возмущений должны превышать характкрный кинетический масштаб - ионный ларморов-ский радиус p = u / Ц (где ui - тепловая скорость ионов), а характерные частоты быть много меньше частоты циклотронного вращения ионов (со □ Г2Д Эти параметры можно оценить по соотношениям Ц [sec-1] □ О.ЩиТ] ,

p [m]= 4.6 • 10%/W[kev]/B[nT],

Здесь УА = В / VМг,Р - альвеновская скорость, которую можно оценить по соотношению V. [km/s] = 22 ки=к1 и к^ - волновой вектор вдоль и поперек В0.

Здесь V = VТ / М - ионно-звуковая

скорость, определяемая температурой электронов Te и массой ионов Mi. ММЗ напоминает колебания воздуха внутри труб органа, и энергия ММЗ распростра-

щ[m/s]=4.4•10^Л/W[këv], W - энергия частицы. При таких низких частотах проводимость плазмы вдоль внешнего магнитного поля В0 очень велика, что позволяет считать продольное электрическое поле отсутствующим —» 0 ^ а силовые линии - эквипотенциальными. При этом в отсутствие столкновений движения не слишком горячей плазмы и магнитных силовых линий как бы вморожены друг в друга, Е + [V х В ] = 0, т. е. поперечное электрическое поле Е неразрывно связано со

скоростью движения плазмы V.

В холодной плазме, где отношение теплового и магнитного давлений

р = 2Иор / <<1, могут распространяться два типа МГД волн - быстрые магнито-звуковые (БМЗ) и альвеновские. БМЗ волны аналогичны обычным звуковым волнам, только возмущение в них передается не сжатием газа, а сжатием магнитного поля, и роль газового давления P играет давление магнитного поля B2/2дo. Как и акустические волны, БМЗ распространяется практически изотропно. Дисперсионное уравнение для БМЗ в наиболее важном для ситуации в околоземной плазме пределе, когда поперечный масштаб волны много меньше ее продольного масштаба к1 » к^, имеет вид:

(1)

В плазме конечного давления (Р~1) возникает мода медленного магнитного звука (ММЗ). В том же пределе к1 » к^ дисперсионное уравнение для ММЗ имеет вид:

(2)

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

а = к Vs

вофазе. Таким образом, ММЗ волны обладают свойством диамагнитности.

Если БМЗ и ММЗ волны во многом аналогичны акустическим волнам в газе, то альвеновские волны принципиально отличаются от обычных волновых процессов в оптике, акустике или сейсмологии. Специфика магнитосферы как колебательной системы заключается в наличии в спектре ее собственных колебаний частотных интервалов, соответствующих альвеновским колебаниям (альвеновскому континууму) [Alperovich and Fedorov, 2007]. Наглядно альвеновский континуум можно представлять как арфу с набором непрерывно изменяющихся струн. Их спектр даже в ограниченной системе непрерывен, а собственные функции сингулярны. Альвенов-ские волны обладают и другими интересными особенностями - одномерному характеру распространения вдоль поля В0 и отсутствию геометрического затухания, возможности переноса нестационарного продольного тока вдоль силовых линий, отсутствию незатухающих решений даже в бездиссипативной системе, необратимой конверсии БМЗ возмущений в альвенов-ские волны в неоднородной среде.

В МГД возмущениях продольная компонента электрического поля пренебрежимо мала E ^ 0 из-за высокой проводимости плазмы вдоль внешнего магнитного поля. Поперечное к полю B0 = B0ez (ez - единичный вектор вдоль внешнего поля) 2D электрическое поле E± МГД волны можно представить как сумму потенциальной и вихревой компонент E±=-У±Ф + [V±x eг , где

V± = (dx, - поперечный Лапласиан.

Потенциалы Ф и ¥ характеризуют поле альвеновской и БМЗ мод.

БМЗ волна, подобно звуковым волнам, стремится накапливаться в области с минимумом фазовой скорости. Таким области с минимумом альвеновской скорости оказываются естественными волноводами для БМЗ волн. В неоднородной по Х плазме, для БМЗ волны возникает точка отражение (поворота), в которой эффективное волновое число к2(х) = с2 / У2 (х) - к2 - к2г обращается в 0. Областью прозрачности для БМЗ волны является область, с > см , где сСМ = (к2у + к2)У2(х). в БМЗ моде поперечная магнитная компонента безвихревая, [У± х В ] = 0, и отсутствует продольный ток _/г = 0 .

В альвеновской волне поперечные компоненты электрического и магнитные полей характеризуются потенциалом □, т. е. Е = -У±Ф и ^В = фгхУ1]5гФ. Альве-

новская волна переносит поперечный .и и продольный 7г токи, которые определяются соотношениями . 1 = 2 АУА 1д( Е ± и

Л = 2а^ 1® , где 2а =1/ Моуа - волновая альвеновская проводимость,

УА = В0 / ^ЦоР - альвеновская скорость, □

- плотность плазмы. В однородной плазме дисперсионное уравнение для альвенов-ской волны имеет вид с = са = к2УА . Характерное волновое число альвеновской волны кА = с / УА . Электрическое поле

альвеновского возмущения Е , распространяющегося в однородной плазме вдоль магнитного поля, описывается волновым уравнением:

(3)

(5 2 - v/a?) Ex (t, z) = Ио dtf;\t, z)

В правой части фигурирует ток внешне- отношению к В0 и бездивергентна

го драйвера , возбуждающего альве- у±В1 =0, при этом пр°д°льная компо-

новскую волну. Магнитная компонента нента исчезает В = 0 . поля альвеновской волны поперечна по

Рис. 13. Модель «плазменного ящика», моделирующая магнитосферный альвеновский резонатор, который образуется силовыми линиями геомагнитного поля, ограниченными сопряженными ионосферами ( г = ±Ь ). Прямые силовые линии направлены по оси Z, а радиальная неоднородность альвеновской скорости по оси Х моделируется неоднородностью

плотности плазмы Р( х). В направлении Y система однородна

Альвеновская волна, в отличие от БМЗ, не имеет точки отражения. Однако частичное отражение может происходить на границе плазменной системы (ионосфере) или от резких градиентов альвеновской скорости, что приводит к образованию альвеновских резонаторов. Основные представления о магнитосферном альве-новском резонаторе, который образуется силовыми линиями геомагнитного поля, ограниченными сопряженными ионосферами ( г = ±Ь ), были первоначально сформулированы в рамках одномерной модели т. н. «плазменного ящика» (рис. 13). В этой модели прямые силовые линии направлены по оси Z, а радиальная неоднородность альвеновской скорости по оси Х моделируется неоднородностью плотно-

Неоднородность системы проявляется в наличии члена с (к2А)' =дхк2л(х) . Уравнение (4) имеет две характерных точки: син-

сти плазмы Р(х). В направлении Y, соответствующем азимутальному направлению в магнитосфере, система однородна. В рамках модели плазменного ящика альве-новский континуум образуется непрерывным спектром мод магнитосферного резонатора с фундаментальным периодом Та ( х ) = 4 Ь / Га ( х) .

В неоднородной плазме альвеновские и БМЗ волны не могут распространяться независимо друг от друга, а оказываются связанными. Для пространственно-временной гармоники вида

ехр(-ш? + Исуу + гк2г) структура поля связанных БМЗ и альвеновских колебаний описывается уравнением:

(4)

гулярность в точке альвеновского резонанса кл = к2; и точку поворота для БМЗ

волны к\( х) - кА - к2 =0 . Сингулярный ха-

д2Ъг (k2J дЪ

дх2 k2A - k]

дх

+(kA - k] - k]) = о

рактер поля МГД волны в области альве-новского резонанса, полученный в модели плазменного ящика, сохраняется и в плазменных конфигурациях с многомерной неоднородностью [Alperovich and Fedorov, 2007]. Сингулярный характер решения волновых МГД уравнений в неоднородной

ву (х, f) = ь0(/)h~;

Bx (x, f) = b0(f)h-;

где z = x -(f) + iSm, где xa (f) - координата резонансной силовой линии; Sm - полуширина резонансной области. Метрические коэффициенты (hx,опреде-ляются выбором геометрии магнитного поля. Выражения (5) справедливы для самой общей геометрии магнитного поля и обобщают результаты более простой модели плазменного ящика. Например, коэффициент С*1=0 в геометрии с прямыми силовыми линиями или с постоянной кривизной силовых линий. Принципиальной особенностью разложения (5) является наличие сингулярного члена 1/z в компоненте By, ответственного за частотно-зависимое усиление поля волны на резонансной оболочке xa (f).

Качественно решение (5) может быть интерпретировано следующим образом. МГД возмущения БМЗ типа из внешней

Главный член асимптотического разложения (5), описывающий амплитудные и

фазовые характеристики Ву составляющей

Из-за диссипации в системе собственные частоты имеют мнимую добавку:

ша ^ шл - . Декремент затухания у > 0

среде указывает на возможность усиления волны в области альвеновского резонанса ш = шл (х) . Математическое описание

структуры поля в окрестности резонансной оболочки в криволинейном поле можно получить в виде асимптотики:

(5)

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

Тл = {Ч»^ (6)

-I

поля пульсаций в окрестности резонансной оболочки имеет вид:

(7)

связан с шириной резонанса Пт над ионосферой соотношением 5т =-у(дхшА )-1. Исходя из (7) можно качественно предста-

h

z+ C ln(im — z) + ...

h

Ф

im—ln(im — z) + ...

К К

Bv (x, f) = b0(f)

iS_

x - xA (f) + iSm

вить меридиональную структуру поля колебаний в виде суперпозиции сигнала «источника» Ъо(/) и резонансного отклика магнитосферного альвеновского резонатора, т. е. *5т /[х" ха (/) + *5т ] (рис. 14). Сигнал «источника» обусловлен возмущением, переносимым крупномасштабной БМЗ волной, и слабо зависит от координаты. Резонансный же отклик магнитосферного резонатора, связанный с возбуждением альвеновских колебаний, сильно локали-

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

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

ского резонатора

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

частоты дх®л . Обычно в магнитосфере (за исключением узкой области вблизи плаз-мопаузы) частота уменьшается с ростом

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

Широтный профиль фундаментального периода магнитосферного альвеновского

резонатора, т. е. зависимость Т (Ь) от

нормированного радиального расстояния L=r/RE (Яе - радиус Земли), показан на рис. 15. Рассчитанный профиль показывает, что в принципе возможно образование трех областей с резонансным откликом на внешнее возмущение - внутри плазмосфе-ры, на плазмопаузе, и вне плазмосферы (р^та1гои§Ь). Однако при детальном анализе широтного распределения Pc3-5 пульсаций были надежно выявлены субав-роральный и среднеширотный максимумы, которые обусловлены магнитосферными резонансными эффектами в области вне плазмосферы и в плазмосфере, в то время как максимум на плазмопаузе либо не выделялся, либо был весьма слабым. Причина такого кажущегося парадокса, в том, что «каноническое» распределение TA(L) с резким падением на плазмопаузе на самом деле достаточно редкое событие и для его формирования требуется очень резкое па-

дение плотности магнитосферной плазмы с расстоянием: N (г) <х , где ^>8! В результате профиль та^) лишь выполажива-ется в районе плазмопаузы, и резонансный пик размывается и становится слабо заметным. Таким образом, при типичной

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

Рис. 15. Вычисленный широтный профиль фундаментального периода магнитосферного альвеновского резонатора, т. е. зависимость периода собственных тороидальных колебаний TA(L) от нормированного радиального расстояния L=r/RE или высоты г. Вертикальная

стрелка - положение плазмопаузы

Качественно новые особенности магнитосферного альвеновского резонатора начинают проявляться на низких геомагнитных широтах (ниже 30о), где значительная часть силовой линии оказывается погруженной в ионосферную плазму с тяжелыми ионами. «Нагружение» силовых линий тяжелыми ионосферными ионами приводит к: (а) резкому увеличению декремента затухания Рс3 колебаний по мере приближения к экваториальным широтам, и (б) аномальной зависимости резонансного периода Тл (Ф) - росту с уменьшением

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

Объемные и поверхностные моды МГД колебаний

Объемная мода (cavity mode) образуется БМЗ модой, захваченной между магнито-паузой и плазмопаузой или Землей. Часто-

та волноводных мод определяется условием квантования в радиальном направлении кх/,х~ п, где Лх - радиальный масштаб магнитосферного волновода. Оценка основной четвертьволновой моды Т~4Ьх/Уа соответствует частоте порядка первых мГц, если предполагать Lx~5 RE и Ка~800 км/с. В азимутальном направлении распределение плазмы не имеет резких градиентов, и волна может свободно распространяться в ночную магнитосферу. Таким образом, объемная волна является по существу модой магнитосферного проточного волновода.

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

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

В реалистичной неоднородной магнитосфере глобальные поверхностные волны

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

Рис. 16. Иллюстрация магнитосферы как природного мазера для МГД волн, в котором возмущения солнечного ветра возбуждают магнитосферный волновод. Далее они фильтруются и усиливаются в магнитосферном альвеновском резонаторе в резонансной области и высвечиваются к земной поверхности через полупрозрачные зеркала (ионосферы)

Поверхностная волна не является собственным решением МГД уравнений. Она проявляется при отклике на внешнее воздействие системы альвеновских резонаторов с широким спектром, сжатых в очень узкой области. Здесь происходит их интерференция, что приближенно сводится к

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

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

верхностной волны можно оценить по соотношениям [Archer et al., 2024]:

ю2 = 2k2zV2AlV2A2

WSW Tr2 , j/-2 A2

V2 + V

V al + V '

Ysw

я

k - k2

8 k

-{ka)

(8)

При сильном перепаде альвеновских скоростей, УА1 >> УА1, частота поверхностной моды определяется минимальной аль-

веновскои скоростью

ю

SW

ю

Al

. При

наблюдениях поверхностная волна проявляется как затухающий периодический отклик системы с частотой Пз\¥, пространственным радиальным масштабом ~ к^у , и нерезистивным затуханием во времени с относительным декрементом

Уш / тт ~ (kуa) из-за конверсии начального возмущения в альвеновский континуум.

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

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

Глобальные объемные и поверхностные МГД моды играют особую роль в волновой физике [Леонович и Мазур, 2016]. Любая система более эффективно поглощает энергию набегающего турбулентного потока на собственной частоте глобальной моды. По сути, плазменная система с объемной или поверхностной волной работает как резонатор Фабри-Перро, обеспечивающий естественный процесс частотной фильтрации в турбулентной гелиогеофи-зической плазме.

Рис. 17. Основные механизмы генерации МГД волн в околоземном космическом пространстве: флуктуации динамического давления солнечного ветра (слева) и неустойчивость Кельвина-Гельмгольца тангенциального разрыва между потоком плазмы солнечного

ветра и магнитосферой

Источники пульсаций Рс4-5

Источником пульсаций Рс4-5 могут служить флуктуации динамического давления солнечного ветра (рис. 17). Кроме того, сама граница между магнитосферой и обтекающим ее потоком солнечного ветра может стать источником МГД волн при достаточно большом сдвиге скоростей благодаря неустойчивости Кельвина-Гельмгольца [Пудовкин и др., 1976]. Таким образом, неустойчивости тангенциального разрыва на утреннем и вечернем флангах являются одним из основных механизмов генерации МГД волн в околоземном пространстве. Порог возбуждения неустойчивости Кельвина-Гельмгольца зависит от параметров возмущения, но наиболее существенным фактором является скорость обтекания магнитосферы потоком солнечного ветра. Эти представления подтверждаются ростом интенсивности и длительности Рс5 пульсаций с ростом скорости солнечного ветра. Дальнейшая эволюция возбуждаемых колебаний внутри маг-

нитосферы была описана в предыдущих разделах.

Источники пульсаций Рс3

Первичным источником УНЧ волн Рс3 диапазона является турбулентная область перед фронтом стоячей ударной волны, образующейся при набегании потока солнечной плазмы на земную магнитосферу (рис. 18). Возбуждение волн (upstream waves) происходит в форшоковой области в результате циклотронной неустойчивости пучков отраженных протонов на фронте ударной волны при ориентации силовых линий межпланетного магнитного поля (ММП) вдоль нормали к фронту ударной волны (квазипараллельная ударная волна) [Гульельми, 1979]. Вблизи лобовой области ударной волн условие квазипараллельности означает малость конусного угла 9 - угла между направлением ММП и линией Земля-Солнце. Эти соображения подтверждаются надежно установленным эффектом модуляции эффективности возбуждения РсЗ пульсаций ориентацией ММП (т. е. конусного угла □).

Рис. 18. Турбулентная область перед фронтом стоячей ударной волны (форшок) как источник волн Рс3 диапазона в магнитосфере

Частота возбуждаемой ЭМИЦ волны определяется комбинацией дисперсионного соотношения и резонансного условия,

Из-за слабого поля В в форшоке частота возбуждаемых волн попадает в диапазон Рс3. Далее, эти волны сносятся потоком солнечного ветра из форшока к фронту ударной волны и маг-нитослою. Затем ЭМИЦ излучения просачиваются через переходную область, проникают в магнитосферу, и в конечном счете попадают на Землю \Yumoto, 1986]. Эти соображения хорошо подтверждаются свя-

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

Методы наземного мониторинга околоземной плазмы с использованием волн Рс3-5

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

®и (УА/и,)П

логии, для геомагнитных пульсации весьма ориентировочно известны своИства их источников (местонахождение, спектральный состав и т. п.). Тем не менее, существует принципиальная возможность использования регистрируемых на земнои поверхности геомагнитных пульсации Рс3-5 для МГД диагностики - определения свойств магнитосферной и ионосферной плазмы [Menk and Waters, 2013]. Физической основой МГД диагностики является эффект резонансной трансформации МГД волн в магнитосфере. Резонансная частота jа (Ф) определяется для данной широты □ локальными свойствами распределения магнитосферной плазмы вдоль силовой линии, а добротность магнитосферного альвеновского резонатора - диссипатив-ными свойствами ионосферы. Таким образом, уверенное выделение резонансных эффектов и наличие надежных методов экспериментального определения параметров магнитосферного альвеновского резонатора открывают возможность мониторинга плотности околоземной плазмы и проводимости ионосферы по наземным данным [Федоров и др., 1992].

Наиболее эффективным образом удается отстроиться от влияния спектра источника и выделить локальные резонансные особенности с помощью градиентного метода, основанного на прецизионных измерениях поля пульсаций на малой базе (-нескольких десятков км) [Pilipenko and Fedorov, 1993]. Для определения собственной частоты силовых линий, ее радиального градиента и ширины резонансной области наиболее надежным инструментом наземной МГД диагностики магнитосферы служит разновидность градиентно-

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

Генерация низкочастотных волн энергичными частицами

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

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

магнитному полю Ц- (строго говоря, из-за индукционного электрического поля, возбуждаемого вариациями компоненты ^и). Для УНЧ волн, частота которых много меньше гирочастоты ионов (со << Ц ), изменение кинетической энергии частиц W описывается соотношением:

5W

= qV • E + qVD • E + И

dt

8t

(9)

где И = mVj2 / 2B - магнитный момент, q - заряд частицы.

Рис. 19. Схематичная картина дрейфов частиц в неоднородном магнитосферном магнитном поле

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

и2

Ув = V + V , где ^ - градиентный

дрейф, VR =

ir

2QR

- дрейф кривизны,

) - 0 . При этом

рек поля Во. В равновесной конфигурации радиальные градиенты магнитного поля кв=дг 1п В и давления плазмы

кр = дг 1п Р не произвольны, а связаны с

радиусом кривизны поля Я соотношением

А 1п А-Д 1 + ——

протоны дрейфуют на запад, а электроны -на восток. Схематично картина дрейфов частиц в магнитосферном магнитном поле показана на рисунке 19.

Наиболее эффективно обмен энергией между УНЧ волнами и частицами происходит в случае баунс-дрейфового резонанса:

ми, и - тепловые скорости вдоль и попе-

а - mad - Kab = 0

(10)

где = (К + ^я )/ я - усредненная угловая скорость магнитного дрейфа, К -номер баунс-гармоники. Волновой вектор ку связан с азимутальным волновым числом т в дипольном поле соотношением ку = т / ЬЯЕ. Благодаря свойствам симметрии дрейфовый резонанс (К=0) может генерировать нечетные гармоники колебаний силовой линии (в том числе и основную гармонику), а баунс-дрейфовый резо-

нанс К=П1 может генерировать четные гармоники колебаний силовой линии.

Будет ли в результате резонансного взаимодействия происходить перенос энергии от волны к частицам (бесстолкно-вительное затухание волны) или возбуждение волн (развитие неустойчивости), зависит от функции распределения частиц ^(Ж,х) по энергии Ж и пространству. Выражение для скорости роста/затухания волны (инкремент/декремент) имеет вид:

у <х Im ^

со-о*.

dF

i с - mo, - Kc dW i

(11)

где <...> означает усреднение по пространству скоростей,

[^xV.F] kill юл=---1 = ——к- - частота лар-

* QdF / dW 2Q P

моровского дрейфа плазмы, пропорциональная радиальному градиенту функции распределения. Для равновесного макс-велловского распределения dF / dW < 0, и волна испытывает бесстолкновительное затухание П<0. Возбуждение неустойчивости (П>0) может возникать, когда функция распределения имеет инверсное распределение по энергиям (положительный наклон) dF / dW > 0 (распределение типа bump-on-tail). Инверсная функция распределения может естественным образом появиться в магнитосфере из-за дисперсии энергий, инжектированных в суббурю частиц, поскольку частицы с большей энергией дрейфуют быстрее, чем частицы с меньшей энергией.

Другой источник свободной энергии для роста плазменных неустойчивостей -резкие неоднородности давления энергичных протонов кольцевого тока. Дрейфовые кинетические неустойчивости из-за больших градиентов давления плазмы могут развиваться при с < с*. Физическая природа дрейфовых неустойчивости состоит в том, что при с < с* волна обладает отрицательной энергией, в результате чего бес-столкновительная диссипация при dF / dW < 0 приводит к возбуждению неустойчивости. Дрейфовые неустойчивости могут развиваться либо на внутренней, либо внешней кромке кольцевого тока. Скорости распространения колебаний в азимутальном направлении определяются скоростью ларморовского дрейфа плазмы, порядка первых км/c, и определяются знаком и величиной градиента давления горячей компоненты плазмы. Резонансное возбуждение частицами квазимонохроматических колебаний диапазона Рс4-5 имеет место в фазу восстановления магнитной бури, когда неравновесные распределения горячих протонов медленно релаксируют к

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

Хотя теория дрейфовых неустойчиво-стей, вызываемых неоднородностями горячей плазмы, активно развивается в физике плазмы, интерпретация наблюдений волн в магнитосферной плазме потребовала разработки новых теоретических моделей, учитывающих ее характерные особенности: конечное давление плазмы, сопоставимое с давлением магнитного поля ( Р 1); многокомпонентность плазмы -наличие холодных и горячих частиц; и анизотропию = /7^,-1 продольной и поперечной температур. Обнаруженные неустойчивости оказались новыми не только для космической геофизики, но и для физики плазмы в целом. Когда меж-модовое зацепление поперечно-мелкомасштабных колебаний из-за неоднородности плазмы оказывается малым, зацепление различных мод в плазме конечного давления определяется эффектами конечного ларморовского радиуса ионов рг. Возмущения с поперечными длинами

волн, большими ларморовского радиуса, т. е. (к^р )2 «1, распадаются на независимые колебания - сжимаемую альвенов-скую волну и дрейфово-зеркальную моду магнитозвукового типа. В альвеновской моде в плазме конечного давления поддерживается баланс поперечного давления плазмы и магнитного поля p + БСБ0 / 4ж ~0.

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

да

K =-да

теоретических моделях рассматривается либо альвеновские колебания, либо более сложные моды в плазме конечного давления, зацепленные с альвеновской волной. Прямые измерения этого масштаба практически невозможны, т. к. все спутниковые измерения во внешней магнитосфере проводятся вблизи экваториальной плоскости магнитосферы. Тем не менее, по-видимому, существуют колебания Рс5 диапазона с малым продольным масштабом, локализованные вблизи вершины силовой линии, которые не могут быть интерпретированы как альвеновские волны. Предположительно, таким колебаниями являются компрессионные Рс5 волны, т. е. волны сжатия магнитного поля. Традиционно эти волны объясняли как проявление дрейфово-зеркальной моды, которая для возбуждения требует большой анизотропии А и больших значений р. Другим возможным объяснением природы компрессионных Рс5 волн является дрейфово-компрессионная мода, которая может возбуждаться в неоднородной плазме конечного давления при менее жестких условиях [КИттИкт, 2007]. Дрейфово-компрессионные моды похожи на зеркально-дрейфовую моду, т. е. они представляют собой медленные магнитные волны сжатия и обладают диамагнитным свойством: противофазными колебаниями плазменного и магнитного давлений.

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

гулярных возмущений типа Р13 и стимулироваться интенсивными магнитосфер-ными токами вдоль силовых линий. Полная картина возможных механизмов возбуждения УНЧ колебаний на разных фазах магнитных бурь активно разрабатывается, но еще не построена.

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

Знание структуры стоячих колебаний силовых линий в магнитосферном альве-новском резонаторе важно для корректной гидромагнитной диагностики, понимания механизма возбуждения колебаний, определения эффективности резонансного взаимодействия волна-частица, и оценки эффектов взаимодействия волн с ионосферой [Ляцкий и Мальцев, 1983]. Структуру альвеновских колебаний можно определить по наблюдениям на геомагнитно-сопряжённых магнитных станциях. При высокой проводимости ионосферы (что заведомо выполняется в дневных условиях) на ионосферных концах колеблющейся силовой линии образуется узел электрического поля Е и смещения плазмы □. Тогда, как качественно показано на рисунке 20 (слева), фундаментальная мода колебаний с пучностью в вершине силовой линии в магнитосферной экваториальной плоскости вызывают наземные пульсации магнитного поля с синфазными вариациями Н (С-Ю) компоненты. При низкой проводимости ионосферы фундаментальная моды дает противофазное соотношение (рис. 20, справа).

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

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

чае генератора напряжения отклик должен быть больше под высокопроводящей ионосферой.

Для разрешения дихотомии

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

уравнения для альвеновских колебаний (1) при условии отражения от ионосферы (16) следует отношение между наземным магнитным откликом в северном (К) и южном (Б) полушариях:

B

( N )

B

(5 )

ZN) ^ - iz; tang,/.) ) Za - ¿К tan(kAL)

(12)

Как вытекает из соотношения (12), асимметрия магнитных откликов зависит

от соотношения между временным масштабом возмущения и локальным альве-

новским периодом TA. В случае нерезонансного низкочастотного возбуждения, когда ® и (°а (или кАЬ □ 1) отношение

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

Поскольку в реальной ионосфере

^н / ^p ~ const, из соотношения (14) следует, что отношение магнитных откликов в сопряженных точках не зависит от ионосферной проводимости, т. е. реализуется режим генератора тока.

Изменение режима магнитосферного генератора при переходе от нерезонансного к резонансному возбуждению аналогично антиотражательному покрытию в оптике, где входное сопротивление многослойной системы резко изменяется, когда длина волны падающего света кратна ширине слоя. Результаты сопряженных наблюдений различных типов УНЧ-возмущений хорошо согласуются с описываемой концепцией [Pilipenko et al., 2021]. Оказалось, что квазипостоянные нерезонансные возмущения, такие как глобальные магнитосферные системы продольных токов и SC на низких широтах, действительно соответствуют режиму генератора напряжения. Такие же магнитосферные явление, как Pc5 волны и конвективные вихри (Traveling Convection Vortices), следует рассматривать как резонансный от-

Резонансное условие (15) означает, что частота волны с учетом допплеровского сдвига совпадает с циклотронной частотой

между наземными магнитными возмущениями в сопряженных точках:

(13)

ния. В случае возбуждения колебаний на резонансной частоте, т. е. когда О ^ Ол (или kAL ^ п / 2 ), отношение магнитных откликов становится равным:

(14)

клик магнитосферных силовых линий, соответствующий режиму генератора тока. Межполушарная симметрия/асимметрия в геомагнитных наземных откликах проявляется в различных магнитосферно-ионосферных явлениях, крупномасштабных токовых системах, авроральных сияниях, УНЧ волнах и импульсах, т. д.

Pel пульсации (ионно-циклотронные волны)

Пульсации диапазона Рс1 ассоциируются с электромагнитными ионно-циклотронными (ЭМИЦ) волнами, дошедшими из магнитосферы до ионосферы и земной поверхности. ЭМИЦ волны наблюдаются в диапазоне частот от долей Гц до первых Гц, и являются характерной чертой внутренней магнитосферы [Kangas et al., 1998]. Они генерируются вблизи магнитосферной экваториальной плоскости в результате циклотронной неустойчивости при резонансном взаимодействии волн и ионов [Трахтенгерц, Райкрофт, 2011]. Дисперсионное соотношение и резонансное условие для ЭМИЦ волн имеет вид

(15)

ионов Ог.. В неподвижной системе отсчета волны распространяются противоположно

B

( N )

B

( 5 )

-(N ) "н

г( 5 )

н

d(N) y(N) у(5)

f5_.д Н ^Р

B

(5)

?(Ю у(S) ^P ^н

О =

УЛл/1 - ^

1+ О/ Q,

+ (о/ — )(k /k)

co-kV Q =0

скорости движения частиц вдоль внешнего поля Во, поскольку а < П;. Пороговое

условие неустойчивости A¡+1) > . Анизотропия температуры A = T±/% -1, необходимая для развития неустойчивости, может возрасти вследствие бетатрон-ного ускорения при инжекции протонов кольцевого тока во внутреннюю магнитосферу, или из-за скачков динамического давления солнечного ветра. ЭМИЦ волны могут генерироваться как ионами из плазменного слоя на открытых дрейфовых траекториях, так и ионами кольцевого тока на замкнутых дрейфовых траекториях. Инкремент нарастания ЭМИЦ волн становится заметным при условии u¡ U VA. Таким

образом, благоприятные условия для развития этой неустойчивости для горячей анизотропной плазмы с ионами с энергией 10-100 кэВ, возникают либо вблизи плаз-мопаузы, где высока концентрация фоновой холодной плазмы, либо за пределами геосинхронной орбиты L>7, где мало маг-нитосферное магнитное поле. Поперечный масштаб волны предположительно определяется режимом квазипараллельного распространения, когда эффективный поперечный волновой вектор k < к" =ка(со / £1).

Обычно считается, что ЭМИЦ неустойчивость протонов кольцевого тока конвективная, т. е. область неустойчивости в вершине силовой линии работает как усилитель бегущих альвеновских волн [Demekhov, 2007]. Сами волновые пакеты при этом осциллируют между сопряженными ионосферами и испытывают усиление при каждом прохождении через экваториальную область магнитосферы (модель bouncing wave packet). Волновые пакеты, попеременно появляющиеся в сопряженных точках, характерны для дискретных Рс1 пульсаций («жемчужины»). Однако зачастую Pel пульсации наблюдаются в виде непрерывного излучения, что свидетельствует об абсолютном режиме неустойчивости, т.е. работе неустойчивости как генератора излучения. Модельные расчеты показывают, что конвективная неустойчивость может перейти в абсо-

лютную при A > 3/2 . Однако экспериментальных подтверждений условий развития разных режимов неустойчивости не найдено.

В многокомпонентной плазме с примесью тяжелых ионов (He+, O+) в той области, где частота пакета равна перекрестной частоте, может происходить инверсия поляризации, и частичное отражение и поглощение волны. Это может приводить к образованию интервалов поглощения в спектральной структуре ЭМИЦ волн. Пример такой структуры со спектральным пробелом вблизи частоты ионов гелия приведен на рисунке 8.

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

Импульсные ночные колебания Pi2-Pi3

Интерес к Pi2 пульсациям связан с тем, что они возникают при резкой активизации магнитосферной активности - начале взрывной фазы суббури (breakup) и по существу являются «маркером», отмечающим переход медленной фазы роста возмущения ночной магнитосферы во взрывную фазу [Olson, 1999]. Pi2 также могут быть откликом на импульсные потоки плазмы в хвосте магнитосферы (bursty bulk flows). Возбуждение Pi2 пульсаций сопровождается интенсификацией полярных сияний, ионосферными и продольными токами, высыпанием энергичных электронов и другими процессами, происходящими во время взрывной фазы суббури [Keiling and Takahashi, 2011].

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

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

спектрального состава высокоширотных Pi2 пока нет.

Pi2 пульсации на низких широтах интерпретируются как проявление глобальных БМЗ колебаний ночной плазмосферы (рис. 21). Эти БМЗ колебания практически без ослабления проникают через ночную ионосферу к земной поверхности и обуславливают синхронный глобальный отклик на низких широтах. Пример таких колебаний, практически синхронных как в верхней ионосфере, так и на сети низкоширотных наземных станций, приведен на рисунке 10.

Рис. 21. Pi2 пульсации на низких широтах как глобальные БМЗ колебаний ночной плазмо-

сферы

Высказывалось немало гипотез о происхождении Pi3/Ps6 пульсаций - квазипериодической последовательности импульсов с квази-периодом ~10-20 минут на фоне суббуревой магнитной бухты. Их объясняли вариациями параметров межпланетной среды, возникновением периодического режима магнитосферно-ионосферного взаимодействия, квазипериодическим режимом пересоединения в токовом слое хвоста магнитосферы, не-устойчивостями горячей компоненты маг-нитосферной плазмы \Keiling, 2009]. Высказывались предположения, что пульсации Pi3 могут быть связаны с собственными колебаниями магнитного хвоста Земли.

В отличие от других типов УНЧ возмущений (пульсации Pc3-5) пульсации Pi3 вряд ли связаны с резонансом силовых линий в магнитосфере, поскольку их период больше, чем период альвеновской моды на авроральных широтах. Хотя высказывались предположения, что столь большой период может быть связан с вытягиванием

силовых линий геомагнитного поля в хвост магнитосферы во время сильных суббуревых возмущений. Возможным источником возбуждения пульсаций Рi3 может быть развитие баллонной неустойчивости. Действительно, баллонная мода, представляющая собой связанную альве-новскую и ММЗ моды, имеет две ветви. Период «альвеновской» ветви примерно совпадает с альвеновским периодом колебаний силовой линии Га, а поляризация этой ветви смешанная (В2~В ППП «Медленная» ветвь баллонных колебаний имеет гораздо больший период 7>>Гл, а ее поляризация соответствует компрессионной моде (В2»В □.

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

ции, а для некоторых нет - также не решен.

Взаимодействие магнитосферных УНЧ волн с ионосферой

Передача энергии из околоземной среды в земную атмосферу происходит в ионосфере, поэтому описание процессов взаимодействия УНЧ волн с ионосферой имеет ключевое значение для солнечно-земной физики. К тому же наземные наблюдения зачастую являются единственным источником информации об УНЧ волновой активности в магнитосфере. Однако даже для горизонтально однородной ионосферы процесс взаимодействия падающей МГД волны с системой ионосфера-атмосфера-Земля довольно сложен - он включает в себя взаимную трансформацию альвеновской и БМЗ волн в анизотропной ионосфере, суперпозицию падающих и отражённых волн, связь с электрической и магнитной модами в атмосфере [Alperovich andFedorov, 2007].

Качественные представления о прохождении МГД волны к земной поверхности можно получить на основе модели «тонкой» ионосферы, в которой ионосфера представляется тонкой пленкой с интегральными холловской и педерсеновской проводимостями, находящейся на высоте h над земной поверхностью [Ляцкий и Мальцев, 1983]. Эта модель хорошо описывает взаимодействие крупномасштабных МГД

мод с системой ионосфера-атмосфера-земля, для которых скин-длина и длина волны в ионосфере много больше толщины проводящего слоя Е (диапазон Рс3-5 и Pi2-3). Для магнитосферных колебаний ионосфера является частично отражающей и поглощающей стенкой, на которой выполняется граничное условие импедансно-го типа Ву = +/л0ЪрЕх ПрИ г = +Ь .

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

Характерные особенности отражения альвеновских волн от тонкой ионосферы видны из выражения для коэффициента отражения:

R = zp -za sinI

Zp +ZA Sin 1

(16)

где I - наклонение геомагнитного поля (для вертикального поля I = л /2 ). Альве-новские волны хорошо отражаются от ионосферы при сильном контрасте между ионосферной проводимостью Пр и волновой проводимостью Па. Это условие хорошо выполняется для дневной магнито-

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

у 5 1

— = 5m =-ln| R |

( a жп

(17)

где a - поперечный размер неоднородности альвеновской скорости в магнитосфере, п - номер гармоники колебаний.

При прохождении через ионосферу происходит расплывание резонансной области в соответствии с соотношением

8 ^ 8т + Ъ. Наземный отклик в X-компоненте на тороидальные альвенов-ские колебания с амплитудой Б((п) над ионосферой дается формулой:

у

Bx = B(ym) уН exp (-kh)

(18)

Здесь k - горизонтальный волновой вектор возмущения. Соотношение (18) показывает, что мелкомасштабные возмущения с масштабом меньше, чем высота ионосферы (^>>1), экранируются от наземных магнитометров. Для альвеновских волн с очень малыми поперечными масштабами (<2 км) следует дополнительно учитывать

сильное поглощение из-за столкновитель-ной диссипации в ионосфере.

Взаимодействие БМЗ волны с системой ионосфера-атмосфера-Земля принципиально отличатся от взаимодействия альве-новской волны. Соотношение между возмущением горизонтальной магнитной компоненты над ионосферой и наземным откликом имеет вид:

B

(m)

B

= 1 - ip

coh

р = —

К

(19)

где Кс = (^с) - «каулинговская» скорость в ионосфере, определяемая кау-линговской проводимостью

2С = Ър + Ъ2Н / Ър. Параметр р контролирует прохождение БМЗ моды через ионосферу к земной поверхности: при p<<1 ионосфера прозрачна для БМЗ волн (что заведомо выполняется для ночной ионосферы), а ионосферное поглощение проявляется при p~1 в дневные часы. Поворота эллипса поляризации при прохождении БМЗ через ионосферу не происходит.

Для более высокочастотных пульсаций диапазона Рс1 приближение тонкой ионосферы становится неприменимым. Характерной особенностью ионосферы в этом частотном диапазоне является наличие волновода и резонатора, которые могут накапливать энергию электромагнитных волн в диапазоне от долей Гц до нескольких десятков Гц: ионосферный альвенов-ский резонатор (ИАР) и ионосферный волновод для БМЗ волн \Demekhov, 2012]. ИАР образуется за счет отражения альве-новских волн от резкого градиента вертикального профиля альвеновской скорости VA(z) на высотах около 103 км, где длина

волны становится сопоставимой с размером неоднородности, т. е. происходит нарушение приближения геометрической оптики. Нижней границе ИАР служит проводящий Е-слой ионосферы. Наземные проявления ИАР - множественные спектральные полосы, разделенные по частоте на □/-0,3-0,5 Гц, регулярно наблюдаются в ночное время на низких, средних, и даже авроральных широтах. Резонансные свойства ионосферы обеспечивают зависимость свойств пропускания/отражения Pc1 пульсаций от частоты. Таким образом, ионосфера с учетом наличия ИАР может действовать как многодиапазонный полосовой фильтр для магнитосферных ЭМИЦ волн, формируя спектральную структуру излучений Pc1, наблюдаемых на Земле.

Волновод для БМЗ волн образуется в области минимума альвеновской скорости в F-слое ионосферы. В результате магни-тосферная ЭМИЦ волна внутри ионосферы частично преобразуется в волноводную моду, которая переносится на тысячи километров относительно места падения волны на ионосферу. Волноводные моды с частотой ниже критической частоты вол-

новода (~0.5 Гц) не могут распространяться вдоль ионосферы. Горизонтальное распространение в волноводе наиболее эффективно вдоль магнитного меридиана. Таким образом, ионосфера функционирует как огромная «приемная антенна», собирающая магнитосферные ЭМИЦ сигналы с большой области вокруг наземного наблюдателя. Наземные наблюдения с использованием сети магнитометров показали, что амплитуда волны Pc1, распространяющиеся из области инжекции, спадала определенным образом: затухание было очень быстрым (~10 дБ/100км) в центре инжекции, но становилось намного слабее (~2.5 дБ/100км) в области за пределами ~500 км от центра.

Мелкомасштабные УНЧ волны

Спутниковые наблюдения с высоким пространственно-временным разрешением во время различных геофизических процессов показывают, что большая плотность мощности электромагнитных возмущений имеет тенденцию концентрироваться в малых временных и пространственных масштабах [Pilipenko and Enge-bretson, 2002]. Таким образом, импульсные и мелкомасштабные процессы, несмотря на их относительно непродолжительность и малый масштаб, могут давать суще-

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

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

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

Рис. 22. Механизм формирования альвеновской авроры в силовой трубке с AAR (auroral

acceleration region).

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

^ =а>А= к^л , имеем:

- в кинетическом пределе

а = аА^1 1 + ;

а.

- в инерционном пределе

с =

где К = с / аре - инерционная электронная длина, а - электронная плазменная частота.

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

объединены в одно а □ аА^ 1 ± к^р^ , где

Ра ^ ~ обобщенный дисперсион-

ный радиус.

В отличие от альвеновских волн в приближении идеальной МГД, мелкомасштабные структуры с масштабами, сравнимыми с радиусом дисперсии р, обладают собственным продольным электри-

ческим полем

\E,\={klPdf{kJk1)\E1\ .

Ей

следующего вида

Благодаря

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

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

77*

поля -hj, ускорению и высыпанию электронов. Значительное нерезиститвное падение потенциала вдоль силовых линий может создаваться пробочным механизмом (mirror forces). Пробочный механизм обуславливает появление продольного электрического поля Е0 за счет разницы питч-угловых распределений электронов и ионов квази-нейтральной плазмы. В кинетической теории удается рассчитать самосогласованное стационарное распределение потенциала U(s) и продольного тока

№ вдоль силовой линии. В общем случае, j является функционалом от ü(s), т.

е. h =/(G(s)). Однако в широком диапазоне параметров плазмы оказывается справедливым линейное нелокальное соотношение между вытекающим из ионосферы током j\ и падением потенциала □ □ □ □ (формула Knight) [Vogt, 2002]

АФ = Qj\

или j= К АФ

(20)

Продольную проводимость силовой трубки K=1/Q можно оценить как

К □ Ые1 / тае. Благодаря возникающему продольному электрическому полю авро-ральные электроны могут переносить существенные токи между магнитосферой и нижней ионосферой даже при очень малых значениях конуса потерь. Хотя ряд концептуальных проблем продольных электрических полей остается не до конца выясненным, нелокальный закон Ома (20) широко используется в физике аврораль-ной плазмы. По существу, возникновение продольной разности потенциала является кинетическим эффектом, однако в рамках

МГД приближения можно полагать, что существует нелокальная вольтамперная характеристика, следующая из кинетического рассмотрения и связывающая вытекающий ток Л с падением потенциала АФ.

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

Наличие области ускорения авроральных частиц (auroral acceleration region - AAR) -локализованной по высоте области со значительным падением электрического потенциала вдоль силовых линий, является важной особенностью системы магнитосфера-ионосфера на авроральных широтах. Наличие AAR видоизменяет свойства альвеновской волны, в частности приводит к появлению у нее продольного электрического поля Ег □ jn(Q / d), где d - толщина AAR области. Предположительно, этот механизм ответственен за образование «альвеновской авроры» - ускорение электронов и уярчение полярных сияний под действием магнитосферных УНЧ волн.

Наиболее эффективно конверсия энергии альвеновских волн в ускорение аврораль-ных электронов происходит при поперечных масштабах волны k± XA ~1, где параметр ХА = ^JQZa . При типичных параметрах авроральной области масштаб Xa составляет первые десятки км. Если квазистационарные продольные поля приводят к формированию моноэнергетических пучков высыпающихся электронов (т. н. inverted V структуры), то в «альвеновской авроре» ускорение и высыпание захватывает большой диапазон энергий, как схематично показано на рис. 23.

Рис. 23. Энергетические спектры авроральных электронов в области альвеновской авроры (слева) и в области квазистационарных продольных полей (inverted-V структуры)

УНЧ волны как активный фактор космической погоды

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

Геомагнитные УНЧ волны в широком интервале пространственных и временных масштабов традиционно используются в качестве одного из основных средств «мониторинга космической погоды» и «диагностики околоземного пространства» [Troitskaya and Guglielmi, 1967]. Однако, все больше появляется экс-

периментальных свидетельств того, что УНЧ возмущения, заполняющие околоземное пространство, не только индикатор, но и активный фактор воздействия на околоземную среду и на наземные и спутниковые технологические системы. Резкие изменения геомагнитного поля могут приводить к сбоям в работе энергетических и транспортных систем, вызванных геоин-дуцированными токами (ГИТ), которые индуцируются в протяженных заземленных системах. При этом, хотя космическая погода активирует глобальные электромагнитные и плазменные процессы в околоземной среде, однако наибольший риск ГИТ связан не непосредственно с крупномасштабными процессами с огромным выделением энергии, а с более локальными и быстрыми процессами [Пилипенко, 2021]. Оказалось, что электрические сети наиболее чувствительны к помехам от естественных геоэлектрических полей, которые вызываются геомагнитными пульсациями Pi3/Ps6 с периодами порядка первых мин - десятков мин. Именно они вызывают наиболее интенсивные всплески ГИТ с размахом до сотен А в ЛЭП и способствуют возникновению ГИТ в протяженных инженерных системах, таких как железные дороги и

трубопроводы [Belakhovsky et al., 2019]. В качестве примера, на рис. 24 приведено сопоставление ГИТ в ЛЭП на Кольском полуострове с магнитными вариациями на близкой магнитной станции во время магнитной бури 27-28 мая 2017 г. Видно, что возбуждение всплесков ГИТ вызвано воздействием Pi3 пульсаций, наложенных на магнитную бухту, во время развития магнитной бури.

27-28 MAY 2017

21 22 23 00 01 02 03 04 05 06 07

lt

Рис. 24. Возбуждение ГИТ в ЛЭП «Северный Транзит» геомагнитными пульсациями Pi3. Сверху вниз показаны вариации индекса SYM-H, характеризующего интенсивность магнитной бури, магнитограммы горизонтальных магнитных компонент на близкой станции IVA, вариабельность магнитного поля dX/dt и dY/dt на той же станции, и вариации ГИТ на

подстанции «Выходной».

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

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

Динамика УНЧ волн в земной магнитосфере тесно связана с динамикой частиц. Магнитосферные УНЧ волны могут эффективно модулировать потоки частиц, и приводить к модулированному высыпанию в ионосферу, создающему помехи в транс-ионосферной радиосвязи. УНЧ волны могут модулировать полное электронное содержание (ПЭС) ионосферы, созда-

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

УНЧ волны и динамика радиационного пояса

Электродинамика околоземного пространства определяется солнечным ветром и ММП, но ускорение магнитосферных электронов до релятивистских энергий и усиление внешнего радиационного пояса Земли обусловлено внутримагнитосфер-ными процессами. В бесстолкновительной космической плазме поток солнечного ветра не взаимодействует напрямую с захваченными электронами магнитосферы, и электромагнитные поля играют роль посредника, передающего энергию от частиц низкой энергии небольшой группе частиц высокой энергии. Динамика внешнего радиационного пояса определяется сложным действием нескольких механизмов ускорения, переноса и потерь энергичных частиц. Главная фаза сильной магнитной бури сопровождается быстрым падением потока МэВ-электронов, за которым следует медленное адиабатичсекое восстановление до предбуревого уровня. Однако почти в 50% бурь электронные потоки нарастают на 2-3 порядка, далеко за пределы адиабатического уровня. Усиление потоков релятивистских электронов происходит сначала быстро в сердцевине пояса за несколько часов (Ь~3-4), а затем медленно (за дни-недели) на больших магнитных оболочках (Ь>6). Механизм ускорения магнитосфер-ных электронов до релятивистских энергий до сих пор остается нерешенной проблемой космофизики.

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

гии, вызвать сбои в электронике, дифференциальную зарядку внешних и внутренних диэлектриков с последующим пробоем, повреждение панелей солнечных батарей и др. Серьезные аномалии в работе магнитосферных космических аппаратов наблюдались сразу после резкого подъема потока релятивистских магнитосферных электронов, в связи с чем эти электроны даже получили название «электроны-убийцы» (killer electrons). Интересно, что увеличение интенсивности потоков «электронов-убийц» тесно следует за увеличением волновой активности в магнитосфере, характеризуемой «волновым индексом», пропорциональным интегральной спектральной плотности магнитных вариаций в Рс5 диапазоне на наземных станциях [Pilipenko et al., 2017]. Интенсивный рост потоков релятивистских электронов хорошо соответствовал изменениям глобального УНЧ индекса, так что этот индекс гораздо лучше характеризует динамику электронов, чем любой другой геомагнитный индекс. Эта закономерность позволила высказать гипотезу, что именно Рс5 волны ускоряют энергичные электроны до релятивистских энергий на восстановительной фазе магнитных бурь [Mann et al., 2013].

Если быструю фазу ускорения электронов до релятивистских энергий естественно связать с развитием суббури, то медленную фазу связывают с ускорением диффузионного типа и радиальным переносом захваченных частиц к Земле, вызванными УНЧ активностью [Shprits et al., 2008]. УНЧ волны могут быть возможным посредником, передающим энергию от высокоскоростных потоков солнечного ветра электронам магнитосферы. Дрейфовый резонанс между дрейфующими электронами и азимутально распространяющимися волнами возникает, когда азимутальная фазовая скорость D/m совпадает с угловой скоростью электрона Dd (своеобразный «геосерфотрон»), т. е. при с / m = с . Качественно работа магнито-

сферного «геосерфотрона» показана на рисунке 25. Скорость магнитного дрейфа электронов с энергией W и период обра-

щения Td можно оценить по соотношениям

Ff[km/s] □ 8.5Z2r[MeV]g(a) 7j[min] □

44

-g(ce) g(a) = 0.7 + 0.3sina /91 \ 1Г[МеУ]

где □ - питч-угол частицы на экваторе. могут резонансно взаимодействовать с Из этого соотношения следует, что ~1 мэВ волнами в диапазоне Рс5 (5-10 мин). электроны, дрейфующие на оболочке L~6,

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

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

УНЧ волны могут приводить не только к ускорению электронов радиационного пояса, но и вызывать их сброс в атмосферу. ЭМИЦ волны диапазона Рс1 могут резонансно рассеивать по питч-углам электроны релятивистских энергий и быть ответственными за опустошение радиационных пояса Земли. ЭМИЦ излучения также ограничивают время жизни энергичных протонов кольцевого тока и вызывают локализованное протонное сияние за счет резонансного воздействия на ~100 кэВ протоны. Однако волновой пакет типичных пульсаций Рс1 имеет малые масштабы в поперечном направлении, поэтому насколько заметным может оказаться его воздействие на радиационный пояс или кольцевой ток в целом - остается невыясненным.

Возмущения геомагнитного поля антропогенного происхождения

Детальное знание морфологии и природы естественных вариаций геомагнитного

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

- Антропогенные возмущения, вызванные взрывами и стартами ракет.

- Отклики на метеоявления (ураганы, тайфуны, грозы);

- Сейсмо-электромагнитные шумы и импульсы, связанные с процессами в зонах разрушения горной породы во время подготовки землетрясений [Surkov and Hayakawa, 2014].

Все эти процессы связаны с выделением большого количества энергии и затрагивают все геофизические оболочки (литосферу, атмосферу, ионосферу и магнитосферу) и геофизические поля [Pokhotelov et al., 1995]. В годы холодной войны именно появление специфических геомагнитных пульсаций диапазона Рс1 оказалось наиболее надежным индикатором для дистанционного обнаружения высотных ядерных взрывов в атмосфере и ионосфере.

В конкретной ситуации могут реализоваться различные механизмы воздействия взрывов и землетрясений на ионосферу и геомагнитное поле. Распространяющиеся от места взрыва или землетрясения сейсмические волны благодаря индукционному эффекту порождают сопутствующий электромагнитный сигнал [Surkov and Pili-penko, 1997]. Акустическая волна от взрыва или землетрясения, достигающая ионосферных высот, индуцирует токи в проводящем слое ионосферы (Е-слой). Система ионосферных токов оказывается незамкнутой и порождает продольные токи вдоль силовых линий геомагнитного поля на фронте акустического импульса. Импульс продольного тока может быть зарегистрирован на низкоорбитальном спутнике как всплеск возмущения поперечной компоненты геомагнитного поля. Кроме того, индуцируемые акустической волной взрыва или землетрясения в ионосфере токи возбуждают распространяющиеся вдоль ионосферы МГД волны.

Электромагнитные УНЧ эффекты сопровождают и запуски мощных ракет. Продукты сгорания от факела ракетного двигателя подобно снежному плугу расталкивают плазму в нижних слоях ионосферы и генерируют интенсивные токи и гидромагнитные шумы. Часть этих шумов, с частотами более 1 Гц, захватывается в ионосферный БМЗ волновод и распространяется на значительные расстояния от места старта. Такие возмущения диапазона Рс1 наблюдались на сети индукционных магнитометров после запусков межконтинентальных ракет.

Заключение

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

вые системы передачи информации. Это времени за «электромагнитным дыхани-позволит в ближайшее время создать гло- ем» Земли. бальные системы мониторинга в реальном

Библиографический список

1. Гульельми А.В., Троицкая В.А. Геомагнитные пульсации и диагностика магнитосферы. - М.: Наука, 1973. - 208 с.

2. Гульельми А.В. МГД волны в околоземной плазме. - М., Наука, 1979. - 121 c.

Леонович А.С., Мазур В.А., Линейная теория МГД колебаний в магнитосфере. - М.,

Физматлит, 2016.

3. Ляцкий В.Б., Мальцев Ю.П., Магнитосферно-ионосферное взаимодействие. - М., Наука, 1983.

4. Нишида А. Геомагнитный диагноз магнитосферы. - М.: Мир, 1980. - 299 с.

5. Пилипенко В.А. МГД волны в космосе и на Земле: Исторический аспект // Вестник ОНЗ РАН. - 2021. - № 13. - NZ2001. - [Электронный ресурс]. - Режим доступа: https://doi.org/10.2205/2021NZ000371.

6. Пилипенко В.А. Воздействие космической погоды на наземные технологические системы // Солнечно-земная физика. - 2021. - Т. 7, №3. - С. 73-110. - [Электронный ресурс]. - Режим доступа: https://doi.org/10.12737/szf-73202106.

7. Пудовкин М.И., Распопов О.М., Клейменова Н.Г. Возмущения электромагнитного поля Земли. Часть 2. Короткопериодные колебания геомагнитного поля. - Л.: ЛГУ, 1976. -271 с.

8. Пулинец С.А., Ведешин Л.А. Космос как природная плазменная лаборатория // Земля и Вселенная. - 2023. - №6. - С. 35-53. - DOI: 10.7868/S0044394823060038.

9. Трахтенгерц В.Ю. Райкрофт М.Дж., Свистовые и альфвеновские циклотронные мазеры в космосе. - М.: Физматлит, 2011.

10. Федоров Е.Н., Пилипенко В.А., Курнева Н.А., Баранский Л.Н. Магнитотеллуриче-ское зондирование и гидромагнитная диагностика магнитосферы // Исслед. по геомагнетизму, аэрономии и физике Солнца. - 1992. - Вып. 98. - С. 49-79.

11. Яхнин А.Г., Т.А. Яхнина, Н.В. Семенова, Т.А. Попова, А.Г. Демехов, Протонные сияния к экватору от овала как проявление ионно-циклотронной неустойчивости в магнитосфере Земли (краткий обзор) // Геомагнетизм и Аэрономия. - 2018. - Т. 58, № 5. -С. 599-607.

12. Alperovich L.S. Fedorov, E.N. (2007). Hydromagnetic Waves in the Magnetosphere and the Ionosphere // Series: Astrophysics and Space Science Library, Vol. 353, XXIV, 418p., Hardcover, ISBN: 978-1-4020-6636-8.

13. Archer M.O., Pilipenko V.A., Li B., Sorathia K., Nakariakov V.M., Elsden T. Nykyri K. (2024), Magnetopause MHD surface wave theory: progress & challenges. Front. Astron. Space Sci. 11:1407172. DOI: 10.3389/fspas.2024.1407172.

14. Belakhovsky V., Pilipenko V., Engebretson M., Sakharov Ya. Selivanov V. Impulsive disturbances of the geomagnetic field as a cause of induced currents of electric power lines // J. Space Weather and Space Climate. 9. A18. 2019.

15. Demekhov A.G. Recent progress in understanding Pc1 pearl formation // J. Atmos. Solar-Terr. Phys. 69. 1609-1622. 2007.

16. Demekhov A.G., Coupling at the atmosphere-ionosphere-magnetosphere interface and resonant phenomena in the ULF range, Space Science Review (2012) 168:595-609 DOI 10.1007/s11214-011-9832-6.

17. Guglielmi A.V. Pokhotelov O.A., Geoelectromagnetic Waves, Institute of Physics, Bristol and Philadelphia, 1996. - 382 p.

18. Jacobs J.A. Geomagnetic micropulsations. Berlin-Heidelberg-New York, 1970. - 179 p.

19. Kangas J., Guglielmi A., Pokhotelov O. Morphology and physics of short-period magnetic pulsations // Space Science Reviews. - 1998. - Vol. 83. - P. 435-512.

20. Keiling A., Alfven waves and their roles in the dynamics of the Earth's magnetotail: A Review, Space Science Review. - 2009. - № 142. - P. 73-156. - DOI:10.1007/s11214-008-9463-8

21. Keiling A., Takahashi K. Review of Pi2 models // Space Sci Rev. - 2011. - V. 161. P. 63148.

22. Klimushkin D.Yu. How energetic particles construct and destroy poloidal high-m Alfven waves in the magnetosphere // Planet. Space Sci. - 2007. - Vol. 55. - P. 722-730.

23. Mann I., Lee E., Claudepierre S. et al. Discovery of the action of a geophysical synchrotron in the Earth's Van Allen radiation belts // Nature Communication. - 2013. - № 4. - P. 2795.

- URL: https://doi.org/10.1038/ncomms3795.

24. Menk F.W., Waters C.L., Magnetoseismology: Ground based remote sensing of the Earth magnetosphere, Wiley, 2013.

25. Nakariakov V.M., V.A. Pilipenko, B. Heilig, et al. Magnetohydrodynamic oscillations in the solar corona and Earth's magnetosphere: Towards consolidated understanding, Space Science Review. - 2016. - №200. - P. 75-203.

26. Olson J.V. Pi2 pulsations and substorm onsets: A review // J. Geophys. Res. - 1999. -V. 104. - P. 17499-17520.

27. Pilipenko V.A. ULF waves on the ground and in space // J. Atmos. Terr. Phys. - 1990. -Vol. 52, №12. - P. 1193-1209.

28. Pilipenko V.A., Kozyreva O.V., Engebretson M.J., Soloviev A.A. ULF wave power index for the space weather and geophysical applications: A review, Russ. J. Earth Science. - 2017. -Vol. 7, № 2. - ES1004. - DOI:10.2205/2017ES000597.

29. Pilipenko V.A., Engebretson M.J., Hartinger M.D., Fedorov E.N., Coyle S. Electromagnetic fields of magnetospheric disturbances in the conjugate ionospheres: Current/voltage dichotomy, in "Cross-Scale Coupling and Energy Transfer in the Magnetosphere-Ionosphere-Thermosphere System", ed. by T. Nishimura, Elsevier B.V., Amsterdam, 2021, ISBN: 9780128213667.

30. Pilipenko V.A., Fedorov E.N. Magnetotelluric sounding of the crust and hydromagnetic monitoring of the magnetosphere with the use of ULF waves // Ann. Geofisica. - 1993. -Vol. 36, №5-6. - P. 19-33.

31. Pilipenko V., Engebretson M. Ground images at high latitudes of ULF wave processes in the outer magnetosphere // J. Atmosph. Solar-Terrestrial Physics. - 2002. - Vol. 64. - P. 183201.

32. Pokhotelov O.A., Parrot M., Pilipenko V.A., et al. Response of the ionosphere to natural and man-made acoustic sources // Ann. Geophys. - 1995. - Vol. 13. - P. 10197-10210.

33. Saito T. Geomagnetic pulsations // Space Sci.Rev. - 1969. - Vol. 10. - Pp. 319-412.

34. Shprits Y.Y., et al. Review of modeling of losses and sources of relativistic electrons in the outer radiation belt I: Radial transport, J. Atmospheric and Solar-Terrestrial Physics. - 2008.

- Vol. 70. - P. 1679-1693.

35. Solar Wind Sources of Magnetospheric Ultra-Low-Frequency Waves. ed. Engebretson M.J., Takanashi K., Scholer M. Geophys. Monogr. Ser., V. 81, AGU, Washington, D.C., 1994.

36. Surkov V., Pilipenko V. Magnetic effects due to earthquakes and underground explosions: a review// Ann. Geofis. - 1997. - Vol. 40, № 2. - P. 1-13.

37. Surkov V., Hayakawa M. (2014) Ultra and Extremely Low Frequency Electromagnetic Fields. Springer Geophysics, London. DOI: 10.1007/978-4-431-54367-1.

38. Vogt J., Alfven wave coupling in the auroral current circuit, Surveys of Geophysics 23, 335-377, 2002.

39. Walker A.D.M. Plasma Waves in the Magnetosphere, series Physics and Chemistry in Space, Springer Berlin, Heidelberg, 1993. - URL: https://doi.org/10.1007/978-3-642-77867-4.

40. Yumoto K. Generation and propagation mechanisms of low-latitude magnetic pulsation -A review // J. Geophys. Res. - 1986. -V. 60. - P. 79-105.

ULTRA-LOW FREQUENCY WAVES IN SPACE AND ON EARTH

V.A. Pilipenko, Head of Laboratory1, Chief Researcher2 D.D. Pozdnyakova1, Engineer N.V. Savelieva1, Junior Researcher

institute of Physics of the Earth of the Russian Academy of Sciences 2Geophysical Center of the Russian Academy of Sciences (Russia, Moscow)

Abstract. A brief overview of modern ideas about the physical nature of wave variations of the geomagnetic field in the ultra-low frequency (UHF) range from fractions of MHz to the first Hz is given. The noise and fluctuations of the UHF range reflect the turbulent nature of the interaction of solar wind plasma with the Earth's magnetosphere. The excitation of some types of UHF waves is a manifestation of the resonant interaction of waves and particles in near-Earth plasma. The physics of UHF waves is based on the theoretical apparatus of plasma physics and magnetic hydrodynamics, and the experimental research base is the data from observations on satellites and ground stations. The morphological characteristics of the main types of geomagnetic pulsations are given, and then the physical mechanisms of their generation and the possible role in the course of plasma processes in the near-Earth environment are considered. Geomagnetic UHF waves turn out to be not only one of the main diagnostic tools for near-Earth space, but also an active factor affecting both the near-Earth environment and terrestrial and satellite technological systems. The presentation is intended for readers working in related fields (space technology, geophysics) who would like to get a general idea of the physics of UHF waves in near-Earth space.

Keywords: Earth's magnetosphere, UNH.