|San José State University|
& Tornado Alley
the Vertical Temperature Profile of the Atmosphere
Established by the Greenhouse Effect
Short wavelength radiation, visible and ultraviolet light, enters the Earth's atmosphere where some of it is absorbed by the ozone and some reflected back out by clouds. The rest passes through to the surface where part is absorbed and part reflected. That which is reflected and passes back up through the atmosphere where again some is absorbed by ozone. The greenhouse effect occurs when the heated surface of the Earth radiates infrared (long wavelength) radiation which greenhouse gases such as water vapor and carbon dioxide absorb and heat the air which in turn radiates infrared radiation upward and downward where it is reabsorbed and reradiated. The analysis is carried out first in terms of the simplest version of the model. That simplest version involves an atmosphere in horizon planes. A spherical model is considered elsewhere.
Radiation of intensity I traveling through a medium experiences diminishment according to the Beer-Lambert Law
where dz is distance through the medium, ρ is molecular volume density and λ is an absorption coefficient.
The level of density depends upon height and the value of λ depends upon the composition of the air; i.e., the proportional amounts of greenhouse gases such as water vapor, carbon dioxide and ozone in the air and the wavelength of the radiation.
The level of the dimension of the product λρ is inverse length. The values λ and ρ may vary with position z.
The Beer-Lambert Law implies that
The quanty exp(−∫0zλρds) is called the transmission coefficient and the quantity ∫0zλρds is called the optical path length. The analysis is improved, conceptually and analytically, if the optical path length is used as the length variable. Let y be the optical length (or optical height or optical depth) where
The optical path length depends on the wavelength of the radiation.
The crucial thing is that there are monotonic functions relating x and z and y and z.
The atmosphere does not end at some height hmax but trails off with lower and lower density. Thus z may not have have an upper limit but x and y may nevertheless have finite levels.
Since λρ and κρ have dimension of inverse length, the optical path lengths x and y are dimensionless variables.
Let J(z) be the intensity of the short wavelength radiation at height z. The intensity may be measured in various forms but the most convenient is in terms of energy per unit area, watts per square meter. In terms of optical path length the Beer-Lambert Law takes the form
The optical path here is measured from the surface upward. Since there is directionality to the radiation flow, let J- be the downward flowing radiation and J+ the upwardward flowing radiation. The two equations for the intensities are
The sign is right for J- because the intensity decreases as the optical height y decreases from ymax. Therefore as y increases J- also increases; i.e., dJ-/dy is positive.
At the surface, y=0, there is reflectance so
where the reflectance coefficient α is called the albedo.
The solutions for the intensities is simple.
Thus the short wavelength intensities, upward and downward, as functions of optical height y are fully determined.
In a medium which absorbs radiant energy the temperature increases and it radiates energy according to the Stefan-Boltzmann formula σT4, which is in the long wavelength band. The Beer-Lambert Law is replaced by the Schwarzchild Equation
where R(z)=σT4(z). This formula presumes Kirchhoff's Law that the emissivity of a substance is equal to its absorptivity. It also presumes that half of the thermal radiation is in the same direction as the radiation flow.
In terms of the optical path length x the Schwarzchild equation takes the form
The directionality must be taken into account. Let I+ be the radiation flux intensity in the direction of increasing optical height x and I- the flux intensity in the opposite direction.
The two equations for the optical path gradients of radiation intensities are
These relationships are for long wavelength radiation.
For further analysis it is convenient to define two new variables
The reason for these definitions is the equations for the rates of change of these two variables along the optical path take particularly interesting form. First note that from the gradient equations given above
From these conditions it follows that:
Now the matter of energy balance and temperature equilibrium must be considered.
Consider an infinitesimal element of thickness dz and cross-sectional area dA. For the moment consider only the long wavelength radiation intensities, I+ and I-.
The amount of energy absorbed in the element is λρ(I*+I-)dzdA. The amount radiated out is λρRdzdA. This total is from the previous relationships for the optical path gradients
The net inflow of radiative energy to the element is
Had the short wave length radiation intensities been included the above result would have been
That net inflow of energy to the infinitesimal element must increase its temperature according to
where ρ has to be mass density, cp is specific heat at constant pressure and t is time.
At equilibrium dT/dt=0 so it must be that
This means that (I+-I-)+(J+-J-)] has to be constant over the optical path and hence φ=−½(J+-J-), a known quantity. This means that the differential equation for μ can be solved by integration; i.e.,
From the knowledge that φ=−½(J+-J-) and that dφ/dx=−μ+R it follows that
Thus a formal solution has been found for the radiation intensities and temperature as a function of height.
Some further insights into the nature of the solution may be obtained by considering the relationships that exist for the gradients of the J's. Let the ratio of the absorptivities of the two wavelengths of radiation (κ/λ)=(dy/dx) be denoted as ε(x). Then note that
If the atmosphere were well mixed at all levels then ε would be constant. However there are much higher levels of water vapor in the troposphere than the stratosphere. On the other hand, there are higher levels of ozone in the stratosphere than in the troposphere. Nevertheless the case of constant ε is an important special case for mathematical analysis. For this case it would follow that
This holds however only for a vertically uniform atmosphere.
The determining condition for the equilibrium of the system is that the outflow of long wavelength energy at the top of the atmosphere must be equal to the net inflow of short wavelength radiation at the top of the atmosphere. At the top of the atmosphere I-=0 and thus I+=2φ and μ=φ. Thus if the net energy input of short wavelength radiation is S=[J-(xmax)−J+(xmax)] then
One interesting implication of the solution is that there is a temperature discontinuity at the surface. In other words, the surface is at a higher temperature than the air just above it. This is actually a familiar phenomenon as at the beach where the sand may be too hot to walk upon barefooted while the air temperature is not uncomfortable.
The analysis indicates that the air temperature decreases with altitude as a result of the greenhouse effect of the absorption of thermal radiation by water vapor, carbon dioxide and the other greenhouse gases except ozone. Ozone absorbs visible and ultraviolet radiation and the temperature effect of this absorption increases with height. In the real atmosphere temperature decreases with height up to a particular level called the tropopause and above that level the temperature increases with altitude (in the stratosphere).
(To be continued.)
HOME PAGE OF Thayer Watkins