Home / Section index | www.icosaedro.it |
Updated: 2019-06-24
In this article we will show how to read the ephemeris tables and how to locate the celestial bodies by using a magnetic compass and the astrolabe we build in the previous article Celestial Navigation at Home.
The Azimuth is the horizontal angular distance of the celestial body from the North direction. If that North is the geographic North Pole, then we have the Geographic Azimuth; if that North is the Magnetic North, then we have the Magnetic Azimuth. Since we will use a magnetic compass, this latter is the Azimuth we are dealing with.
The Elevation, or Altitude, of the celestial body is its angular distance from the horizon.
Once the our position, the local magnetic variation, the Magnetic Azimuth and the Elevation of the celestial body are known, we can easily locate that celestial body in the sky by using our compass and astrolabe.
Our latitude and longitude can be obtained from several sources: from the world map, from a GPS receiver, or by visiting maps.google.com. In this latter case, the URL tells something like this:
https://www.google.com/maps/@44.4950661,11.3535124,15z
where +44.50 is my latitude and -11.35 is my longitude; note how the sign of the longitude is reversed, but apart from that the longitude can be easily recognized anyway.
By definition, the Magnetic North is the direction of the projection of the Earth magnetic field vector on the horizontal plane. For this reason the needle of the magnetic compass is so thin and flat and must rotate freely over the horizontal plane. The Magnetic North does not coincide with the Geographic North, but instead it deviates more or less from it, either toward East or toward West. The amplitude of this deviation, also named magnetic variation, varies from site to site and varies over time, so it must be checked and updated periodically, at least every few months. The magnetic variation is positive if the needle of the compass deviates toward East respect to the Geographic North, and it is negative if it deviates toward West.
The Earth magnetic field is constantly measured and several mathematical models have been developed in an attempt to predict its evolution over time for any location in the world. For example, on https://www.ngdc.noaa.gov/geomag/calculators/magcalc.shtml we can choose among several mathematical models and get the respective estimation of the magnetic variation. In my case I found:
which finally tells the local magnetic variation I must account for today is +3.1° ± 0.3°. This value should then be re-checked again later, for example after some week or some month.
For more about the Earth magnetic field and the magnetic variation, see Magnetic declination (Wikipedia).
We will use the Nautical Almanac TheNauticalAlmanac.com we already see in the article before, where the positions of the Sun, of the Moon, of several planets and several stars are reported, day by day and hour by hour.
Locate the table of the ephemeris about the Sun for the day of interest, for example 2019-06-24 12:00 UTC:
The data we are interested to are highlighted in the figure above; it is convenient to translate a bit these number to fit our needs as we already shown in the previous article:
Sun_lat = N23° 24.5' = +23.41°
Sun_lon = 359° 24.1' GHA = 359.40° GHA = -0.60°
For the Moon and the planets we may proceed in the same exact way. For example, Saturn 2019-06-24 02:00 UTC:
Here again, we translate these numbers according to our convention obtaining:
Saturn_lat = S21° 50.6' = -21.84°
Saturn_lon = 12° 07.2' GHA = +12.12°
Locate the table for the stars for the days 24, 25 and 26 of June 2019:
Note that there is just only one single table for the entire period of three days. This is because the location of the stars in the sky changes so slowly over time that them are also referred as fixed stars, although nothing is fixed in the Universe. Anyway, the "Dec" column reports the latitude just like we saw for Sun and planets, while the SHA column reports the sidereal hourly angle rather that the GHA (Greenwich hourly angle). The SHA is just like the GHA but the angle is measured respect to the meridian of Aries (or γ) point. The GHA of the star is then the sum of its SHA and the GHA of Aries at the given date and time:
Star_GHA = Star_SHA + Aries_GHA
The result of the formula above must be normalized to the range [0°,360°[, so if it's greater or equal to 360°, subtract 360°. The GHA of the star can then be easily converted to its longitude:
Star_lon = Star_GHA if Star_GHA ≤ 180° Star_lon = 360° - Star_GHA if Star_GHA > 180°
In our example , the current GHA of Aries is reported in the table to the left in the same page:
Aries_GHA = 301° 56.5'
Lets take for example Aldebaran. Its latitude can be read straight from the table, while its longitude must be calculated by adding its SHA (more or less "fixed") to the current GHA of Aries:
Aldebaran_lat = 16° 32.7' = +16.55 Aldebaran_lon = Aldebaran_SHA + Ariete_GHA = 290° 45.0' + 301° 56.5' (replacing numbers in formula) = 290.75° + 301.94° (converting to decimal degrees) = 592.69° (performing addition) = 232.69° (normalizing GHA to the [0°,360°[ range) = -127.31° (normalizing lon to the [-180°,+180°[ range)
Now that our position P and the position of the celestial body S are known, it's a matter to solve the spherical triangle with vertices P, S and N (the North Pole):
Two sides 90°-Slat and 90°-Plat, and the included angle |Slat-Plat| are known. From them we may calculate all the remaining elements of the spherical triangle PSN, that is the opposite side 90°-Selev and the spherical angle in P giving the Geographic Azimut of the star. Care must be paid to the longitudes of S and P, as the Azimuth changes depending from their relative positions (ence the two drawings above). For the details of the calculation, see for example https://en.wikipedia.org/wiki/Spherical_trigonometry. To the Geographic Azimuth we must subtract the local magnetic variation to obtain the Magnetic Azimuth; the elevation should be corrected for the atmospheric refraction.
The details of the calculation are quite boring and are summarized briefly here below:
If Slon > Plon: A = Slon - Plon and star_west = true If Slon ≤ Plon: A = Plon - Slon and star_west = false If A > 180°: A = 360° - A and star_west = NOT star_west c = 90° - Slat b = 90° - Plat a = acos( cos(b) cos(c) + sin(b) sin(c) cos(A) ) s = (a + b + c)/2 C = 2 asin( sqrt( (sin(s-a) sin(s-b)) / (sin(a) sin(b)) ) ) If star_west is true: Saz = 360° - C If star_west is false: Saz = C Sazmag = Saz - Pmagvar Selev = 90° - a
where Saz il the geographic Azimuth of the star, Sazmag is its magnetic Azimuth, and Pmagvar is the local magnetic variation, positive East.
The page Star Relative Location Calculator can do all these calculations for you. Enter the position, the local magnetic variation, and the position of the celestial body (for example, Saturn):
and then press the Compute button to get:
Final result: the day 2019-06-24 02:00 UTC, Saturn will be visible in the sky at 200° magnetic azimuth, 20° over the horizon.
Umberto Salsi | Comments | Contact | Site map | Home / Section index |
Still no comments to this page. Use the Comments link above to add your contribute.