Геомагнетизм и аэрономия, 2021, T. 61, № 3, стр. 309-315

Ионный форшок околоземной ударной волны: теоретическая модель и наблюдательные данные

Г. Н. Кичигин *

Институт солнечно-земной физики СО РАН (ИСЗФ СО РАН)
г. Иркутск, Россия

* E-mail: king@iszf.irk.ru

Поступила в редакцию 20.07.2020
После доработки 02.11.2020
Принята к публикации 28.01.2021

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

Аннотация

Рассматриваются процессы, связанные с образованием перед фронтом околоземной бесстолкновительной ударной волны так называемой форшок-области. В рамках предлагаемой теоретической модели показано, что в рампе фронта ударной волны происходит ускорение ионов до значительных энергий. Выясняются физические условия, при которых эти ускоренные в ударном фронте ионы в итоге оказываются перед фронтом ударной волны и удаляются от него, вследствие чего их иногда относят к категории отраженных ионов. Делается вывод о том, что именно эта популяция отраженных энергичных ионов (чаще всего их называют “продольными пучками” – field-aligned beams) играет основную роль в формировании границы ионного форшока. Основные свойства продольных пучков, следующие из рассматриваемой модели, сравниваются с наблюдательными данными, полученными с помощью космических аппаратов.

1. ВВЕДЕНИЕ

С помощью космических аппаратов (КА) сейчас твердо установлено [Paschmann et al., 1981; Bonifazi and Moreno, 1981; Eastwood et al., 2005; Oka et al., 2005; Balogh and Treumann, 2013], что перед фронтом околоземной бесстолкновительной ударной волны (БУВ) существует так называемая форшок-область (foreshock region). Наблюдения показывают, что существует две границы, которые отделяют от невозмущенной плазмы солнечного ветра (СВ), соответственно, области электронного и ионного форшоков. Схематически в плоскости эклиптики картина форшока представлена на рис. 1 в случае, когда угол между нормалью к фронту БУВ в подсолнечной точке и вектором межпланетного магнитного поля (ММП) θBn близок к 90°, т.е. силовые линии ММП почти перпендикулярны линии Солнце–Земля. На рисунке 1 граница электронного форшока практически совпадает с касательной к фронту БУВ линии ММП (на рис. 1 – это точка “a”). По обе стороны относительно точки “а” наблюдается две области ионного форшока (на рис. 1 – заштрихованы). Одна из областей располагается на утренней стороне магнитосферы (начало – в точке “b”), другая – в вечернем секторе (начало – в точке “с”), как это показано на рис. 1. Граница ионного форшока на утренней стороне магнитосферы – это линия (жирная), начинающаяся в точке “b” фронта БУВ. Для наиболее типичного случая, когда θBn ≈ 45°, наибольший интерес представляет утренний ионный форшок, так как он со стороны Солнца занимает значительную часть области перед ударной волной.

Рис. 1.

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

Мы ограничимся рассмотрением процессов, происходящих в утреннем ионном форшоке. Здесь, как и в вечернем форшоке, регистрируются колебания и волны с достаточно большой амплитудой, а также наблюдаются ионы с энергиями от единиц до сотен кэВ [Paschmann et al., 1981; Bonifazi and Moreno, 1981; Fuselier, 1994; Eastwood et al., 2005; Oka et al., 2005; Balogh and Treumann, 2013; Burgess et al., 2012]. По виду энергетического распределения ионы в форшоке условно подразделяют на 5 групп [Oka et al., 2005]: (1) – пучки, движущиеся вдоль границы форшока, которые называют или “отраженными ионами”, или “продольными пучками” (FAB), (2) – диффузные, (3) – промежуточные (с признаками первых двух), (4) – гировращающиеся, (5) – ионы, сгруппированные в пакеты с одинаковой гирофазой. Из этих пяти групп только происхождение гировращающихся ионов имеет достаточно простое объяснение – это ионы, отраженные от скачка потенциала, существующего во фронте БУВ. Эти отраженные ионы совершают вращение в ММП перед фронтом БУВ и создают подножье БУВ. По поводу механизма образования всех остальных групп энергичных ионов существует множество предположений и моделей, но до сих пор общепринятого мнения о механизмах их возникновения нет. Весьма актуальная проблема – происхождение продольных пучков в форшоке – вызывает наибольший интерес, ей посвящено много публикаций, в том числе и данная работа.

2. СВОЙСТВА ПРОДОЛЬНЫХ ПУЧКОВ, СЛЕДУЮЩИЕ ИЗ НАБЛЮДЕНИЙ

Свойства пучков ионов, движущихся вдоль границы ионного форшока перед фронтом квазиперпендикулярной ударной волны (45° < θBn < 90°) достаточно хорошо изучены [Paschmann et al., 1981; Bonifazi and Moreno, 1981; Fuselier, 1994; Eastwood et al., 2005; Oka et al., 2005; Burgess et al., 2012; Balogh and Treumann, 2013]. Мы их будем называть “продольные пучки”. Приведем здесь наблюдательные данные о типичных свойствах и параметрах этой популяции энергичных ионов: (1) – для типичного диапазона параметров околоземной БУВ продольные пучки характеризуются относительно узким спектром ионов по энергии и наблюдаются преимущественно в областях ударного фронта, где величина угла θBn лежит в интервале от 60° до 85°, т.е. в тех областях, где БУВ считается квазиперпендикулярной; (2) – cвое движение они начинают из некоторой локальной области фронта квазиперпендикулярной БУВ и продолжают далее остронаправленно двигаться вдоль границы ионного форшока; (3) – измеренная скорость движения продольного пучка колеблется в пределах 3–5 U, где U – скорость СВ; (4) – типичный энергетический спектр протонов пучка имеет максимум в районе 5 кэВ или больше, температура их 100–700 эВ, полная ширина спектра до 30–40 кэВ; (5) концентрация продольных протонных пучков составляет менее одного процента от концентрации протонов в плазме СВ.

3. ТЕОРЕТИЧЕСКАЯ МОДЕЛЬ ПРОИСХОЖДЕНИЯ ПРОДОЛЬНЫХ ПУЧКОВ

Все перечисленные выше наблюдательные свойства продольных пучков убедительно объясняются в работах [Кичигин, 2009, 2018]. В этих работах предложена простая модель ударного фронта, учитывающая величину перепада потенциала и его пространственный размер в рампе фронта БУВ. Эта модель, схематически представлена на рис. 2 и используется в данной работе. Хотя принятая нами модель довольно проста, тем не менее, она претендует на описание структуры околоземной БУВ, но главное достоинство модели – это то, что на ее основе получены почти все основные свойства продольных пучков. Одно из допущений модели – это одномерное приближение для каждого небольшого участка фронта вдоль поверхности ударной волны. Такой упрощенный подход обосновывается тем, что пространственный размер фронта ударной волны бесконечно мал по сравнению с радиусом кривизны околоземной БУВ, имеющей примерно форму параболоида. В модели также предполагается, что в пространственной области фронта, где расположен перепад потенциала, влияние магнитного поля на динамику ионов несущественно.

Рис. 2.

Модель структуры ударной волны, представленная в системе отсчета волны. Передняя плоскость рампа – это плоскость y0z, n – вектор нормали к рампу. Поток плазмы набегает перпендикулярно к плоскости y0z со скоростью u. Вектор магнитного поля В направлен под углом θ к нормали n и имеет компоненты Bx и Bz. Перепад потенциала сосредоточен в рампе шириной d, а электрическое поле на отрезке 0 ≤ xd направлено навстречу набегающему потоку плазмы и линейно нарастает от нуля до амплитудного значения Е0. Конвекционное стационарное электрическое поле Ey = uBz/с направлено вдоль оси y (перпендикулярно плоскости рисунка).

В принятой модели, в системе отсчета волны поток плазмы набегает по нормали к плоскости фронта (плоскость y0z) со скоростью u. Вектор магнитного поля В направлен под углом θ к нормали и имеет компоненты Bx и Bz. Перепад потенциала сосредоточен в рампе шириной d, а электрическое поле на отрезке 0 ≤ xd направлено навстречу набегающему потоку плазмы и линейно нарастает от нуля до амплитудного значения Е0. Конвекционное стационарное электрическое поле Ey = uBz/с направлено вдоль оси y.

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

(1)
${{d{{{v}}_{x}}} \mathord{\left/ {\vphantom {{d{{{v}}_{x}}} {dt}}} \right. \kern-0em} {dt}} = {{q{{Е}_{х}}} \mathord{\left/ {\vphantom {{q{{Е}_{х}}} m}} \right. \kern-0em} m} + {{q{{{v}}_{y}}{{B}_{z}}} \mathord{\left/ {\vphantom {{q{{{v}}_{y}}{{B}_{z}}} {mc,}}} \right. \kern-0em} {mc,}}$
(2)
${{d{{{v}}_{y}}} \mathord{\left/ {\vphantom {{d{{{v}}_{y}}} {dt}}} \right. \kern-0em} {dt}} = {{q{{E}_{y}}} \mathord{\left/ {\vphantom {{q{{E}_{y}}} m}} \right. \kern-0em} m} + {{q({{{v}}_{z}}{{B}_{x}} - {{{v}}_{x}}{{B}_{z}})} \mathord{\left/ {\vphantom {{q({{{v}}_{z}}{{B}_{x}} - {{{v}}_{x}}{{B}_{z}})} {mc,}}} \right. \kern-0em} {mc,}}$
(3)
${{d{{{v}}_{z}}} \mathord{\left/ {\vphantom {{d{{{v}}_{z}}} {dt}}} \right. \kern-0em} {dt}} = {{ - q{{{v}}_{y}}{{B}_{x}}} \mathord{\left/ {\vphantom {{ - q{{{v}}_{y}}{{B}_{x}}} {mc,}}} \right. \kern-0em} {mc,}}$
в которых компоненты поля во всем пространстве
${{B}_{z}} = B{\text{sin}}\theta ,\,\,\,\,{{B}_{x}} = - B{\text{cos}}\theta ,\,\,\,\,{{E}_{y}} = {{uB{\text{sin}}\theta } \mathord{\left/ {\vphantom {{uB{\text{sin}}\theta } {c,}}} \right. \kern-0em} {c,}}$
где с – скорость света.

Для строго перпендикулярной МЗУВ (θ = π/2) уравнения (1)(3), описывающие динамику ионов в рампе в плоскости х0y, имеют аналитическое решение:

$х = ({{{v}}_{x}}_{0} - {{{v}}_{d}}){{{\text{sin}}({{\Omega }_{1}}{{\omega }_{Н}}t)} \mathord{\left/ {\vphantom {{{\text{sin}}({{\Omega }_{1}}{{\omega }_{Н}}t)} {{{\Omega }_{1}}}}} \right. \kern-0em} {{{\Omega }_{1}}}} + {{{v}}_{d}}{{\omega }_{Н}}t,$
${{{v}}_{x}} = ({{{v}}_{{x0}}} - {{{v}}_{d}}){\kern 1pt} {\text{cos(}}{{\Omega }_{1}}{{\omega }_{Н}}t{\text{)}} + {{{v}}_{d}},\,\,\,\,{{{v}}_{y}} = {{\omega }_{Н}}(ut - x),$
где vx0 – начальная скорость ионов, попадающих в рамп, ωН ионная циклотронная частота, Ω = = сЕ0(2eφA/m)–1/2/В, Ω1= (1 + Ω2)1/2, vd = ${u \mathord{\left/ {\vphantom {u {\Omega _{1}^{2}}}} \right. \kern-0em} {\Omega _{1}^{2}}}.$ Как видно из этих формул, частица в рампе в данном случае совершает движение с постоянной скоростью vd, на которое накладываются колебания с частотой ωН1.

Для косых МЗУВ (0 < θ < π/2) в частном, но чрезвычайно важном для рассматриваемой проблемы случае, а именно, при описании динамики захваченных в рампе частиц, аналитическое решение тоже удалось найти. Это оказалось возможным вследствие того, что для захваченных в рампе ионов величина скорости vx $ \ll $ u. Учитывая этот факт, в нулевом приближении в уравнении (2) можно положить vx = 0 и тогда система уравнений (2)–(3) имеет точное решение:

${{{v}}_{y}} = u{\text{tg}}\theta {\text{sin}}{{\tau }_{1}},\,\,\,\,{{{v}}_{z}} = u{\text{tg}}\theta (1 - {\text{cos}}{{\tau }_{1}}),$
где τ1 = ωНt cos θ. Эти формулы описывают колебательные движения ионов в плоскости y0z с частотой ωН cos θ. В следующем приближении, учитывая принятые нами неравенства vx $ \ll $ u, Ω $ \gg $ 1, Ω $ \gg $ cos θ, получим соотношения для координаты х и компоненты скорости vx:

$\begin{gathered} х = ({{{v}}_{{x0}}} - {{u{\text{si}}{{{\text{n}}}^{2}}\theta } \mathord{\left/ {\vphantom {{u{\text{si}}{{{\text{n}}}^{2}}\theta } {{{\Omega }^{2}}}}} \right. \kern-0em} {{{\Omega }^{2}}}}){{{\text{sin}}(\Omega {{\omega }_{Н}}t)} \mathord{\left/ {\vphantom {{{\text{sin}}(\Omega {{\omega }_{Н}}t)} {(\Omega {{\omega }_{Н}})}}} \right. \kern-0em} {(\Omega {{\omega }_{Н}})}} + \\ + \,\,{{u{\text{tg}}\theta {\text{sin}}\theta {\text{sin}}{{\tau }_{1}}} \mathord{\left/ {\vphantom {{u{\text{tg}}\theta {\text{sin}}\theta {\text{sin}}{{\tau }_{1}}} {({{\omega }_{Н}}{{\Omega }^{2}}),}}} \right. \kern-0em} {({{\omega }_{Н}}{{\Omega }^{2}}),}} \\ \end{gathered} $
${{{v}}_{x}} = ({{{v}}_{{x0}}} - {{u{\text{si}}{{{\text{n}}}^{2}}\theta } \mathord{\left/ {\vphantom {{u{\text{si}}{{{\text{n}}}^{2}}\theta } {{{\Omega }^{2}}}}} \right. \kern-0em} {{{\Omega }^{2}}}}){\text{cos}}(\Omega {{\omega }_{Н}}t) + {{u{\text{si}}{{{\text{n}}}^{2}}\theta {\text{cos}}{{\tau }_{1}}} \mathord{\left/ {\vphantom {{u{\text{si}}{{{\text{n}}}^{2}}\theta {\text{cos}}{{\tau }_{1}}} {{{\Omega }^{2}}.}}} \right. \kern-0em} {{{\Omega }^{2}}.}}$

Из этих решений следует, что на колебания иона с частотой ωН cos θ накладываются колебания с более высокой частотой (Ω $ \gg $ cos θ).

Еще один случай – это динамика ионов в косой МЗУВ при Ω $ \ll $ 1, в котором удалось найти приближенное аналитическое решение уравнений (1)–(3), которое имеет вид

$х \approx {{G{\text{sin(}}{{\Omega }_{2}}{{\omega }_{Н}}t{\text{)}}} \mathord{\left/ {\vphantom {{G{\text{sin(}}{{\Omega }_{2}}{{\omega }_{Н}}t{\text{)}}} {({{\omega }_{Н}}{{\Omega }_{2}})}}} \right. \kern-0em} {({{\omega }_{Н}}{{\Omega }_{2}})}},$
где Ω2= Ω cos θ, G = (uvx0)sin2θ + vx0.

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

1. В системе отсчета, движущейся вместе с фронтом ударной волны, все ионы, набегающие на рамп косых МЗУВ с безразмерной скоростью s0 (нормированной на скорость u), условно можно разделить на следующие категории: (a) пролетные – это те ионы, которые сходу пересекают рамп и больше никогда в него не возвращаются; (b) отраженные это ионы, которые после простых или сложных движений в рампе, выходят из него и удаляются от рампа вверх по течению; (c) вращающиеся это ионы, формирующие перед рампом подножье; (d) захваченные частицы это ионы, которые в пределах рампа медленно перемещаются относительно него, поэтому длительное, но конечное время находятся в рампе. Как вращающиеся, так и захваченные ионы присутствуют во фронте МЗУВ временно. В итоге, и те и другие оказываются в строго поперечных МЗУВ только пролетными, а в косых МЗУВ либо пролетными, либо отраженными.

2. Найден критический угол θ*, который разделяет плоскость s0θBn, на которой возможно существование косых МЗУВ, на две области. Пример такой кривой взят из работы [Кичигин, 2018] и приведен на рис. 3. В одной области (с углами 90° ≥ θBn > θ*) все набегающие на рамп МЗУВ ионы становятся пролетными, а во второй (с углами 0 < θBn < θ*) все набегающие ионы становятся отраженными. Величина критического угла θ*, разделяющего набегающие на рамп ионы на пролетные и отраженные, зависит как от скорости налетающих на рамп частиц s0, так и от пространственного размера рампа ударного фронта БУВ, который контролируется параметром D = = МАψ(cpi)/d, где МА – альфвеновское число Маха, ψ – безразмерный перепад потенциала в рампе, нормированный на энергию набегающих на рамп ионов mu2/2, ωpi = (4πn0е2/т)1/2 – плазменная ионная частота (n0 – плотность невозмущенной ударной волной плазмы).

Рис. 3.

Зависимость предельного угла θ* от скорости падающих на рамп ионов s0 (в системе отсчета, движущейся вместе с фронтом БУВ). Приведенная на рис. кривая разделяет все налетающие на рамп ионы на пролетные (для значений θ и s0, находящихся выше кривой) и отраженные (для значений θ и s0, находящихся ниже кривой). Параметр D = 1.

3. При больших значениях параметра D (D $ \gg $ 1, что соответствует малым размерам рампа: d $ \ll $ cpi) для основной части налетающих на рамп ионов, имеющих скорость, близкую к скорости волны в системе отсчета БУВ (s0 ≈ 1), критический угол имеет значение θ* 52°. Когда параметр D становится меньше единицы (D < 1, что соответствует размеру рампа d > cpi) для значений критического угла θ* > 60° отраженные ионы полностью отсутствуют (см., например, рис. 3). При дальнейшем уменьшении параметра D для ионов, имеющих малую скорость (в системе отсчета БУВ s0 ≈ 0), максимальное значение критического угла не превышает 45°, а для основной части ионов, набегающих на рамп (в системе отсчета БУВ s0 ≈ 1), критический угол θ* уменьшается и стремится к нулю при D → 0. (Такая тенденция просматривается и на рис. 3).

4. В косых МЗУВ небольшая часть ионов, вследствие малой начальной скорости (в системе отсчета БУВ s0 ≈ 0), в процессе медленного разворота в рампе оказываются временно захваченными. Эти ионы за счет серфотронного ускорения в рампе приобретают энергию вплоть до 2D2mu2 (при θ* ≈ π/2), с такой энергией ионы выходят из рампа навстречу набегающему потоку плазмы под малым углом ≈(π/2 θ*)/2 к плоскости фронта БУВ, образуя популяцию отраженных частиц. В околоземной ударной волне предельная энергия отраженных ионов (2D 2mu2) при θ* ≈ π/2 может достигать значений порядка 1 МэВ. Предельная энергия отраженных ионов уменьшается в зависимости от угла θ* как 1/cos2θ* (при условии π/4 < θ* < π/2). В работах (Кичигин, 2009; Кичигин 2018) полагается, что эти отраженные ионы формируя границу ионного форшока, движутся остронаправленно вдоль этой границы и именно они представляют популяцию так называемых продольных пучков (field-aligned beam – FAB).

5. Вследствие серфотронного ускорения в рампе энергия захваченных ионов имеет резкий пик в точке, в которой угол θBn равен критическому (θBn = θ*). В ближайшей окрестности критического угла θ*, с одной стороны, где θBn > θ*, все захваченные в ударном фронте ионы окажутся в итоге пролетными, имеющими максимальную энергию ≈D 2mu2, с другой стороны, где θBn < θ* все захваченные ионы окажутся в итоге отраженными, максимальная энергия которых ≈2D 2mu2. Таким образом, если рассматривать случай околоземной БУВ, то около точки θBn = θ* на профиле ударного фронта как пролетные, так и отраженные ионы будут иметь максимальные энергии. Именно из этой точки (точка “b” на рис. 1) энергичные отраженные ионы начинают свое остронаправленное движение, образуя наблюдаемую на КА эжекцию продольных пучков, которые образуют границу ионного форшока.

4. СРАВНЕНИЕ СВОЙСТВ ИОННОГО ФОРШОКА, СЛЕДУЮЩИХ ИЗ ТЕОРЕТИЧЕСКОЙ МОДЕЛИ И ИЗ НАБЛЮДАТЕЛЬНЫХ ДАННЫХ, ПОЛУЧЕННЫХ С ПОМОЩЬЮ КОСМИЧЕСКИХ АППАРАТОВ

В данной работе мы хотим обратить внимание на то, что в рамках идеологии, приведенной в работах [Кичигин, 2009, 2018], достаточно убедительно объясняются как физические основы происхождения продольных пучков, так и наблюдательные данные, полученные при исследовании ионного форшока околоземной БУВ. Начнем с наблюдений, полученных в работе [Croooker et al., 1981], в которой с помощью приборов, установленных на спутниках ISEE 1 и ISEE 3 , в окрестности околоземной ударной волны в августе–сентябре 1978 г. измерялся энергетический спектр ионов и направление ММП. В результате анализа спутниковых данных авторы работы [Croooker et al., 1981] пришли к следующим выводам: (1) – источник энергичных ионов находится не в магнитослое и не на магнитопаузе, а вероятнее всего – во фронте БУВ; (2)– энергичные ионы в магнитослое, а также популяции промежуточных и диффузных ионов перед фронтом БУВ присутствуют (отсутствуют) только тогда, когда угол θBn < 60° (θBn > 60°). На рисунке 4 (это рис. 7 из работы [Croooker et al., 1981]) приведены картинки, на которых изображены области магнитослоя, где, согласно мнению авторов, присутствуют энергичные ионы. Наличие энергичных ионов в указанных на рис. 4 областях магнитослоя вполне обьяснимо в рамках нашей модели: энергичные ионы из заштрихованных областей, указанных на рис. 1, вполне закономерно могут оказаться в заштрихованных областях, отмеченных на рис. 4. Далее, если предположить, что ширина фронта околоземной БУВ, которую исследовали авторы работы [Croooker et al., 1981], была порядка cpi (параметр D ~ 1), что вполне приемлемо, то результаты работы [Croooker et al., 1981] можно уверенно интерпретировать в рамках нашей модели, согласно которой для параметра D ~ 1 критический угол θ*, больше которого перед фронтом БУВ отсутствуют продольные пучки (отраженные ионы), примерно равен 60° (см. рис. 3).

Рис. 4.

Вид магнитослоя в плоскости эклиптики при различных углах θBn. Заштрихованные области – это области, заполненные энергичными ионами. Показанный на рисунке угол, расположенный на границе заштрихованной области (локальный в данной точке ударной волны угол θBn) равен 60°. (Данный рисунок – это рис. 7 из работы [Croooker et al., 1981]).

Критический угол 60° фигурирует также в работах [Oka et al., 2005; Kis et al., 2007]. Авторы работы [Oka et al., 2005] в октябре 1995 г. с помощью аппаратуры, расположенной на спутнике Geotail, исследовали свойства продольных пучков в области ионного форшока, расположенного на утренней стороне околоземной ударной волны. Они обнаружили, что если перемещаться от области квазипараллельного участка околоземной БУВ (где θBn < 45°) вдоль ударного фронта в направлении к подсолнечной точке (угол θBn при этом растет), то за точкой ударного фронта околоземной БУВ, где θBn = 60°, величина потока продольных пучков резко падает. Цитата из абстракта работы [Oka et al., 2005]: “We have then compared the FAB flux normalized by the solar wind flux, F, with the shock angle θBn and have found that F falls off rapidly just above θBn > 60°, but it maintains significant level (F = 0.01%) up to 75°”. Это адекватно объяснятся нашей моделью, в рамках которой такой поток при θBn > 60° должен отсутствовать. (Присутствие мизерного количества продольных пучков при θBn > 60°, наличие которых отмечают авторы работы [Oka et al., 2005], можно объяснить рассеянием пучков на колебаниях и волнах).

В работе [Kis et al., 2007] с использованием данных аппаратуры спутника Cluster (Cluster spacecraft) исследовалось рассеяние продольных пучков перед ударным фронтом околоземной БУВ. Во-первых, авторы работы [Kis et al., 2007] обнаружили, что продольные пучки движутся остронаправленно, а поперечный разброс по скоростям может расти за счет рассеяния их на колебаниях плазмы. Далее, они пишут: “Taking into consideration the direction of the interplanetary magnetic field, the solar wind and the beam velocities and the position and shape of the bow shock, we have calculated the trajectory and the point of origin of the FAB on the bow shock surface. The results show that the FAB originates on the quasi-perpendicular side of the bow shock at a position where the θBn = 60° This is also in good agreement with earlier observations…”. Изложенный в приведенной цитате вывод полностью вписывается в нашу концепцию.

5. ЭФФЕКТ ТРОИЦКОЙ-БОЛЬШАКОВОЙ В СВЕТЕ ПРЕДЛОЖЕННОЙ МОДЕЛИ

Еще один пример, в котором можно использовать нашу модель, касается эффекта Троицкой-Большаковой [Большакова и Троицкая, 1968; Большакова и др., 1987; Potapov et al., 2018]. Троицкая и Большакова обнаружили влияние ориентации ММП на амплитудный режим дневных устойчивых геомагнитных пульсаций Рс3. Они установили, что амплитуда регистрируемых на Земле пульсаций Рс3 падает до нуля, когда отклонение вектора ММП от линии Солнце–Земля равно 90°. В настоящее время по поводу происхождения пульсаций Рс3 установлено: (1) – имеется почти стопроцентная корреляция между регистрацией на Земле пульсаций Рс3, и наличием волн в пространственной области перед ударным фронтом околоземной БУВ, примыкающей к подсолнечной точке; (2) – волны в этой области возникают за счет неустойчивости потока быстрых ионов, движущихся перед фронтом околоземной БУВ [Guglielmi, 1974; Потапов, 1974]. Другими словами, в отсутствие быстрых пучков ионов перед фронтом околоземной БУВ волны в СВ не генерируются, а, следовательно, пульсации Рс3 на Земле не регистрируются.

Эффект Троицкой-Большаковой, т.е. исчезновение пульсаций Рс3 при θBn = 90°, вполне объясним в рамках нашей модели. Вначале отметим, что для строго перпендикулярной МЗУВ (θBn = 90°) перед ее фронтом временно находятся только вращающиеся ионы, создающие подножье, а отраженных ионов (т.е. убегающих от рампа вниз по потоку) здесь не может быть в принципе. Если волна косая (θBn < 90°), то, как это следует из нашей модели, быстрые ионы перед фронтом БУВ будут отсутствовать для углов θBn, больше критического, т.е. в интервале θ* < θBn < 90°. Для околоземной БУВ при θBn = 90° в подсолнечной точке и для значений параметра D ~ 1 область отсутствия продольных пучков ионов в окрестности подсолнечной точки будет максимальной при величине критического угла θ* = 60°. Таким образом, из нашей модели следует, что если в подсолнечной точке θBn = 90°, то перед околоземной ударной волной, в конусе углов 60° < θBn < 90° быстрые ионы и, следовательно, колебания будут отсутствовать, поэтому, согласно эффекту Троицкой–Большаковой, пульсации Рс3 на Земле регистрироваться не будут.

6. ЗАКЛЮЧЕНИЕ

Как величина минимального критического угла, примерно равного 60°, так и другие свойства продольных пучков, следующие из нашей модели, уверенно согласуется с результатами наблюдений, приведенных в работах [Crooker et al., 1981; Oka et al., 2005; Kis et al., 2007]. Также, в рамках нашей модели находит достаточно простое объяснение эффект Троицкой–Большаковой [Большакова и Троицкая, 1968].

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

  1. Большакова О.В., Троицкая В.А. Связь между направлением межпланетного магнитного поля с режимом устойчивых колебаний // ДАН. Т. 180(№ 2). С. 343–346. 1968.

  2. Большакова О.В., Троицкая В.А., Русакова Т.Б. Контроль межпланентным магнитным полем пульсаций Рс3 в солнечном ветре и на Земле // Геомагнетизм и аэрономия. Т. 27. № 2. С. 345–346. 1987.

  3. Кичигин Г.Н. О происхождении энергичных частиц в области форшока околоземной ударной волны // Письма в астроном. журн. Т. 35. № 4. С. 295–303. 2009.

  4. Кичигин Г.Н. Динамика ионов во фронте магнитозвуковых ударных волн // Солнечно-земная физика. Т. 4. № 4. С. 23–31. 2018. https://doi.org/10.12737/szf-44201803

  5. Потапов А.С. Возбуждение геомагнитных пульсаций Рс3 перед фронтом околоземной ударной волны пучком отраженных ионов // Исследования по геомагнетизму, аэрономии и физике Солнца. Т. 34. С. 3–12. 1974.

  6. Balogh A., Treumann R.A. Physics of Collisionless Shocks. Springer Science Business Media. N.Y. 2013. https://doi.org/10.1007/978-1-4614-6099-2

  7. Bonifazi C., Moreno G. Reflected and Diffuse Ions Backstreaming From the Earth’s Bow Shock 1. Basic Properties // J. Geophys. Res. V. 86. P. 4397–4404. 1981a.

  8. Bonifazi C., Moreno G. Reflected and Diffuse Ions Backstreaming From the Earth’s Bow Shock 2. Origin // J. Geophys.Res. V. 86. P. 4405–4412. 1981b.

  9. – Burgess D., Möbius E, Scholer M. Ion acceleration at the Earth’s bow shock // Space Sci. Rev. V. 173. P. 5–47. 2012. https://doi.org/10.1007/s11214-012-9901-5

  10. Crooker N.U., Eastman T.E., Frank L.A., Smith E.J., Russel C.T. Energetic Magnetosheath Ions and the Interplanetary Magnetic Field Orientation // J. Geophys. Res. V. 86. № A6. P. 4455–4460. 1981.

  11. Eastwood J.P., Lucek E.A., Mazelle C., Meziane K., Narita Y., Pickett J., Treumann R.A. Foreshock // Space Sci. Rev. V. 118. P. 41–94. 2005. https://doi.org/10.1007/s11214-005-3824-3

  12. Guglielmi A.V. Diagnostics of the magnetosphere and interplanetary medium by means of pulsations // Space Sci. Rev. V. 16. P. 331–345. 1974.

  13. Kis A., Scholer M., Klecker B., Kucharek H., Lucek E.A., R`eme H. Scattering of field-aligned beam ions upstream of Earth’s bow shock // Ann. Geophys. V. 25. P. 785–799. 2007. https://www.ann-geophys.net/25/785/2007

  14. Oka M., Terasawa T., Saito Y., Mukai T. Field-aligned beam observations at the quasi-perpendicular bow shock: Generation and shock angle dependence // J. Geophys. Res. V. 110. A05101-116. 2005. https://doi.org/10.1029/2004JA010688

  15. Paschmann G., Sckopke I., Papamastorakis J., Asbridge J.R., Bame S.J., Gosling J.T. Characteristics of reflected and diffuse ions upstream from the Earth’s bow shock // J. Geophys. Res. V. 86. P. 4355–4364. 1981.https://doi.org/10.1029/JA086iA06p0435

  16. Potapov A.S., Polyushkina T.N., Guglielmi A.V. Troitskaya–Bolshakova effect as a manifestation of the solar wind wave turbulence // Planet. Space Sci. V. 151. P. 78–84. 2018.

Дополнительные материалы отсутствуют.