Известия РАН. Серия физическая, 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
Аннотация
По данным наземных измерений космических лучей на мировой сети станций исследованы Форбуш-эффект и геомагнитная буря в апреле 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.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], включающую только кольцевой ток и токи на магнитопаузе.
В табл. 1 представлены наблюдаемые (Dst набл) и рассчитанные (Dst расч) значения Dst-индекса на разные фазах развития магнитной бури, а также вклад кольцевого тока (Dst DR) и токов на магнитопаузе (Dst DCF) в полную величину Dst-индекса и силу тока этих токовых систем. Рассчитанные параметры кольцевого тока и тока на магнитопаузе отражают суммарные вклады нескольких токовых систем. Поэтому правильнее говорить об эффективных радиусах токовых систем, дающих вклады в Dst-индекс подобно кольцевому току и току на магнитопаузе.
Таблица 1.
Дата | 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-токов.
ЗАКЛЮЧЕНИЕ
Спектры вариаций КЛ в период ФЭ не являются степенными в широкой области жесткостей, только выше ~10 ГВ жесткостные спектры вариаций КЛ близки к степенным. При аппроксимации спектров вариаций в период ФЭ степенной функцией от жесткости частиц при жесткости более 10 ГВ на фазе максимальной модуляции интенсивности КЛ наблюдается наибольший показатель спектра. В период повышенной амплитуды двунаправленной питч-угловой анизотропии КЛ в апреле 1990 г. Земля находилась в петлеобразной структуре ММП. Показана возможность по данным наземных измерений КЛ на мировой сети станций проводить количественную оценку радиуса и силу кольцевого тока, а также рассчитывать расстояние до подсолнечной точки и силу тока на магнитопаузе в периоды геомагнитных возмущений.
Работа выполнена при финансовой поддержке Минобрнауки России. Результаты получены с использованием оборудования Центра коллективного пользования “Ангара” http://ckp-rf.ru/ckp/3056/ и Уникальной научной установки “Российская национальная наземная сеть станций космических лучей” (Сеть СКЛ) [12].
Список литературы
Antalova A. // Contrib. Astron. Obs. Skalnate Pleso. 1991. V. 21. P. 205.
Wild J. P., Smerd, S.F., Weiss A.A. // Annu. Rev. Astron. Astrophys. 1963. V. 1. P. 291.
Gopalswamy N. // In: The Sun and the Heliosphere as an integrated system. Dordrecht: Springer, 2004. V. 317. P. 201.
http://center.stelab.nagoya-u.ac.jp/WDCCR.
http://satdat.ngdc.noaa.gov./sem/goes/data/new_avg.
Richardson I.G., Dvornikov V.M., Sdobnov V.E., Cane H.V. // JGR. 2000. V. 105. Art. No. 12579.
Кравцова М.В., Сдобнов В.Е. // Изв. РАН. Сер. физ. 2017. Т. 81. № 2. С. 196; Kravtsova M.V., Sdobnov V.E. // Bull. Russ. Acad. Sci. Phys. 2017. V. 81. No. 2. P. 177.
Дворников В.М., Кравцова М.В., Сдобнов В.Е. // Геомагн. и аэроном. 2013. Т. 53. № 4. С. 457; Dvornikov V.M., Kravtsova M.V., Sdobnov V.E. // Geomagn. Aeron. 2013. V. 53. No. 4. P. 430.
Кичигин Г.Н., Сдобнов В.Е. // Геомагн. и аэроном. 2017. Т. 57. № 2. С. 149; Kichigin G.N., Sdobnov V.E. // Geomagn. Aeron. 2017. V. 57. No. 2. P. 132.
Шадрина Л.П., Мамрукова В.П., Плотников И.Я. // Геомагн. и аэроном. 1996. Т. 36. № 3. С. 169.
Shadrina L.P., Starodubtsev S.A., Plotnikov I.Ya. // IJGA. 2001. V. 2. No. 3. P.
https://ckp-rf.ru/usu/433536.
Дополнительные материалы отсутствуют.
Инструменты
Известия РАН. Серия физическая