http://SaturnianCosmology.Org/ mirrored file For complete access to all the files of this collection see http://SaturnianCosmology.org/search.php ========================================================== Eclipse Banner Blue Bar Visibility of Transits of Mercury and Venus Fred Espenak Blue Bar Calculating the Sun's Altitude and Azimuth To determine whether a transit is visible from a specific geographic location, it is simply a matter of calculating the Sun's altitude and azimuth during each phase of the transit. The calculations can be performed on any pocket calculator having trig functions (SIN, COS, TAN). Armed with the latitude and longitude of the location, the transit catalogs provide all the additional information needed to make the calculations. The altitude 'a' and azimuth 'A' of the Sun during any phase of a transit depends on the time and the observer's geographic coordinates. Neglecting the effects of atmospheric refraction and planetary parallax, 'a' and 'A' are calculated as follows: h = 15 * (GST + t - ra ) + l a = ArcSin [ Sin d Sin f + Cos d Cos h Cos f ] A = ArcTan [ - (Cos d Sin h) / (Sin d Cos f - Cos d Cos h Sin f) ] where: h = Hour Angle of the Sun (in degrees) a = Altitude (in degrees) A = Azimuth (in degrees) GST = Greenwich Sidereal Time at 00:00 UT t = Universal Time ra = Right Ascension of the Sun (in hours) d = Declination of the Sun (in degrees) l = Observer's Longitude (East +, West -) f = Observer's Latitude (North +, South -) The Venus catalog lists the following record for the transit of 2004 June 08: Transit Contact Times (UT) ------------------------------------- Minimum Sun Sun Transit Date I II Greatest III IV Sep. RA Dec GST Series h:m h:m h:m h:m h:m " h ° h 2004 Jun 08 05:13 05:33 08:20 11:07 11:26 626.9 5.121 22.89 17.137 3 Using the preceeding formulas and the data above, we will now determine whether the Sun will be above the horizon at greatest transit as seen from Washington DC. The geographic coordinates of Washington DC are: Latitude: f = 38°53´N = +38.9° Longitude: l = 077°02´W = -077.0° From the transit catalog record, we have: Time of Greatest Transit: t = 08:20 = 8.333 Greenwich Sidereal Time at 00:00 UT: GST = 17.137 Right Ascension of the Sun: ra = 5.121 Declination of the Sun: d = 22.89 Thus: Hour Angle of the Sun: h = 15 * (GST + t - ra ) + l = 15 * (17.137 + 8.333 - 5.121) + -077.0 = 15 * (20.346) + -077.0 h = 228.2° Altitude of Sun: a = ArcSin [Sin d Sin f + Cos d Cos h Cos f] = ArcSin [Sin(22.89) Sin(38.9) + Cos(22.89) Cos(228.2) Cos(38.9)] = ArcSin [0.386 * 0.628 + 0.921 * -0.666 * 0.778] = ArcSin [0.242 + -0.477] = ArcSin [-0.235] = -13.6° With an altitude of -13°, the Sun will be below the horizon and will not be visible at greatest transit as seen from Washington DC. Repeating the calculation using the time of contact III (at 11:07 UT), reveals that the Sun's altitude is then 14°. Thus Washington DC will see the end of the transit. The Sun's azimuth can be calculated similarly using the above formula. Geocentric contact times have been used in the above calculations in order to greatly simplify the problem. This approximation may introduce of an error in the actual contact times of up to 3 minutes for Mercury and up to 10 minutes for Venus. This translates into a possible error in the calculated altitudes of a up to 1° for Mercury and 2.5° for Venus. Thus, the proceedure has more that enough accuracy for most historical inquiries. As an aid to historical research, four Excel 97 spreadsheet files have been prepared which perform the above calculations automatically. You simply enter the location name, latitude and longitude. Each of the tables then calculates the altitude of the Sun at that location for every contact and for every transit in the table. The four tables are similar but cover different time periods for Mercury and Venus: / // Transits of Mercury: 1601 CE - 1900 CE / (Excel File - 84 KB) / Transits of Mercury: 1901 CE - 2300 CE / (Excel File - 100 KB) / Transits of Venus: 2000 BCE - 1000 CE / (Excel File - 80 KB) / Transits of Venus: 1000 CE - 4000 CE / (Excel File - 84 KB) Please note that these files will not open properly unless you have Excel 97 (or newer) installed on your computer. When each of the spreadsheets is downloaded, it will be opened automatically in Excel where you will be able to enter the coordinates of any geographic location to calculate the transit circumstances. The spreadsheets are protected so that you can not accidently delete or edit any information required by the calculations. Only the name and coordinates of the geographic location (in the green box of each spreadsheet) may be modified. Blue Bar Transit Predictions Transit predictions are based on algoritmhms and elements published in "Transits" by Jean Meeus (Willmann-Bell, 1989). The value for /delta-T <../../SEhelp/deltaT.html>/ was determined as follows: 1) pre-1600: delta T was calculated from empirical expressions by Stephenson and Houlden [1986] 2) 1600-present: delta T was obtained from published observations 3) future: delta-T was extrapolated from current values All transit calculations are by Fred Espenak, and he assumes full responsibility for their accuracy. Special thanks to eclipse chaser *Michael Gill* for enthusiastically and rigorously beta-testing this page and for catching typographic errors. (2003 March) Permission is freely granted to reproduce this data when accompanied by an acknowledgment: "Transit Predictions by Fred Espenak, NASA/GSFC" Blue Bar Return to / Catalog of Mercury Transits / Return to / Catalog of Venus Transits / Return to / Transits of the Sun <../transit.html>/ Return to /Eclipse Home Page <../../eclipse.html>/ Blue Bar NASA logo Transit Predictions & WebMaster: Fred Espenak e-mail: espenak at gsfc.nasa.gov Planetary Systems Branch - Code 693 NASA/Goddard Space Flight Center, Greenbelt, Maryland 20771 USA /Last revised: 2003 Mar 24 - F. Espenak/