Солнечно-земная Физика

Наблюдения с самолетов и спутников с использованием оптической аппаратуры показало, что кроме аврорального овала, где выделяется основная энергия, связанная с вторжением заряженных частиц, существует еще широкая полоса свечения, которая окаймляет авроральный овал с экватора. Это свечение, по сравнению с овалом, отличается относительной однородностью, но границы его достаточно четкие [Anger and Lui, 1973; Lui et al., 1973]. Измерения проводились с помощью сканирующего фотометра в эмиссиях l =391.4 и 557.7 нм. Экваториальная граница диффузного свечения меняется с широтой, поднимаясь к полудню в более высокие широты. По данным самолетных наблюдений [Buchau et al., 1972] ширина области диффузного свечения может доходить до 4° широты, что подтверждается также наблюдениями на спутнике [Slater et al., 1980].

Высота экваториальной границы диффузного свечения, по данным параллактическим наблюдениям, составляет около 200 км [Алексеев и др., 1975], что больше, чем для дискретных форм сияний, но значительно меньше, чем для среднеширотных красных дуг.

Хотя детальное исследование морфологии диффузного свечения связано с наблюдениями на спутниках, но наличие такого свечения было обнаружено ранее с помощью наземной спектральной аппаратуры [Евлашин, 1961; Rees et al., 1961], где бвло показано, что в вечерние часы наблюдается красное диффузное свечение, которое располагается экваториальнее дискретных форм полярных сияний. Обычно это свечение сопровождается водородным излучением, которое затухает при появлении ярких дискретных форм. Высыпание протонов зарегистрировано и по спутниковым данным [Lui et al., 1977], но их энергетический вклад был существенно ниже, чем для электронов.

Обычно экваториальная граница диффузного свечения относительно гладкая. Однако во время сильных геомагнитных возмущений на ней могут появляться крупномасштабные волнообразные структуры. Пример трех последовательных регистраций таких структур в вечерние часы в южном полушарии спутником DMSP приведен на рис.10 [Lui et al., 1982].

Рис.10. Полоса свечения в вечернем секторе 26 июля 1978 года по данным спутника DMSP [Lui et al., 1982]. Точками показана область диффузного свечения, зачерненными полосами - дискретные формы сияний.

Пространственное совпадение областей диффузного свечения и высыпания подтверждается одновременными измернениями этих границ. В [Slater et al., 1980] сопоставлялось положение границы диффузного свечения в эмиссии l 630 нм по данным наземных наблюдений сканирующим фотометром с наблюдениями электронных потоков на спутниках при прохождении ими над районом наземнных измерений. За границу диффузного свечения принималась область, где интенсивность свечения l 630 нм на 100 R превышала фон ночного неба. Показано, что экваториальная граница свечения, которая была обычно довольно резкой, хорошо совпадает с границей повышения потока энергии вторгающихся электронов до 10^-2эрг/(cм²*с*стер). Зависимость между этими границами приведена на рис.12. Для выбранных пороговых значений положения границ практически совпадают. Хорошее соответствие между границами бесструктурных электронных вторжений и диффузного свечения обнаружено также в [Lui et al., 1977].

Рис.12. Зависимость между положениями экваториальных границ диффузного свечения Fў(I) и диффузного высыпания Fў(F) [Slater et al., 1980].

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

Рис.13. Полярные сияния и высыпания частиц в вечернем секторе 11 октября 1974 г. по данным спутника DMSP [Meng, 1976]. а - взаимное положение дискретных форм (зачерненные участки) и диффузного свечения (точки). Сплошной линией показана траектория спутника со шкалой географических широт, штриховкой по краям - область засветки Солнцем. б - вариации потока энергии (эрг / см Ч с Ч ср.) вдоль траектории пролета спутника.

Так как границы диффузного свечения и диффузного высыпания практически совпадают [Slater et al., 1980], то появляется возможность сопоставлять эти данные. В [Гальперин и др., 1977; Вальчук и др., 1979] по наблюдениям на спутниках серии “Ореол” границы диффузного высыпания определялись по резкому повышению интенсивности потока энергии. Положение экваториальной границы диффузного свечения (или высыпания) смещается в более низкие широты при повышении уровня магнитной активности. Аналогичные результаты по спутниковым наблюдениям были получены также в [Slater et al., 1980; Gussenhoven et al., 1981]. На рис.14, взятому из [Николаенко и др., 1983], приведено изменение положения границ диффузного высыпания в вечерние и полуночные часы в зависимости от магнитной активности по данным [Гальперин и др., 1977] и по наблюдениям полученным на 4 года позднее [Gussenhoven et al., 1981]. Из рисунка видно, что разные авторы на разных спутниках получили практически совпадающие результаты.

Рис.14. Зависимость положения экваториальной границы диффузных высыпаний от магнитной активности в вечерние (а) и полуночные (б) часы. Крестики и сплошные линии - данные [Гальперин и др., 1977], кружки и пунктир - данные [Gussenhoven et al., 1981].

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

Первая простая аппроксимация экваториальной границы овала в зависимости от местного времени и уровне магнитной активности, который, как и в [Feldstein and Starkov, 1967], задавался Q-индексом, дана в [Старков, 1969]. Изменения положения экваториальной границы для Q > 0 c точностью ~ 1° широты описывались формулой

J_e = 18 + 0.9 × Q + 5.1× cos (t - 12° ), (1)

Более подробное математическое. описание границ овала проведено в [Holzworth and Meng, 1975]. Положение экваториальной и приполюсной границ описывалось тремя гармониками, амплитуды и фазы которых зависели от магнитной активности. Было показано, что соответствующим смещением центра координат можно существенно уменьшить амплитуду первой гармоники до величины, меньшей 0.1° широты. Так как амплитуды второй и третьей гармоник не превышали одного градуса широты, то приблизительно с такой точностью границы овала можно аппроксимировать кругом при соответствующем сдвиге начала координат. Однако в [Holzwort and Meng, 1975] не указаны координаты этих новых центров. Кроме того, в качестве магнитной активности тоже использовался дискретный Q-индекс на ночной стороне в зоне сияний, публикация которого в настоящее время отсутствует.

В связи с этим в [ Старков, 1994б ] приведены аппроксимационные формулы, подобные используемым в [Holzworth and Meng, 1975], но у которых коэффициенты зависят от AL. Кроме того, в аналогичной форме описана экваториальная граница диффузного свечения. Положение границ задавалось в виде

J_e = A₀+ A₁cos [15(t + a₁ )] + A₂cos [15(2t + a₂ )] + A₃cos [15(3t + a₃ )],      (2)

где J - коширота границы в исправленных геомагнитных координатах; A_i- амплитуды, выраженные в градусах широты; t - время в часах местного геомагнитного времени; a _i- фазы в том же часовом измерении.

На первом этапе по данным о границах овала, приведенным в [Feldstein and Starkov, 1967; Старков и Фельдштейн, 1968], определялись коэффициенты A_iи a_i для разных значений Q-индексов. Значения Q пересчитывались в AL по формулам регрессии [Cтарков, 1994а] и строились зависимости A_i= j (lg [AL]) и a_i = f (lg [AL]). Эти кривые, в свою очередь, аппроксимировались полиномами третьей степени вида

A_i( a _i) = a₀+ a₁g [AL] + alg[AL] + a₂lg[AL]. (3)

Для диффузного свечения исходные данные взяты из работы [Gussenhoven et al., 1983], где для всех часов местного времени приведены зависимости
F¢ = b₀ +b₁.K_p
Так как b и bзначительно изменялись от часа к часу и точность их определения была различна, то они предварительно сглаживались по формуле трапеции с весовыми множителями, равными коэффициентам корреляции для данного часа, приведенными в [Gussenhoven et al., 1983]. Значения Kр в свою очередь переводились в величины AL-индекса по соответствующим формулам регрессии [Cтарков, 1994а]. Все вычисления А_i и a_i проводились методом наименьших квадратов.

Характер зависимости a_i oт lg(AL) был различен, но в большинстве случаев полином третьей степени достаточно хорошо описывал их поведение. Наибольший разброс наблюдается для фаз, поэтому перед аппроксимацией они предварительно сглаживались по формуле трапеции.

Рис.15. Примеры зависимости коэффициентов A и a от уровня магнитной активности: сплошные кривые - исходные данные ; штриховка - результаты аппроксимации полиномом. а - амплитуда нулевой гармоники для экваториальной границы овала, б - амплитуда второй гармоники для экваториальной границы диффузного свечения, в - фаза третьей гармоники для приполюсной границы овала, г - фаза второй гармоники для экваториальной границы овала [Старков, 1994б].

На рис.15 приведены четыре примера таких аппроксимаций, последний случай есть пример наиболее плохой аппроксимации из-за резких и нерегулярных изменений фазы.
Коэффициенты aдля всех Aи a даны в таблице 1. Для определения границ необходимо определить коэффициенты A и a для данного уровня магнитной активности по формуле (3) и, подставив эти значения в формулу (2), вычислить положение соответствующей границы в выбранном интервале времени для данного значения AL.

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

	А0	А1	А	А2	А2	А3	А3
а0 а1 а2 а3	-0.07 24.54 -12.53 2.15	-10.06 19.83 -9.33 1.24	-6.61 10.17 -5.80 1.19	-4.44 7.47 -3.01 0.25	6.37 -1.10 0.34 -0.38	-3.77 7.90 -4.73 0.91	-4.48 10.16 -5.87 0.9

Экваториальная граница овала

½ a 1.61 -9.58 -2.22 -12.07 -23.98 -6.56 -20.07 ½

½ a 23.21 17.78 1.50 17.49 42.79 11.44 36.67 ½

½ a -10.97 -7.20 -0.58 -7.96 -26.96 -6.73 -24.20 ½

½ a 2.03 0.96 0.08 1.15 5.56 1.31 5.11 ½

Экваториальная граница диффузного свечения

½ a 3.44 -2.41 -1.68 -0.74 8.69 -2.12 8.61 ½

½ a 29.77 7.89 -2.48 3.94 -20.73 3.24 -5.34 ½

½ a -16.38 -4.32 1.58 -3.09 13.03 -1.67 -1.36 ½

½ a 3.35 0.87 -0.28 0.72 -2.14 0.37 0.76 ½

Среднеквадратическая ошибка аппроксимации исходных данных, взятых из [Feldstein and Starkov, 1967], составляет ~ 0.3° широты, и только при больших AL она повышается до 0.4° . Ошибка после аппроксимации коэффициентов (2) полиномами (3) увеличивается на 0.1 - 0.2° . Ошибка определялась сравнением полученных значений границ с данными для границ, приведенными в [Feldstein and Starkov, 1967] для шести моментов местного геомагнитного времени. Ошибка оценки границы диффузного свечения в среднем на 0.1° выше.

Используя формулу (2), можно рассчитать размеры овала и диффузного свечения. Площадь поверхности, ограниченной кривой в полярных координатах при интегрировании от 0 до 2p , равна 1/2 интеграла от r²dj . В нашем случае имеем:

^{J /2 dj = Adj +
A cos(kj +a_k ) dj +}

^{A A cos(kj + a_k ) cos(pj +a_p) dj ,} (4)

где k ¹ p, а J - длина дуги, выраженная в градусах широты. В результате интегрирования получаем:

S = p A (5)

Основной вклад в сумму вносит A, которая в 5 - 20 раз превышает амплитуду A. Последняя, в свою очередь обычно в несколько раз больше амплитуд следующих гармоник.

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

Рис.16. Вариации поверхности полярной шапки (кружки), овала (треугольники) и экваториального диффузного свечения (крестики) в зависимости от уровея магнитной активности [Старков, 1994б].

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

Обычно А₁> А₂ и А₃. Тогда границы свечения можно представить в виде круга с радиусом A, центр которого смещен относительно геомагнитного полюса и на который накладываются сравнительно небольшие гармонические колебания более высокого порядка, связанные с Aи A. Первое представление границ овалов в виде кругов, центры которых смещены по линии “полдень-полночь”, по данным о положении дневных и ночных сияний дано в [Старков и Фельдштейн, 1967а]. Детально этот вопрос исследован в [Holzworth and Meng, 1975], где показано, что можно найти такое положение нового начала координат, при котором амплитуда первой гармоники существенно уменьшается. Практически операция смещения центра эквивалентна приравниванию нулю амплитуды первой гармоники. Действительно, как показали расчеты, расстояние нового центра от полюса равно амплитуде первой гармоники А₁, а азимут его соответствует фазе a₁ .

Выполненные расчеты коэффициентов разложения для смещенного центра показали резкое уменьшение амплитуды первой гармоники. Так для экваториальной границы овала она изменяется от 0.01 до 0.09° со средним значением 0.03° , для полярной границы овала - от 0.02 до 0.12° со средним значением 0.06° . Для границы диффузного свечения значение A1' несколько больше и составляет ~ 0.07° для малых и средних уровней магнитной возмущенности и растет с увеличением AL-индекса, достигая 0.25° при AL » -600 нТл. Это может быть связано с ошибками исходных данных. В целом же амплитуды первой гармоники после смещения центра всегда существенно меньше амплитуд Aи A

Наличие A1'¹ 0 после соответствующего смещения центра свидетельствует о том, что для более точного описания границ аврорального свечения, по-видимому, необходимо учитывать гармоники более высоких порядков.

Рис.17. Аппроксимация кругами границ аврорального свечения. Положение центров (а, б) и радиусы кругов (в) в зависимости от уровня магнитной активности: полярная шапка (кружки), экваториальная граница овала (треугольники), экваториальная граница диффузного свечения (крестики) [Старков, 1994б].

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

r= .

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

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

J ¢ =A+Acos[15(2t+a)+A cos[15(3t+a )

где А- радиус круга и d =J -A.
Тогда s² = d . Заменим суммирование интегрированием, то есть перейдем к бесконечно малым шагам суммирования. В этом случае n равно интервалу интегрирования, то есть 2p . Интегралы от произведения косинусов для пределов 0 - 2p , как и в случае вычисления поверхностей, равны нулю. Остаются только интегралы от A cos (kt + a_k ), поэтому окончательно получаем

     (6)

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

Рис.18. Суточные изменения границ аврорального свечения при разном уровне магнитной активности. а - для AL=15, б - для AL = 160, в -для AL = 460 нТл. Сплошные кривые - результаты расчетов по формулам (2) и (3); штрриховые - аппроксимации кругами со смещенным центром. Кружки, треугольники и крестики - экспериментальные данные о положении приполюсной, экваториальной границ овала [Feldstein and Starkov, 1967] и экваториальной границы диффузного свечения [Gussenhoven et al., 1983] соответственно.

На рис.18 приведены временные вариации положения границ аврорального свечения, полученные по формулам (2) и (3), а также также аппроксимации границ кругами со смещенным центром для трех уровней магнитной активности. Точками нанесены экспериментальные данные [Feldstein and Starkov, 1967; Gussenhoven et al., 1983]. Результаты расчетов по формулам (2) и (3) хорошо согласуются с эксперименом. Для аппроксимации кругами заметное расхождение наблюдается только при высокой магнитной активности.

Примеры положения границ аврорального овала и диффузного свечения для разных значений AL-индекса, рассчитанные по формулам (2) и (3), приведены на рис.19. .

Рис. 19. Положения аврорального овала и диффузного свечения при разном уровне магнитной активности.

При спокойных условиях (рис.19а) овал сияний представляет собой тонкое кольцо, ширина которого 2 - 3° широты, центр которого смещен к полуночным часам. Эта ситуация, согласно [Старков, 1994а], соответствует Q-индексу, равному нулю. Диффузное свечение занимает значительно большую поверхность (см. также рис.16), но в связи с тем, что, как будет показано ниже, интенсивность свечения в овале приблизительно на порядок выше, чем в диффузном свечении, то основное энерговыделение идет в овале, хотя при спокойных условиях оно сравнимо с энергетикой диффузного свечения. Граница свечения по форме близка к кругу, центр которого тоже смещен к экватору на полуночном меридиане. Однако за счет того, что это смещение меньше, чем для овала, то наибольшая ширина диффузного свечения наблюдается в районе полудня.

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

При дальнейшем увеличении магнитной активности (рис.19 д, е) овал принимает форму, похожую на среднестатистический овал, полученный Фельдштейном [Фельдштейн, 1963]. Максимальная ширина диффузного свечения наблюдается в предполуденные часы. При дальнейшем увеличении магнитной активности (рис.19 ж) форма овала начинает искажаться. В дневной части наиболее узкий участок наблюдается в районе 08 MLT, появляется некоторое расширение в 14 часов. Возможно, это связано с тем, что максимальная магнитная активность, для которой, согласно [Feldstein and Starkov, 1967; Старков и Фельдштейн, 1968], определялись положения границ овала, была равна Q = 7, что соответствует AL = 570 нТл, то есть величины коэффициентов в формулах (2) и (3) интерполировались для значений AL, для которых уже не было контролирующих экспериментальных данных.

Надо отметить, что для многих расчетов для определения экваториальной границы овала можно пользоваться простой формулой (1). В этом случае граница овала автоматически аппроксимируется кругом со смещенным центром. Фактически первые два члена в формуле (1), а именно (18 + 0.9× Q), должны быть равны коэффициентам А формулы (2) для соответствующего AL-индекса. На рис.20 приведена зависимость (18 + 0.9× Q) от А, штриховая прямая на рисунке соответствует линии, проходящей под углом в 45° . Видно, что только для Q = 0 имеется расхождение, а в остальных случаях эти величины практически равны. Сопоставление с данными спутников показало, что рассчитанные границы аврорального овала хорошо совпадают с прямыми наблюдениями, см. например[Meng et al., 1977]. Для проверки точности аппроксимации границ диффузного свечения очень удобным оказался финский радар системы CUTLASS. Эта система состоит из двух радаров, наблюдающих отражение из общей области пространства. Один радар находится в Исландии - другой на юге Финляндии. Она похожа на радары STARE , но работает на более низкой частоте ~ 10 Мгц, поэтому лучи этих радаров могут заметно загибаться в ионизированной среде и находить области с нулевыми ракурсными углами.

Лучи финского радара направлены приблизительно по меридиану и ортогональны к авроральному электроджету, дугам сияний и границе диффузного свечения. При относительно спокойных магнитных условиях в вечерние и ночные часы луч этого радара подходит к границе диффузного свечения под ракурсными углами, близкими к нулю, т. е. здесь существуют оптимальные условия для набюдения отражений. В дальнейшем за счет загибания достаточно низкочастотного луча в ионизированной среде происходит ухудшение ракурсных условий и радар должен видеть только экваториальную кромку диффузного свечения. Согласно[Uspensky et al., 1999], в вечерние и ночные часы при небольших уровнях магнитной активности финский радар часто видит узкий сигнал из области, близкой к границе диффузного свечения.

На рис.21 приведены данные для 17 ноября 1996 года, AL-индекс для вечерних часов составлял ~ -100 нТл. Для этого уровня по формулам (2) и (3) были рассчитаны границы овала и диффузного свечения. На верхней панели рис.21 граница диффузного свечения нанесена штрихами, кружками показано положение радиоотражений [Uspensky et al.,1999]. Из рисунка видно, что начиная с 15.30 UT радиоотражения располагаются в районе экваториальной границы и медленно смещаются к экватору по мере приближения к полуночи за счет суточного дрейфа. На нижней пвнели рис.21 показано распределений по высоте, откуда видно, что отражения на границе диффузного свечения располагаются на высотах слоя Е. В период до 18.00 UT оптические наблюдения были невозможны из-за солнечной засветки, тем не менее радарные данные четко фиксируют границу диффузного свечения, которая совпадает с рассчетной. Камера всего неба на ст. Муонио начала работать в 18.00 UT, но из-за низкой чувствительности вначале не регистрировала границу диффузного свечения. Когда после 19.00 UT на аскафильмах можно видеть границу диффузного свечения, то она совпала как с радиоотражениями, так и с рассчетнми данными. В 19.30 UT началась суббуря и радар, наряду с другими отражениями продолжал фиксировать границу диффузного свечения. Таким образом, можно считать, что и для экваториальной границы диффузного свечения формулы (2) и (3) дают достаточно хорошее совпадение с экспериментальными данными.

5. ПЛАНЕТАРНАЯ МОДЕЛЬ ИНТЕНСИВНОСТИ СВЕЧЕНИЯ НА РАЗНЫХ ВЫСОТАХ

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

В качестве модели высыпаний были взяты данные [Spiro et al., 1982], где по набюдениям спутников AE-C и AE-D приведено пространственно-временное распределение потока энергии авроральных электронов в эрг. см^-2.с^-1 (W) и средней энергии в кэВ (Е₀). Всего было проведено около 31 тысяч измерений спектров авроральных электронов по данным 16 энергетических каналов в диапазоне 0,2 - 27 кэВ. Данные для четырех уровней магнитной возмущенности, определяемой AE-индкесом, были собраны в специальные таблицы с разрешением в 1 час по местному времени и с разрешением по широте в 1° для диапазона F ¢ =60-80° и в 2° для F ¢ = 50-60° и 80-88° . Полученные средние значения дополнительно сглаживались по соседним точкам.
Предполагалось, что авроральное свечение вызывается вторгающимися электронами. Методика расчетов свечения основана на исследовании прохождения авроральных электронов через атмосферу Земли, что позволяло оценить потери их кинетической энергии на упругое рассеяние, ионизацию с образованием вторичных электронов и неупругие соударения, приводящие к возбуждению вращательных, колебательных и электронных уровней атмосферных газов с последующим излучением. Расчет проводился методом Монте-Карло, что позволяло достаточно корректно, в отличии от расчетоа методом непрерывных потерь, исследовать процессы переноса электронов малых энергий (Е £ 1 кеВ). При расчетах учитывалось изменения величины дипольного поля Земли в процессе прохождения электронов через атмосферу. Сечения рассеяния электронов атмосферными газами задавалось в аналитической форме в соответствии с экспериментальными данными. Расчет интенсивности свечения проводился в мидели нейтральной атмосферы CIRA-72. Подробно методика расчета описана в [Иванов и Осипов, 1981; Юрова и Иванов, 1989; Иванов и Сергиенко, 1992; Sergienko and Ivanov, 1993].

Для исследования распределения свечения были выбраны в области ближнего ультрафиолета полосы первой отрицательной системы азота 1NGN (B²S - X²S ). Так как время жизни колебательных уровней B²S состояния (~ 10^-7с) намного меньше характерного времени процессов столкновительной дезактивации, то вариации свечения данных полос хорошо следуют за вариациями потоков авроральных электронов и пропорциональны скорости ионообразования N₂. В ближней инфракрасной области исследовалось пространственное распределение наиболее мощных молекулярных полос этой части спектра, а именно, для первой положительной системы N₂(B³P_g - AS ) и системы Мейнела GMN(A²P_u - XS ). Энергетический спектр задавался в экспоненциальной форме

dN = W/E× exp(-E/E) dE (7)

Значения W и Е₀ брались из таблиц, приведенных в [Spiro et al, 1982].

Был проведен расчет свечения атмосферы в вертикальном столбе от бесконечности до разных высотных уровней от 90 до 300 км. Для 90 км при тех средних энергиях авроральных электронов, которые приведены в [Spiro et al., 1982], интенсивность свечения соответствует полному свечению полярных сияний. Для других высотных уровней свечение соответствует свечению в вертикальном столбе от указанного уровня до бесконечности. За бесконечность принята высота 800 км, где плотность нейтральной атмосферы настолько мала, что авроральные частицы на этих высотах практически не теряют энергию. В качестве модели границ аврорального свечения были использованы математические аппроксимации, согласно [Старков, 1994б]. Пересчет AE-индекса, который использовался в [Spiro et al.,1982] в качестве характеристикм магнитной активности, в AL проводился согласно [Старков, 1994а].

Для сопоставления характеристик высыпания с границами аврорального свечения были построены меридиональные профили W, Eи J(H), где J(H) - рассчитанная интенсивность свечения в релеях для разных высот, а также границы овала и диффузного свечения. На рис.22 представлены 6 таких профилей для максимальной магнитной активности (AE > 600 нТл) для ночных и дневных часов. Характерной особенностью всех профилей является хорошее подобие конфигурации свечения для существенно разных высотных уровней, разумеется, при соответствующем изменении интенсивности. Ход свечения на разных высотах подобен изменению потока энергии.

В ночные часы наблюдается хорошее совпадение границ овала с максимумом потока энергии и свечения. Особенно четко это видно в 00 MLT. Надо учитывать, что данные о свечении и потоке энергии приведены в логарифмическом масштабе. Линейный масштаб не позволяет наглядно охватить весь диапазон изменений, но область максимума свечения в этом случае прорабатывается более резко. Границы ширины максимума свечения по половине интенсивности (зачерненные кружки на рис.22), как правило, для ночных часов лежат внутри овала, то есть полуширина максимума свечения меньше ширины овала. При этом экваториальная граница полуширины максимума свечения всегда располагается приполюсней экваториальной границы овала, а приполюсная граница иногда может располагаться в более высоких широтах. Такая ситуация наблюдается только в 21 и 22 MLT.

В дневные часы интенсивность свечения в 5-10 раз ниже, чем в ночные, но максимум свечения также совпадает с положением овала. Только в полуденные часы максимум свечения несколько смещен к полюсу, см. рис.15. Полуширина максимума интенсивности свечения обычно несколько шире овала.

Ход средней энергии авроральных электронов в дневные часы имеет своеобразный двухступенчатый характер. От полюса до 80° энергия частиц пракитчески постоянна, низка и равна ~ 0.4 кэВ. После 80° наблюдается резкое повышение средней энергии до 1.0 - 1.5 кэВ, и эта величина с некоторыми вариациями сохраняется до приполюсной границы овала, после чего энергия опять начинает быстро расти и достигает максимума порядка 3.0 - 3.5 кэВ уже экваториальнее овала. Наиболее четко эти вариации наблюдаются в 09, 10 и 12 MLT. На ночной стороне первый подъем средней энергии иногда наблюдается уже на 85° .

Область диффузного свечения в ночные и дневные часы характеризуются разными параметрами авроральных частиц. В ночные часы диффузное свечение связано как с понижением интенсивности свечения, так и с уменьшением средней энергии. Именно в диффузном свечении наблюдается минимум средней энергии, прчем иногда довольно глубокий, например в 22 MLT. Такая же ситуация наблюдается в 03 и 21 MLT. Этот результат хорошо согласуется с [Lui et al., 1977; Meng, 1978; Feldstein and Galperin, 1985], где показано, что в экваториальном диффузном свечении преобладают мягкие высыпания и красное свечение. Наиболее мягкие электроны набюдаются не на экваториальной границе диффузного свечения, а ближе к авроральному овалу. По мере приближения к экваториальной границе свечения средняя энергия электронов довольно быстро растет, причем рост продолжается и за границей диффузного свечения. Вероятно, этот эффект является результатом осреднения энергии мягких электронов диффузного свечения и жестких электронов радиационных поясов. Если внутри овала жесткость частиц пропорциональна потоку энергии, то для диффузного свечения наблюдается тенденция обратной зависимости.

В дневные часы, наоборот, диффузное свечение связано с повышением средней энергии, причем максимум наблюдается внутри диффузного свечения (9 и 12 MLT) или располагается за его экваториальной границей (14 и 15 MLT). Полученные результаты совпадают с фотометрическими наблюдениями в каспе, когда экваториальнее овала, где преобладают красные сияния, наблюдается повышенное свечение зеленой линии [Воробьев и др., 1984].

Интенсивность диффузного свечения на экваториальной границе в ночные часы составляет 100 R для l 127.8 нм, кроме разреза в 22 MLT, когда эта граница попала на минимум свечения. В среднем по шести разрезам свечение на экваториальной границе составляет (100± 70) R, что хорошо совпадает с другими экспериментальными данными, см. например, [Feldstein and Galperin]. На дневной стороне интенсивность свечения на экваториальной границе несколько ниже и составляет (80± 30) R.

За приполюсной границей овала интенсивность свечения падает, однако средняя энергия остается относительно высокой, после чего наблюдается резкий спад как на дневной, так и на ночной стороне до уровня 0.4 кэВ и ниже. Вероятно, эту довольно резкую границу, хорошо заметную на рис.22, можно принять за границу полярной шапки, а интервал широт от приполюсной границы овала до резкого спада Е отнести к диффузному свечению с приполюсной стороны овала. Для высокой магнитной активности радиус полярной шапки составляет ~ 8° , что совпадает с экспериментальными данными [Eather and Mende, 1971].

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

На рис.22 отчетливо видно, что вариации свечения повторяют изменения потока энергии и не зависят от средней энергии частиц. Для выявления зависимости интенсивности свечения от средней энергии была использована характеристика J/W, где J - интенсивность свечения, а W - вариации потока энергии. Она не зависит от местного времени, так как использовалась одна и та же модель нейтральной атмосферы, что приемлимо для исследования характеристик свечения в темное время при отсутствии дополнительной ионизации за счет ультрафиолетового излучения Солнца, когда возможны наземные оптические наблюдения. Расчет проводился для максимального величины магнитной активности для высотного уровня 90 км. Полученные результаты представлены на рис.23. Для всей светящейся толщи атмосферы при увеличении энергии частиц это отношение увеличивается с ростом средней энергии и при изменении Е⁰ от 0.5 до 3.0 кэВ увеличивается почти в два раза. Вариации интенсивности хорошо описываются формулой:

J = 10(1.45 + 0.89 lg E) × W,

где Е измеряется в кэВ, а W в эрг. см^-2 с^-1. Сплошная линия на рис.23 вычислена по приведенной формуле. Она позволяет по интенсивности свечения и одной из характеристик потока частиц определить второй параметр.

Планетарное распределение свечения l 427.8 нм для четырех значений магнитной активности и трех высотных уровней приведено на рис.24. Градации интенсивности свечения различаются на порядок.

Для высотного уровня 90 км видно хорошее совпадение области максимального свечения с границами аврорального овала, особенно для высокой магнитной активности (AE > 600 нТл и 600 > AE > 300нТл). Экваториальная граница максимального свечения практически точно совпадает с границей овала. Приполюсная граница свечения в ночные часы располагается несколько выше соответствующей границы овала.

Интересной особенностью является “разрыв” в дневные часы, который наблюдается при АЕ > 600 нТл, но отсутствует при 600 > АЕ > 300 нТл, т. е. при максимальной магнитной активности средняя интенсивность свечения в дневные часы ниже, чем при более низкой возмущенности, но область достаточно высокой интенсивности свечения (вторая градация интенсивности) шире. Таким образом, при максимальной магнитной активности яркость свечения в дневном секторе понижается но сама область свечения расширяется. При слабой магнитной активности (300 > АЕ > 100 нТл) и при максимальной магнитной активности для высотного уровня 150 км разрыв в дневные часы появляется снова. При понижении магнитной активности и повышении высотного уровня интенсивность свечения падает, но во всех случаях область максимального свечения наблюдаеся внутри аврорального овала и границы этой области, как правило, совпадают с границами овала.

Экваториальная граница диффузного свечения в свою очередь оконтуривает область с интенсивностью свечения на порядок ниже, хотя сопадение и менее хорошее, чем для овала. С уменьшением магнитной активности граница диффузного свечения и свечение с J @ 0.5 kR сжимаются синхронно. При AE > 600 нТл диффузная экваториальная полоса и свечение с J @ 0.5 kR имеют одинаковую асимметрию с расширением в предполуденные часы. Для слабой магнитной активности (300 > AE > 100 нТл) в дневные часы область свечения l 427.8 нм существенно шире полосы диффузного свечения. В остальные часы наблюдается довольно хорошее соответствие интенсивности свечения и положения границ. В полярной шапке свечение слабое, наблюдается тенденция некоторого повышения интенсивности с уменьшением магнитной активности.

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

Для высотного уровня 300 км (см. рис.24 в) интенсивность свечения уже низка и даже при максимальной магнитной активности не превышает 100 R, но область максимального свечения, как и раньше, имеет овальную форму, схожную с авроральным овалом. Для слабой магнитной активности и спокойных условий экваториальная граница овала довольно хорошо совпадает с границей свечения последней градации яркости, приполюсная граница свечения располагается в области более высоких широт. При высокой магнитной активности для этого уровня экваториальная и приполюсная границы свечения с J ~ 50 R (максимального для этой высоты) в ночном секторе хорошо совпадают с границами овала.

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

Были рассчитаны аналогичные пространственно-временные распределения свечения для наиболее сильных полос систем 1PGN (l 750.5 нм) и GMN (l 1108.7 нм). Общий характер распределения свечения сходен с 1NGN, но интенсивность свечения этих эмиссий выше, чем для l 427.8 нм. При спокойных магнитных условиях область повышенного субаврорального свечения в ночные часы для 1PGN и GMN тоже заметно шире, чем для l 427.8 нм.

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

---

---

назад		вперед		оглавление		литература
		другие обзоры

---

Для связи:
lll@srd.sinp.msu.ru
пароль: "сизиф почти не виден"