Magellan's calculations are based on equations Talwani published in 1964. He takes a 2-d polygon and calculates the contribution from each side of the polygon to the magnetic anomaly at the origin of the coordinate system. The polygon is assumed to extend infinitely along the y-axis.

The calculated anomaly has two components, one vertical and one horizontal. Note that this anomaly is the measured magnetic field after the total magnetic field, at the time of formation of the lithosphere, has been removed. These components are given, respectively, by the following equations

$\displaystyle V = 2 (J_x Q - J_z P)$ (1)

$\displaystyle H = 2 (J_x P + J_z Q)$ (2)

where

$\displaystyle P = \frac{z^{2}_{21}}{z^{2}_{21} + x^{2}_{12}} (\theta_1 - \theta_2) + \frac{z_{21} x_{12}}{z^{2}_{21} + x^{2}_{12}} log(\frac{r_2}{r_1})$ (3)

$\displaystyle Q = \frac{z_{21} x_{12}}{z^{2}_{21} + x^{2}_{12}} (\theta_1 - \theta_2) - \frac{z^{2}_{21}}{z^{2}_{21} + x^{2}_{12}}log(\frac{r_2}{r_1})$ (4)

and additionally


$\displaystyle x_{12} = x_1-x_2,$   $\displaystyle z_{21} = z_2-z_1,$  
$\displaystyle r_1 = \sqrt{x^2_{1} + z^2_{1}},$   $\displaystyle r_2 = \sqrt{x^2_{2} + z^2_{2}}$  

Here the $ \vec{J}$ is the magnetization vector where its horizontal and vertical components are, respectively,


$\displaystyle J_x = Jcos(I)cos(D-B)$   $\displaystyle J_z = Jsin(I)$  

The magnetic anomalies that are measured arise due to magnetization of the oceanic floor and this magnetization is by induction, that is $ \vec{J} = k\vec{F}$ where k is the magnetic susceptibility and $ \vec{F}$ is the total magnetic field. Here I is the inclination of the total magnetic field, D is the declination and B is the strike of spreading direction clockwise from north.

The total anomaly vector is a projection of its horizontal and vertical components along the direction of $ \vec{F}$

$\displaystyle \vec{T} = Vsin(I) + Hcos(I)cos(C-D)$ (5)

Magellan's algorithm is based on these equations and assumes that the total field is either in the direction of Earth's current magnetic field or in the reverse direction. For instance if we are calculating the contribution from a block that is reversly magnetized then we multiply $ \vec{T}$ with a minus sign which results in a negative contribution from that block.


Ásdís Benediktsdóttir 2009-02-21