1. Top-of-atmosphere
Top-of-atmosphere
(TOA) graphs of radiation flux from earth show ‘notches’ centered at the
nominal wavenumbers of long wavelength infrared (LWIR) active gas molecules,
mostly CO2, which do not condense in the atmosphere. Figure 1 appears
to be a typical mid-latitude TOA graph. Radiation flux intensity is plotted vs.
wavenumber at any one geometric location (line-of-sight). Wavenumber is simply
the number of radiation wavelengths per centimeter. Wavenumber, cm-1
and wavelength in microns are easily converted: Divide 10,000 by either to get
the other.
Part of the power (flux, energy-rate) radiated from the surface goes directly to space through the ‘atmospheric window’ in the approximate wavenumber range 770-1230 (13-8.13 microns) where no gas molecules (except ozone, O3) significantly absorb radiation. Over the rest of the total range of significant terrestrial radiation (wavelength 5-2000 microns, wavenumber 2000-5 cm-1) the radiation is absorbed by Long-Wave-InfraRed (LWIR) active gases (which, including water vapor (WV), are misleadingly called greenhouse gases (ghg)).
Height of the black line trace is the flux (energy-rate) as a function of wave number at the TOA at a single line-of-sight. This line accounts for the power that radiates from the surface plus that from latent and convective heat from the surface and solar radiation absorbed by clouds and WV in the atmosphere. The vertical distance on the graph from a curve representing radiation from the earth’s surface to the black line indicates the amount at each wave number that the radiation at TOA (energy rate) is less than the equivalent energy rate radiating from the surface. Energy is conserved so the rate is what changes, i.e. the flow of energy is slowed down as it passes through the atmosphere. The slowing is proportional to the logarithm of the number of molecules of the responsible ghg.
2. Atmospheric structure
At the scale of
atoms, the atmosphere can be visualized as molecules bouncing elastically (no
energy loss) off each other in empty space. At sea level conditions, the time
between collisions for the air molecules (molecule diameter about 4 Angstroms;
4E-10 m) is extremely short, less than 0.0002 microseconds [4]. Among other
properties, these collisions are the basis for thermal conduction in the gas.
Therefore, electromagnetic radiation (EMR) energy absorbed by ghg molecules is
immediately (within a few nanoseconds) shared with surrounding molecules both
ghg and non-ghg. The sharing is thermal conduction in the gas. The process of
absorbing radiation energy and sharing it with surrounding molecules is called
thermalization.
Logic mandates that the elapsed time between when a molecule absorbs a photon and when it emits one must be more than zero or there would be no indication the photon had been absorbed. This elapsed time, called decay time, has been calculated for CO2 and averages from about 10 microseconds at sea level up to about 1.1 second at TOA [2, 3]. Decay time is very long compared to the relaxation time of a nanosecond or so which is the time required for thermalization to take place. I suspect that, early in thermalization, nearby ghg molecules are stimulated to emit photons that might sum to energy-rate different from that absorbed.
Power in Watts is defined as energy-rate by W = J/s where
J is energy in Joules and s is time in seconds. One Joule for one second is one
W by definition.
3. Steady-State
Power Balance
As a starting point in the analysis, the flow of power is
determined partly using estimates from other sources, some rational
refinements, and an approximation of the power loss. A graphic that was copied
from [15] is useful as a qualitative description of the various power
flows. It is shown in Fig 1.1 with some revisions. Some power flows are more
explicitly defined than at the NASA site and the values presented there are
replaced by symbols (alpha characters) with values determined herein. All power
magnitudes are in watts per square meter.
4. Calculation
of the Magnitudes of Power (Energy Flow Rate)
Measured values of power from the sun vary only slightly
depending on the solar cycle and slight other long-term variations in the
‘solar constant’. The value used here is 1365 W/m^2 of insolation. The power
from the sun is intercepted by the planet according to the cross-section area
but is distributed according to the surface area so this is divided by 4 to
account for the ratio of surface area to cross-section area of a sphere.
Q = 1365/4 = 341.25 W/m^2
Part of this incident power is reflected. The fraction of
the incident power that is reflected is often assumed to be determined by the
earth’s albedo which, from what appears to be a credible source [13], is about
0.297. However, albedo is primarily diffuse reflection and does not include
Fresnel-type low-incident-angle specular reflection [14] from the 71% of the
planet that is covered by water. The fraction of the total power reflected from
the planet is slightly higher than that determined by albedo alone and is about
0.308 (my calculation).
R = Q * 0.308 = 105.105 W/m^2
The amount of power reflected from the clouds and
atmosphere was determined as required to maintain the same ratio of cloud + atmosphere
reflection to surface reflection as is reported by Kiehl and Trenberth [15].
A = R/ (1+23/79) = 81.4 W/m^2
Part of the remaining power is absorbed by clouds and
atmosphere. I selected a value of 48% for this to be consistent with the
absorbed power rate assessment by Kiel & Trenberth. An additional factor of
0.6 accounts for the observation that clouds cover about 60% of the planet.
B = (Q-A) * 0.48 * 0.6 = 74.84 W/m^2
The incoming power that gets through the atmosphere is
simply that which is not reflected or absorbed by the clouds and atmosphere.
C = Q - A - B = 185.01 W/m^2
The part of this that is reflected by the surface is simply
that part of the total reflection that was not reflected by the clouds and
atmosphere.
D = R - A = 23.71 W/m^2.
The power that is absorbed by the surface is the part of C
(the entering power that got through the atmosphere) that is not reflected.
E = C - D = 161.31 W/m^2.
The power leaving the surface as convection (thermal),
evaporation (latent) and radiation from the surface that directly leaves the
planet are all as presented on the K & T chart except G which is based on
average global precipitation.
F = 17 W/m^2
G = 78 W/m^2
J = 40 W/m^2
The radiation from the surface is gray-body radiation as
calculated by the Stephan-Boltzmann equation. Most (71%) of the surface of the
planet is covered by sea water with an emissivity of about 0.995 at water
temperature. Part is from snow with emissivity of about 0.99 (clean) and the
rest is from land that is mostly about 0.99 but with a few small local areas
that are substantially lower. The over-all average emissivity from the surface
is taken to be 0.98. The average global temperature is taken to be 288 K as
typically reported. The S-B constant, σ, is 5.6697E-08 W/m^2 /T^4. (Temperature
in degrees Kelvin)
U = 0.98 * 5.6697E-8 * 288^4 = 382.26 W/m^2
Clouds are fine particles of liquid water or ice and thus
also radiate according to the S-B equation. Clouds cover about 60% of the
planet and have an average emissivity [16] of about 0.5. Their average
temperature was determined to be 258 K in an entirely different analysis that
produced an average global temperature of 288K. The thin air and low
temperature at high altitude means that there is very little water vapor so radiation
up from the clouds nearly all gets directly to space. The small part that
doesn’t is ignored resulting in:
P = 0.6 * 0.5 * 5.6697E-8 * 258^4 = 75.36 W/m^2
The same power flux as from the top of the clouds exists
also from the bottom of the clouds. The down flux encounters substantial
absorption prior to reaching the ground. K & T report 40 W/m2 or 40/396 =
10.1 % of the radiation making it from the surface through the ghg to space. A
calculation using data from Barrett [17] for 60% cloud cover produces 10.8 %.
Since about half of the source radiation is stopped by clouds in going from the
surface to space and not stopped by clouds when going from clouds to the
surface, the fraction getting to the surface is about 21% of that going from
the surface directly to space.
K = 0.21 * P = 15.83 W/m^2.
The total power being radiated from the planet must be the
amount from the sun minus the amount reflected.
M = Q - R = 236.10 W/m^2.
The power radiated from the atmosphere must be the total
for the planet minus that just from the clouds to space minus that which goes
directly from the surface to space.
N = M - P - J = 120.78 W/m^2.
Energy must be conserved so the rate entering the
atmosphere must be slowed. The
energy-rate (power) is slowed by about 33 ms due to ‘random walk’ of photons up
through the atmosphere. This equates to a power loss of
H = (M – J)/ (1-0.033) –(M-J) = (236 – 40)/ (1.0.033) = 6.692
W/m^2
The net power that is returned to the surface is
I = U + G + F + B -M – H -K = 293.3 W/m^2.
Figure 1.1: Power flow in the atmosphere.
5. Power assessment
for the atmosphere:
1. (75 W/m^2) The power entering the
atmosphere from the sun (absorbed by WV in the 1.38- and 1.97-micron bands,
2. (78 W/m^2) 3. latent heat from
hydrologic cycle,
3. (17 W/m^2) convection,
4. (382 W/m^2) radiation from the
surface at emissivity of 0.98 and T = 288 K,
5. (˗293.3 W/m^2), back radiation to
surface,
6. (˗16 W/m^2) radiation from clouds
to ground thru atmospheric window,
7. (˗40 W/m^2) radiation from surface
to space thru atmospheric window.
This all adds up to 202.7 W/m^2 of unreduced power. The photon
energy travels from ghg molecule to ghg molecule in a ‘random walk’ up thru the
atmosphere. The average path length of the ‘random walk’ is approximated by
calculation [10] to be a million 10-meter steps. This path is traversed at the
speed of light in the atmosphere in about 33 ms. This slows the energy-rate to
0.967 * 202.7 = 196 W/m^2 which matches the power leaving the atmosphere. Water
vapor increase increases the path length of the random walk and contributes to the
increase to the misleadingly named GHE and the increase in average global
temperature.
6. Atmospheric water vapor
A common observation which shows that WV increase
contributes to global warming is that cloudless nights cool faster and farther when
absolute WV content of the atmosphere is lower. Clear nights cool faster and
farther in the desert than where it is humid.
Below the tropopause (below about 8 to 16 km depending mostly on latitude; more at the equator) WV decreases from an average of about 8,000 ppmv (0.8%, 4% near the equator) at ground level, because of the low temperature (about negative 50 °C), to 203 ppmv at most, at the tropopause (assume about 12 km). Saturation vapor pressure of ice at ˗50 C/total pressure at 12 km = 3.94 Pa/19400 Pa = 0.000203 = 203 ppmv (at saturation). Thus, at the troposphere and above, radiation to space is mostly from CO2 and other IR active molecules that do not condense in the atmosphere. Increased CO2 in the extremely thin air there actually counters planet warming.
In addition to the decline in WV ppmv due to temperature
decline, is the decline resulting from pressure decline with altitude of
19.4/101.1 = 0.191. The total WV molecule population gradient from surface to
tropopause is thus about 8000/203 * 101.1/19.4 = 206 to 1.
7. Thermalization and redirection
Thermalization causes some of the energy absorbed by CO2
molecules to be redirected by gaseous thermal conduction to WV molecules. The
Schwarzschild equation [12] accounts for this with “as radiation passes through
an isothermal layer, its monochromatic intensity exponentially approaches that
of blackbody radiation corresponding to the temperature of the layer.” Much of
the radiant energy along a line-of-sight is ‘redirected’ with respect to wave
number so CO2 molecules become transparent to it.
8. Water Vapor Factor
Water vapor is a transparent gas
that, molecule for molecule, is more effective at absorb/emit of
earth-temperature infrared radiation (IR) than carbon dioxide (at low altitude,
absorb lines for CO2 overlap). From Jan 1988 thru Dec 2024 NASA/RSS
accurately measured and reported monthly the global average WV as Total
Precipitable Water (TPW). The anomaly data are reported at [11]. The nominal
value is about 29 kg/m^2 so the trend from Jan 1988 thru Dec 2024 is about 1.5
% per decade.
Given that at ground level
average global WV is about 0.8% or 8,000 ppmv (parts per million by volume)
(Ref 34, 35 of [10]), the increase in WV molecules in 3.6 decades is about 0.015
* 8000 * 3.6 = 432 ppmv. From Mauna Loa data at [9] the CO2 increase in that time period is 423 - 349 = 74 ppm. Per
ideal gas laws, ppm = ppmv. With that, at ground level, WV molecules have been
increasing 432/74 = 5.8 times faster than CO2
molecules. Thus, at ground level, regardless of the initial source of warming, WV
molecules have been increasing about 5.8 times faster than CO2 molecules. The idea that CO2
starts the increase is ludicrous considering the comparatively large
fluctuations in surface temperature.
Radiation from WV molecules can
be in any direction but, because of the steep decline with altitude of the
population gradient of WV molecules, the distance traveled by a photon before
it encounters another WV molecule is greater towards space than towards earth
so the prevailing direction of IR flux is towards space. This is shown on a Top
of Atmosphere (TOA) graph of radiation flux vs wavenumber by the jagged line
below about wavenumber 600. Because of the characteristic absorb/emit signature
of every gas, no other gas can significantly absorb or emit radiation in the
wavenumber range occupied by WV. The line is jagged because radiation that
reaches TOA/space is from WV molecules at different temperatures/altitudes.
9. Top-of-atmosphere
with black body overlay
Fig 1.5 is a TOA graph with
overlaid constant temperature curves for black body radiation at noted
temperatures. For any specified atmosphere, these are also constant elevation
curves. The temperature and associated altitude for standard atmosphere are
shown in the upper right corner of Fig 1.5.
Figure 1.5: Typical TOA
radiant emission. (U Chicago version of MODTRAN [8]
At about 2 km and higher, the
outward directed radiation from WV can make it all the way to space. Below the
tropopause, much of the energy absorbed by CO2 and
other IR active molecules is redirected with respect to wave number via thermalization
to WV molecules. This mitigates any warming from increased CO2 in the troposphere (or any other IR active gas that does
not condense at earth temperatures).
The end result is that CO2 does not cause significant climate change, sequestering
it is an expensive mistake and the Green New Deal would have no significant
effect on climate.
10.
Radiance calculated by MODTRAN6
MODTRAN6 [5] is
a computer program developed for the Airforce Research Laboratory which
(besides other things) can calculate the radiation flux at selected elevations
in the atmosphere for specified constituents and conditions. It contains
default values for several environments including the tropics and the 1976
Standard Atmosphere. Values for WV change rate and atmospheric temperature vary
with altitude for different latitudes and seasonal conditions as shown in
MODTRAN6 documentation [7].
Most of the photons emitted by the WV molecules are at wavelengths different from the comparatively narrow band that CO2 molecules can absorb so they are absorbed by other WV molecules. Effectively, below the tropopause, much of the terrestrial thermal radiation energy and other energy absorbed by CO2 (and other non-condensing ghg) is thermalized, redirected to, and radiated to space from WV.
The WV content of the atmosphere diminishes rapidly as the temperature decreases with increasing altitude. Above the tropopause it has declined to a level where emission from WV ceases to dominate and emission from CO2 molecules becomes significant. The result is most of the residual energy not emitted to space by WV at low altitude is, at high altitude, redirected back to the wavenumber range 600-740 cm-1 and emitted to space by CO2. The ’redirection’ is not geometric because all wavenumbers refer to photons at essentially the same line-of-sight.
Energy redirection helps explain why average global temperature tracks average global WV and not CO2 [6].
References:
1. NASA/GISS TOA
graph https://www.giss.nasa.gov/research/briefs/2010_schmidt_05/
2. Average elapsed time to emit a photon https://sealevel.info/Happer_UNC_2014-09-08/Another_question.html
3. Average elapsed time to emit a photon http://rabett.blogspot.com/2013/04/this-is-where-eli-came-in.html
4. Time between
gas molecule collisions http://hyperphysics.phy-astr.gsu.edu/hbase/Kinetic/frecol.html
5. MODTRAN6
calculator http://modtran.spectral.com/modtran_home#plot
6. Climate
change drivers http://globalclimatedrivers2.blogspot.com
7. MODTRAN6
defaults http://modtran.spectral.com/modtran_faq
8. MODTRAN
calculator http://climatemodels.uchicago.edu/modtran/
9. Mauna
Loa data for CO2: https://www.co2.earth/monthly-co2
10. Water Vapor
vs CO2 for Planet Warming, Sect 1 of: https://watervaporandwarming.blogspot.com
11. NASA/RSS
measurements thru Dec 2023 of Total Precipitable Water anomalies https://data.remss.com/vapor/monthly_1deg/old_files/tpw_v07r02_198801_202312.time_series.txt
12. Wallace and
Hobbs: http://www.gnss-x.ac.cn/docs/Atmospheric%20Science%20An%20Introductory%20Survey%20(John%20M.%20Wallace,%20Peter%20V.%20Hobbs)%20(z-lib.org).pdf
13. Earth’s
albedo: https://ui.adsabs.harvard.edu/abs/2001GeoRL..28.1671G/abstract
14. Specular
reflection from water: https://en.wikipedia.org/wiki/Reflectance#Reflectivity
15, Kiel &
Trenberth graphic: https://chriscolose.wordpress.com/2008/12/10/an-update-to-kiehl-and-trenberth-1997/
16. Cloud
emissivity: https://www.arm.gov/publications/proceedings/conf08/extended_abs/shippert_tr.pdf
17. Barrett,
greenhouse molecule function: https://www.warwickhughes.com/papers/barrett_ee05.pdf


