Analysis: Global Average Longwave Radiation Emissions

Note: This page is highly mathematical.

Longwave (LW) thermal radiation emissions from a planet’s surface play an important role in planetary heat flows and the global energy balance. The overall energy balance operates at the level of global totals or, equivalently, global-average values.

With that in mind, let us obtain useful formulas for the average of the radiative flux of both surface longwave radiation emissions (\SLR) and outgoing longwave radiation emissions (\OLR).

Average Surface Longwave Radiation Emissions (SLR)

The Stefan-Boltzmann (SB) law tells us that the surface temperature, \tsurf, and the radiative flux of surface LW thermal radiation emissions, \slr, are related by:

(1)   \begin{equation*}  \SLR = \emis \sigma \tsurf^4 \end{equation*}

where \sigma is the Stefan-Boltzmann constant and \emis is the emissivity, a value between 0 and 1, typically fairly close to 1, which depends the the characteristics of the material doing the emitting.

The surface temperature varies in time, and also depending on where one is on the planet. Emissivity varies also, mostly by location, but sometimes over time.

Climate data is usually averaged over some time window. For purposes of looking at the global energy balance, it is also useful to average over location, i.e, latitude and longitude.

I will denote the average of a quantity X over a time window as \ex{X}_t, and the average of X over latitude and longitude as \ex{X}_\loc. (Formulas for calculating those averages are available here.) I will denote that result of averaging over both a time window and latitude and longitude as \ex{X}. You can calculate \ex{X} by averaging over time and longitude and latitude in either order, and you’ll end up with the same result. In other words, .

(2)   \begin{equation*}  \ex{X} = \ex{\ex{X}_t}_\loc = \ex{\ex{X}_\loc}_t \end{equation*}

If two quantities are equal, then the average values of those quantities will also be equal. So, we can average equation the SB equation to obtain:

(3)   \begin{equation*}  \ex{\SLR} = \sigma \ex{\emis \tsurf^4} \end{equation*}

Because \sigma is a constant, it’s okay to bring it outside of the average. But because both \tsurf and \emis vary, they can’t be average separately.

I like to see how things look when expressed using the concept of a weighted average. I’ll use the notation \exw{X}{W} to denote the average of X, weighted by W. This is calculated as:

(4)   \begin{equation*} \exw{X}{W} = \frac{\ex{X W}}{\ex{W}} \end{equation*}

This allows \ex{X W} to be expressed as {\ex{W}}\exw{X}{W}. Sometimes I find this helpful in thinking about what an expression means. But, your mileage may vary on this issue. The formulation with or without weighted averages is equally valid.

The consequence of this definition is that \ex{X}{W} still provides information about “typical” values of X, but it treats values of X as more important when W is large and less important when W is small.

Given that definition, we can rewrite the equation for \ex{\SLR} as:

(5)   \begin{equation*}  \ex{\SLR} = \sigma \exw{\emis}{\tsurf^4} \ex{\tsurf^4} =  \sigma \glob{\emis} \ex{\tsurf^4}  \end{equation*}

where the global emissivity, \glob{\emis}, or standardly-weighted average emissivity, \avg{\emis}, is defined as

(6)   \begin{equation*} \glob{\emis} = \avg{\emis} = \exw{\emis}{\tsurf^4}  \end{equation*}

Because \glob{\emis} or \avg{\emis} is a weighted average of \emis, it is constrained to be between between the minimum and maximum values of \emis.

Alternatively, this could be expressed as:

(7)   \begin{equation*} \glob{\emis} = \avg{\emis} = \gammaemis \cdot \ex{\emis} \end{equation*}

where \gammaemis is defined as

(8)   \begin{equation*} \gammaemis = \frac{\exw{\emis}{\tsurf^4}}{\ex{\emis}} = \frac{\ex{\emis \tsurf^4}}{\ex{\emis} \ex{\tsurf^4}} \end{equation*}

The factor \gammaemis may be thought of as the temperature-variation emissivity adjustment factor. If \emis varies only within a small range, then \gammaemis will necessarily be close to 1.

It is convenient to be able to express the relationship between average surface emissions,\ex{\slr}, and average surface temperature, \ex{\tsurf}. We can do that by defining the quantity, \tveb:

(9)   \begin{equation*} \tveb = \frac{\ex{\tsurf^4}}{\ex{\tsurf}^4} \end{equation*}

Then, the equation for \slr becomes:

(10)   \begin{equation*}  \ex{\SLR} = \tveb \cdot \glob{\emis} \cdot \sigma \ex{\tsurf}^4  = \tveb \gammaemis \cdot \ex{\emis}  \cdot \sigma \ex{\tsurf}^4 \end{equation*}

What is interesting to notice about this is that, in expressing the relationship between emissions and temperature in terms of the average temperature, \ex{\tsurf}, it is as if we have had to replace the emissivity with the emissivity multiplied by \tveb. So, \tveb acts as a factor to effectively modify the emissivity. Does it increase or decrease the emissivity?

It can be proven mathematically (via an application of Holder’s inequality) that it is always true that:

(11)   \begin{equation*} \tveb \ge 1 \end{equation*}

Thus, it’s clear that \tveb effectively increases emissivity. So, I’ll refer to \tveb as the temperature-variation emissivity boost factor.

Note that while the “adjustment factor” \gammaemis will never adjust the emissivity to be outside the range 0 to 1, the “boost factor” \tveb could potentially lead to an effective emissivity, \tveb\cdot\glob{\emis} greater than 1.

What we have achieved here is a transformation of the Stefan-Boltzmann equation from an equation that applies at an individual location and time to a similar equation which applies to suitably averaged values, \ex{\tsurf}, \avg{\emis}, and \ex{\SLR}. In the process, it has been necessary to introduce a new parameter, \ex{\tveb}, which effectively boosts the emissivity. This reflects the fact that a system in which the temperature varies from the average value emits more than a system with a uniform temperature. This is true whether those variations occur over time or between locations around the globe.

Average Outgoing Longwave Radiation Emissions (OLR)

The Greenhouse effect is defined quantitatively as:

(12)   \begin{equation*} \GHE = \SLR - \OLR \end{equation*}

I can also be useful to work with the dimensionless normalized Greenhouse effect:

(13)   \begin{equation*} \nghe = \lea = \frac{\GHE}{\SLR} \end{equation*}

The symbol \nghe is used for this quantity in IPCC reports, but for simplicity, I typically use the symbol \lea.

The quantity \lea can be thought of either as the normalized Greenhouse effect or as the longwave effective absorptance (LEA) of the atmosphere. The standard definition of absorptance is the fraction of radiation power incident on an absorbing medium that is absorbed. In comparison, the LEA, \lea, is the the fraction of longwave thermal radiation power emitted by the surface that, as viewed from space, is absorbed by the atmosphere without having been replaced by an equivalent amount of power from atmospheric longwave emissions. So, LEA is like absorptance, but it accounts for both absorption and emissions.

Note that 0 \leq \lea < 1. For a planet with a longwave-transparent atmosphere (with nothing in it to absorb, emit, reflect, or scatter LW radiation) the value of \lea would necessarily have 0.

Given the definition of \lea, it follows that:

(14)   \begin{equation*} \OLR = (1-\lea) \cdot \SLR \end{equation*}

(15)   \begin{equation*} \GHE = \lea \cdot \SLR \end{equation*}

Taking the average of the latter equation:

(16)   \begin{equation*} \ex{\GHE} = \ex{\lea \cdot \SLR} = \frac{\ex{\lea \cdot \SLR}}{\ex{\SLR}} \ex{\SLR} \end{equation*}

This can be written as:

(17)   \begin{equation*} \ex{\GHE} =\glob{\lea} \cdot \ex{\SLR} = \gammalea \cdot \ex{\lea} \cdot  \ex{\SLR} \end{equation*}

where the global normalized GHE or global LEA, \glob{\lea}, the standard-weighted average \lea, \avg{\lea}, and emission-variability LEA adjustment factor, \gammalea, are defined as:

(18)   \begin{equation*} \glob{\lea} = \avg{\lea} = \exw{\lea}{\SLR} =  \gammalea \cdot \ex{\lea} \end{equation*}

(19)   \begin{equation*} \gammalea = \frac{\exw{\lea}{\SLR}}{\ex{\lea}} = \frac{\ex{\lea\cdot\SLR}}{\ex{\lea}\cdot\ex{\SLR}} \end{equation*}

Because \glob{\lea}=\avg{\lea} is a weighted average of \lea, it constrained, to be between the minimum and maximum values of \lea. If \lea varies only within a small range, then \gammalea will necessarily be close to 1.

Combining our equation for \ex{OLR} with the prior equation for \ex{SLR} leads to:

(20)   \begin{equation*}  \ex{\OLR} = \glob{\lea} \cdot \tveb \cdot \glob{\emis} \cdot \sigma \ex{\tsurf}^4  = \tveb \gammalea  \gammaemis \cdot \ex{\lea} \cdot \ex{\emis}  \cdot \sigma \ex{\tsurf}^4 \end{equation*}

Conclusion

This analysis provides formulas for the global average radiative flux of longwave radiation thermal emissions from the surface, \ex{\SLR}, and the global average radiative flux of longwave radiation thermal emissions outgoing from the top of the atmosphere, \ex{\OLR}

These are exact results, involving no approximations or assumptions.

The results rely only on the definitions of the terms involved. Beyond those definitions, the resulting formulas are the exact result of mathematical transformations and tautologies.

The well-established Stefan-Boltzmann law for thermal emissions is the only principle of physics that has been used. However, the SB law can be considered to be simply a definition of the emissivity, \emis. So, in that sense, this analysis has been purely mathematical, with full mathematical rigor.