Lunar Limb Profile

Eclipse contact times, magnitude and duration of totality all depend on the angular diameters and relative velocities of the Moon and Sun. Unfortunately, these calculations are limited in accuracy by the departure of the Moon's limb from a perfectly circular figure. The Moon's surface exhibits a rather dramatic topography, which manifests itself as an irregular limb when seen in profile. Most eclipse calculations assume some mean radius that averages high mountain peaks and low valleys along the Moon's rugged limb. Such an approximation is acceptable for many applications, but if higher accuracy is needed, the Moon's actual limb profile must be considered. Fortunately, an extensive body of knowledge exists on this subject in the form of Watts' limb charts [Watts, 1963]. These data are the product of a photographic survey of the marginal zone of the Moon and give limb profile heights with respect to an adopted smooth reference surface (or datum). Analyses of lunar occultations of stars by Van Flandern [1970] and Morrison [1979] have shown that the average cross-section of Watts' datum is slightly elliptical rather than circular. Furthermore, the implicit center of the datum (i.e., the center of figure) is displaced from the Moon's center of mass. In a follow-up analysis of 66,000 occultations, Morrison and Appleby [1981] have found that the radius of the datum appears to vary with libration. These variations produce systematic errors in Watts' original limb profile heights that attain 0.4 arc seconds at some position angles. Thus, corrections to Watts's limb data are necessary to ensure that the reference datum is a sphere with its center at the center of mass.

The Watts charts have been digitized by Her Majesty's Nautical Almanac Office in Herstmonceaux, UK, and transformed to grid-profile format at the U. S. Naval Observatory. In this computer readable form, the Watts limb charts lend themselves to the generation of limb profiles for any lunar libration. Ellipticity and libration corrections may be applied to refer the profile to the Moon's center of mass. Such a profile can then be used to correct eclipse predictions which have been generated using a mean lunar limb.

Along the path, the Moon's topocentric libration (physical + optical) in longitude ranges from l = -3.1° to l = -4.6°. Thus, a limb profile with the appropriate libration is required in any detailed analysis of contact times, central durations, etc.. But a profile with an intermediate value is useful for planning purposes and may even be adequate for most applications. The lunar limb profile presented in Figure 8 includes corrections for center of mass and ellipticity [Morrison and Appleby, 1981]. It is generated for 13:00 UT, which corresponds to western Zambia near the border with Angola. The Moon's topocentric libration is l = -4.39°, and the topocentric semi-diameters of the Sun and Moon are 944.3 and 988.0 arc seconds, respectively. The Moon's angular velocity with respect to the Sun is 0.364 arc seconds per second.

The radial scale of the limb profile in Figure 8 (at bottom) is greatly exaggerated so that the true limb's departure from the mean lunar limb is readily apparent. The mean limb with respect to the center of figure of Watts' original data is shown (dashed) along with the mean limb with respect to the center of mass (solid). Note that all the predictions presented in this publication are calculated with respect to the latter limb unless otherwise noted. Position angles of various lunar features can be read using the protractor marks along the Moon's mean limb (center of mass). The position angles of second and third contact are clearly marked along with the north pole of the Moon's axis of rotation and the observer's zenith at mid- totality. The dashed line with arrows at either end identifies the contact points on the limb corresponding to the northern and southern limits of the path. To the upper left of the profile are the Sun's topocentric coordinates at maximum eclipse. They include the right ascension R.A., declination Dec., semi-diameter S.D. and horizontal parallax H.P.. The corresponding topocentric coordinates for the Moon are to the upper right. Below and left of the profile are the geographic coordinates of the center line at 13:00 UT while the times of the four eclipse contacts at that location appear to the lower right. Directly below the profile are the local circumstances at maximum eclipse. They include the Sun's altitude and azimuth, the path width, and central duration. The position angle of the path's northern/southern limit axis is PA(N.Limit) and the angular velocity of the Moon with respect to the Sun is A.Vel.(M:S). At the bottom left are a number of parameters used in the predictions, and the topocentric lunar librations appear at the lower right.

In investigations where accurate contact times are needed, the lunar limb profile can be used to correct the nominal or mean limb predictions. For any given position angle, there will be a high mountain (annular eclipses) or a low valley (total eclipses) in the vicinity that ultimately determines the true instant of contact. The difference, in time, between the Sun's position when tangent to the contact point on the mean limb and tangent to the highest mountain (annular) or lowest valley (total) at actual contact is the desired correction to the predicted contact time. On the exaggerated radial scale of Figure 8, the Sun's limb can be represented as an epicyclic curve that is tangent to the mean lunar limb at the point of contact and departs from the limb by h through:

h = S (m - 1) (1 - cosC)      [8]

where:

h = departure of Sun's limb from mean lunar limb
S = Sun's semi-diameter
m = eclipse magnitude
C = angle from the point of contact

Herald [1983] has taken advantage of this geometry to develop a graphical procedure for estimating correction times over a range of position angles. Briefly, a displacement curve of the Sun's limb is constructed on a transparent overlay by way of equation [8]. For a given position angle, the solar limb overlay is moved radially from the mean lunar limb contact point until it is tangent to the lowest lunar profile feature in the vicinity. The solar limb's distance d (arc seconds) from the mean lunar limb is then converted to a time correction by:

[Equation 9]

where:

= correction to contact time (seconds)
d = distance of Solar limb from Moon's mean limb (arc sec)
v = angular velocity of the Moon with respect to the Sun (arc sec/sec)
X = center line position angle of the contact
C = angle from the point of contact

This operation may be used for predicting the formation and location of Baily's beads. When calculations are performed over a large range of position angles, a contact time correction curve can then be constructed.

Since the limb profile data are available in digital form, an analytical solution to the problem is possible that is quite straightforward and robust. Curves of corrections to the times of second and third contact for most position angles have been computer generated and are plotted in Figure 8. The circular protractor scale at the center represents the nominal contact time using a mean lunar limb. The departure of the contact correction curves from this scale graphically illustrates the time correction to the mean predictions for any position angle as a result of the Moon's true limb profile. Time corrections external to the circular scale are added to the mean contact time; time corrections internal to the protractor are subtracted from the mean contact time. The magnitude of the time correction at a given position angle is measured using any of the four radial scales plotted at each cardinal point.

For example, Table 15 gives the following data for Lusaka, Zambia:

Second Contact = 13:09:19.3 UT      P2=118°
Third Contact = 13:12:32.8 UT      P3=247°

Using Figure 8, the measured time corrections and the resulting contact times are:

C2 = +4.0 seconds  Second Contact = 13:09:19.3 + 4.0s  = 13:09:23.3 UT
C3 = -1.2 seconds  Third Contact = 13:12:32.8 - 1.2s  = 13:12:31.6 UT

The above corrected values are within 0.2 seconds of a rigorous calculation using the actual limb profile.

Lunar limb profile diagrams for a number of other positions/times along the path of totality are available via a special web site of supplemental material for the total solar eclipse of 2001 June 21.


Next section

Table of Contents