Известия РАН. Серия физическая, 2021, T. 85, № 11, стр. 1623-1626

Форбуш-эффект и геомагнитная буря в апреле 1990 г.

М. В. Кравцова 1*, С. В. Олемской 1, В. Е. Сдобнов 1

1 Федеральное государственное бюджетное учреждение науки Ордена Трудового Красного Знамени Институт солнечно-земной физики Сибирского отделения Российской академии наук
Иркутск, Россия

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

Поступила в редакцию 21.06.2021
После доработки 05.07.2021
Принята к публикации 28.07.2021

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

Аннотация

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

ВВЕДЕНИЕ

В период умеренной солнечной активности [1] 9 апреля 1990 г. произошла геомагнитная буря с внезапным началом ~08:42 UT, вызванная приходом ударной волны от солнечной вспышки в активной области 6007 (координаты 24° с.ш., 28° в.д.). Вспышка сопровождалась радиовсплеском II типа, что указывает на наличие коронального выброса массы [2, 3]. К сожалению, в данный период имеются только отрывочные данные наблюдений межпланетного магнитного поля (ММП), но можно отметить, что скорость солнечного ветра (СВ) на орбите Земли составляла ~500 км/с и выше, модуль ММП и его Bz-составляющая в 19:00 UT 10 апреля доходили до ~28 и ~–25 нТл, соответственно. В этот период на мировой сети станций космических лучей (КЛ) наблюдался Форбуш-эффект (ФЭ). Следует отметить, что максимальное понижение интенсивности КЛ произошло почти на сутки раньше максимальной депрессии геомагнитного поля. Так максимальное понижение скорости счета нейтронного монитора (НМ) в Иркутске составило ~–8% в 21:00 UT 9 апреля, а минимальное значение Dst-индекса составило –280 нТл 10 апреля 1990 г. в 19:00 UT.

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

ДАННЫЕ И МЕТОД ИССЛЕДОВАНИЯ

Для анализа использовались данные мировой сети станций КЛ (42 НМ) [4], исправленные на давление и усредненные за часовые интервалы, а также данные космического аппарата GOES-7 (протоны в семи энергетических интервалах 0.8–4, 4–9, 9–15, 15–40, 40–80, 80–165 и 165–500 МэВ) [5]. Амплитуды модуляции отсчитывались от фонового уровня 8 апреля 1990 г. Выбор спокойного периода обусловлен тем, что в этот период электромагнитная обстановка в межпланетном пространстве и геомагнитная обстановка были спокойными, а спектр галактических КЛ наименее модулирован.

Анализ выполнен с применением метода спектрографической глобальной съемки (СГС) [6, 7]. Результаты интерпретированы в рамках модели модуляции КЛ регулярными электромагнитными полями гелиосферы [8] и осесимметричной модели ограниченной магнитосферы Земли [9].

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

На рис. 1 приведены амплитуды вариаций нейтронной компоненты КЛ (ΔI/I) на ст. Иркутск (Rc = 3.64 ГВ), амплитуды вариаций КЛ с жесткостью 10 ГВ на границе магнитосферы (ΔJ/J), амплитуды первой (А1) и второй (А2) гармоник питч-угловой анизотропии для частиц с жесткостью 4 ГВ, а также изменения жесткости геомагнитного обрезания (ЖГО) (ΔR) в Иркутске совместно с Dst-индексом 6–15 апреля 1990 г. Амплитуда понижения интенсивности КЛ на ст. Иркутск (см. рис. 1а) 9 апреля в 21:00 UT составила ~ –8%, а амплитуда эффекта в это время для частиц с жесткостью 10 ГВ на границе магнитосферы составила ~–15%. Максимальные значения А1 (более 20%) наблюдались 7, 8 и 10 апреля, а А2 (~15%) – 9 апреля (см. рис. 1б). В рамках модели модуляции КЛ регулярными электромагнитными полями гелиосферы [8] возрастание амплитуд A1 и A2 (двунаправленной) наблюдается при изменениях напряженности ММП. Перед началом ФЭ амплитуды А1 наблюдались в отдельные часы 06:00–08:00 UT, а также перед понижением интенсивности КЛ в конце 9–начале 10 апреля. Появление А2 указывает на пребывание Земли в магнитном поле петлеобразной структуры. В это время повышенный поток частиц наблюдался из направлений с координатами в солнечно-эклиптической геоцентрической системе ψ = ~110°, λ = ~20° и ψ = ~270°, λ = = ~–20°.

Рис. 1.

Вариации нейтронной компоненты КЛ на ст. Иркутск (сплошная кривая) совместно с вариациями изотропной составляющей интенсивности первичных КЛ с жесткостью 10 ГВ (штриховая кривая) (а). Амплитуды первой A1 (сплошная кривая) и второй A2 (штриховая кривая) гармоник питч-углового распределения КЛ с жесткостью 4 ГВ (б); Временной ход изменений ЖГО в Иркутске (Rc = 3.64 ГВ) (сплошная кривая) совместно с Dst-индексом (штриховая кривая) (в).

Изменения ЖГО в Иркутске в период магнитной бури, сопровождавшей ФЭ, доходили до ~‒1.1 ГВ (см. рис. 1в). Максимальное понижение ЖГО и Dst-индекса произошло 10 апреля в 15:00 UT и в 19:00 UT, соответственно. Коэффициент корреляции между изменениями ЖГО в Иркутске и Dst -индексом составил ~0.82. ФЭ на Земле произошел на ~18 ч раньше, чем максимальное понижение ЖГО и на ~22 ч раньше, чем депрессия Dst-индекса. Причиной большой разницы начала ФЭ и геомагнитный бури согласно [10, 11] является геометрический фактор пересечения Землей переднего фронта крупномасштабных возмущений СВ.

Спектры вариаций КЛ в период рассматриваемого ФЭ не являются степенными в широкой области жесткостей. Только в диапазоне жесткостей выше ~10 ГВ эти спектры по жесткости близки к степенным. Получены средние значения показателя степени γ на разных стадиях события при аппроксимации жесткостных спектров вариаций степенной функцией от жесткости частиц в диапазоне жесткостей выше 10 ГВ. В моменты максимальной модуляции КЛ спектры вариаций мягче (γ = 0.93), чем на фазах спада (γ = 0.80), и восстановления (γ = 0.87) интенсивности КЛ.

На основе рассчитанных планетарных изменении ЖГО в период магнитной бури в апреле 1990 г. в рамках осесимметричной модели ограниченной магнитосферы, учитывающей токи на магнитопаузе и кольцевой ток [9], мы рассчитали расстояние до подсолнечной точки и радиус кольцевого тока, силу кольцевого тока (IDR), тока на магнитопаузе (IDCF), а также соотношение между вкладами от этих токов в Dst-индекс на разных фазах магнитной бури.

На рис. 2 представлены изменения ЖГО в отдельные моменты на разных фазах геомагнитной бури в апреле 1990 г. в зависимости от ЖГО. Видно, что на фазах наиболее интенсивного развития магнитной бури и на фазах ее затухания результаты расчетов изменений ЖГО по данным наблюдений интенсивности КЛ на мировой сети станций хорошо описываются осесимметричной моделью ограниченной магнитосферы [9], включающую только кольцевой ток и токи на магнитопаузе.

Рис. 2.

Зависимость изменения ЖГО от ЖГО на разных фазах геомагнитной бури в апреле 1990 г. Сплошная линия – результаты, полученные из данных мировой сети станций КЛ методом СГС, штрих-пунктирная линия – расчет по осесимметричной модели ограниченной магнитосферы с кольцевым током. Штриховая линия – вклад в изменения ЖГО от кольцевого тока, полученного по осесимметричной модели ограниченной магнитосферы.

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

Таблица 1.  

Вклад DR и DCF-токов в развитии магнитной бури в апреле 1990 г. на разных стадиях эффекта

Дата Dstнабл, нТл Dstрасч, нТл DstDR, нТл DstDCF, нТл IDR, 106 A IDCF, 106 A
10.04.1990 г.
07:00 UT
–65 –38 –171 133 14.9 4.9
10.04.1990 г.
20:00 UT
–257 –254 –386 132 15.4 3.1
13.04.1990 г.
23:00 UT
–93 –118 –236 118 13.5 3.1

В табл. 2 приведены средние радиусы кольцевого тока (rс) в единицах радиуса Земли и средние расстояния до подсолнечной точки (rm) в рамках осесимметричной модели ограниченной магнитосферы для разных фаз развития геомагнитной бури. Минимальные радиусы кольцевого тока наблюдаются в главной фазе геомагнитной бури с последующим возрастанием на фазе восстановления геомагнитного возмущения. Поскольку в осесимметричной модели ограниченной магнитосферы не учитываются вклад многих токовых систем, а рассчитанные значения Dst-индекса в некоторые моменты близки к наблюдаемым, можно предположить, что в эти периоды магнитной бури основной вклад в ее интенсивность вносят симметричная составляющая DR и DCF-токов.

Таблица 2.  

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

10 апреля 1990 г. 13 апреля 1990 г.
начальная фаза в 07:00 UT главная фаза в 20:00 UT фаза восстановления в 23:00 UT
rс rm rс rm rс rm
~4.8 ~8.4 ~3.8 ~8.5 ~4.2 ~8.7

ЗАКЛЮЧЕНИЕ

Спектры вариаций КЛ в период ФЭ не являются степенными в широкой области жесткостей, только выше ~10 ГВ жесткостные спектры вариаций КЛ близки к степенным. При аппроксимации спектров вариаций в период ФЭ степенной функцией от жесткости частиц при жесткости более 10 ГВ на фазе максимальной модуляции интенсивности КЛ наблюдается наибольший показатель спектра. В период повышенной амплитуды двунаправленной питч-угловой анизотропии КЛ в апреле 1990 г. Земля находилась в петлеобразной структуре ММП. Показана возможность по данным наземных измерений КЛ на мировой сети станций проводить количественную оценку радиуса и силу кольцевого тока, а также рассчитывать расстояние до подсолнечной точки и силу тока на магнитопаузе в периоды геомагнитных возмущений.

Работа выполнена при финансовой поддержке Минобрнауки России. Результаты получены с использованием оборудования Центра коллективного пользования “Ангара” http://ckp-rf.ru/ckp/3056/ и Уникальной научной установки “Российская национальная наземная сеть станций космических лучей” (Сеть СКЛ) [12].

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

  1. Antalova A. // Contrib. Astron. Obs. Skalnate Pleso. 1991. V. 21. P. 205.

  2. Wild J. P., Smerd, S.F., Weiss A.A. // Annu. Rev. Astron. Astrophys. 1963. V. 1. P. 291.

  3. Gopalswamy N. // In: The Sun and the Heliosphere as an integrated system. Dordrecht: Springer, 2004. V. 317. P. 201.

  4. http://center.stelab.nagoya-u.ac.jp/WDCCR.

  5. http://satdat.ngdc.noaa.gov./sem/goes/data/new_avg.

  6. Richardson I.G., Dvornikov V.M., Sdobnov V.E., Cane H.V. // JGR. 2000. V. 105. Art. No. 12579.

  7. Кравцова М.В., Сдобнов В.Е. // Изв. РАН. Сер. физ. 2017. Т. 81. № 2. С. 196; Kravtsova M.V., Sdobnov V.E. // Bull. Russ. Acad. Sci. Phys. 2017. V. 81. No. 2. P. 177.

  8. Дворников В.М., Кравцова М.В., Сдобнов В.Е. // Геомагн. и аэроном. 2013. Т. 53. № 4. С. 457; Dvornikov V.M., Kravtsova M.V., Sdobnov V.E. // Geomagn. Aeron. 2013. V. 53. No. 4. P. 430.

  9. Кичигин Г.Н., Сдобнов В.Е. // Геомагн. и аэроном. 2017. Т. 57. № 2. С. 149; Kichigin G.N., Sdobnov V.E. // Geomagn. Aeron. 2017. V. 57. No. 2. P. 132.

  10. Шадрина Л.П., Мамрукова В.П., Плотников И.Я. // Геомагн. и аэроном. 1996. Т. 36. № 3. С. 169.

  11. Shadrina L.P., Starodubtsev S.A., Plotnikov I.Ya. // IJGA. 2001. V. 2. No. 3. P.

  12. https://ckp-rf.ru/usu/433536.

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