Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Алания М.В. Квазипериодические вариации космических лучей

.pdf
Скачиваний:
23
Добавлен:
23.10.2023
Размер:
5.11 Mб
Скачать

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

Из рис. 5 следует, что в области г,:Д0,1% для станции Норикура (^нейтрошшя =0,29% , ажесткпя = 0,21 % см. таблицу. 2 из §2, а также работу [15]) суточная вариация имеет в основном статистический характер как для жесткой, так и для нейтронной компоненты.

При гх ~ 0,2% доля статистической вариации состав­ ляет 15—20% для нейтронной и 10—15% Для жесткой ком­ понент. При гх —0,3% доля статистической вариации стано­ вится уже очень малой (<Л%, т. е., другими словами, в этом

случае вероятность того, что наблюдаемая вариация статис­ тического происхождения меньше 1%).

Рис. 6. а) Дифференциальное частотное распределение максимума фаз первой гармсикки суточной вариации нейтронной компоненты интенсив­ ности космических лучей по данным станции Норикура за период МГГ. б) Интегральное частотное распределение максимума фаз первой гармо­ ники суточной вариации нейтронной компоненты интенсивности космиче­

ских лучей по данным станции Норикура за период МГГ

21

Перейдем к рассмотрению частотных распределений фаз. Мы нашли их для нейтронной компоненты станции Норпку-

ра. На рис. 6а показаны частотные

распределения

фаз

при

изменениях

гг от

0,1%

до 1,5%.

На

рис. 66 представлены

частотные

распределения

для всех rt (i — 1,2- ■• -529) (сплош­

ная кривая)

и для

0,2% (пунктирная кривая).

1Тз

рис.

6а видно,

что при

/^< 0,2% трудно

выделить

максимум

частотного распределения фаз. Этот результат находится в полном согласии с рис. 5, из которого следует, что при

г1^0,2% большой вклад должна вносить суточная вариа­

ция статистического характера с неопределенной фазой. Най­ денное выше частотное распределение, первой гармоники

W(rv а), дает возможность определить математическое ожи­ дание амплитуды первой гармоники, как случайной величи­

ны [25, 26], т. е.

 

 

В \

r1W(rl,<y)dr1

 

со

 

 

В \

W (rv a) d,\

 

I'

 

 

О

 

 

3.1 г*

 

/фехр

drx

 

с т -

У

___

 

 

3 '

2 ~ 2 '

3.1 rj

 

Hexр

dr±

 

а 2

 

 

Посколыку вторая гармоника также определяется из тех же данных xv х2, • • -х2р по формулам (1.3.2) распределение, по­ лученное для первой гармоники, справедливо и для второй

гармоники. Таким

образом, при вычислении

амплитуд

пер­

вой

и второй гармоник ожидается

значение

амплитуды за

счет

флуктуации

статистического

характера

а

(где

rh<----^

к = 1, 2).

Согласно нашим вычислениям (§ 2) для двухчасовых данных нейтронной компоненты для станции Клаймакс и

22

Норпкура — о«0,28% 4-0,25% и ожидаемые амплитуды пер­

ст вой и второй гармоник суточных колебаний будут — «0,15% .

§ 4 К вопросу об определении полусуточной вариации космических лучей

При зондировании межпланетных условий солнечно-су­ точными вариациями коом,нчеС'КИХ лучей во многих случаях недостаточно ограничиваться только лишь первыми гармони­ ками Фурье-аналпза. Привлечение второй гармоники не толь­ ко обогащает информацию, но и отвечает на такие вопросы, на которые невозможно ответить исчерпывающе, пользуясь только первыми гармониками. Существуют различные приемы определения второй гармоники. Естественно, что наиболее полную информацию можно получить при сферическом ана­ лизе [4, 5] путем привлечения мировой сети станций, про­ изводя глобальную съемку космических лучей. Но бывают и такие случаи, когда не требуется такой строгий подходи та­ кие громоздкие вычисления, поэтому рассматриваются поч­ ти без ущерба различные модификации второй гармоники солнечно-суточных вариаций. Одним из таких приближений является разделение полусуточной гармоники на дневную и ночную части [27]. При вычислении полусуточной вариа­ ции космических лучей обычным методом подразумевают идентичность полусуточной вариации в обеих половинах су­ ток. В действительности же они могут отражать разные ус­ ловия в различных половинах космического пространства. Обычным и общепринятым способом выделения полусуточ­ ной вариации является метод гармонического анализа экспе­ риментальных данных по формуле (1.3.1) Бесселя [23].

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

Таким образом, для определнпя воторой гармоники ме­

тодом гармонического

анализа, расчеты надо проводить по

следующим формулам:

 

аг

1

ш

Р/2

№ '

 

i= l

23

(1.4.1)

Полусуточной гармонике отождествляют член из фор­ мулы (1.4.1) .соответствующий К = 2.

При К = 2 формулу

(1.3.1) можно переписать и как сум­

му двух независимых

составляющих

а.

 

Pi2

 

/= 1 ;= р + 1

и если обозначим:

Р

р

1=1

i = i

(1.4.3)

 

< = Р + 1

 

 

то (1.4.2)

можно

переписать:

»

 

 

_/

.

 

 

С1о — С12

“р С1‘>

 

 

Ь2 = ^2 "Ь К'

Определив

а2, К,

аг2 и Ь"2 мы за каждые сутки получаем две

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

Таким образом, для определения второй гармоники ме­

тодом простого гармонического анализа, надо проводить рас­ четы по следующим формулам:

а

i ,

24

2p

a,

-Its

,

I cos ----i

 

P

 

 

 

(1.4.5).'

r,sin

4ts

,

s i n ----i

 

P

 

 

4ts

} al+bf,a-

ь:

sin ---- i, где i\

 

P

 

a2 и o.j, b'2 и b"2—коэффициенты Фурье, соответствующие ссответственно первой (а[ и £д) и второй (а'2 и К) половинам сутокПользуясь формулами (1.4.5) можно определять вторые гармо­ ники для любой половины суток, например, для дневной и ноч­ ной стороны и др.

§ 5. Барометрический эффект нейтронной компоненты космических лучей

При исследованиях модуляционных эффектов космичес­ ких лучей, эффектов, вызванных явлениями в космическом' пространстве, если наблюдения ведутся прибором изнут­ ри атмосферы Земли, требуется освобождение от влияние

атмосферного

происхождения

и магнитного

поля Земл№

[28,

29,

30].

 

 

 

 

Поскольку во многих случаях многие исследователи ог­

раничиваются

исключением

барометрического

эффекта в

данных наблюдения космических лучей, в особенности в- данных нейтронной компоненты (следует подчеркнуть, что эта компонента имеет ничтожно малый температурный эф­ фект [28, 30, 31]), мы несколько подробно остановимся наосвещенип методических вопросов, связанных с барометри­ ческим эффектом нейтронной компоненты. Здесь нет необ­ ходимости описывать процессы, обуславливающие погло-- щенпе и ядерные взаимодействия космических лучей с ато­ мами воздуха. Для этой цели вполне достаточно будет сос­ латься на довольно богатую литературу [28, 30 и др.].

5.1. Ш и р о т н а я з а в и с и м о с т ь

б а р о м е т р и ­

ч е с к о г о к о э ф ф и ц и е н т а

Частицы космических лучей, попадая

в земную атмос­

феру, после взаимодействия с атомами воздуха создают це­

25»

лый каскад вторичных частиц. В огромном потоке вторич­ ных частиц особо выделяется нейтронная компонента.

Интенсивность нейтронной компоненты космических лу­ чен уменьшается при прохождении нейтронов через слой ат­ мосферы Земли, за счет неупругого и упругого взаимодей­ ствия с атомами воздуха. Эффекты иеупругого п упругого взаимодействия нейтронов с атомами воздуха зависят от энергий вторичных нейтронов п теоретически ожидается уве­

личение интегрального пробега для поглощения L' (обрат­

ная величина значения барометрического коэффициента) с увеличением энергий нейтронов. Поскольку станции косми­ ческих лучей мировой сети расположены в различных пунк­

тах

Земного шара

и имеют разные

геомагнитные пороги

жесткости обрезания

R

(в пределах л; 0,5 ч-176Г),

в работе

Л.

И. Дормана, О.

Г.

Рогава, Л. X.

Шаташвнлп

[32] бы­

ла поставлена задача выяснить как влияет изменение мини­ мального порога жесткости на барометрический эффект дан­ ной станции, т. е. определить планетарное распределение ■барометрического коэффициента нейтронной компоненты кос­ мических лучей. Аналогичная работа была выполнена ч [33] по данным 25 станций мировой сети вблизи минимума солнечной активности за период 1962—1963 гг.

Работа [321 должна вызвать интерес и по той простои причине, поскольку эта работа выполнена по данным мини­ мума солнечной активности (1964—1965 гг период МГСС), когда модулирующее воздействие солнечного ветра на кос­ мические лучи было относительно минимальным и интен­ сивность вторичной нейтронной компоненты в фиксирован­ ном пункте наблюдения возникла в основном за счет изме­ нения барометрического давления атмосферы Земли. В ми­ нимуме солнечной активности легко установить, каково вли­ яние чисто геофизических факторов па планетарное распре­ деление интенсивности космических лучей.

Ясно, что прежде всего представляет интерес выяснить планетарное распределение барометрических коэффициентов р нейтронной компоненты космических лучей. В работе [32], эта задача выполнена по данным периода МГСС.

В [32, 34] барометрические коэффициенты вычислялись для интервалов наблюдений, когда происходили заметные изменения атмосферного давления п соответствующие изме­ нения интенсивности нейтронной компоненты по данным ми­

ровой сети за период МГСС (1954-—1965

гг), в соответст­

вии с рекомендациями, изложенными в

[35]. Проводилась

линейная корреляция логарифма интенсивности нейтронной

компоненты со значениями атмосферного давления

h. Бы­

ли использованы данные нейтронной компоненты 23

стан-

•26

шпг. 1. Туле,

2. Инувик, 3.

Оулу, 4. Колледж,

5. Апатиты,

6. Гусбей,

7.

Дпп-Рпвер, 8.

Оттава, 9. г. Вашингтон,

10.

Киль, 11.

Москва, 12. Халле,

13. Клаймакс,

14. Иркутск,

15. Мюнхен,

16. Цугспптце,

17.

Тбилиси, 18. Алма-Ата,

19.

Уанкайо, 20. Ахмадабад, 21. Лаэ, 22. Кодайканал, 23. Кордо­

ва.

На рис. 7 представлена кривая шпротной зависимости (по 17 пунктам наблюдений) барометрических коэффициен­ тов усредненных по группам наземных станций с близкими порогами обрезания по данным 1964—1965 гг. Детальный анализ показывает, что на долготах 50°±20°Е Апатиты,

.Москва, Тбилиси, Алма-Ата, Ахмадабад получаются несколь­ ко заниженные значения барометрических коэффициентов. Несколько аномальное значение барометрического коэффици­ ента получается также для южно-американской станппн Кордова (жесткость обрезания 11,6 Бв) (3=(—0,69± 0,01) %/мб. На рис. 7 значение барометрического коэффициента для станции Кордова выделено звездочкой и при выводе

кривой широтной зависимости не

 

предусматривалось.

Возможно

~ Л ./т$

завышенное значение барометри­

ческого коэффициента для Кор­

 

дова связано с несколько

ано­

 

мальной широтной зависимостью

 

интенсивности космических

лу­

 

чей на долготе Бразилии-Арген­

 

тины (36].

 

 

 

Быть может занижения зна­

M 6V)

чения барометрических коэффи­

 

циентов нейтронной

компоненты-

рис. 7. Широтная зависимость баро-

космпческих лучей связаны С не-^метрических коэффициентов нейтрон­

ной западно-восточной аснммет-

и0" компоненты космических лучей,

рней магнитного поля

Земли, на

°сРВДне..ных для групп станций с

которсе указывали

Сарабай и

близкими порогами обрезания, по даи-

HhIV 19з4_,дз5 гг.

Пай [37]. Эта асимметрия, если

 

она существуют, должна давать

несколько иные, чем даваемые

в [32, 38] жесткости

обрезания для космических лучей п соот­

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

Не исключена возможность, что указанные отклонения барометрических коэффициентов обусловлены нестандарт-

27

ностыо аппаратуры, которая использовалась в период МГСС.

на мировой сети.

показывает, что в

целом (за

исключени­

Исследование

ем этих аномальных случаев)

для наземных станций в сред-

Д|3

г

 

в то время,

 

 

 

Л|3

нем—-- — 1%/Ьв,

как для высокогорных-^- «

Р

 

 

 

МГСС. Сравнение

приведенных

нами

~ 0,5%/Бв в период

барометрических

коэффициентов с результатами

[33]

пока­

зывает,

что [32]

во

многих случаях

значения |3

для

эпохи

вблизи

минимума

 

солнечной

активности

несколько

выше,

чем данные широтной зависимости, которые приведены в на­ шей работе.

Таким образом, как в [33], так и в [32] обнаружено уменьшение барометрического эффекта с увеличением энер­ гий вторичных нейтронов, что хорошо согласуется с теорией.

5. 11- л е т н я я з а в и с и м о с т ь б а р о м е т р и ч е с к о г о

ко э ф ф и ц п е и т а

В§ 52 наказано, что барометрический эффект нейтрон­ ной компоненты космических лучей зависит от жесткости частиц. Естественно, напрашивается вопрос, наблюдается ли

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

В ряде работ [33, 36, 39, 40, 41] была обнаружена тен­ денция к 11-летнему изменению барометрического коэф­ фициента интенсивности нейтронной компоненты космичес­ ких лучей, однако в связи с тем, что на наш взгляд не су­ ществует единой методики для определения барометричес­ кого коэффициента, результаты, полученные разными авто­ рами, естественно могут отличаться друг от друга. Как по­ казывает анализ, эти результаты могут различаться лишь

в 1,5 2 раза. Поэтому исследование в этом направлении

мы считаем естественным и оправданным. Есть основания полагать, что при вычислении барометрических коэффициен­ тов целесообразно там, где это возможно пользоваться фор­ мулой множественной корреляции [31].

Формула множественной корреляции [31] в нашем слу­ чае представляется как

ГMnN

NM'hM

■ЧпЛ'

Р =

 

(1.5.1)

1 -г ? ,Ш

 

28

где /i — атмосферное давление, N и М — интенсивности нейтронной компоненты космических лучей, соответственно неисправленный и исправленный на барометрический эф­ фект, а1пЛ, и ак — то же, что и среднеквадратичные откло­ нения параметров 1пЛ, н h.

rh\aN— коэффициент корреляции между изменениями ат­ мосферного давления h и логарифма интенсивности нейт­ ронной компоненты космических лучей N (неисправленной на барометрический эффект), rNM Y— коэффициент корреляции между изменениями интенсивности нейтронной компоненты

космических лучей (не исправленной

N и исправленной М

на барометрический эффект),

гш

— коэффициент корреля­

ции между

изменениями атмосферного давления h и N

ин­

тенсивности

 

нейтронной

компоненты

космических

лучей

(исправленной на барометрический эффект).

 

 

 

 

Для двух станций А и В с близкими порогами обреза­

ния первичных космических

лучей

(А—Черчилль и В—Отта­

ва)

вычислены

серин

барометрических

коэффициентов

 

(где

/ =

1,2,3...

участки определения барометрических коэффициен­

тов в разные периоды

солнечной активности, i =

А.В-названпя

станции и К =

1,2,3...

номера последовательных

приближений

для данного участка у) следующим образом:

 

 

 

 

Для

периода i был

вычислен

барометрический

коэффи­

циент |3£и

(по данным станции Черчилль) по формуле

множественной корреляции

[31],

используя в

качестве

па­

раметра

М

(характеризующей вариации первичного спект­

ра)

исправленные на

средний

барометрический коэффици­

ент данные станции В. После этого первичные данные стан­

ции А были исправлены на давление с

барометрическим

коэффициентом рАл и уже исправленные данные были

исполь­

зованы в качестве

параметра /И при

расчете

для

станции

В (Оттава). Далее перзичные данные

станции В были

исправ­

лены на давление

с барометрическим

коэффициентом

и

исправленные данные были использованы в качестве параметра М при расчете {3^<2 и т. д.

Таким образом, были вычислены серии $Bik и §>Bk и конеч­ ные значения определены при условии, что

IP W -

0,005%/М б ,

(1-5.2)\

29

I ^ 0,005 % /мб.

Очевидно, условие (1.5.2) является приемлемым при расчетах барометрического коэффициента по данным нейт­

ронных мониторов

типа МГГ (Симпсона).

Результаты вычисления барометрических коэффициен­

тов приведены

на

рис. 8а (Черчилль) и 86 (Оттава). На

этих рисунках

по осям ординат отложены значения баромет-

Рис. 8 а, б. Среднегодовые барометрические коэффициенты ин­ тенсивности нейтронной компоненты космических лучей в зависи­ мости от эпохи 11—летного цикла солнечной активности по данным

станций Черчнль (а), Оттава (в)

рнческого коэффициента, а по осям абсцисс — годы. Из этих рисунков видно, что по данным обеих станций Черчилль и Оттава наблюдается изменение барометрического коэффи­

циента ~ (6-р-7)/%мб в зависимости от 11-летней циклич­ ности солнечной активности.

5. 3 Учет роли р - ме з онов п л о к а л ь н ы х

п р о д у к т о в з а г р я з н е н и я

Выше (в § 5.1) было показано, что существует широт­

ный эффект барометрического

коэффициента нейтронной

компоненты космических лучей ~ (1

1,5)%/Бв и 11-летнее изме­

нение—1 %/ на год. Можно утверждать, что найденная по экспери­ ментальным данным яркая зависимость барометрического эффекта от порога жесткости обрезания однозначно дока­ зана, в то время, как вопрос малого эффекта за счет П-лет- него цикла (~1 % на год) все еще требует дополнительных дока­

зательств. Действительно, барометрический коэффициент также может зависеть от средней интенсивности нейтронов

30

Соседние файлы в папке книги из ГПНТБ