Extraterrestrial Irradiance

Total Solar Irradiance and the Solar Constant

Extraterrestrial irradiance refers to solar irradiance outside the earth’s atmosphere. Total Solar Irradiance (TSI) is the amount of incoming solar electromagnetic radiation per unit area at the top of the Earth’s atmosphere (TOA).   The solar constant is the mean TSI.

The following R code downloads TSI measurements across a range of experiments and satellite data collection efforts to present a composite time series of TSI spanning the period 1976 to 2013.  The objective is to confirm the numeric value of the solar constant and to define the amount of variation in solar irradiance entering the Earth’s atmosphere.

All data is sourced from the World Radiation Center at the Physikalisch-Meteorologische Observatorium Davos (PMOD/WRC).  The code uses unpublished data from the VIRGO Experiment on the cooperative ESA/NASA Mission known as SoHO.1

The data visualization shows daily TSI with an 81-day moving average.  Variations in TSI are seen to have an 11-year cycle, which corresponds to variations in sun spot activity.2  The presence of an 11-year solar cycle is sometimes used to explain climate variation or periodic seasonal departures from normal.

The data shows that cycle variation has only a minimal impact on the measured  amount of solar radiation entering the Earth’s atmosphere.  The mean irradiance at the top of the atmosphere, or the solar constant, is calculated to be 1365.88 W/m2. However, the amount of variation around the mean is well contained to 5 W/m2.  A solar cycle is present, indeed, but the amplitude is very small relative to the total irradiance entering the Earth’s atmosphere. Notwithstanding, when this small variation is magnified over the total surface area of the Earth’s disc facing the Sun, then the amount of variation is significant (5 W/m2  *  1.29e14 meters).

Sunlight and the Solar Spectrum

The solar constant includes all types of solar radiation, not just visible light.  The following code decomposes TSI into its spectral components to determine the distribution or the relative weights of different wavelengths.  The code utilizes data from the Renewable Resource Data Center (RReDC) at NREL to define the spectral distribution at the top of the Earth’s atmosphere.  Also presented are two standard reference distributions of sunlight at the Earth’s surface, one for global horizontal irradiance (GHI) and one for direct normal irradiance (DNI).  The reference distributions are found in ASTM G-173-03 and define standard test conditions for solar photovoltaic (PV) and concentrator equipment.

Plot131c

The profile of the measured solar spectrum closely resembles that of a blackbody emitter calculated previously (see Blackbody Radiation).  The Y-axis values on the blackbody chart are  slightly different than those above since the current chart normalizes the solar radiation by wavelength to present spectral intensity.  

Again we see that the peak of the solar spectrum is dominated by the visible spectrum between 400 and 700 nanometers (nm = 1 billionth of a meter = 1000 * \mu m).  We have already shown that TOA irradiance is 1,366 W/m2 using direct satellite observations, and this number is roughly confirmed by taking the integral of the AM0 sunlight curve, which yields 1,356 W/m2.[ There are two possible reasons why the WRC data and the ASTM reference curves do not agree: Either the WRC measured data is biased for some reason (e.g. sensor measurement or inconsistent data treatment across satellite experiments) or the ASTM reference curve model is biased (e.g. incorrect blackbody temperature, model error, etc.).  Either way, the discrepancy is small and represents less 1% of the WRC estimate)]. The new data also reveals that the energy distribution of sunlight entering the atmosphere is 10% ultraviolet light, 40% visible light, and 50% infrared light. 

As sunlight travels through the atmosphere, chemicals interact with the sunlight and absorb certain wavelength bands, transforming radiant energy into heat energy. The best known example of atmospheric filtering is the stripping of high-energy, ultraviolet light by ozone molecules in the upper atmosphere. A more active component of the filtering process is water vapor, which has absorption bands at numerous wavelengths, as shown by the dramatic dips in the GHI curve.

At ground level, the energy composition of sunlight decreases to roughly 1000 W/m2 using the ASTM curves. The reference curve data also reveals that the energy fractions are now 3% ultraviolet, 44% visible light, and 57% infrared.  In practice, ground level energy intensity and the solar spectrum will differ from standard test conditions based on location (e.g. latitude), elevation, time of day, and atmospheric conditions.

Solar Orbit and Extraterrestrial Irradiance

The mean distance between the Sun and Earth is 1 astronomical unit (149,597,871 km).  The distance varies as the Earth follows an elliptical orbit, moving from the perihelion in January (closest to the Sun) to the aphelion in July (farthest from the sun).3  At the average distance, light travels from the Sun to the Earth in 8 minutes and 19 seconds.    

Johannes Kepler is credited as the first to validate the heliocentric model of the Earth rotating around the Sun.4  Specifically, Kepler showed that the Earth’s orbit is ellipse with an eccentricity of 0.17.   The following diagram shows the relation of the Sun and earth with an exaggerated eccentricity for clarity.

SolarOrbitDates

 

The orbital plot shows that the four seasons are determined by perpendicular axes: the solar equinox occurs when the Earth’s tilt is at a maximum angle to the Sun, and the solar solstice occurs when the Earth’s tilt is at a minimum angle to the Sun.  The Earth is closest to the Sun at the periapsis in January and furthest from the Sun at the apoasis in July.  As a consequence, the Southern hemisphere enjoys more solar energy during its summer than the Northern hemisphere.

The distance between the Sun and Earth is defined by basic trigonometry and the function for an ellipse.  The same function determines the intensity of TOA radiation as a function of time.  The orbital distance can be estimated using the following equation:5

(1)   \begin{equation*} D = \frac{a(1-e^2)}{(1+e\cos\Theta)} \end{equation*}

where
D is the distance from the Sun;
a is the semi-major axis of the ellipse,
e is the eccentricity, and \Theta is the angle around from perihelion. 

We know a equals 1 AU or 149,598,261 km, and e equals 0.01671123. The tricky part here is linking \Theta with time (or the day of year). A simplifying approximation is to assume that the Earth orbits at a constant angular rate.6 Hence, time in days since perihelion is given by time = \Theta * 365.25/360 since there are  365.25 days in the year and it takes 360 degrees to orbit once.   \Theta is found by defining time as the number days since perihelion,  rearranging terms, and converting degrees to radians:

(2)   \begin{equation*} \Theta = model~date - perhihelion~date * \frac{360}{365.256363004} * \frac{\Pi}{180} \end{equation*}

Solving for D as a function time yields the following plot. Plot135

Extraterrestrial Irradiance

TOA irradiance will vary throughout the year with Earth’s elliptical orbit.  The variation effect can be represented empirically with the following equation:7

(3)   \begin{equation*} E_{a} = E_{SC} * \left(\frac{D_{avg}}{D}\right)^2 \end{equation*}

where
E_{SC} is the solar constant;
D_{avg} is the mean sun-earth distance; and
D is the actual sun-earth distance as a function of the day of the year.

Plot136

Optical Air Mass

AirMass_EarthImage

Air mass is a relative measure of the optical length of the atmosphere. AM0 refers to the entry point into the atmosphere, where the air mass equals zero and solar irradiance equals the solar constant.  At sea level, when the sun is directly overhead  and the zenith angle equals 0,  the air mass is equal to one or AM1.  As the zenith angle increases, the direct path of sunlight through the atmosphere grows longer. Air mass approaches its maximum value (AM38) as the zenith angle approaches 90°.  At this angle, the path length of sunlight through the atmosphere is roughly 38 times longer than when the zenith angle is zero.

Estimating the impact of atmospheric filtering on terrestrial solar irradiance is a challenge since the atmosphere is highly variable in nature.  The following sections discuss interpolative and atmospheric models of air mass, their role in solar resource assessment, and the impact on the performance of solar power systems.

B.Horn

Back | Next

Show 7 footnotes

  1. ftp://ftp.pmowrc.ch/pub/data/irradiance/composite/
  2. Willson, R. C.; Hudson, H. S. (1991). “The Sun’s Luminosity Over a Complete Solar Cycle”. Nature 351 (6321): 42–4.
  3. “Equinoxes, Solstices, Perihelion, and Aphelion, 2000–2020”. US Naval Observatory. 31 January 2008.
  4. Johannes Kepler (1609), New Astronomy Based upon Causes, or Celestial Physics, Treated by Means of Commentaries on the Motion of the Star Mars, from Observations of Tycho Brahe.
  5. Weisstein, Eric (2007) “Periapsis”, Wolfram Research.
  6. Kepler figured out that planets orbit a little bit faster when they are near the Sun than when they are away from it. As a result, relating theta to time since perihelion is a bit trickier. For more information see Weisstein, Eric (2007) “Anomalies for Planetary Motion” Wolfram Research and Weisstein, Eric (2007) “Kepler’s Equation” Wolfram Research.
  7. Stoffel et. al, Best Practices Handbook for the Collection and Use of Solar Resource Data, 2010, NREL, p.4

Leave a Reply

Your email address will not be published. Required fields are marked *

18 − twelve =