солнечный пятно сейсмический земля
Таблица 1 Параметр сходства категории магнитуда в четырех моделях
|
Категория |
Число случаев |
M160(L,R)) |
M160(RA,L,R) |
M170(RA,L,P)) |
M170(RES,L,R) |
|
|
A40-Mb=4,0 |
1362 |
-18,593 |
-12,149 |
-14,462 |
-20,037 |
|
|
A41-Mb=4,1 |
1580 |
-8,208 |
-6,109 |
-3,077 |
-6,813 |
|
|
A42-Mb=4,2 |
1796 |
-1,797 |
-0,407 |
2,698 |
0,637 |
|
|
A43-Mb=4,3 |
2224 |
4,889 |
5,069 |
8,640 |
7,037 |
|
|
A44-Mb=4,4 |
2744 |
-4,445 |
-3,647 |
-1,136 |
-4,510 |
|
|
A45-Mb=4,5 |
3358 |
-0,364 |
0,486 |
3,546 |
1,890 |
|
|
A46-Mb=4,6 |
4119 |
8,102 |
7,114 |
2,359 |
8,161 |
|
|
A47-Mb=4,7 |
4768 |
11,303 |
7,139 |
1,431 |
10,887 |
|
|
A48-Mb=4,8 |
4954 |
13,381 |
8,852 |
2,838 |
12,764 |
|
|
A49-Mb=4,9 |
5008 |
17,503 |
10,799 |
6,597 |
15,614 |
|
|
A50-Mb=5 |
4904 |
26,141 |
19,409 |
11,462 |
22,451 |
|
|
A51-Mb=5,1 |
4582 |
29,802 |
20,595 |
14,006 |
26,773 |
|
|
A52-Mb=5,2 |
4134 |
25,101 |
18,918 |
13,014 |
23,767 |
|
|
A53-Mb=5,3 |
3563 |
30,307 |
22,368 |
17,016 |
29,035 |
|
|
A54-Mb=5,4 |
3010 |
28,031 |
20,712 |
13,964 |
24,969 |
|
|
A55-Mb=5,5 |
2367 |
24,238 |
17,492 |
10,680 |
22,462 |
|
|
A56-Mb=5,6 |
1940 |
32,491 |
23,640 |
21,515 |
31,337 |
|
|
A57-Mb=5,7 |
1460 |
31,689 |
23,066 |
22,448 |
31,225 |
|
|
A58-Mb=5,8 |
1179 |
28,16 |
22,760 |
20,935 |
28,725 |
|
|
A59-Mb=5,9 |
864 |
24,738 |
19,841 |
18,139 |
25,450 |
|
|
A60-Mb=6,0 |
656 |
23,495 |
20,666 |
17,798 |
25,074 |
|
|
A61-Mb=6,1 |
453 |
29,394 |
25,123 |
22,496 |
29,190 |
|
|
A62-Mb=6,2 |
319 |
23,934 |
20,695 |
21,355 |
25,934 |
|
|
A63-Mb=6,3 |
202 |
35,651 |
28,722 |
28,901 |
36,412 |
|
|
A64-Mb=6,4 |
137 |
42,81 |
37,909 |
39,716 |
45,038 |
|
|
A65-Mb=6,5 |
87 |
55,802 |
47,917 |
51,895 |
62,048 |
|
|
A66-Mb=6,6-7,0 |
68 |
63,169 |
53,843 |
56,396 |
68,781 |
|
|
B40-Mb=4,0 |
446 |
28,361 |
18,736 |
23,423 |
30,855 |
|
|
B41-Mb=4,1 |
660 |
18,857 |
12,704 |
14,817 |
21,036 |
|
|
B42-Mb=4,2 |
835 |
9,046 |
6,722 |
6,265 |
9,861 |
|
|
B43-Mb=4,3 |
955 |
-6,859 |
-2,669 |
-8,075 |
-8,597 |
|
|
B44-Mb=4,4 |
1099 |
-14,986 |
-8,869 |
-10,093 |
-13,509 |
|
|
B45-Mb=4,5 |
1223 |
-1,494 |
-0,076 |
2,080 |
0,992 |
|
|
B46-Mb=4,6 |
1455 |
-1,678 |
-1,113 |
1,722 |
-0,022 |
|
|
B47-Mb=4,7 |
1612 |
-1,317 |
0,130 |
1,591 |
-1,349 |
|
|
B48-Mb=4,8 |
1817 |
3,926 |
4,190 |
5,574 |
3,183 |
|
|
B49-Mb=4,9 |
1636 |
3,783 |
3,855 |
2,818 |
5,508 |
|
|
B50-Mb=5 |
1428 |
15,492 |
11,992 |
11,574 |
16,131 |
|
|
B51-Mb=5,1 |
1206 |
13,093 |
10,956 |
9,329 |
12,782 |
|
|
B52-Mb=5,2 |
936 |
22,413 |
19,267 |
16,803 |
22,561 |
|
|
B53-Mb=5,3 |
617 |
24,482 |
20,240 |
17,872 |
22,979 |
|
|
B54-Mb=5,4 |
422 |
21,194 |
18,464 |
15,913 |
18,494 |
|
|
B55-Mb=5,5 |
261 |
30,487 |
23,114 |
26,050 |
31,170 |
|
|
B56-Mb=5,6 |
180 |
33,408 |
28,321 |
30,616 |
37,841 |
|
|
B57-Mb=5,7 |
93 |
56,994 |
47,736 |
52,339 |
60,019 |
|
|
B58-Mb=5,8 |
73 |
61,514 |
52,517 |
55,427 |
68,742 |
|
|
B59-Mb=5,9-6,4 |
69 |
62,605 |
54,655 |
59,161 |
68,168 |
|
|
C40-Mb=4,0 |
130 |
53,018 |
42,336 |
51,118 |
60,331 |
|
|
C41-Mb=4,1 |
259 |
35,327 |
27,456 |
34,395 |
40,927 |
|
|
C42-Mb=4,2 |
362 |
22,208 |
14,780 |
20,365 |
24,192 |
|
|
C43-Mb=4,3 |
436 |
21,621 |
15,521 |
18,421 |
24,474 |
|
|
C44-Mb=4,4 |
482 |
19,875 |
14,973 |
11,957 |
19,482 |
|
|
C45-Mb=4,5 |
468 |
11,685 |
9,127 |
8,575 |
10,782 |
|
|
C46-Mb=4,6 |
515 |
2,741 |
3,117 |
-0,391 |
1,491 |
|
|
C47-Mb=4,7 |
501 |
5,517 |
5,687 |
5,406 |
6,989 |
|
|
C48-Mb=4,8 |
450 |
10,42 |
9,536 |
8,836 |
12,503 |
|
|
C49-Mb=4,9 |
447 |
15,73 |
13,178 |
14,786 |
17,397 |
|
|
C50-Mb=5 |
356 |
20,558 |
16,937 |
15,841 |
22,289 |
|
|
C51-Mb=5,1 |
293 |
22,202 |
18,100 |
17,483 |
22,638 |
|
|
C52-Mb=5,2 |
166 |
35,702 |
30,512 |
31,979 |
38,506 |
|
|
C53-Mb=5,3-6,0 |
105 |
48,043 |
37,886 |
44,190 |
52,017 |
|
|
Среднее значение |
1374 |
20,995 |
17,063 |
16,874 |
22,052 |
Рис. 1
Ранее было показано /13-16/, что матрица информативности (2) может быть использована для выявления и визуализации когнитивных функциональных зависимостей в фрагментированных и зашумленных данных большой размерности в задачах прогнозирования климата и сейсмической активности. Кратко поясним суть этого метода. Матрица информативностей (2) рассчитывается на основе системной теории информации /19/ непосредственно на основе эмпирических данных и представляет собой таблицу, в которой столбцы соответствуют обобщенным образам классов, т.е. будущим состояниям моделируемой системы, строки - значениям факторов, влияющих на эту систему, а на пересечениях строк и столбцов находится количество информации, которое содержится в факте действия значения фактора, соответствующего строке, на переход системы в состояние, соответствующее столбцу. Максимальное количество информации, которое может быть в значении фактора, определяется числом будущих состояний моделируемой системы. Модуль количества информации отражает силу влияния значения фактора, а знак - направление этого влияния, т.е. то, способствует он или препятствует наступлению данного состояния. Если последовательности классов и значений факторов образуют порядковые шкалы или шкалы отношений, т.е. соответственно, на них определены отношения «больше-меньше» или, кроме того, единица измерения, начало отсчета и арифметические операции, то матрица информативностей допускает наглядную графическую визуализацию, традиционного для функций типа, когда значения факторов рассматриваются в качестве значений аргумента, а классы, о наступлении которых в этих значениях факторов содержится максимальное количество информации - в качестве значений функции. Другие классы, менее обусловленные данным значением фактора, а также те, наступлению которых это значение препятствует в большей или меньшей степени, также могут отображаться соответствующими цветами, и это также может представлять интерес, т.к. позволяет задействовать мощные способности человека к анализу изображений. Когнитивные функции, представляемые в форме матрицы информативностей, соответствуют очень общему виду функциональной зависимости: многозначной функции многих аргументов, т.к. каждое значение фактора влияет на все состояния моделируемого объекта, и каждое его состояние обусловлено всеми значениями факторов /20/.
На рис. 2 представлены фрагменты матрицы информативностей моделей M170(RA,L,Р) и M170(В,L). Визуализация этих фрагментов осуществлялась на основе программы InfVisual v. 2.2 (авторы: Д. Бандык и Е.В. Луценко). Из представленных на рис. 2 данных следует, что существует влияние солнечных пятен на сейсмические события, но это влияние не подобно влиянию индукции магнитного поля. Отметим, что в работе /21/ было установлено подобие влияния долготы Сатурна и вариаций магнитного поля на сейсмические события.
Рис. 2 Зависимость частоты категорий однократных (А), двукратных (В) и трехкратных (С) сейсмических событий с магнитудой от параметра RADAILY и вариаций вертикальной компоненты индукции магнитного поля Земли на трех станциях в 1963-2006 гг по данным /10/ в модели М170
На рис. 3. представлены фрагменты матрицы информативностей моделей M160(L,Р365) и M160(L,R). Из представленных на рис. 3 данных следует, что существует влияние солнечных пятен на сейсмические события через посредство комплексов Р3 Урана и Нептуна, и что это влияние подобно влиянию индукции магнитного поля и долготы Сатруна.
Рис. 3 Зависимость частоты категорий однократных (А), двукратных (В) и трехкратных (С) сейсмических событий с магнитудой от параметра P3 вычисленного для Урана и Нептуна при Q0i =365, долготы Сатурна и вариаций вертикальной компоненты индукции магнитного поля Земли на станции SOD в 1963-2006 гг по данным /10/ в модели М160
Можно предположить, что основную роль во влиянии на сейсмические события играет возбуждение магнитосферы небесных тел солнечным ветром. Однако, как было установлено, влияние Юпитера, обладающего самой большой по объему магнитосферой среди планет, не подобно влиянию Урана и Нептуна. На рис. 4. представлены фрагменты матрицы информативностей моделей M160(L,Р365) и M160(L,R). Из представленных на рис. 4 данных следует, что существует влияние солнечных пятен на сейсмические события через посредство комплекса Р3 Юпитера, но оно не подобно влиянию магнтного поля земли, Уран, Нептуна и Сатурна.
Рис. 4 Зависимость частоты категорий однократных (А), двукратных (В) и трехкратных (С) сейсмических событий с магнитудой от параметра P3 вычисленного для Юпитера при Q=1; Q0 =68 и Q0 =365, и долготы Юпитера
В рамках семантических информационных моделей можно решить вопрос о влиянии взаимных аспектов небесных тел, которые определяются как разность долготы двух небесных тел. Исследовалось влияние на сейсмические события разности долготы Солнца и четырех небесных тел - Юпитера, Сатурна, Урана и Нептуна. В результате установлено, что заметного влияния аспектов Солнца с Юпитером, Сатурном, Ураном и Нептуном на сесмические события нет, хотя система «Эйдос-астра» позволяет отфильтровать несколько пиков для каждого из аспектов - см. рис. 5.
Исследовалось влияния разности долготы четырех небесных тел - Юпитера, Сатурна, Урана и Нептуна между собой. Наиболее значимы оказались взаимные аспекты (в порядки значимости):
1) Уран/Нептун - влияние на общее число землетрясений, на суммарную магнитуду, на среднюю энергию и на средний объем;
2) Сатурн/Нептун - влияние на общее число землетрясений, на суммарную магнитуду, на среднюю магнитуду;
3) Сатурн/Уран - влияние на среднюю магнитуду;
4) Юпитер/Сатурн - влияние на общее число землетрясений и на суммарную магнитуду.
У перечисленных аспектов ярко выраженное влияние, похожее на влияние магнитного поля Земли - см. рис. 5 и 2.
Рис. 5 Зависимость частоты категорий однократных (А), двукратных (В) и трехкратных (С) сейсмических событий с магнитудой от взаимных аспектов Солнца, Юпитера, Сатурна, Урана и Нептуна
На данных для аспектов Юпитера на рис. 5 хорошо просматриваются всплески, соответствующие некоторым аспектам Юпитер/Уран и Юпитер/Нептун. Полосы, обнаруженные на данных для долготы Юпитера на рис. 4, видимо, соответствуют этим аспектам.
Влияние солнечных пятен на магнитное поле Земли
Как известно, число солнечных пятен связано с изменением полярности солнечного магнитного диполя /1, 6/. Поэтому, можно ожидать, что есть прямая связь между параметром RADAILY и вариациями индукции магнитного поля земли. На рис. 6. приведены коэффициенты корреляции параметра RADAILY с вертикальной компонентой индукции магнитного поля на протяжении 16032 дней по измерениям на 23 станциях, перечисленных выше. Средний коэффициент корреляции по всем станциям составляет около -0,00409, что, на первый взгляд, свидетельствует об отсутствии связи между числом пятен на Солнце и вариациями магнитного поля Земли. С другой стороны, данные /2-5/ и другие указывают на наличие такой связи. Для объяснения этого парадокса была разработана специальная методика вычисления геомагнитных индексов, которые изменяются подобно параметру RADAILY или числу Вольфа /21/. Такие индексы составляются на основе часового, трехчасового, суточного, месячного или годового изменения энергии магнитного поля. В частности, в монографии /21/ были построены месячные и годичные средние комбинации вариаций магнитного поля, которые имеют коэфиициент корреляции с числом пятен на Солнце 0,654 и 0,884 соответственно.
Начиная с 1 января 1964 года IAGA было рекомендовано не использовать вертикальную компоненту индукции магнитного поля при вычислении индексов геомагнитной активности, поскольку влияние местных аномалий может быть весьма велико для вертикальной компоненты поля. Это представляется несколько странным, так как в средних и высоких широтах вертикальная компонента поля превосходит горизонтальную составляющую.
Рис. 6
Можно предположить, что некоторые выборки их рядов данных RADAILY будут иметь более высокий коэффициент корреляции с вертикальной компонентой индукции магнитного поля, нежели весь ряд. Таким образом, были образованы новые ряды данных, удовлетворяющие условию RADAILY>C, а также годичный ряд и ряд 74 дня - таблица 2. Укороченные ряды действительно имеют более высокий коэффициент корреляции на некоторых станциях, однако средний по всем 23 станциям коэффициент корреляции имеет максимум 0,06 для годичного ряда, что значительно меньше величины 0,884, полученной авторами /21/ для корреляции среднегодичных величин геомагнитного индекса и числа солнечных пятен.
Полученные результаты означают, что влияние солнечных пятен на магнитное поле Земли является опосредованным. Согласно /21-22/ в механизм этого влияния включены такие сложные системы, как земная магнитосфера и ионосфера, а также солнечная корона.
С другой стороны, в работах /13, 16/ и других было показано, что на магнитное поле Земли существенное влияние оказывают небесные тела Солнечной системы. Особая роль в этом влиянии принадлежит Урану и Нептуну. Гипотеза опосредованного влияния солнечных пятен на земную магнитосферу находит свое подтверждение в рамках модели /18/, в которой параметры (5) вычисляются с учетом влияния числа солнечных пятен. Физический смысл выражений (5) довольно очевиден: опосредованное влияние солнечных пятен включает медленный интегральный процесс и быстрый процесс, который распространяется со скоростью света. При этом предполагается, что солнечные пятна индуцируют дипольное излучение путем воздействия на магнитосферу небесного тела.
Таблица 2 Изменение коэффициента корреляции параметра RADAILY с вариацией вертикальной компоненты индукции магнитного поля на шести станциях и среднее значение по 23 станциям для различной выборки данных.
|
RADAILY |
BzGNA |
BzGUA |
BzIRT |
BzTHL |
BzDRV |
BzLRV |
AVERAGE23 |
|
|
N=16032 |
-0,172 |
-0,186 |
-0,145 |
0,121 |
0,143 |
0,152 |
-0,00408696 |
|
|
N(RA>9)=13890 |
-0,187 |
-0,209 |
-0,152 |
0,119 |
0,134 |
0,124 |
-0,01105652 |
|
|
N(RA>50)=8459 |
-0,238 |
-0,257 |
-0,189 |
0,133 |
0,0904 |
0,129 |
-0,02345044 |
|
|
N=5162(RA>90) |
-0,227 |
-0,248 |
-0,211 |
0,127 |
0,00867 |
0,0934 |
-0,03587739 |
|
|
N=3749(RA>110) |
-0,249 |
-0,242 |
-0,213 |
0,113 |
0,00505 |
0,096 |
-0,0374187 |
|
|
N=365 |
-0,0578 |
0,146 |
0,295 |
-0,126 |
-0,129 |
-0,212 |
0,06005909 |
|
|
N=74 |
-0,152 |
-0,00877 |
-0,146 |
0,172 |
0,521 |
0,227 |
-0,02464864 |
На рис. 7 представлены абсолютные значения коэффициентов корреляции комплексов (5) в случае Урана и Нептуна с вертикальной компонентой индукции магнитного поля Земли. Отметим, что в отличие от параметра RADAILY - см. рис. 6, комплексы (5) имеют высокий, вплоть до 0,98, коэффициент корреляции с вертикальной компонентой индукции. Интегральный параметр, фигурирующий в выражениях (5), принят постоянным и равным H0i =365. Это значение не было оптимизировано для указанных небесных тел.