Declination
In astronomy, declination (abbrev. dec or δ) is one of the two coordinates of the equatorial coordinate system, the other being either right ascension or hour angle. Declination in astronomy is comparable to geographic latitude, but projected onto the celestial sphere. Declination is measured in degrees north and south of the celestial equator. Points north of the celestial equator have positive declinations, while those to the south have negative declinations.
- An object on the celestial equator has a declination of 0°.
- An object at the celestial north pole has a declination of +90°.
- An object at the celestial south pole has a declination of −90°.
The sign is customarily included even if it is positive. Any unit of angle can be used for declination, but it is often expressed in degrees, minutes, and seconds of arc.
As seen from locations in the Earth's Northern Hemisphere, celestial objects with declinations greater than 90° − φ (φ = observer's latitude) are always above the horizon. This similarly occurs in the Southern Hemisphere for objects with declinations less than, i.e. more negative than, -90° − φ. Such stars appear to circle daily around the celestial pole without dipping below the horizon, and are therefore called circumpolar stars. An extreme example is the pole star which has a declination near to +90°, so it is circumpolar as seen from anywhere in the Northern Hemisphere except very close to the equator. The Sun's declination varies with the seasons (see below). As seen from arctic or antarctic latitudes, the Sun is circumpolar near the local summer solstice, leading to the phenomenon of it being above the horizon at midnight, which is called midnight sun.
When an object is directly overhead its declination is almost always within 0.01 degree of the observer's latitude; it would be exactly equal except for two complications. The first complication applies to all celestial objects: the object's declination equals the observer's astronomic latitude, but the term "latitude" ordinarily means geodetic latitude, which is the latitude on maps and GPS devices. The difference (the vertical deflection) usually doesn't exceed a few thousandths of a degree. For practical purposes the second complication only applies to solar system objects: "declination" is ordinarily measured at the center of the earth, which isn't quite spherical, so a line from the center of the earth to the object is not quite perpendicular to the Earth's surface. It turns out that when the moon is directly overhead its geocentric declination can differ from the observer's astronomic latitude by up to 0.005 degree. The importance of this complication is inversely proportional to the object's distance from the earth, so for most purposes it's not a concern for the sun and planets.
Contents |
[edit] Stars
A star lies in a nearly constant direction as viewed from Earth, with its declination roughly constant from year to year, but right ascension and declination do both change gradually due to precession of the equinoxes, proper motion, and annual parallax.
The declinations of all solar system objects change much more quickly than those of stars.
[edit] Declination of the Sun as seen from Earth
The declination of the Sun, δ☉, is the angle between the rays of the Sun and the plane of the Earth's equator. The Earth's axial tilt (called the obliquity of the ecliptic by astronomers) is the angle between the Earth's axis and a line perpendicular to the Earth's orbit. The Earth's axial tilt changes gradually over thousands of years, but its current value is about ε = 23°26'. Because this axial tilt is nearly constant, solar declination (δ☉) varies with the seasons and its period is one year.
At the solstices, the angle between the rays of the Sun and the plane of the Earth's equator reaches its maximum value of 23°26'. Therefore δ☉ = +23°26' at the northern summer solstice and δ☉ = −23°26' at the southern summer solstice.
At the moment of each equinox, the center of the Sun appears to pass through the celestial equator, and δ☉ is 0°.
The Sun's declination at any given moment is calculated by:
Where EL is the ecliptic longitude. (On some calculator keyboards, and elsewhere, arcsin is written as sin−1.) Since the Earth's orbital eccentricity is small, its orbit can be approximated as a circle which causes up to 1 degree of error. The circle approximation means the EL would be 90 degrees ahead of the solstices in Earth's orbit (at the equinoxes), so that sin(EL) can be written as sin(90+NDS)=cos(NDS) where NDS is the number of days after the December solstice. By also using the approximation that arcsin[sin(d)*cos(NDS)] is close to d*cos(NDS), the following frequently used formula is obtained:
where N is the day of the year beginning with N=0 at midnight Coordinated Universal Time as January 1 begins (i.e. the days part of the ordinal date -1). The number 10, in (N+10), is the approximate number of days after the December solstice to January 1. This equation overestimates the declination near the September equinox by up to +1.5 degrees. The sine function approximation by itself leads to an error of up to 0.26 degrees and has been discouraged for use in solar energy applications. The 1971 Spencer formula[1] (based on a fourier series) is also discouraged for having an error of up to 0.28 degrees.[2] An additional error of up to 0.5 degrees can occur in all equations around the equinoxes if not using a decimal place when selecting N to adjust for the time after Coordinated Universal Time midnight for the beginning of that day. So the above equation can have up to 2.0 degrees of error, about 4 times the Sun's angular width, depending on how it's used.
The declination can be more accurately calculated by not making the two approximations, using the parameters of the Earth's orbit to more accurately estimate EL:[3]
which can be simplified by evaluating constants to:
N is the number of days since midnight Coordinated Universal Time as January 1 begins (i.e. the days part of the ordinal date -1) and can include decimals to adjust for local times later or earlier in the day. The number 2, in (N-2), is the approximate number of days after January 1 to the Earth's perihelion. The number 0.0167 is the current value of the eccentricity of the Earth's orbit. The eccentricity varies very slowly over time, but for dates fairly close to the present, it can be considered to be constant. The largest errors in this equation are less than +/- 0.2 degrees, but are less than +/- 0.03 degrees for a given year if the number 10 is adjusted up or down in fractional days as determined by how far the previous year's December solstice occurred before or after noon on December 22. These accuracies are compared to NOAA's advanced calculations[4][5] which are based on the 1999 Jean Meeus algorithm that is accurate to within 0.01 degree.[6]
(The above formula is related to a reasonably simple and accurate calculation of the Equation of Time, which is described here.)
More complicated algorithms[7] [8] correct for changes to the ecliptic longitude by using terms in addition to the 1st-order eccentricity correction above. They also correct the 23.44-degree obliquity which changes very slightly with time. Corrections may also include the effects of the moon in offsetting the Earth's position from the center of the pair's orbit around the Sun. After obtaining the declination relative to the center of the Earth, a further correction for parallax is applied, which depends on the observer's distance away from the center of the Earth. This correction is less than 0.0025 degrees. The error in calculating the position of the center of the Sun can be less than 0.00015 degrees. For comparison, the Sun's width is about 0.5 degrees. The declination calculations do not include the effects of the refraction of light in the atmosphere, which causes the apparent angle of elevation of the Sun as seen by an observer to be higher than the actual angle of elevation, especially at low Sun elevations. For example, when the Sun is at an elevation of 10 degrees, it appears to be at 10.1 degrees. The Sun's declination can be used, along with its right ascension, to calculate its azimuth and also its true elevation, which can then be corrected for refraction to give its apparent position.[9][5]
[edit] See also
- Celestial coordinate system
- Ecliptic
- Geographic coordinate system
- Inclination
- Lunar standstill
- Setting circles
- Euler angles
- Equation of time
- Right ascension
[edit] References
- ^ J. W. Spencer (1971). Fourier series representation of the position of the sun. http://www.mail-archive.com/sundial@uni-koeln.de/msg01050.html.
- ^ Sproul, Alistair B.. "Derivation of the solar geometric relationships using vector analysis". http://www.physics.arizona.edu/~cronin/Solar/References/Irradiance%20Models%20and%20Data/SPR07.pdf. Retrieved 28 February 2012.
- ^ "SunAlign". http://www.green-life-innovators.org/tiki-index.php?page=Sunalign. Retrieved 28 February 2012.
- ^ "NOAA Solar Calculator". Earth System Research Laboratory. http://www.esrl.noaa.gov/gmd/grad/solcalc. Retrieved 28 February 2012.
- ^ a b "Solar Calculation Details". Earth System Research Laboratory. http://www.esrl.noaa.gov/gmd/grad/solcalc/calcdetails.html. Retrieved 28 February 2012.
- ^ "Astronomical Algorithms". http://www.jgiesen.de/elevaz/basics/meeus.htm. Retrieved 28 February 2012.
- ^ "Computing the Solar Vector". Solar Energy 70 (5): 431–441. 2001. http://www.assembla.com/spaces/sun_follower/documents/d0PIA0oe0r347_eJe5avMc/download/10022519534118276.pdf.
- ^ Ibrahim Reda and Afshin Andreas. "Solar Position Algorithm for Solar Radiation Applications". http://www.nrel.gov/docs/fy08osti/34302.pdf. Retrieved 28 February 2012.
- ^ "Atmospheric Refraction Approximation". National Oceanic and Atmospheric Administration. http://www.srrb.noaa.gov/highlights/sunrise/atmosrefr.gif. Retrieved 28 February 2012.
[edit] External links
- NOAA's very accurate declination and sun position calculator (code can be viewed in the Javascript)
- Table of the Declination of the Sun: Mean Value for the Four Years of a Leap-Year Cycle (source unknown)
- Declination function for Excel, CAD or your other programs. The Sun API is free and extremely accurate. For Windows computers.
- How to compute planetary positions by Paul Schlyter.
![\delta_\odot = \arcsin \left [ \sin \left ( -23.44^\circ \right ) \cdot \sin \left ( EL \right ) \right ]](http://upload.wikimedia.org/wikipedia/en/math/0/6/e/06ee042f52c3bd529b2fbc676bdd4c6f.png)
![\delta_\odot = - 23.44^\circ \cdot \cos \left [ \frac{360^\circ}{365} \cdot \left ( N + 10 \right ) \right ]](http://upload.wikimedia.org/wikipedia/en/math/d/f/8/df812187e5aa53eff8d6f7cd5c7873ec.png)
![\delta_\odot = \arcsin \left [ \sin \left ( -23.44^\circ \right ) \cdot \cos \left ( \frac{360^\circ}{365.24} \left (N + 10 \right ) + \frac{360^\circ}{\pi} \cdot 0.0167 \sin \left ( \frac{360^\circ}{365.24} \left ( N - 2 \right ) \right ) \right ) \right ]](http://upload.wikimedia.org/wikipedia/en/math/0/8/8/088be3cec34d60ecf62d7c3f9b173115.png)
![\delta_\odot = - \arcsin \left [ 0.39779 \cos \left ( 0.98565 \left (N + 10 \right ) + 1.914 \sin \left ( 0.98565 \left ( N - 2 \right ) \right ) \right ) \right ]](http://upload.wikimedia.org/wikipedia/en/math/7/b/f/7bf1af04dd08acbc82427369e0d54355.png)
