# 3.10: A Geophysical Example

Assume that planet Earth is spherical and that it has a little magnet or current loop at its centre. By “little” I mean small compared with the radius of the Earth. Suppose that, at a large distance from the magnet or current loop the geometry of the magnetic field is the same as that of an electric field at a large distance from a simple dipole. That is to say, the equation to the lines of force is \(r=a\sin^2 \theta \)(equation 3.7.13), and the differential equation to the lines of force is \(\frac{dr}{dθ} = \frac{2r}{tan \theta}\) (equation 3.7.12).

Show that the angle of dip \(D\) at geomagnetic latitude \(L\) is given by

\[\tan D = 2\tan L \]

The geometry is shown in figure \(III\).16.

The result is a simple one, and there is probably a simpler way of getting it than the one I tried. Let me know (jtatum@uvic.ca) if you find a simpler way. In the meantime, here is my solution.

I am going to try to find the slope \(m_1\) of the tangent to Earth (i.e. of the horizon) and the slope \(m_2\) of the line of force. Then the angle *D* between them will be given by the equation (which I am hoping is well known from coordinate geometry!)

\[\tan D = \frac{m_1 -m_2}{1+m_1m_2}.\]

The first is easy:

\[m_1=\tan (90^\circ + \theta ) = -\frac{1}{\tan \theta}\label{3.10.3}\]

For \(m_2\) we want to find the slope of the line of force, whose equation is given in polar coordinates? So, how do you find the slope of a curve whose equation is given in polar coordinates? We can do it like this:

\[\begin{align} x&=r\cos \theta , \\ y&=r\sin \theta , \\ dx&=\cos \theta dr - r\sin \theta d\theta , \\ dy&=\sin \theta dr + r\cos \theta d\theta . \\ \end{align}\]

From these, we obtain

\[\frac{dy}{dx}=\frac{\sin \theta \frac{dr}{d\theta}+r\cos \theta }{\cos \theta \frac{dr}{d\theta}-r\sin \theta }.\label{3.10.8}\]

In our particular case, we have \(\frac{dr}{d\theta}=\frac{2r}{\tan \theta}\) (equation 3.7.12), so if we substitute this into Equation \ref{3.10.8} we soon obtain

\[m_2 = \frac{3\sin \theta \cos \theta}{3\cos^2 \theta -1}.\label{3.10.9}\]

Now put Equations \ref{3.10.3} and \ref{3.10.9} into equation 3.10.2, and, after a little algebra, we soon obtain

\[\tan D = \frac{2}{\tan \theta}=2\tan L . \]

Here is another question. The magnetic field is generally given the symbol \(B\). Show that the strength of the magnetic field \(B(L)\) at geomagnetic latitude \(L\) is given by

\[B(L) = B(0)\sqrt{1+3\sin^2 L } ,\]

where \(B(0)\) is the strength of the field at the equator. This means that it is twice as strong at the magnetic poles as at the equator.

Start with equation 3.7.2, which gives the electric field at a distant point on the equator of an electric dipole. That equation was \(E=\frac{p}{4\pi\epsilon_0 y^3}\). In this case we are dealing with a magnetic field and a magnetic dipole, so we’ll replace the electric field \(E\) with a magnetic field \(B\). Also \(p/(4π\epsilon_0)\) is a combination of electrical quantities, and since we are interested only in the geometry (i.e. on how \(B\) varies from equation to pole, let’s just write \(p/(4\pi\epsilon_0 )\) as \(k\). And we’ll take the radius of Earth to be \(R\), so that equation 3.7.2 gives for the magnetic field at the surface of Earth on the equator as

\[B(0) = \frac{k}{R^3}.\]

In a similar vein, equations 3.7.10 for the radial and transverse components of the field at geomagnetic latitude \(L\) (which is 90º − \(θ\)) become

\[B_r ( L) = \frac{2k\sin L}{R^3} \quad \text{and} \quad B_{\theta}(L)=\frac{k\cos L}{R^3}.\]

And since \(B=\sqrt{B_r^2 + B_{\theta}^2}\) the result immediately follows.