and
longitude , using the hour angle and solar declination (where is latitude of subsolar point, and is relative longitude of subsolar point) The average annual solar radiation arriving at the top of the Earth's atmosphere is about 1361W/m2. This represents the power per unit area of solar irradiance across the spherical surface surrounding the Sun with a radius equal to the distance to the Earth (1
AU). This means that the approximately circular disc of the Earth, as viewed from the Sun, receives a roughly stable 1361W/m2 at all times. The area of this circular disc is , in which is the radius of the Earth. Because
the Earth is approximately spherical, it has total area 4 \pi r^2, meaning that the solar radiation arriving at the top of the atmosphere, averaged over the entire surface of the Earth, is simply divided by four to get 340W/m2. In other words, averaged over the year and the day, the Earth's atmosphere receives 340W/m2 from the Sun. This figure is important in
radiative forcing.
Derivation The distribution of solar radiation at the top of the atmosphere is determined by
Earth's sphericity and orbital parameters. This applies to any unidirectional beam incident to a rotating sphere. Insolation is essential for
numerical weather prediction and understanding
seasons and
climatic change. Application to
ice ages is known as
Milankovitch cycles. Distribution is based on a fundamental identity from
spherical trigonometry, the
spherical law of cosines: \cos(c) = \cos(a) \cos(b) + \sin(a) \sin(b) \cos(C) where , and are arc lengths, in radians, of the sides of a spherical triangle. is the angle in the vertex opposite the side which has arc length . Applied to the calculation of
solar zenith angle , the following applies to the spherical law of cosines: \begin{align} C &= h \\ c &= \Theta \\ a &= \tfrac{1}{2}\pi-\varphi \\ b &= \tfrac{1}{2}\pi-\delta \\ \cos(\Theta) &= \sin(\varphi) \sin(\delta) + \cos(\varphi) \cos(\delta) \cos(h) \end{align} This equation can be also derived from a more general formula: \begin{align} \cos(\Theta) = \sin(\varphi) \sin(\delta) \cos(\beta) &+ \sin(\delta) \cos(\varphi) \sin(\beta) \cos(\gamma) + \cos(\varphi) \cos(\delta) \cos(\beta) \cos(h) \\ &- \cos(\delta) \sin(\varphi) \sin(\beta) \cos(\gamma) \cos(h) - \cos(\delta) \sin(\beta) \sin(\gamma) \sin(h) \end{align} where is an angle from the horizontal and is the
solar azimuth angle. The derivation of the cosine of solar zenith angle, \cos(\Theta), based on vector analysis instead of spherical trigonometry is also available in the article about
solar azimuth angle. The separation of Earth from the Sun can be denoted and the mean distance can be denoted , approximately The
solar constant is denoted . The solar flux density (insolation) onto a plane tangent to the sphere of the Earth, but above the bulk of the atmosphere (elevation 100 km or greater) is: Q = \begin{cases} S_o \frac{R_o^2}{R_E^2}\cos(\Theta) & \cos(\Theta) > 0 \\ 0 & \cos(\Theta) \le 0 \end{cases} The average of over a day is the average of over one rotation, or the
hour angle progressing from to : \overline{Q}^\text{day} = -\frac{1}{2\pi} {\int_\pi^{-\pi}Q\,dh} Let be the hour angle when becomes positive. This could occur at sunrise when \Theta = \tfrac{1}{2} \pi, or for as a solution of \sin(\varphi) \sin(\delta) + \cos(\varphi) \cos(\delta) \cos(h_o) = 0 or \cos(h_o) = -\tan(\varphi)\tan(\delta) If , then the sun does not set and the sun is already risen at , so . If , the sun does not rise and \overline{Q}^\text{day} = 0. \frac{R_o^2}{R_E^2} is nearly constant over the course of a day, and can be taken outside the integral \begin{align} \int_\pi^{-\pi}Q\,dh &= \int_{h_o}^{-h_o}Q\,dh \\[5pt] &= S_o\frac{R_o^2}{R_E^2}\int_{h_o}^{-h_o}\cos(\Theta)\, dh \\[5pt] &= S_o\frac{R_o^2}{R_E^2}\Bigg[ h \sin(\varphi)\sin(\delta) + \cos(\varphi)\cos(\delta)\sin(h) \Bigg]_{h=h_o}^{h=-h_o} \\[5pt] &= -2 S_o\frac{R_o^2}{R_E^2}\left[ h_o \sin(\varphi) \sin(\delta) + \cos(\varphi) \cos(\delta) \sin(h_o) \right] \end{align} Therefore: \overline{Q}^{\text{day}} = \frac{S_o}{\pi}\frac{R_o^2}{R_E^2}\left[ h_o \sin(\varphi) \sin(\delta) + \cos(\varphi) \cos(\delta) \sin(h_o) \right] Let
θ be the conventional polar angle describing a planetary
orbit. Let
θ = 0 at the
March equinox. The
declination δ as a function of orbital position is \delta = \varepsilon \sin(\theta) where is the
obliquity. (Note: The correct formula, valid for any axial tilt, is \sin(\delta) = \sin(\varepsilon) \sin(\theta).) The conventional
longitude of perihelion ϖ is defined relative to the March equinox, so for the elliptical orbit: R_E = \frac{R_o(1-e^2)}{1 + e\cos(\theta - \varpi)} or \frac{R_o}{R_E} = \frac{1 + e\cos(\theta - \varpi)}{1-e^2} With knowledge of
ϖ,
ε and
e from astrodynamical calculations and So from a consensus of observations or theory, \overline{Q}^\text{day}can be calculated for any latitude
φ and
θ. Because of the elliptical orbit, and as a consequence of
Kepler's second law,
θ does not progress uniformly with time. Nevertheless,
θ = 0° is exactly the time of the March equinox,
θ = 90° is exactly the time of the June solstice,
θ = 180° is exactly the time of the September equinox and
θ = 270° is exactly the time of the December solstice. A simplified equation for irradiance on a given day is: Q \approx S_0 \left (1 + 0.034 \cos \left (2 \pi \frac{n}{365.25} \right ) \right ) where
n is a number of a day of the year.
Variation Total solar irradiance (TSI) changes slowly on decadal and longer timescales. The variation during
solar cycle 21 was about 0.1% (peak-to-peak). In contrast to older reconstructions, most recent TSI reconstructions point to an increase of only about 0.05% to 0.1% between the 17th century
Maunder Minimum and the present. However, current understanding based on various lines of evidence suggests that the lower values for the secular trend are more probable. Ultraviolet irradiance (EUV) varies by approximately 1.5 percent from solar maxima to minima, for 200 to 300 nm wavelengths. However, a proxy study estimated that UV has increased by 3.0% since the Maunder Minimum. Some variations in insolation are not due to solar changes but rather due to the Earth moving between its
perihelion and aphelion, or changes in the latitudinal distribution of radiation. These orbital changes or
Milankovitch cycles have caused radiance variations of as much as 25% (locally; global average changes are much smaller) over long periods. The most recent significant event was an axial tilt of 24° during boreal summer near the
Holocene climatic optimum. Obtaining a time series for a \overline{Q}^{\mathrm{day}} for a particular time of year, and particular latitude, is a useful application in the theory of Milankovitch cycles. For example, at the summer solstice, the declination δ is equal to the obliquity
ε. The distance from the Sun is \frac{R_o}{R_E} = 1 + e\cos(\theta - \varpi) = 1 + e\cos\left(\frac{\pi}{2} - \varpi\right) = 1 + e \sin(\varpi) For this summer solstice calculation, the role of the elliptical orbit is entirely contained within the important product e \sin(\varpi), the
precession index, whose variation dominates the variations in insolation at 65°N when eccentricity is large. For the next 100,000 years, with variations in eccentricity being relatively small, variations in obliquity dominate.
Measurement The space-based TSI record comprises measurements from more than ten radiometers and spans three solar cycles. All modern TSI satellite instruments employ
active cavity electrical substitution radiometry. This technique measures the electrical heating needed to maintain an absorptive blackened cavity in thermal equilibrium with the incident sunlight which passes through a precision
aperture of calibrated area. The aperture is modulated via a
shutter. Accuracy uncertainties of 2 per century. For various reasons, the sources do not always agree. The Solar Radiation and Climate Experiment/Total Irradiance Measurement (
SORCE/TIM) TSI values are lower than prior measurements by the Earth Radiometer Budget Experiment (ERBE) on the
Earth Radiation Budget Satellite (ERBS), VIRGO on the
Solar Heliospheric Observatory (SoHO) and the ACRIM instruments on the
Solar Maximum Mission (SMM),
Upper Atmosphere Research Satellite (UARS) and
ACRIMSAT. Pre-launch ground calibrations relied on component rather than system-level measurements since irradiance standards at the time lacked sufficient absolute accuracies. Recommendations to resolve the instrument discrepancies include validating optical measurement accuracy by comparing ground-based instruments to laboratory references, such as those at
National Institute of Science and Technology (NIST); NIST validation of aperture area calibrations uses spares from each instrument; and applying
diffraction corrections from the view-limiting aperture. and . Estimates from space-based measurements range +3–7W/m2. SORCE/TIM's lower TSI value reduces this discrepancy by 1W/m2. This difference between the new lower TIM value and earlier TSI measurements corresponds to a climate forcing of −0.8W/m2, which is comparable to the energy imbalance. == On Earth's surface ==