# Moon Phases and Eclipses

by Phill Edwards

## Moon Phases

To calculate the Moon phase, first calculate the longitude of the Earth λE and the longitude of the Moon λM for the date. Then calculate the lunar solar elongation D = λE - λM. The Moon phase can then be determined from the phase angle as shown in the table.

Phase Angle Phase
D = 0° New Moon
0° < D < 90° Waxing Crescent
D = 90° First Quarter
90° < D < 180° Waxing Gibbous
D = 180° Full Moon
180° < D < 270° Waning Gibbous
D = 270° Last Quarter
00° < D < 180° Waning Crescent

## Eclipses

For a Solar eclipse to occur it needs to be close to a New Moon. New and Full Moons where the Sun, Earth and Moon are aligned are called syzygy. There is not an eclipse at every New Moon because the Moon's orbit is inclined at about 5° to the ecliptic. So, for a Solar eclipse to occur the lunar latitude at the syzygy βsyz needs to be within a certain range.
The angular radius of the Sun ρS = RS/rE
The angular radius of the Moon ρM = RS/rM
The maximum parallax of the Moon δM = RE/rM
An eclipse, either total, annular or partial, occurs if |βsyz| < δM + ρM + ρS
A total eclipse occurs if ρM > ρS and |βsyz| < δM + ρM - ρS
An annular eclipse occurs if ρS > ρM and |βsyz| < δM + ρS - ρM

### Lunar Parallax

The Moon is sufficiently large that the Moon as seen from a hypothetical viewer at the centre of the Earth (geocentric) will appear in a different postion to an observer on the Earth's surface (topocentric). In fact observers at different latitudes on the Earth's surface will see the Moon at a different position. This effect is due to parallax. There is also a similar, but much smaller parallax offset for the Sun. The Solar parallax is considered a derived astronomical constant with a value of 8".794. Only the sine of this angle is required for calculations.

Given the following parameters:

• λ is the geocentric ecliptical longitude of the Moon
• β is the geocentric ecliptic latitude of the Moon
• s is the semi diameter of the Moon
• φ is the observer's geographic latitude
• ρ and ψ define the observer's geocentric latitude as described in Coordinates
• ε is the obliquity of the ecliptic
• θ is the local sidereal time
• π is the equatorial horizontal parallax of the Sun
We need to calculate the topocentric values λ', β' and s'.

Then: $S=\rho sin\psi$ $C=\rho cos\psi$ $N=cos\lambda cos\beta -Csin\pi sin\theta$

$tan\mathrm{\lambda \text{'}}=\frac{sin\lambda cos\beta -sin\pi \left(Ssin\epsilon +Ccos\epsilon sin\theta \right)}{N}$

$tan\mathrm{\beta \text{'}}=\frac{cos\mathrm{\lambda \text{'}}\left(sin\beta -sin\pi \left(Scos\epsilon -Csin\epsilon sin\theta \mathrm{\right)\right)}}{N}$

$sin\mathrm{s\text{'}}=\frac{cos\mathrm{\lambda \text{'}}cos\mathrm{\beta \text{'}}sins}{N}$

### Eclipses Revisited

Terms:

• λS the solar longitude
• λM the lunar longitude
• D = λS - λM the lunar solar elongation
• β the lunar latitude
• rS the Earth Sun distance
• rM the Earth Moon distance
• rMS the Moon Sun distance
• rUM the Moon umbral point distance

At the new moon syzygy D = 0
By the cosine rule rMS2 = rS2 + rM2 - 2rSrMcosβ
${r}_{\mathrm{UM}}=\frac{{R}_{M}{r}_{\mathrm{MS}}}{{R}_{S}-{R}_{M}}$

## Besselian Elements

Besselian Elements were introduced in 1824 by the Prussian astronomer and meathematician Frederick Bessel. They are a good method of calculating the characteristics and locations of an eclipse which are still valid. The principle is to define a fundamental plane, which is a plane which passes through the Earth's centre and is normal to the Sun-Moon centre line. The x axis point East, the y axis point North and the z axis point towards the Sun. There are six Besselian Elements which are:

• x, y the x and y positions of the centre of the shadow cone in units of the Earth's equatiorial radius
• d the declination of the shadow axis
• μ the hour angle of the shadow axis
• L1 the radius of the penumbral shadow on the fundamental plane
• L2 the radius of the umbral shadow on the fundamental plane
• f1 the angle which the penumbral shadow cone makes with the shadow axis
• f2 the angle which the umbral shadow cone makes with the shadow axis

The Besselian Elements are normally calculated accurately at five hourly intervals and then a least squares interpolation is performed. This may not be necessary using modern computers. The elements are calculated for the nearest hour to the maximum eclipse and for the two hours before and after the greatest eclipse.

Let r be the Earth-Moon vector and R be the Sun-Earth vector and k = (0, 0, 1) is the unit vector pointing North.
The unit vector K which is normal to the fundamental plane, pointing in the direction of right ascension a and declination b is:

$K=\frac{R-r}{|R-r|}=\left(cosdcosa,cosdsina,sind\right)$

The other two unit vectors in the fundamental plane are the unit vectors created from the cross products:

$I=\frac{Kxk}{|Kxk|}=\left(-sina,cosa,0\right)$ $J=\frac{IxK}{|IxK|}=\left(-cosasind,-sinasind,cosd\right)$

The Besselian Elements x, y and d plus the z coordinate are:

$x=r.I$ $y=r.J$ $z=r.K$ $sind=K.k$

The hour of the shadow axis μ in terms of the right ascension a and Grenwich Sidereal Time Ts is:

$\mu ={T}_{s}-a$

The penumbral and umbral cone angles f1 and f2 are:

$sin\mathrm{f1}=\frac{{R}_{S}+{R}_{M}}{|R-r|}$ $sin\mathrm{f2}=\frac{{R}_{S}-{R}_{M}}{|R-r|}$

The penumbral and umbral cone vertices c1 and c2 are:

$\mathrm{c1}=z+\frac{{R}_{M}}{{R}_{E}}cosec\mathrm{f1}$ $\mathrm{c2}=z-\frac{{R}_{M}}{{R}_{E}}cosec\mathrm{f2}$

The penumbral and umbral cone radii in the fundamental plane L1 and L2 are:

$\mathrm{L1}=\mathrm{c1}tan\mathrm{f1}$
$\mathrm{L2}=\mathrm{c2}tan\mathrm{f2}$

## Besselian to Geocentric Coordinates

Given the Earth's equatorial radius a = 6378.1km and the Earth's polar radius b = 6356.8. Then the flattening f and the ellipticity e are defined as:

$f=\frac{a-b}{b}$ $e=\sqrt{1-\frac{{b}^{2}}{{a}^{2}}}$

Given an observer's longitude λ, latitude φ and distance from the centre of the Earth ρ then the observer's Geocentric position ρG is:

${\rho }_{G}=\rho \left(cos\phi cos\lambda ,cos\phi sin\lambda ,sin\phi \right)$

The observer's location in the fundamental plane is:

${\rho }_{F}=\left(\xi ,\eta ,\zeta \right)$

The x, y and z axis rotation matrices are R1, R2 and R3.

Then the relationship between observer coordinates is:

${\rho }_{F}={R}_{1}\left(\frac{\pi }{2}-d\right){R}_{3}\left(\frac{\pi }{2}-\mu \right){\rho }_{G}$