Модель линейной регрессии для параметров индукции магнитного поля имеет вид:
(4)
Здесь S - станция наблюдения. В таблице 4 приведены географические координаты магнитных обсерваторий, данные которых были использованы в настоящем исследовании. Для каждой станции из таблицы 4 вычислялись коэффициенты линейной регрессии, согласно (4), по которым восстанавливались значения зарядов небесных тел.
Таблица 4 Коды и географические координаты магнитных обсерваторий
|
Station |
Dumont d'Urville |
Hermanus |
Gnangara |
Fuquene |
Guam |
Alibag |
Honolulu |
Kakioka |
|
|
IAGA Code |
DRV |
HER |
GNA |
FUQ |
GUA |
ABG |
HON |
KAK |
|
|
Lat |
-66,667 |
-34,425 |
-31,8 |
5,47 |
13,59 |
18,638 |
21,32 |
36,232 |
|
|
Long |
140,009 |
19,227 |
118 |
286,265 |
144,89 |
72,874 |
204 |
140,188 |
|
|
Station |
Chambon-la-Foret |
Fursten-feldbruck |
Dourbes |
Hartland |
Niemegk |
Patrony |
Meanook |
Lerwick |
|
|
IAGA Code |
CLF |
FUR |
DOU |
HAD |
NGK |
IRT |
MEA |
LER |
|
|
Lat |
48,017 |
48,17 |
50,1 |
51 |
52,072 |
52,167 |
54,615 |
60,133 |
|
|
Long |
2,269 |
11,3 |
4,8 |
355,519 |
12,675 |
104,47 |
246,655 |
358,817 |
|
|
Station |
Baker Lake |
Resolute Bay |
Qaanaaq (Thule) |
Sodankyla |
Alma Ata |
Addis Ababa |
Antananarivo |
San Juan |
|
|
IAGA Code |
BLC |
RES |
THL |
SOD |
AAA |
AAE |
TAN |
SJG |
|
|
Lat |
64,333 |
74,69 |
77,483 |
67,367 |
43,25 |
9,03 |
-18,917 |
18,117 |
|
|
Long |
263,969 |
265,105 |
290,833 |
26,633 |
76,917 |
38,765 |
47,552 |
293,85 |
|
|
Station |
MBour |
Faraday Islands |
Apia |
Bangui |
Tucson |
Byrd Station 2 |
Meanook |
Cape Wellen (Uelen) |
|
|
IAGA Code |
MBO |
AIA |
API |
BNG |
TUC |
BYR |
MEA |
CWE |
|
|
Lat |
14,384 |
-65,245 |
-13,807 |
4,333 |
32,17 |
66,163 |
54,615 |
66,163 |
|
|
Long |
343,033 |
295,742 |
188,225 |
18,867 |
249,27 |
240,483 |
246,653 |
190,165 |
|
|
Station |
Nurmijarvi |
Valentia |
Wingst |
College |
Kanozan |
Niemegk |
Memambetsu |
Leirvogur |
|
|
IAGA Code |
NUR |
VAL |
WNG |
CMO |
KNZ |
NGK |
MMB |
LRV |
|
|
Lat |
60,508 |
51,933 |
53,743 |
64,87 |
35,256 |
52,072 |
43,91 |
64,183 |
|
|
Long |
24,655 |
349,75 |
9,073 |
212,14 |
139,956 |
12,675 |
144,189 |
338,5 |
Оценка зарядов получается неоднозначной - см. таблицу 5, поэтому для каждого небесного тела можно выделить гипотетическую величину, используя среднее значение по всем станциям. В свою очередь, из средних значений можно выбрать наиболее значимое, в результате были получены следующие порядки величин (Кулон): Солнце - 1016 ; Луна - 1010 ; Меркурий - 1013; Венера - 1013; Марс - 1013; Юпитер - 1014; Сатурн - 1015; Уран - 1017; Нептун - 1018.
Эти данные показывают, что небесные тела должны обладать гигантскими зарядами, чтобы обеспечить наблюдаемую вариацию индукции магнитного поля Земли. В этой связи заметим, что Земля обладает отрицательным зарядом, распределение которого создает в приземном слое атмосферы электрическое поле напряженностью около 100 В/м. По теореме Гаусса этому полю соответствует суммарный отрицательный заряд около 451615 кулон, что на 5-10 порядков меньше, чем приведенные оценки.
Проблема большого заряда легко разрешается, если предположить, что носителями электричества выступают не ионы и электроны, а крупные образования - кластеры, кристаллы, капли и т.п. Максимальный заряд, который может удержать небесное тело, определяется из баланса силы гравитационного притяжения и силы Кулона отталкивания одноименных зарядов в виде
(5)
Здесь электрическая и гравитационная постоянная соответственно, - общее число носителей электричества. Отсюда находим заряд небесного тела
(6)
Таблица 5 Коэффициенты моделей линейной регрессии вариаций вертикальной компоненты индукции магнитного поля и соответствующий им средний заряд
|
STATION |
KAK |
MMB |
CLF |
FUR |
HAD |
Q, Кулон |
|
|
R= |
0,998 |
0,998 |
0,999 |
0,999 |
0,998 |
||
|
Constant |
3,797575 |
3,918072 |
6,015683 |
5,439575 |
5,836794 |
||
|
SUN1 |
2,156975 |
2,320956 |
3,876403 |
3,479264 |
3,718508 |
6,41E+16 |
|
|
SUN2 |
0,008794 |
0,009831 |
0,015294 |
0,011703 |
0,009992 |
2,29E+14 |
|
|
SUN3 |
0,003658 |
0,005581 |
0,002942 |
0,003197 |
0,002389 |
7,32E+13 |
|
|
MOON1 |
0,000193 |
0,000131 |
-8,67E-05 |
-9,19E-05 |
-0,00003 |
4,75E+11 |
|
|
MOON2 |
2,81E-05 |
1,97E-05 |
-1,03E-05 |
-1,16E-05 |
-2,22E-06 |
9,73E+10 |
|
|
MOON3 |
1,56E-05 |
1,36E-05 |
-3,11E-07 |
1,07E-06 |
-2,62E-07 |
1,22E+11 |
|
|
MERCURY1 |
-0,00398 |
-0,00472 |
-0,00237 |
-0,00278 |
-0,00252 |
-6,74E+13 |
|
|
MERCURY2 |
-0,00241 |
-0,00166 |
-0,00459 |
-0,00262 |
-0,003 |
-5,88E+13 |
|
|
MERCURY3 |
-0,01361 |
-0,0141 |
-0,01045 |
-0,00924 |
-0,01026 |
-2,38E+14 |
|
|
VENUS1 |
-0,00167 |
-0,00179 |
-0,0012 |
-0,00089 |
-0,00195 |
-3,09E+13 |
|
|
VENUS2 |
0,00435 |
0,004903 |
0,00205 |
0,002908 |
0,004819 |
7,84E+13 |
|
|
VENUS3 |
0,003433 |
0,002081 |
0,003881 |
0,003125 |
0,002356 |
6,13E+13 |
|
|
MARS1 |
-0,00529 |
-0,00523 |
-0,00371 |
-0,00239 |
-0,00285 |
-8,02E+13 |
|
|
MARS2 |
0,001369 |
0,002572 |
0,004842 |
0,004797 |
0,00016 |
5,66E+13 |
|
|
MARS3 |
0,005525 |
0,003436 |
-0,00411 |
-0,00534 |
-0,00086 |
-5,54E+12 |
|
|
JUPITER1 |
0,163169 |
0,131408 |
0,042344 |
0,039258 |
0,034583 |
1,69E+15 |
|
|
JUPITER2 |
-0,01562 |
-0,01259 |
0,05185 |
0,033114 |
0,033647 |
3,72E+14 |
|
|
JUPITER3 |
-0,06074 |
-0,05493 |
0,017375 |
0,004247 |
0,005156 |
-3,66E+14 |
|
|
SATURN1 |
0,900206 |
0,678317 |
0,538817 |
0,294714 |
0,760256 |
1,31E+16 |
|
|
SATURN2 |
0,906956 |
0,900878 |
-0,4955 |
-0,28283 |
-0,38613 |
2,65E+15 |
|
|
SATURN3 |
-0,1365 |
4,97E-05 |
0,498336 |
0,272967 |
0,986864 |
6,68E+15 |
|
|
URANUS1 |
-17,214 |
-13,9422 |
-24,1043 |
-22,4837 |
-26,4272 |
-4,29E+17 |
|
|
URANUS2 |
10,36017 |
7,186717 |
5,329775 |
3,175997 |
8,928747 |
1,44E+17 |
|
|
URANUS3 |
16,00536 |
14,50029 |
-9,95766 |
-9,43402 |
-10,3066 |
3,33E+15 |
|
|
NEPTUNE1 |
54,83559 |
49,98008 |
94,94126 |
87,9602 |
96,60819 |
1,58E+18 |
|
|
NEPTUNE2 |
4,4182 |
18,91763 |
17,80636 |
16,79089 |
3,133958 |
2,52E+17 |
|
|
NEPTUNE3 |
60,22368 |
55,4146 |
63,46015 |
65,44635 |
59,93677 |
1,25E+18 |
Если бы Земля состояла только из носителей электричества, тогда, согласно (6), ее максимальный заряд составил бы 5,15*1014 кулон, что по порядку величины совпадает с оценкой заряда Юпитера, приведенной выше. Реально же общая масса носителей электричества нашей планеты, вычисленная по уравнению (6), составляет всего около 4500 тонн, поэтому ее суммарный заряд не столь велик, как у планет-гигантов.
При нарушении условия равновесия (5) слишком большой заряд не сможет удержаться на поверхности небесного тела, а элементарные носители в виде кристаллов или даже крупных глыб будут выброшены на орбиту, где образуют скопление типа колец Сатурна, Урана или Нептуна. Кратковременные всплески излучения, характерные для электрических разрядов, а также необычные образования типа спиц, обнаруженные в кольце В, свидетельствуют о наличии значительного электрического заряда у глыб льда, образующих кольца Сатурна. Если предположить, что величины суммарных зарядов колец Сатурна, Урана и Нептуна согласуются с полученными данными, тогда механизм возмущения магнитного поля Земли можно полностью объяснить на основе модели (1). Кроме того, на основе этой модели можно объяснить суточные колебания напряженности электрического поля в приземном слое, достигающие 30В/м, а также возникновение грозовых облаков.
Излучение Вавилова-Черенкова при сверхсветовом относительном движении небесных тел
В работе /1/ были построены модели линейной регрессии для данных по магнитному полю /3/ в зависимости только от параметров дипольного излучения Урана и Нептуна. Действительно, как следует из полученных данных, вклад Урана и Нептуна в вариации магнитного поля Земли значительно превышает вклады других небесных тел.
Столь сильное влияние Нептуна, возможно, объясняется тем, что в системе отсчета, связанной с поверхностью вращающейся Земли, относительная скорость Нептуна, периодически превышает скорость света - рис. 4. При такой относительной скорости движения необходимо учитывать релятивистские поправки в уравнениях электродинамики и механики, тогда как в современных астрономических вычислениях в пределах Солнечной системы используются уравнения классической небесной механики.
Дадим качественную оценку электродинамических эффектов, связанных с субсветовым и сверхсветовым движением Нептуна, основанную на преобразованиях Лоренца для поля равномерно движущегося заряда и потенциалах Лиенара-Вихерта /6-7/. Предположим, что заряд движется с до световой скоростью. В этом случае с учетом релятивистских эффектов уравнение (1) имеет вид /6/:
(7)
Здесь - угол между направлением движения и радиус-вектором небесного тела. Согласно второму уравнению (7), электрическое поле содержит особенности (стремится к бесконечности) в точках , которые, очевидно, являются источниками возмущений электрического и магнитного поля Земли. Согласно уравнениям (7) в этом случае наведенные электрические и магнитные поля неограниченно возрастают по величине, что на практике эквивалентно электромагнитному импульсу.
В настоящее время нет теории, позволяющей описать движение Нептуна со сверхсветовой скоростью относительно поверхности Земли, чтобы оценить влияние его электромагнитных и гравитационных полей на параметры земной ионосферы, магнитосферы и скрость прецессии земной оси. Так, в учебном пособии /6/ утверждается, что вращающейся системой координат можно пользоваться лишь до расстояний, равных , поскольку для больших расстояний такую систему отсчета невозможно осуществить материальными телами. В случае Земли это составляет 27,5566 астрономических единиц, тогда как орбиты Нептуна и Плутона пролегают несколько дальше этого расстояния (например, на 22.07.2010 расстояние от Земли до указанных небесных тел составило 29,12 а.е. и 30,96 а.е. соответственно). С другой стороны, астрономические наблюдения за удаленными объектами Солнечной системы, звездами и галактиками, которые движутся относительно поверхности Земли со сверхсветовой скоростью, осуществляются во вращающейся евклидовой системе координат, которая простирается сколь угодно далеко. Следовательно, можно использовать вращающуюся евклидову систему координат с учетом возможных релятивистских поправок.
Отметим, что задача о поведении электромагнитных полей при сверхсветовой относительной скорости движения заряда была решена Хевисайдом (1888), де Кудром (1900) и Зоммерфельдом (1904) еще до открытия Эйнштейном релятивистской формы уравнений электродинамики в 1905 г /7-8/. В последнее время наблюдается повышенный интерес к такого рода задачам, что связано с обнаружением релятивистских струй, испускаемых ядрами галактик /9/, а также с различными модификациями эффекта излучения Вавилова-Черенкова /8/. Как известно, электромагнитные потенциалы заряда, движущегося со сверхсветовой скоростью, определены в конусе, вершина которого совпадает с зарядом, а угол раствора определяется из уравнения . Поверхность конуса, ограничивающего решение, является огибающей волн Вавилова-Черенкова. Эти волны можно наблюдать в сплошных средах с большим показателем преломления, типа воды или свинцового стекла, в которых скорость видимого света меньше скорости света в вакууме, при движении заряженных частиц со сверхсветовой скоростью. Излучение Вавилова-Черенкова распространяется пол углом к скорости движения, а сам угол определяется из уравнения /7-8/.
Будем исходить из уравнений потенциалов электромагнитного поля движущегося заряда в форме Лоренца
(8)
В частном случае движения с постоянной скоростью векторный потенциал выражается через скалярный потенциал в виде /6/
Рассмотрим автомодельное решение первого уравнения (8) в однородном случае, когда плотность заряда равна нулю. Предположим, что заряд движется вдоль оси OX. Введем цилиндрическую систему координат с осью, совпадающей с траекторией движения, положим
(9)
Подставляя зависимость (9) в первое уравнение (8), находим
(10)
Разрешая уравнение (10), окончательно получим
(11)
Используя решение (11), найдем компоненты электрического и магнитного поля в виде
(12)
Выражения (12) описывают излучение заряженной частицы, движущейся со сверхсветовой скоростью. Полное решение задачи о поле заряда, движущегося со сверхсветовой скоростью, является суммой решений, первое из которых соответствует частному решению первого уравнения (8) для точечного заряда, т.е. определяется на основе потенциалов Лиенара-Вихерта /7/, а второе, описывающее поле излучения, задается в виде (12).
Неопределенную константу, фигурирующую в выражениях (12) можно доопределить, предположив, что в лобовой точке сферического заряда конечного радиуса радиальная компонента электрического поля обращается в ноль, т.е. . Отсюда находим
Здесь - радиус сферического заряда. Следовательно, имеем
(13)
Из полученных выражений (12-13) следует, что электромагнитное поле точечного заряда, движущегося со сверхсветовой скоростью, определено в конусе, вершина которого совпадает с зарядом, а угол раствора определяется из уравнения .
Согласно последнему уравнению (12), магнитное поле излучения заряда, движущегося со сверхсветовой скоростью во вращающейся системе координат, не зависит от расстояния. Действительно, в общем случае выражение (1) необходимо модифицировать с учетом влияния релятивистских эффектов. Используя выражения (7) и (13), получим
Подставляя сюда выражение электрического поля (13), находим, что магнитное поле излучения заряда содержит слагаемое, не зависящее от расстояния. Это слагаемое появляется лишь при относительном сверхсветовом движении. Если, например, заряд Нептуна излучает волны Вавилова-Черенкова в какой-либо проводящей среде - ионосфере, морской воде, ядре или мантии, то эти волны отразятся на величине магнитного поля Земли. Согласно (13), интенсивность этих волн зависит от величины электростатического потенциала на поверхности небесного тела.