|San José State University|
& Tornado Alley
the Latitudinal Temperature Profile
and Latitudinal Heat Flux
In meteorological studies of the transfer of heat energy from the equatorial regions toward the polar regions the problem is characterized as one of establishing the flow of heat energy from regions of a heat surplus to regions of a heat deficit. The regions of heat surplus are those in which the absorption of radiant energy from the Sun exceeds the thermal radiation. The heat deficit regions where more energy is radiated away into space than comes in from the sun. This characterization of the problem implies that there are some regions that are intrinsically heat surplus areas and others that are heat deficit areas. This gets the causalities mixed up. The temperature profile is established simultaneous with the heat flows. One does not determine the other; they are simultaneously and mutually determined. The models examined below are ones in which this simultaneous determination is made explicit.
Let φ be latitude and F(φ) be the incoming solar radiation energy (insolation) at φ. The albedo (reflectance) α could be a function of latitude and/or temperature T. The temperature T is generally a function of latitude, T(φ).
The energy radiated from the surface is equal to εσT4, where σ is the Stefan-Boltzmann constant and ε is an emissivity factor representing the possibility that the surface is not a perfect blackbody radiator.
A balance of incoming and outgoing radiative energy at the surface requires that
The computation of the insolation F(φ) as a function of latitude is complicated for a planet whose axis of rotation is not perpendicular to the plane of its orbit, but for a planet with a zero angle of inclination the value is simply μcos(φ) where μ is the solar constant. For a constant albedo then
which has the shape shown below.
In materials heat diffuses through molecular interaction from areas of high temperature to areas of low temperature. In a bar of material of cross-sectional area A the amount of heat-flow S at a position x is proportional to the temperature gradient but in the opposite direction; i.e., it is given by
where k is a constant, the coefficient of heat transmission. This type of heat transmission equation could apply to the atmosphere through other processes of diffusion such as eddy diffusion.
To work out the equation that must be satisfied by the temperature gradient and the heat fluxes consider a ring-wall of atmosphere between latitude φ−Δφ/2 and latitude φ+Δφ/2. Let a be the radius of the Earth and for convenience let that radius extend to the midpoint of the atmosphere which is of depth H. The ring-wall then has the following dimensions;
These dimensions then result in:
It is desirable to use a linear transmission coefficient k which would be such that the heat flux through an area A is
The heat flux through the wall side at φ−Δφ/2 is then
The net gain in heat within the volume of the ring-wall is then
The expression within the bracket may be replaced by
When the radiation flows are included the net gain in heat to the ring-wall is
Now consider the volume of the ring-wall. Let cv and ρ be the heat capacity and the mass density of the atmosphere. The volume of the ring-wall is (2πacos(φ))(aΔφ)H. The rate of increase of heat energy in the ring-wall is
This net rate of increase of heat energy in the ring-wall has to equal the net flux of heat energy through the surface of the ring-wall. There is a common factor of 2πΔφ which can be factored leaving as the requirement for the heat balance
where ∂T/∂t is the derivative of temperature T with respect to time t.
Dividing through by the coefficient of ∂T/∂t gives
At equilibrium where ∂T/∂t=0
In principle this equation can be solved subject to the boundary conditions that ∂T/∂φ=0 at φ=0 and φ=π/2.
(To be continued.)
HOME PAGE OF Thayer Watkins