3. The Sun’s Position

In order to understand how to collect energy from the sun, one must first be able to predict the location of the sun relative to the collection device. In this chapter we develop the necessary equations by use of a unique vector approach. This approach will be used in subsequent chapters to develop the equations for the sun’s position relative to a fixed or tracking solar collector, (Chapter 4) and the special case of a sun-tracking mirror reflecting sunlight onto a fixed point (Chapter 10). Once developed, the sun position expressions of this chapter are used to demonstrate how to determine the location of shadows and the design of simple sundials. In outline form, our development looks like this:

 

Although many intermediate steps of derivation used to obtain the equations described in this chapter have been omitted, it is hoped that there are adequate comments between steps to encourage the student to perform the derivation, thereby enhancing understanding of the materials presented. Brief notes on the transformation of vector coordinates are included as Section 3.5 and a summary of sign conventions for all of the angles used in this chapter is given in Table 3.3 at the end of this chapter. Figures defining each angle and an equation to calculate it are also included.

One objective in writing this chapter has been to present adequate analytical expressions so that the solar designer is able to develop simple computer algorithms for predicting relative sun and collector positions for exact design conditions and locations. This will eliminate the need to depend on charts and tables and simplified equations.

3.1 Earth-Sun Angles

The earth revolves around the sun every 365.25 days in an elliptical orbit, with a mean earth-sun distance of 1.496 x 1011 m (92.9 x 106 miles) defined as one astronomical unit (1 AU). This plane of this orbit is called the ecliptic plane. The earth's orbit reaches a maximum distance from the sun, or aphelion, of 1.52 × 1011 m (94.4 × 106 miles) on about the third day of July. The minimum earth-sun distance, the perihelion, occurs on about January 2nd, when the earth is 1.47 × 1011 m (91.3 × 106 miles) from the sun.  Figure 3.1 depicts these variations in relation to the Northern Hemisphere seasons.

Figure 3.1 The ecliptic plane showing variations in the earth-sun distance and the equinoxes and solstices. The dates and day numbers shown are for 1981 and may vary by 1 or 2 days.

The earth rotates about its own polar axis, inclined to the ecliptic plane by 23.45 degrees, in approximately 24-hour cycles. The direction in which the polar axis points is fixed in space and is aligned with the North Star (Polaris) to within about 45 minutes of arc (13 mrad). The earth’s rotation about its polar axis produces our days and nights; the tilt of this axis relative to the ecliptic plane produces our seasons as the earth revolves about the sun.

3.1.1 Time

We measure the passage of time by measuring the rotation of the earth about its axis. The base for time (and longitude) measurement is the meridian that passes through Greenwich, England and both poles. It is known as the Prime Meridian.

Today, the primary world time scale, Universal Time (previously called Greenwich Mean Time), is still measured at the Prime Meridian. This is a 24-hour time system, based on mean time, according to which the length of a day is 24 hours and midnight is 0 hours.

Mean time is based on the length of an average day. A mean second is l/86,400 of the average time between one complete transit of the sun, averaged over the entire year. In fact, the length of any one specific day, measured by the complete transit of the sun, can vary by up to 30 seconds during the year. The variable day length is due to four factors listed in order of decreasing importance (Jesperson and Fitz-Randolph, 1977):

Standard Time Zones - Since it is conventional to have 12:00 noon be approximately in the middle of the day regardless of the longitude, a system of time zones has been developed. See (Blaise 2000) for an interesting story about how this unification developed. These are geographic regions, approximately 15 degrees of longitude wide, centered about a meridian along which local standard time equals mean solar time. Prior to about 1880, different cities (and even train stations) had their own time standards, most based on the sun being due south at 12:00 noon.

Time is now generally measured about standard time zone meridians. These meridians are located every 15 degrees from the Prime Meridian so that local time changes in 1-hour increments from one standard time zone meridian to the next. The standard time zone meridians east of Greenwich have times later than Greenwich time, and the meridians to the west have earlier times.

Ideally, the meridians 7 degrees on either side of the standard time zone meridian should define the time zone. However, boundaries separating time zones are not meridians but politically determined borders following rivers, county, state or national boundaries, or just arbitrary paths. Countries such as Spain choose to be on `European Time (15o E) when their longitudes are well within the adjacent Standard Time Meridian (0o). Figure 3.2 shows these time zone boundaries within the United States and gives the standard time zone meridians (called longitudes of solar time).

Figure 3.2 Time zone boundaries within the United States. From Jesperson and Fitz-Randolph 1977.

Daylight Savings Time -To complicate matters further in trying to correlate clock time with the movement of the sun, a concept known as daylight savings time was initiated in the United States in the spring of 1918 to "save fuel and promote other economies in a country at war" (Jesperson and Fitz-Randolph, 1977). According to this concept, the standard time is advanced by 1 hour, usually from 2:00 AM on the first Sunday in April until 2:00 AM on the last Sunday in October. Although various attempts have been made to apply this concept uniformly within the country, it is suggested that the designer check locally to ascertain the commitment to this concept at any specific solar site.

Sidereal Time - So far, and for the remainder of this text, all reference to time is to mean time, a time system based on the assumption that a day (86,400 seconds) is the average interval between two successive times when a given point on the earth faces the sun. In astronomy or orbital mechanics, however, the concept of sidereal time is often used. This time system is based on the sidereal day, which is the length of time for the earth to make one complete rotation about its axis.

The mean day is about 4 minutes longer than the sidereal day because the earth, during the time it is making one revolution about its axis, has moved some distance in its orbit around the sun. To be exact, the sidereal day contains 23 hours, 56 minutes, and 4.09053 seconds of mean time. Since, by definition, there are 86,400 sidereal seconds in a sidereal day, the sidereal second is slightly shorter than the mean solar second is. To be specific: 1 mean second = 1.002737909 sidereal seconds. A detailed discussion of this and other time definitions is contained in another work (Anonymous, 1981, Section B).

3.1.2 The Hour Angle  

To describe the earth's rotation about its polar axis, we use the concept of the hour angle . As shown in Figure 3.3, the hour angle is the angular distance between the meridian of the observer and the meridian whose plane contains the sun. The hour angle is zero at solar noon (when the sun reaches its highest point in the sky). At this time the sun is said to be ‘due south’ (or ‘due north’, in the Southern Hemisphere) since the meridian plane of the observer contains the sun. The hour angle increases by 15 degrees every hour.

 

Figure 3.3 The hour angle . This angle is defined as the angle between the meridian parallel to sun rays and the meridian containing the observer.

Solar Time - Solar time is based on the 24-hour clock, with 12:00 as the time that the sun is exactly due south. The concept of solar time is used in predicting the direction of sunrays relative to a point on the earth. Solar time is location (longitude) dependent and is generally different from local clock time, which is defined by politically defined time zones and other approximations.

Solar time is used extensively in this text to define the rotation of the earth relative to the sun. An expression to calculate the hour angle from solar time is

                         (3.1)

where ts is the solar time in hours.

EXAMPLE: When it is 3 hours after solar noon, solar time is 15:00 and the hour angle has a value of 45 degrees. When it is 2 hours and 20 minutes before solar noon, solar time is 9:40 and the hour angle is 325 degrees (or 35 degrees).

The difference between solar time and local clock time can approach 2 hours at various locations and times in the United States, For most solar design purposes, clock time is of little concern, and it is appropriate to present data in terms of solar time. Some situations, however, such as energy demand correlations, system performance correlations, determination of true south, and tracking algorithms require an accurate knowledge of the difference between solar time and the local clock time.

Knowledge of solar time and Universal Time has traditionally been important to ship navigators. They would set their chronometers to an accurately adjusted tower clock visible as they left port. This was crucial for accurate navigation. At sea a ship's latitude could be easily ascertained by determining the maximum altitude angle of the sun or the altitude angle of Polaris at night. However, determining the ship’s longitude was more difficult and required that an accurate clock be carried onboard. If the correct time at Greenwich, England (or any other known location) was known, then the longitude of the ship could be found by measuring the solar time onboard the ship (through sun sightings) and subtracting from it the time at Greenwich. Since the earth rotates through 360 degrees of longitude every 24 hours, the ship then has traveled 1 degree of longitude away from the Prime Meridian (which passes through Greenwich) for every 4 minutes of time difference. An interesting story about developing accurate longitude measurements may be found in Sobel, 1999.

Equation of Time - The difference between mean solar time and true solar time on a given date is shown in Figure 3.4. This difference is called the equation of time (EOT). Since solar time is based on the sun being due south at 12:00 noon on any specific day, the accumulated difference between mean solar time and true solar time can approach 17 minutes either ahead of or behind the mean, with an annual cycle.

The level of accuracy required in determining the equation of time will depend on whether the designer is doing system performance or developing tracking equations. An approximation for calculating the equation of time in minutes is given by Woolf (1968) and is accurate to within about 30 seconds during daylight hours.

                    (3.2)

where the angle x is defined as a function of the day number N

 

                     (3.3)

with the day number, N being the number of days since January 1. Table 3.1 has been prepared as an aid in rapid determination of values of N from calendar dates.

Table 3.1 Date-to-Day Number Conversion

Month

Day Number, N

Notes

January

d

February

d + 31

March

d + 59

Add 1 if leap year

April

d + 90

Add 1 if leap year

May

d + 120

Add 1 if leap year

June

d + 151

Add 1 if leap year

July

d + 181

Add 1 if leap year

August

d + 212

Add 1 if leap year

September

d + 243

Add 1 if leap year

October

d + 273

Add 1 if leap year

November

d + 304

Add 1 if leap year

December

d + 334

Add 1 if leap year

Days of Special Solar Interest

Solar Event

Date

Day Number, N

Vernal equinox

March 21

80

Summer solstice

June 21

172

Autumnal equinox

September 23

266

Winter solstice

December 21

355

NOTES:

  1. d is the day of the month
  2. Leap years are 2000, 2004, 2008 etc.
  3. Solstice and equinox dates may vary by a day or two. Also, add 1 to the solstice and equinox day number for leap years.

Figure 3.4 The equation of time (EOT). This is the difference between the local apparent solar time and the local mean solar time.

EXAMPLE: February 11 is the 42nd day of the year, therefore N = 42 and x is equal to 40.41 degrees, and the equation of time as calculated above is 14.35 minutes. This compares with a very accurately calculated value of -14.29 minutes reported elsewhere (Anonymous, 198l). This means that on this date, there is a difference between the mean time and the solar time of a little over 14 minutes or that the sun is "slow" relative to the clock by that amount.

To satisfy the control needs of concentrating collectors, a more accurate determination of the hour angle is often needed. An approximation of the equation of time claimed to have an average error of 0.63 seconds and a maximum absolute error of 2.0 seconds is presented below as Equation (3.4) taken from Lamm (1981). The resulting value is in minutes and is positive when the apparent solar time is ahead of mean solar time and negative when the apparent solar time is behind the mean solar time:

                     (3.4)

Here n is the number of days into a leap year cycle with n = 1 being January 1 of each leap year, and n =1461 corresponding to December 31 of the 4th year of the leap year cycle. The coefficients Ak and Bk are given in Table 3.2 below. Arguments for the cosine and sine functions are in degrees. 

Table 3.2 Coefficients for Equation (3.4)

k

Ak (hr)

Bk (hr)

0

2.0870 × 10-4

0

1

9.2869 × 10-3

1.2229 ´ 10-1

2

5.2258 ´ 10-2

1.5698 ´ 10-1

3

1.3077 ´ l0-3

5.1 602 ´ 10-3

4

2.1867 ´ l0-3

2.9823 ´ 10-3

5

1.5100 ´ 10-4

2.3463 ´ 10-4

Time Conversion - The conversion between solar time and clock time requires knowledge of the location, the day of the year, and the local standards to which local clocks are set. Conversion between solar time, ts and local clock time (LCT) (in 24-hour rather than AM/ PM format) takes the form

                     (3.5)

where EOT is the equation of time in minutes and LC is a longitude correction defined as follows:

                     (3.6)

and the parameter D in Equation (3.5) is equal to 1 (hour) if the location is in a region where daylight savings time is currently in effect, or zero otherwise.

EXAMPLE: Let us find the clock time for solar noon at a location in Los Angeles, having a longitude of 118.3 degrees on February 11. Since Los Angeles is on Pacific Standard Time and not on daylight savings time on this date, the local clock time will be:

 

3.1.3 The Declination Angle  

The plane that includes the earth 's equator is called the equatorial plane. If a line is drawn between the center of the earth and the sun, the angle between this line and the earth's equatorial plane is called the declination angle , as depicted in Figure 3.5. At the time of year when the northern part of the earth's rotational axis is inclined toward the sun, the earth 's equatorial plane is inclined 23.45 degrees to the earth-sun line. At this time (about June 21), we observe that the noontime sun is at its highest point in the sky and the declination angle = +23.45 degrees. We call this condition the summer solstice, and it marks the beginning of summer in the Northern Hemisphere.

As the earth continues its yearly orbit about the sun, a point is reached about 3 months later where a line from the earth to the sun lies on the equatorial plane. At this point an observer on the equator would observe that the sun was directly overhead at noontime. This condition is called an equinox since anywhere on the earth, the time during which the sun is visible (daytime) is exactly 12 hours and the time when it is not visible (nighttime) is 12 hours. There are two such conditions during a year; the autumnal equinox on about September 23, marking the start of the fall; and the vernal equinox on about March 22, marking the beginning of spring. At the equinoxes, the declination angle is zero.

Figure 3.5 The declination angle . The earth is shown in the summer solstice position when = +23.45 degrees. Note the definition of the tropics as the intersection of the earth-sun line with the surface of the earth at the solstices and the definition of the Arctic and Antarctic circles by extreme parallel sun rays.

The winter solstice occurs on about December 22 and marks the point where the equatorial plane is tilted relative to the earth-sun line such that the northern hemisphere is tilted away from the sun. We say that the noontime sun is at its "lowest point" in the sky, meaning that the declination angle is at its most negative value (i.e., = -23.45 degrees). By convention, winter declination angles are negative.

Accurate knowledge of the declination angle is important in navigation and astronomy. Very accurate values are published annually in tabulated form in an ephemeris; an example being (Anonymous, 198l). For most solar design purposes, however, an approximation accurate to within about 1 degree is adequate. One such approximation for the declination angle is

                     (3.7)

where the argument of the cosine here is in degrees and N is the day number defined for Equation (3.3) The annual variation of the declination angle is shown in Figure 3.5.

3.1.4 Latitude Angle  

The latitude angle is the angle between a line drawn from a point on the earth’s surface to the center of the earth, and the earth’s equatorial plane. The intersection of the equatorial plane with the surface of the earth forms the equator and is designated as 0 degrees latitude. The earth’s axis of rotation intersects the earth’s surface at 90 degrees latitude (North Pole) and -90 degrees latitude (South Pole). Any location on the surface of the earth then can be defined by the intersection of a longitude angle and a latitude angle.

Other latitude angles of interest are the Tropic of Cancer (+23.45 degrees latitude) and the Tropic of Capricorn (- 23.45 degrees latitude). These represent the maximum tilts of the north and south poles toward the sun. The other two latitudes of interest are the Arctic circle (66.55 degrees latitude) and Antarctic circle (-66.5 degrees latitude) representing the intersection of a perpendicular to the earth-sun line when the south and north poles are at their maximum tilts toward the sun. As will be seen below, the tropics represent the highest latitudes where the sun is directly overhead at solar noon, and the Arctic and Antarctic circles, the lowest latitudes where there are 24 hours of daylight or darkness. All of these events occur at either the summer or winter solstices.

3.2 Observer-Sun Angles

When we observe the sun from an arbitrary position on the earth, we are interested in defining the sun position relative to a coordinate system based at the point of observation, not at the center of the earth. The conventional earth-surface based coordinates are a vertical line (straight up) and a horizontal plane containing a north-south line and an east-west line. The position of the sun relative to these coordinates can be described by two angles; the solar altitude angle and the solar zenith angle defined below. Since the sun appears not as a point in the sky, but as a disc of finite size, all angles discussed in the following sections are measured to the center of that disc, that is, relative to the "central ray" from the sun.

3.2.1 Solar Altitude, Zenith, and Azimuth Angles

The solar altitude angle is defined as the angle between the central ray from the sun, and a horizontal plane containing the observer, as shown in Figure 3.6. As an alternative, the sun’s altitude may be described in terms of the solar zenith angle which is simply the complement of the solar altitude angle or

                     (3.8)

The other angle defining the position of the sun is the solar azimuth angle (A). It is the angle, measured clockwise on the horizontal plane, from the north-pointing coordinate axis to the projection of the sun’s central ray.

Figure 3.6 Earth surface coordinate system for observer at Q showing the solar azimuth angle Α, the solar altitude angle and the solar zenith angle for a central sun ray along direction vector S. Also shown are unit vectors i, j, k along their respective axes.

The reader should be warned that there are other conventions for the solar azimuth angle in use in the solar literature. One of the more common conventions is to measure the azimuth angle from the south-pointing coordinate rather than from the north-pointing coordinate. Another is to consider the counterclockwise direction positive rather than clockwise. The information in Table 3.3 at the end of this chapter will be an aid in recognizing these differences when necessary.

It is of the greatest importance in solar energy systems design, to be able to calculate the solar altitude and azimuth angles at any time for any location on the earth. This can be done using the three angles defined in Section 3.1 above; latitude , hour angle , and declination . If the reader is not interested in the details of this derivation, they are invited to skip directly to the results; Equations (3.17), (3.18) and (3.19).

For this derivation, we will define a sun-pointing vector at the surface of the earth and then mathematically translate it to the center of the earth with a different coordinate system. Using Figure 3.6 as a guide, define a unit direction vector S pointing toward the sun from the observer location Q:

                     (3.9)

 

where i, j, and k are unit vectors along the z, e, and n axes respectively. The direction cosines of S relative to the z, e, and n axes are Sz, Se and Sn, respectively. These may be written in terms of solar altitude and azimuth as

                     (3.10)

Similarly, a direction vector pointing to the sun can be described at the center of the earth as shown in

Figure 3.7. If the origin of a new set of coordinates is defined at the earth’s center, the m axis can be a line from the origin intersecting the equator at the point where the meridian of the observer at Q crosses. The e axis is perpendicular to the m axis and is also in the equatorial plane. The third orthogonal axis p may then be aligned with the earth’s axis of rotation. A new direction vector pointing to the sun may be described in terms of its direction cosines m , S'e and S'p relative to the m, e, and p axes, respectively. Writing these in terms of the declination and hour angles, we have

                     (3.11)