Thursday, May 27, 2010

Water Vapor Continuum

It may be that we can determine the absorption length of long-wave radiation in dry air with 330 ppm CO2 simply by dividing the absorption length pure CO2 by 0.00033. The absorption length of 15-μm radiation in pure CO2 is 10 mm (see CO2 Spectrum), so its absorption length in 330 ppm by volume CO2 would, according to this calculation, be 30 m.

The absorption spectrum of water vapor is not so simple. We have already presented our readers with the absorption spectrum of liquid water. But we cannot simply divide this spectrum by the fraction of water present in air and so obtain the spectrum for water vapor. Water vapor is water dissolved in air. The spectrum of liquid water applies to the Earth's oceans, and to the tiny droplets of water in clouds, but it does not apply to water vapor, as we saw in Vapor Versus Liquid. Water vapor absorbs several orders of magnitude less long-wave radiation than we would expect from dividing the liquid water spectrum by the water vapor concentration.

In liquid water, there are hydrogen bonds between each water molecule and its neighbors. These bonds absorb long-wave radiation. Hydrogen bonds are far less common in water vapor because the molecules are separated. Hydrogen bonds occur in water vapor when pairs of molecules join together to form dimers. The hydrogen bonds in these dimers absorb long-wave radiation. The chance of a water molecule finding a partner to form a dimer increases with concentration. At the same time, the number of water molecules per unit volume increases with concentration. The result is that the absorbing power of water vapor increases as the square of the concentration. Absorption by water vapor is also temperature-dependent because dimers are broken apart by kinetic energy, so colder dimers live longer and have more opportunity to absorb long-wave radiation.

There are also weak long-wave absorbing bonds between water and nitrogen, and water and oxygen. The absorption of long-wave radiation by water vapor is an on-going subject of research at many institutes. The relationship between molecules in water vapor is called the water vapor continuum. In Water vapor absorption coefficients in the 8−13-μm spectral region: a critical review (1990), Grant et al review existing measurements of the water vapor continuum spectrum. In Recent Developments in the Water Vapor Continuum (1999), Mlawer et al review existing calculations of the spectrum. The following graph of water vapor continuum absorption agrees well with measurements such as these.



Instead of plotting the spectrum versus wavelength, this graph plots the spectrum versus the inverse of the wavelength, in units of cm−1. The spectrum itself takes the form of a quantity in units of cm2 molecule−1 atm−1, plotted for various temperatures between 220 K and 330 K. In our next post, we will convert this graph into one of absorption length versus wavelength, and show how we can calculate the absorption of long-wave radiation by water vapor in the atmosphere.

PS. Thanks to Peter Newnam for his help researching this post.

UPDATE: See CO2 Continuum for the spectrum of diluted CO2.

Wednesday, May 19, 2010

CO2 Doubling Temperature

One claim made by climatologists is that each doubling in the concentration of CO2 will cause the same increase in the Earth's surface temperature. So far as we can tell, the basis for this claim is that the absorption of radiation by greenhouse gases is "logarithmic". See Nelson, Earthguide and Wikipedia for examples. We never understood this "logarithmic" argument. We have never seen any calculations to back it up.

The greenhouse effect in our Planetary Greenhouse example is the result of logarithmic absorption by a greenhouse gas. The atmosphere contains a hypothetical gas that is opaque to some, but not all, long-wave radiation. As we showed in our previous post, successive doublings in the concentration of this gas do not cause the same increase in surface temperature. Indeed, the graph produced by our computer program shows that each successive doubling causes half as much warming as the previous doubling.

Thus logarithmic absorption alone is not sufficient justification for the existence of a "CO2 Doubling Temperature".

Monday, May 17, 2010

Planetary Greenhouse, Part Two

We continue our Planetary Greenhouse discussion. The atmosphere around our imaginary planet contains nitrogen and oxygen, which we choose because these diatomic gases are transparent to both short-wave and long-wave radiation. But the atmosphere also contains an imaginary greenhouse gas that is transparent to wavelengths shorter than λT, but absorbs longer wavelengths with constant absorption length, LT.

We wrote a computer program to determine the total heat radiated per square meter by the planet and its atmosphere. Before we can run the program, we must choose λT and LT.

Consider the absorption spectrum of pure CO2 at roughly 300 K and 100 kPa, which we present here. Between 8 μm and 11 μm the gas is transparent. The graph of radiated power versus wavelength for a black body at 300 K, which we present here, shows a peak at around 10 μm, right in the middle of CO2's transparent band. Meanwhile, the gas is at its most opaque for wavelengths greater than 13 μm. The absorption length drops to as little as 1 cm between 13μm and 17 μm, but is close to 1 m for 17 μm and higher. Let us choose λT = 13 μm and LT = 1 m.

Suppose our greenhouse gas makes up 330 ppm by volume of the atmosphere around our imaginary planet. Suppose the pressure at the planet surface is 100 kPa and the temperature is in the neighborhood of 300 K. The absorption length of the atmosphere at at the surface will be 1 m / 330 ppm = 3 km for wavelengths greater than λT. In The Upper Gas we showed how the absorption length at the surface is related to the pressure in the tropopause. With absorption length 3 km at the surface, the tropopause pressure will be 30% of the pressure at the surface, or 30 kPa. If we double the concentration of the greenhouse gas from 330 ppm to 660 ppm, the absorption length at the surface decreases to 1.5 km, and our tropopause must rise to an altitude where the pressure is only 15 kPa.

Now we can run our program, and we do so. We calculate the sum of the power radiated directly into space by the planet surface and by the troposphere. We plot this total radiated power versus the temperature of the planet surface, and we do so for a range of tropopause pressures.



The graph legend gives tropopause pressure in kPa. Each halving of pressure corresponds to a doubling in greenhouse gas concentration. The 100 kPa graph represents the case where the concentration is so low that the troposphere is at the planet surface. The planet radiates like a black body directly into space. The 1.6 kPa graph represents a concentration so high that the the troposphere radiates hardly any heat into space. The atmosphere acts to block all radiation at wavelengths greater than λT.

For thermal equilibrium in our Planetary Greenhouse, the total power radiated into space, EP + ET, must be equal to the total power received from the sun, ES. In Solar Heat we showed how the average power reaching the Earth from the Sun was 350 W/m2, and we calculated that the surface of a black planet with no atmosphere would be 280 K. In the 100 kPa graph above, we see that the surface temperature, TP, crosses 280 K when the total radiated power is 350 W/m2.

Let us suppose that the total radiated power in our Planetary Greenhouse must be 350 W/m2. With a tropopause pressure of 50 kPa, the planet surface must be at 299 K. If we double the concentration of our greenhouse gas, the planet warms by 13 K. Subsequent doublings in concentration raise the temperature by an additional 7 K, 5 K, 3K, and 1 K.

Wednesday, May 12, 2010

Planetary Greenhouse

We now enhance our Extreme Greenhouse to produce a more interesting model. We incorporate the decrease in atmospheric pressure with altitude. We use our understanding of atmospheric convection to determine how atmospheric temperature changes with pressure. We enhance the properties of our greenhouse gas by making it transparent to some long-wave radiation. The result is a simple model of a Planetary Greenhouse, as shown below.



In our new model, the atmosphere continues to infinite altitude, getting thinner and thinner as we go. The lowest layer of the atmosphere is the troposphere. Above the troposphere is the tropopause. The tropopause is the layer that radiates heat into space. Above the tropopause, the atmosphere is so thin that it is transparent.

For now, we ignore day and night, and we ignore clouds. We ignore heat flowing through the solid mass of the planet. Short-wave radiation arriving from the sun passes straight through the atmosphere and is of uniform and constant power ES. The temperature and pressure at the planet surface are Tp and Pp respectively. A near-adiabatic convection cycle carries heat up through the troposphere to the tropopause, where it is radiated into space. The temperature and pressure of the tropopause are TT and PT respectively.

In our Extreme Greenhouse, the atmosphere was entirely opaque to all long-wave radiation. In our new Planetary Greenhouse, the atmosphere is transparent to all radiation with wavelength less than some threshold, λT. At greater wavelengths, a trace greenhouse gas in the atmosphere absorbs radiation. The greater the concentration of this greenhouse gas, the shorter the absorption length of the atmosphere for wavelengths greater than λT. As we saw in The Upper Gas, the shorter the absorption length, the lower the tropopause pressure. As we saw in Atmospheric Convection, the lower the tropopause pressure, the lower its temperature.

As in our Extreme Greenhouse, we assume that the planet radiates like a black body. The radiation it emits below λT passes out into space. This radiated power is EP. As we described in Black Bodies, we can use Plank's Law to determine the power radiated at all wavelengths, so we can calculate the power radiated below λT as a function of TP.

The atmospheric convection cycle allows us to calculate the temperature of the tropopause given only TP and the ratio PT/PP. The tropopause radiates power ET into space. Using Plank's Law and radiative symmetry, we can calculate the power radiated by the tropopause as a function of TT. The atmosphere emits radiation only at wavelengths greater than λT.

For thermal equilibrium, we must have the total power radiated outward by the body and the tropopause equal to the total power arriving from the sun. Thus we must have ES = EP + ET. And so we arrive at a robust model of the planetary greenhouse effect when the greenhouse gas has a step-like absorption spectrum. Given only λT and PP, we can plot a graph of EP + ET versus temperature. This graph crosses ES at the equilibrium value of TP. We wrote a computer program to calculate EP + ET versus TP for a range of tropopause pressures. You will find the program here.

In our next post, we will choose useful values for λT and ES, and we will see how doubling the concentration of our greenhouse gas affects the temperature of the planet.

Wednesday, May 5, 2010

Dissipation by Convection

Now we return to our Extreme Greenhouse thought experiment. In our previous post on the subject, Work by Convection, we showed that atmospheric convection does work. We wondered where the energy for that work comes from, and where it goes to.

In our Atmospheric Convection post, the air taking part in our convection cycle warmed by 15 K while in contact with the warm surface of the Body, and cooled by 10 K while in the tropopause radiating its heat into space. The warming and cooling took place at constant pressure. The amount of heat required to warm 1 kg of a gas by 1 K at constant pressure is the heat capacity at constant pressure of the gas. Air is a close to an ideal gas. Its heat capacity at constant pressure is close to 1 kJ/kg/K at all the temperatures and pressures we are likely to encounter in the atmosphere of our Extreme Greenhouse. Each 1 kg of air absorbs 15 kJ of heat from the Body and radiates 10 kJ of heat into space. There are 5 kJ unaccounted for. These 5 kJ are converted into work.

Each kilogram of air taking part in the cycle rushes upwards with such force that it will do 5 kJ of work for us and still reach the tropopause. As it rushes upwards, it mixes with other air. Its momentum is diverted into whirlwinds. The viscosity of air is low, but it is not zero. The air moves upwards and sideways and in circles, rubbing against itself until the viscous friction it encounters consumes all the work available from convection. Viscosity turns the 5 kJ into heat.

As a result, the air in the tropopause will be warmer by 5 K than we expected with perfect, adiabatic expansion. Given that the air in the tropopause is over 100 K cooler than the air at the surface of the Body, this 5 K makes little difference to the size of our greenhouse effect. But it does mean that the air must cool in the tropopause by 15 K instead of 10 K, or else our convection cycle will not be closed. All 15 kJ of heat absorbed by the air from the Body will be radiated into space at the tropopause.

We know from our previous calculations that the atmosphere of our Extreme Greenhouse must transport 350 W/m2 from the Body surface to the tropopause (the boundary between the Filter Gas and the Upper Gas). Of these 350 W/m2, roughly 230 W/m2 will be carried directly by the heat capacity of the atmosphere, but the remaining 120 W/m2 will first be turned into work, later converted into heat by viscous action, and finally radiated into space.

The atmosphere of our Extreme Greenhouse has 120 W/m2 available to cause weather. This work is available only because the upper atmosphere is radiating heat into space. If there were no greenhouse gases in our Extreme Greenhouse, there would be no heat flow to the upper atmosphere, and no convection cycle to generate weather.