–PAGE_BREAK–Магнитное поле земли
Английский ученый Уильям Гильберт, придворный врач королевы Елизаветы, в 1600 г. впервые показал, что Земля является магнитом, ось которого не совпадает с осью вращения Земли. Следовательно, вокруг Земли, как и около любого магнита, существует магнитное поле. В 1635 г. Геллибранд обнаружил, что поле земного магнита медленно меняется, а Эдмунд Галлей провел первую в мире магнитную съемку океанов и создал первые мировые магнитные карты (1702 г.). В 1835 г. Гаусс провел сферический гармонический анализ магнитного поля Земли. Он создал первую в мире магнитную обсерваторию в Гёттингене.
О распределении силовых линий магнитного дипольного поля и о магнитных полюсах наклонения Пс, Пю можно судить по рисунку.
Составляющие геомагнитного поля определены следующим образом. В любой точке О вектор напряженности магнитного поля В может быть разложен на составляющие, как это показано на рисунке. Можно выбрать в качестве составляющих абсолютную величину полного вектора В (модуль) и два угла: Dи I. Угол Dобразован направлением на север и горизонтальной составляющей вектора В, т. е. Н; I– это угол между В и Н, Угол D считается положительным, если Н отклоняется к востоку, а I положительно при отклонении В вниз от горизонтальной плоскости. Величина D называется магнитным склонением, а I– наклонением. Вертикальная плоскость, которая проходит через Н, именуется местной магнитной меридиональной плоскостью.
Используется также разложение В на северную (X) и восточную (Y) составляющие вектора Н. Третьей служит вертикальная составляющая Z, которая считается положительной, если В направлено вниз. Напряженности B, H, Z, X, Yизмеряются в гауссах (Гс) или гаммах (g). 1g=10-5Гс. Углы D и I измеряются в дуговых градусах и минутах. Все приведенные семь величин В, Н, D, I, X, У, Zназываются магнитными элементами. Соотношения между ними ясны из рисунка.
H=B cos I, Z=B sin I=H tg I,
X=H cos D, Y=H sin D,
X2+Y2=H2 X2+Y2+Z2=H2+Z2=B2
Ясно, что для полного описания вектора В достаточно иметь три независимых элемента. По ним могут быть рассчитаны все остальные.
Обычная стрелка магнитного компаса уравновешивается, вращаясь горизонтально на вертикальной оси. В северной полусфере Земли почти везде северный полюс магнитной стрелки направлен вниз (т. е. Iположительно), а в южном полушарии Iотрицательно, поскольку вниз направлен южный полюс стрелки. Линия, которая разделяет области положительного и отрицательного I, называется магнитным экватором или экватором наклонения. Естественно, что на ней I=0, т. е. магнитная стрелка в любой точке на этой кривой располагается горизонтально.
На полюсах магнитного наклонения горизонтальная компонента полного вектора В исчезает и магнитная стрелка устанавливается вертикально. Эти точки еще называют полюсами наклонения. Таких точек в принципе может быть несколько. Две основные из них обычно называются магнитными полюсами Земли. Они расположены в Арктике и в Антарктиде. Координаты их 75°,6 с. ш., 101° з. д. и 66°,3 ю.ш., 141° в. д. Местоположение магнитных полюсов не является постоянным. Приведенные выше координаты относятся к эпохе 1965 г.
Чтобы определить азимут[3] вектора Н, нужно выбрать некоторое нулевое направление, от которого можно отсчитывать магнитное склонение D. За такое направление принято направление на северный географический полюс. Таким образом, D определяется относительно условного направления, поскольку ось вращения Земли не связана непосредственно с конфигурацией геомагнитного поля. То же относится и к элементам Х и Y. Поэтому D, X, Y называют относительными магнитными элементами, тогда как H, Z и I именуются собственными магнитными элементами.
Несколько слов о магнитных картах. Обычно через каждые 5 лет распределение магнитного поля на поверхности Земли представляется магнитными картами трех или более магнитных элементов. На каждой из таких карт проводятся изолинии, вдоль которых данный элемент имеет постоянную величину. Линии равного склонения Dназываются изогонами, наклонения I– изоклинами, величины полной силы В – изодинамическими линиями или изодинами. Изомагнитные линии элементов H, Z, Х и Yназываются соответственно изолиниями горизонтальной, вертикальной, северной или восточной компонент.
Направление оси магнитного диполя практически не меняется с 1829 г. При этом магнитный момент диполя систематически уменьшался. Его уменьшение может быть аппроксимировано выражением
m=(15,77-0,003951t)×1025 Гс×см3,
где t— время в годах, отсчитываемое вперед пли назад от 1900 г. н. э. По этой формуле можно рассчитать, что если уменьшение магнитного момента будет продолжаться с такой же скоростью, то к 3991 г. магнитный момент станет равным нулю.
Мы будем постоянно иметь дело с геомагнитными силовыми линиями, а также различного рода координатами.
Геомагнитные дипольные координаты — это дополнение к широте q’ и восточной долготе j’. Они определяются относительно полярной оси и нулевого меридиана. Если точка Р имеет географические координаты qи j, то геомагнитные координаты могут быть вычислены по следующим формулам:
cosq’=-cosqcosq0 — sinqsinq0 cos(j-j),
sinj’=sinq×sin(j-j) cosecq’.
Магнитное склонение дипольного поля Y– это угол, образованный магнитным и географическим меридианами в точке Р. Он определяется из выражения
sin(–y)= sinq(sin(j-j0)/sinq’)
Существуют таблицы, которые содержат геомагнитные координаты сетки точек, расположенных через ровные угловые интервалы в географических координатах qи j. Имеются также сетки географических и геомагнитных координат. По этим сеткам можно легко найти геомагнитные координаты любой точки с известными географическими координатами, и наоборот.
Обратный переход от геомагнитных координат к географическим можно произвести по формулам
cosq=cosq’ ×cosq0 – sinq’ ×sinqcosj’
Если рассматривать только дипольную часть геомагнитного поля в любой точке Р с геомагнитными координатами q’ и j’, то потенциал V1, описываемый членами первого порядка, равен V1= –m(cosq/r2) Tак как V1не зависит от долготы, то восточная компонента дипольного поля В равна нулю. Северная Я и вертикальная Z составляющие поля получаются равными
H=m(sinq’/r3)=H0(a/r)3sinq’,
Z=2m(cosq’/r3)=Z0(a/r)3cosq’; Z0=2H0
где Z0и Н0– максимальные значения Z и Hна геоцентрической сфере радиуса а, содержащей точку Р. Hсоответствует полю на геомагнитном экваторе, а Z0– на северном полюсе. На южном полюсе Z= –Z.
Наклонение Iи магнитную широту l’ можно определить из следующих уравнений:
tgI=(Z/H)2ctgq’, tgl’=1/2tgI.
Каждая силовая линия дипольного поля лежит в плоскости геомагнитного меридиана. Ее уравнение
r=re×sin2q’
где re– радиальное расстояние, на котором данная силовая линия пересекает плоскость геомагнитного экватора, с величиной поля равной m/re3Величину re, можно принять за параметр, определяющий силовую линию.
продолжение
–PAGE_BREAK–