Idealized greenhouse model

The temperatures of a planet's surface and atmosphere are governed by a delicate balancing of their energy flows. The idealized greenhouse model is based on the fact that certain gases in the Earth's atmosphere, including carbon dioxide and water vapour, are transparent to the high-frequency solar radiation, but are much more opaque to the lower frequency infrared radiation leaving Earth's surface. Thus heat is easily let in, but is partially trapped by these gases as it tries to leave. Rather than get hotter and hotter, Kirchhoff's law of thermal radiation says that the gases of the atmosphere also have to re-emit the infrared energy that they absorb, and they do so, also at long infrared wavelengths, both upwards into space as well as downwards back towards the Earth's surface. In the long-term, the planet's thermal inertia is surmounted and a new thermal equilibrium is reached when all energy arriving on the planet is leaving again at the same rate. In this steady-state model, the greenhouse gases cause the surface of the planet to be warmer than it would be without them, in order for a balanced amount of heat energy to finally be radiated out into space from the top of the atmosphere.

Essential features of this model where first published by Svante Arrhenius in 1896. It has since become a common introductory "textbook model" of the radiative heat transfer physics underlying Earth's energy balance and the greenhouse effect. The planet is idealized by the model as being functionally "layered" with regard to a sequence of simplified energy flows, but dimensionless (i.e. a zero-dimensional model) in terms of its mathematical space. The layers include a surface with constant temperature Ts and an atmospheric layer with constant temperature Ta. For diagrammatic clarity, a gap can be depicted between the atmosphere and the surface. Alternatively, Ts could be interpreted as a temperature representative of the surface and the lower atmosphere, and Ta could be interpreted as the temperature of the upper atmosphere, also called the skin temperature. In order to justify that Ta and Ts remain constant over the planet, strong oceanic and atmospheric currents can be imagined to provide plentiful lateral mixing. Furthermore, the temperatures are understood to be multi-decadal averages such that any daily or seasonal cycles are insignificant.

Simplified energy flows
The model will find the values of Ts and Ta that will allow the outgoing radiative power, escaping the top of the atmosphere, to be equal to the absorbed radiative power of sunlight. When applied to a planet like Earth, the outgoing radiation will be longwave and the sunlight will be shortwave. These two streams of radiation will have distinct emission and absorption characteristics. In the idealized model, we assume the atmosphere is completely transparent to sunlight. The planetary albedo αP is the fraction of the incoming solar flux that is reflected back to space (since the atmosphere is assumed totally transparent to solar radiation, it does not matter whether this albedo is imagined to be caused by reflection at the surface of the planet or at the top of the atmosphere or a mixture). The flux density of the incoming solar radiation is specified by the solar constant S0. For application to planet Earth, appropriate values are S0=1366 W m−2 and αP=0.30. Accounting for the fact that the surface area of a sphere is 4 times the area of its intercept (its shadow), the average incoming radiation is S0/4.

For longwave radiation, the surface of the Earth is assumed to have an emissivity of 1 (i.e. it is a black body in the infrared, which is realistic). The surface emits a radiative flux density F according to the Stefan–Boltzmann law:



F=\sigma T^4 $$

where σ is the Stefan–Boltzmann constant. A key to understanding the greenhouse effect is Kirchhoff's law of thermal radiation. At any given wavelength the absorptivity of the atmosphere will be equal to the emissivity. Radiation from the surface could be in a slightly different portion of the infrared spectrum than the radiation emitted by the atmosphere. The model assumes that the average emissivity (absorptivity) is identical for either of these streams of infrared radiation, as they interact with the atmosphere. Thus, for longwave radiation, one symbol ε denotes both the emissivity and absorptivity of the atmosphere, for any stream of infrared radiation.



The infrared flux density out of the top of the atmosphere is computed as:

F\uparrow =\epsilon \sigma T_a^4 + (1-\epsilon) \sigma T_s^4 $$ In the last term, ε represents the fraction of upward longwave radiation from the surface that is absorbed, the absorptivity of the atmosphere. The remaining fraction (1-ε) is transmitted to space through an atmospheric window. In the first term on the right, ε is the emissivity of the atmosphere, the adjustment of the Stefan–Boltzmann law to account for the fact that the atmosphere is not optically thick. Thus ε plays the role of neatly blending, or averaging, the two streams of radiation in the calculation of the outward flux density.

The energy balance solution
Zero net radiation leaving the top of the atmosphere requires:

-\frac{1}{4}S_0(1-\alpha_p)+\epsilon \sigma T_a^4 + (1-\epsilon) \sigma T_s^4= 0 $$ Zero net radiation entering the surface requires:

\frac{1}{4}S_0(1-\alpha_p)+\epsilon \sigma T_a^4 - \sigma T_s^4 = 0 $$ Energy equilibrium of the atmosphere can be either derived from the two above equilibrium conditions, or independently deduced:

2 \epsilon \sigma T_a^4 - \epsilon \sigma T_s^4 = 0 $$ Note the important factor of 2, resulting from the fact that the atmosphere radiates both upward and downward. Thus the ratio of Ta to Ts is independent of ε:
 * $$ T_a = { T_s \over 2^{1/4} }

= { T_s \over 1.189 } $$ Thus Ta can be expressed in terms of Ts, and a solution is obtained for Ts in terms of the model input parameters:

\frac{1}{4}S_0(1-\alpha_p)=\left( 1-\frac{\epsilon}{2} \right) \sigma T_s^4 $$ or

T_s=\left[ \frac{S_0(1-\alpha_p)}{4\sigma} \frac{1}{1-{\epsilon \over 2}} \right]^{1/4} $$ The solution can also be expressed in terms of the effective emission temperature Te, which is the temperature that characterizes the outgoing infrared flux density F, as if the radiator were a perfect radiator obeying F=σTe4. This is easy to conceptualize in the context of the model. Te is also the solution for Ts, for the case of ε=0, or no atmosphere:

T_e \equiv \left[ \frac{S_0(1-\alpha_p)}{4\sigma} \right]^{1/4} $$ With the definition of Te:

T_s= T_e \left[ \frac{1}{1-{\epsilon \over 2}} \right]^{1/4} $$ For a perfect greenhouse, with no radiation escaping from the surface, or ε=1:

T_s= T_e 2^{1/4} = 1.189 T_e \qquad T_a=T_e $$

Application to Earth
Using the parameters defined above to be appropriate for Earth,
 * $$ T_e = 255 ~\mathrm{K} = -18 ~\mathrm{C} $$

For ε=1:
 * $$ T_s = 303 ~\mathrm{K} = 30 ~\mathrm{C} $$

For ε=0.78,
 * $$ T_s = 288.3 ~\mathrm{K} \qquad T_a = 242.5 ~\mathrm{K} $$.

This value of Ts happens to be close to the published 287.2 K of the average global "surface temperature" based on measurements. ε=0.78 implies 22% of the surface radiation escapes directly to space, consistent with the statement of 15% to 30% escaping in the greenhouse effect.

The radiative forcing for doubling carbon dioxide is 3.71 W m−2, in a simple parameterization. This is also the value endorsed by the IPCC. From the equation for $$F\uparrow$$,
 * $$ \Delta F\uparrow = \Delta\epsilon \left( \sigma T_a^4 -\sigma T_s^4 \right) $$

Using the values of Ts and Ta for ε=0.78 allows for $$ \Delta F\uparrow $$= -3.71 W m−2 with Δε=.019. Thus a change of ε from 0.78 to 0.80 is consistent with the radiative forcing from a doubling of carbon dioxide. For ε=0.80,
 * $$ T_s = 289.5 ~\mathrm{K} $$

Thus this model predicts a global warming of ΔTs = 1.2 K for a doubling of carbon dioxide. A typical prediction from a GCM is 3 K surface warming, primarily because the GCM allows for positive feedback, notably from increased water vapor. A simple surrogate for including this feedback process is to posit an additional increase of Δε=.02, for a total Δε=.04, to approximate the effect of the increase in water vapor that would be associated with an increase in temperature. This idealized model then predicts a global warming of ΔTs = 2.4 K for a doubling of carbon dioxide, roughly consistent with the IPCC.

Extensions
The one-level atmospheric model can be readily extended to a multiple-layer atmosphere. In this case the equations for the temperatures become a series of coupled equations. These simple energy-balance models always predict a decreasing temperature away from the surface, and all levels increase in temperature as "greenhouse gases are added". Neither of these effects are fully realistic: in the real atmosphere temperatures increase above the tropopause, and temperatures in that layer are predicted (and observed) to decrease as GHG's are added. This is directly related to the non-greyness of the real atmosphere.

An interactive version of a model with 2 atmospheric layers, and which accounts for convection, is available online.