Модуляция космических лучей в период минимума 24 солнечного

advertisement
Солнечно-земная физика. Вып. 20. (2012) С. 36–39
УДК 523.745, 539.1633
ОСОБЕННОСТИ МОДУЛЯЦИИ КОСМИЧЕСКИХ ЛУЧЕЙ
В МИНИМУМЕ 24-го СОЛНЕЧНОГО ЦИКЛА
В.М. Дворников†, М.В. Кравцова, В.Е. Сдобнов
PECULIARITIES OF COSMIC-RAY MODULATIONS DURING THE MINIMUM
OF THE SOLAR CYCLE 24
V.M. Dvornikov†, M.V. Kravtsova, V.E. Sdobnov
По данным наблюдений космических лучей (КЛ) на мировой сети станций и космических аппаратах показано, что
определяющим фактором увеличения потоков галактических КЛ на орбите Земли в период минимума 24-го солнечного
цикла является уменьшение потерь энергии при движении частиц в регулярных электромагнитных полях гелиосферы.
Cosmic-ray observations from the worldwide network of stations and from space vehicles have been used to demonstrate that
the crucial factor of the increase in galactic CR fluxes at the Earth’s orbit during the solar cycle 24 minimum is the reduction of
energy losses when particles travel in regular electromagnetic fields of the heliosphere.
_A_STAR/]. Так как магнитные поля на Солнце являются источником гелиосферного магнитного поля, то оно тоже ослабло и на орбите Земли в 2009 г.
неоднократно опускалось ниже 4 нТл. Это на ~1 нТл
ниже, чем в предшествовавшие минимумы активности [http://omniweb.gsfc.nasa.gov/]. В связи с этим
наблюдаемые потоки галактических КЛ достигли
самого высокого уровня за весь период, в течение
которого проводились регулярные измерения космического излучения [Свиржевская и др., 2010].
Так, например, по данным наблюдений на ст. Иркутск (435 м над уровнем моря) интенсивность
нейтронной компоненты КЛ в октябре 2009 г. повысилась на ~2 % по сравнению с апрелем 2008 г. Повидимому, природа данного явления аналогична
природе эффекта, описанного в работе [Базилевская
и др., 2011], где показано, что в минимуме между
23-м и 24-м циклами солнечной активности модуляция потоков КЛ была минимальной за весь период
прямых наблюдений КЛ. В настоящей статье по
наблюдениям интенсивности КЛ на мировой сети
станций и на космическом аппарате GOES на основе
предположения, что определяющим фактором в модуляции КЛ является изменение их энергии при
движении в регулярных электромагнитных полях
гелиосферы, проведены расчеты дифференциального жесткостного спектра КЛ на орбите Земли в апреле 2008 и октябре 2009 г.
Введение
Для описания модуляционных эффектов космических лучей (КЛ) в гелиосфере принципиально
возможны два подхода. Согласно первому из них
модуляция КЛ обусловлена изменением плотности
частиц за счет их выноса мелкомасштабными магнитными неоднородностями, вмороженными в
плазму солнечного ветра (СВ), и диффузией вдоль
спирального межпланетного магнитного поля
(ММП) [Parker, 1963]. Амплитуда модуляции в линейном приближении пропорциональна скорости
СВ, размеру области модуляции и обратно пропорциональна коэффициенту диффузии. Анизотропия в
распределении частиц по направлениям прихода
мала и в основном имеет токовую структуру [Крымский, 1969].
В соответствии со вторым подходом предполагается, что ММП в большей степени структурировано,
чем турбулизовано. Определяющим фактором в модуляции КЛ является изменение их энергии при
движении в регулярных электромагнитных полях
гелиосферы. Изменение же энергии, в свою очередь,
определяется величиной и природой электрических
полей межпланетной среды и временем взаимодействия частиц с этими полями [Dvornikov, Sdobnov,
2002], которое в значительной мере определяется
структурами ММП типа магнитных ловушек. Анизотропия КЛ может достигать значительной величины
и, наряду с первой, иметь более высокие гармоники
[Dvornikov et al., 1983; Richardson et al., 2000].
Минимум между циклами солнечной активности
23 и 24 был необычно продолжительным и глубоким [ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUN
SPOT_NUMBERS/]. Этот период, названный в работе
[Базилевская и др., 2011] минимумом цикла 24, характеризуется минимальным за последние несколько
солнечных циклов сглаженным числом солнечных
пятен Rz. В работе [Базилевская и др., 2011] отмечается, что магнитное поле солнечных пятен систематически уменьшается с 1992 г. независимо от фазы
цикла [Livingston, Penn, 2009]. Магнитное поле
Солнца как звезды тоже уменьшилось. В мае 2009 г.
оно составляло 3.6 мкТл, тогда как в предшествующие минимумы активности было ~6–7 мкТл
[ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUN_AS
Модуляция космических лучей регулярными
электромагнитными полями межпланетной среды
Для описания модуляционных эффектов в рамках второго подхода использовалось кинетическое
уравнение в дрейфовом приближении без учета эффектов выноса КЛ мелкомасштабными магнитными
неоднородностями, вмороженными в плазму СВ, и
диффузией частиц вдоль спирального ММП. Следует отметить, что данное приближение не исключает
наличия рассеяния частиц на магнитных неоднородностях, но предполагает, что этот процесс не является марковским. В этом случае возможна ситуация,
когда при ненулевом среднеквадратичном отклонении частиц в угловых и пространственных координатах их средние отклонения будут близки нулю.
При таком характере рассеяния необходимо учиты36
В.М. Дворников, М.В. Кравцова, В.Е. Сдобнов
вать эффекты, накапливающиеся при многократном
вращении частиц по ларморовской окружности, т. е.
использовать дрейфовое приближение, пренебрегая
эффектами конвекции и диффузии. В данном случае
изменение интенсивности КЛ в гелиосфере происходит вследствие изменения энергии в регулярных
электромагнитных полях межпланетной среды в
соответствии с теоремой Лиувилля, которая в дрейфовом приближении имеет следующий вид [Сивухин, 1963]:
F
F
F
F
r
 p||
 p
 0,
t
r
p||
p
vR

2
 vμ cos   u sin   uB 
f
2
2
 cos λ cos  sin  (  1) 
r
1 vR
2
2 f

sin  cos  (1  cos  )(1   ) 
2
2 B0 r0


r  v|| h 
B
2
EB 
1 cp v
2 ze B
(1)
cp
1
p||  ze Eh    ph  || Eh  (hh) 
2
B
cp p   B

 ||   
h  (h)h  h  (h)h  .
2 ze B  B

1
p   v p|| h 
2
cp 
B

  h  E  Eh 
 Eh(h)  h  
2B 
B


cp|| p v||  B
где μ=cosθ, θ – питч-угол частицы в ММП;
 r sin λ  – угол между радиусом и си  arctan 

 u 
ловой линией ММП.
При изменении структуры поля соответствующим образом изменятся и выражения перед частными производными в данном уравнении и, следовательно, изменится вид функции распределения частиц. Таким образом, в предлагаемом подходе к
проблеме модуляции КЛ вариации функции распределения частиц отражают изменения в структуре
электромагнитных полей межпланетной среды.
При движении в электромагнитном поле (5), (6)
заряженные частицы космического излучения,
наряду с движением вдоль поля и дрейфом под действием силы zeE, будут дрейфовать в направлении,
обратном направлению электрического поля, и в
результате этого терять часть своей кинетической
энергии. Дрейф в указанном направлении обусловлен двумя причинами: во-первых, кривизной силовых линий (центробежный дрейф) и, во-вторых,
наличием нормальной по отношению к силовой линии компоненты градиента поля (магнитный или
градиентный дрейф).
Изменение жесткости частицы R за счет данного
процесса описывается уравнением [Дворников, Матюхин, 1976]

cp
B
 h(h  h)  h    v|| || h  (h)h.
B 
ze B

(2)
(3)
(4)

h  (h)h  h  (h )h  .
2 ze B  B

p|| , v|| и p , v – составляющие импульса и скорости частиц, параллельные и перпендикулярные вектору магнитного поля; ze – их заряд; с – скорость
света; h – единичный вектор в направлении поля В;
r – пространственный вектор; t – время.
Отметим, что функция распределения F относится не к распределению самих частиц, а к распределению их ведущих центров.
При стационарном и пространственно-однородном
СВ радиальная и азимутальная компоненты ММП
имеют вид [Parker, 1963]



dR 1 u

R 1  cos 4  1  2 .
dt 2 r
(8)
Для решения уравнения (7) произведем замену
переменных в уравнении (8).
В связи с тем, что изменение энергии происходит
вследствие центробежного и градиентного дрейфа
по координате  антипараллельно индуцированному
электрическому полю, в (8) можно перейти от переменной t к переменной :
2
r 
Br (r )  Br (r0 )  0  ;
r
r
Bφ (r , λ)=Br (r )
sin λ,
u
1v
cos  (1  cos 2  )(1   2 ) 
2r
f 1 u
f
 
(1  cos 4  )(1   2 ) p
 0,
 2 r
p
где
c
(7)
(5)
где u – скорость СВ,  – угловая скорость вращения
Солнца, а индуцированное электрическое поле
dR
1 u
λ=  R (1  cos 4 ψ)(1+μ 2 ).
dλ
2 r
1
E   u  B.
(6)
c
С учетом (2)–(4) и (5), (6) из (1) получим уравнение (в переменных r, p, ), описывающее влияние
регулярных полей СВ на функцию распределения
КЛ в межпланетном пространстве F(r, p, ):

1 vR
2
sin  cos  
2 B0 r0
(10)
(1  cos  )(1   ).
2
2
С учетом (10) уравнение (9) примет вид
37
(9)
Особенности модуляции космических лучей в минимуме 24-го солнечного цикла
pdp
c 2 p 2  ε 02

ze B0 r02
c
sin λdλ.
(11)
Проинтегрировав (11) по λ, получим
ε 
ze B0 r02
c
(1  cos  E )  zeU ,
(12)
где U=ΩB 0 r 0 2 (1–cosλ E ) – потенциал электрического поля (6) на гелиошироте E.
Выражение (12) показывает, насколько изменится энергия частицы при дрейфе от полюса до плоскости гелиоэкватора.
Из уравнения (12) следует, что потери энергии
частиц при потенциальном электрическом поле не
зависят от жесткости частиц и при характерной
напряженности ММП (5 нТл) составляют 200 МэВ.
В том случае, когда ММП имеет структуру, отличную от спиральной (например, при наличии петлеобразных структур или магнитных облаков, т. е. магнитных ловушек), значения ∆ε будут зависеть от напряженности магнитных полей данных структур и скорости СВ и почти на порядок превосходить 200 МэВ
[Dvornikov, Sdobnov, 2002].
В соответствии с предлагаемой концепцией, если
жесткостной спектр в Галактике описывается степенной функцией от жесткости частиц R, то из
условия J(R)/p 2 =const следует, что жесткостной
спектр в гелиосфере будет иметь следующий вид:
(ε 2  ε02 )
3/ 2
 
ε  ε
J ( R)  A
2
2
ε (ε  ε)  ε 0  T  ε 0
Рис. 1. Расчетные спектры КЛ при различных значениях Δε.
ции отсчитывались от спокойного уровня апреля
2008 г.
Для определения зависимости  от R методом
минимизации функционала

S    J (R

k
k
A
3/ 2
 
ε  ε
2
2 T  ε

ε (ε  ε)  ε 0 
0
2
2
0
γ



2
решалось уравнение (13) путем подстановки наблюдаемых значений интенсивности частиц в левую
часть уравнения для каждого из используемых интервалов энергий с целью определения неизвестной величины , соответствующей данному энергетическому
интервалу. В результате были определены зависимость
 от жесткости частиц, соответствующей среднему
значению энергии в используемых интервалах энергии, т. е. (R), и жесткостной спектр протонов на
орбите Земли в апреле 2008 г. и октябре 2009 г.
На рис. 2 представлены часовые значения изменения энергии КЛ в электромагнитных полях гелиосферы в апреле 2008 г. и в октябре 2009 г.
γ
,
(ε  ε )
(13)
где  – полная энергия частиц с жесткостью R;  –
ее изменения в электромагнитных полях гелиосферы; 0 – энергия покоя; T0 – кинетическая энергия
частиц, при которой интенсивность КЛ соответствующей жесткости в Галактике равна А; А и  –
спектральные индексы галактического спектра.
Полученное выражение удовлетворяет кинетическому уравнению в дрейфовом приближении, т. е.
теореме Лиувилля. Это легко проверяется подстановкой выражения для жесткостного дифференциального спектра (13), с учетом (12), в уравнение
(7). На рис. 1 представлены расчетные спектры
КЛ при различных значениях Δε, соответствующих наблюдаемым значениям модуля ММП и
скорости СВ, полученные с использованием выражения (13).
Результаты и выводы
В настоящей работе использованы данные наблюдений интенсивности протонов на спутнике GOES-11
в энергетических диапазонах 4–9, 9–15, 15–40, 40–80,
80–165 и 165–500 МэВ [http://spidr.ngdc.gov/spidr/
index.html]. При более высоких энергиях использовалась информация о вариациях первичного спектра КЛ,
полученная методом спектрографической глобальной
съемки (СГС) [Dvornikov, Sdobnov, 1997; 2002] по
наземным измерениям на мировой сети нейтронных мониторов (44 станции). Амплитуды модуля-
Рис. 2. Изменения энергии КЛ в электромагнитных
полях гелиосферы в апреле 2008 г. и в октябре 2009 г.
38
В.М. Дворников, М.В. Кравцова, В.Е. Сдобнов
СПИСОК ЛИТЕРАТУРЫ
Базилевская Г.А., Крайнев М.Б., Махмутов В.С. и др.
Особенности вариаций космических лучей в фазе минимума между 23 и 24 солнечными циклами // Известия
РАН. Сер. физ. 2011. Т. 75, № 6. С. 831–834.
Дворников В.М., Матюхин Ю.Г. Энергетические потери космических лучей при движении в регулярном поле
солнечного ветра // Изв. АН СССР. Сер. физ. 1976. Т. 39,
№ 3. С. 624–626.
Крымский Г.Ф. Модуляция космических лучей в межпланетном пространстве. М.: Наука, 1969. 152 с.
Свиржевская А.К., Базилевская Г.А., Крайнев М.Б. и
др. Особенности минимума 24-го солнечного цикла в солнечных и гелиосферных характеристиках и в интенсивности галактических космических лучей // Научная сессия
МИФИ-2010: Сборник научных трудов. Т. 4. М.: НИЯУ
МИФИ, 2010. С. 95–98.
Сивухин Д.В. Дрейфовая теория движения заряженных частиц в электромагнитных полях // Вопросы теории
плазмы. М.: Госатомиздат, 1963. Т. 1. 287 с.
Dvornikov V.M., Sdobnov V.E. Time variations of the
cosmic ray distribution function during a solar event of
September 29, 1989 // J. Geophys. Res. 1997. V. 102, N A11.
P. 24209–24219.
Dvornikov V.M., Sdobnov V.E. Variations in the rigidity
spectrum and anisotropy of cosmic rays at the period of Forbush effect on 12–15 July 1982 // International Journal of
Geomagnetism and Aeronomy. 2002. V. 3, N 3. P. 217–226.
Dvornikov V.M., Sdobnov V.E., Sergeev A.V. Analysis of
cosmic ray pitch-angle anisotropy during the Forbush-effect in
June 1972 by the method of spectrographic global survey //
Proc. 18th ICRC. Bangalor, 1983. V. 3. P. 249–252.
Glesson L.J., Axford W.I. Solar modulation of galactic
cosmic rays // Astrophys. J. 1968. V. 154. P. 111.
Livingston W., Penn M. Are sunspots different during this
solar minimum? // Eos. Trans. Amer. Geophys. Union. 2009.
V. 90. P. 257–264. doi:10.1029/2009EO300001.
Parker E.N. Interplanetary dynamical processes. New
York; London: Inerscience Publishers a Division of John
Wiley and Sons, 1963. 362 с.
Richardson I.G., Dvornikov V.M., Sdobnov V.E., Cane H.V.
Bidirectional particle flows at cosmic ray and lower (1 MeV)
energies and their association with interplanetary coronal mass
ejections // J. Geophys. Res. 2000. V. 105, N A6. Р. 12579.
Usoskin I.G., Alanko-Yuotari K., Kovaltsov G.A., Mursula K. Heliospheric modulation of cosmic rays: Monthly reconstruction for 1951–2004 // J. Geophys. Res. 2005. V. 110.
A12108. doi:10.1029/2005JA011250.
ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_
NUMBERS/
ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUN_AS_A
_STAR/
http://omniweb.gsfc.nasa.gov/
http://spidr.ngdc.gov/spidr/index.html
Рис. 3. Дифференциальные спектры протонов на орбите Земли в апреле 2008 г. и октябре 2009 г.
Изменения энергии КЛ в электромагнитных полях гелиосферы в апреле 2008 г. в среднем составили
~0.50 ГэВ, а в октябре 2009 г. – ~0.37 ГэВ. Эти величины близки значениям модуляционного параметра для периодов минимума солнечной активности, полученным в работе [Usoskin et al., 2005] на
основе модели силового поля (force field model)
[Glesson, Axford, 1968]. Результаты расчетов жесткостного спектра протонов показаны на рис. 3, где
штриховая кривая – жесткостной спектр протонов в
апреле 2008 г., сплошная – в октябре 2009, треугольники и кружки – данные наблюдений.
Видно, что результаты расчетов хорошо соответствуют данным наблюдений мировой сети станций
космических лучей, а изменения энергии КЛ в электромагнитных полях гелиосферы в апреле 2008 г. и
октябре 2009 г. не противоречат значениям модуляционного параметра для периодов минимума солнечной активности, полученным на основе модели
силового поля.
Из вышеизложенного можно сделать вывод о
том, что определяющим фактором увеличения потоков галактических КЛ на орбите Земли в октябре
2009 г. по сравнению с апрелем 2008 г. является
уменьшение потерь энергии при движении частиц в
регулярных электромагнитных полях гелиосферы.
Данная работа поддержана программой Президиума РАН «Фундаментальные свойства материи и
астрофизика» в рамках проекта «Космические лучи
в гелиосферных процессах по наземным и стратосферным наблюдениям».
Институт солнечно-земной физики СО РАН, Иркутск, Россия
39
Download