Radiative transfer#

Radiation modelling is a notorious challenge due to its highly coupled nature. Our model’s goals must be defined so that when approximations are made to the problem, the sources of discrepancy can be understood. Throughout the book, most of the effort in the chapters has been focused around \(\textrm{CO}_2\). Therefore, to continue this focus, the developed radiation model is focused on modelling the greenhouse effect caused by molecules re-radiating longwave radiation into the atmosphere. We take the only source of radiation to be black body radiation emitted by the earth, assuming it is a perfect black body, at some temperature \(T\) and wavenumber \(\nu\), the spectral radiance is

(41)#\[B(T, \nu) = 2hc^2 \nu^3 \frac{1}{e^{hc\nu/(k_bT)} - 1}.\]

Note

The hitran convention for wavenumber is used where it is denoted as \(\nu\) rather than \(\tilde{\nu}\) in many texts.

Additionally, because the effects of incoming solar radiation are not considered, which performs a substantial amount of heating in the stratosphere with the absorption of UV by \(\textrm{O_3}\), limiting the model to only extend up to 10km enables the model to both avoid warming due to shortwave and relative gas concentration changes as discussed in Standard Atmosphere Model.

A simple atmosphere model is a single block of some height with an associated optical depth. An observer looks down towards the Earth’s surface, and radiation is emitted perpendicular to the surface. It travels through the block and arrives at the observer’s eye. The Beer-Lambert Law describes this

(42)#\[-\frac{d}{d OD_\nu} I(OD_\nu, \hat{n}, \nu) = I(OD_\nu, \hat{n}, \nu).\]

The parameter \(OD\) is the optical depth of the medium. This equation describes a radiation model that only models the amount of light, which is not scattered by the medium. This is solved by a simple exponential decay, which has a decay length inversely proportional to the optical depth of the medium.

(43)#\[I = I_0 \exp (-OD)\]

This is an insufficient model for our use case; firstly, a key behaviour of greenhouse gases is that they re-radiate isotropically, and secondly, they have differing optical depth depending on wavenumber.

We iterate upon the Beer-Lambert law to include the isotropic re-radiation with the Schwarzschild equation. The Schwarzschild equation adds a black body radiation term to the Beer-Lambert law, which is dependent on the optical depth of the medium. Further, it is assumed that the incident radiation is incident on the medium with a direction defined by a unit vector \(\hat{n}\)

(44)#\[\frac{d}{d \tau_nu} I(OD_{\nu}, \hat{n}, \nu)= -\frac{1}{\cos\theta} [I(OD_{\nu}, \hat{n}, \nu) - B(T(OD_{\nu}), \nu)].\]

Where \(OD = k(\nu, p, T)u\), \(u\) is the column number density of absorbing species and \(k(\nu, p, T)\) is the absorption coefficient per molecule per unit length. The re-radiation term is the plank function \(B(T(OD_\nu), \nu)\). Additionally, the \(\cos \theta\) term comes from the distance through the medium of optical depth \(\delta OD\), where, if the path of the incoming radiation is at some angle \(\theta\) to the perpendicular from the surface, the distance the radiation travels is further \(\delta OD / \cos \theta\), as seen in the figure below from[Pie10]. This model is taken from sections of chapter 4 in \(principles of planetary climate\)[].

If it is assumed that the incoming and outgoing radiation isotropic as it is emitted black body radiation, then the directional dependence of the equation is removed. The Solution can then be integrated over a hemisphere,and two separate ordinary differential equations for upward and downward propagating spectral radiance can be performed, assuming that \(\cos \bar{\theta} = 1/2\)

(45)#\[\begin{align} \int_0^{\pi/2}I(T(OD_{\nu}),\nu) d\Omega &= \int_0^{\pi/2} -\frac{1}{\cos\theta} [I(OD_{\nu}, \hat{n}, \nu) - B(T(OD_{\nu}), \nu)] d\Omega, \\ \frac{d}{d OD_{\nu}} I_+ &= - I_+ + \pi B(T_{\nu}, \nu). \end{align}\]

\(I_+\) is all the radiation in the upward direction from the block of atmosphere. In the downward direction \(I_-\), the equation is similar

(46)#\[\begin{align} \int_{\pi/2}^{\pi}I(T(OD_{\nu}), \nu) d\Omega &= \int_{\pi/2}^{\pi} -\frac{1}{\cos\theta} [I(OD_{\nu}, \hat{n}, \nu) - B(T(OD_{\nu}), \nu)] d\Omega, \\ \frac{d}{d OD_{\nu}} I_- &= I_- - \pi B(T_{\nu}, \nu). \end{align}\]

As the temperature of the radiation block is a function of the atmosphere’s optical depth, Ideally, it should be calculated rather than taken from the ISA model. However, this isn’t easy. Often, atmospheres are not in just radiative equilibrium because they are comprised of gases. Convection plays a key role in energy transfer as well. Thus blocks of the atmosphere may infect heat and not be in thermal equilibrium with their surroundings. A model which includes convection is much better approximated by an adiabatic atmosphere where warm blocks of air rise, expand, and cool adiabatically. This, is typically the assumption for the lower atmosphere, troposphere[Pie10]. For higher altitudes, convection is less important and radiative equilibrium provides a better approximation.

Despite it not being a good approximation for the lower atmosphere, radiative equilibrium provides a more straightforward path to a solution. With this boundary condition two equations for the upwards and downward spectral radiance can be formed:

(47)#\[\begin{align} 0 = \frac{d}{d OD_{\nu}}I_+ - I_- &= -(I_+ + I_-) + 2 \pi B(T(OD_{\nu}, \nu), \\ \frac{d}{d OD_{\nu}}I_+ + I_- &= -(I_+ - I_-). \end{align}\]

As an initial start point, a grey gas atmosphere is considered. This is where optical depth is independent of wavenumber.

Grey Gas Model#

With wavenumber independence, we can consider the total flux given by the integral over the integral over wavenumber and the solid angle of the plank function, which gives the Stephan-Boltzmann law

(48)#\[\sigma T^4 = \int \int d \Omega d \nu B(T,\nu).\]

With the solution to eq.label now,

(49)#\[I_+ - I_- = C,\]

Thus

(50)#\[I_+ + I_- = -c OD + c'\]

Using the boundary condition that the outgoing longwave radiation(OLR) at the top of the atmosphere is a constant \(I_{OLR}\)

(51)#\[I_+ + I_- = c' - c OD_{\infty} = I_{OLR},\]

Then

(52)#\[\begin{align} c &= I_{OLR}, \\ c' &= I_{OLR}(OD_{\infty} + 1). \end{align}\]

Finally, substituting into eq.(51)

(53)#\[T = \left(\frac{(1 + OD_{\infty} - OD)I_{\infty}}{2 \sigma} \right)^{1/4}.\]

This gives us the expression for the Grey Gas at zero altitude. A comparison between this and the planet’s surface temperature is plotted below.

import matplotlib.pyplot as plt
import matplotlib as mpl
import numpy as np
from simpletrans import isa
import scipy.constants as cst
from dur_utils import colours

mpl.style.use("../CDS.mplstyle")
max_alt = 10000  # Alt strat stop
SBC = cst.physical_constants["Stefan-Boltzmann constant"][0]
olr = SBC * 2 * isa.get_temperature(max_alt) ** 4
points = 100
altitudes = np.linspace(0, max_alt, points, dtype=float)
taus_inf = np.logspace(-2, 1, 8)
isa_temp_profile = isa.get_temperature(altitudes)


def ground_gas_temp(od_inf: float) -> float:
    """
    Returns the ground temperature for a particular optical thickness of the
    whole atmosphere tau_infinity.
    Args:
        od_inf: The total optical Thickness of the atmosphere.

    Returns:
        t_ground

    """
    t_ground = ((1 + od_inf) * olr / (2 * SBC)) ** 0.25
    return t_ground


def ground_temp(od_inf: float) -> float:
    """
    Returns The surface temperature of the planet assuming radiative equilibruim as a function of
    Args:
        od_inf:

    Returns:

    """
    t_ground = ((2 + od_inf) * olr / (2 * SBC)) ** 0.25
    return t_ground

plt.figure(figsize=(10,6))
plt.semilogx(
    taus_inf,
    ground_gas_temp(taus_inf),
    c=colours.durham.ink,
    label="Grey gas temperature at ground level",
)
plt.semilogx(
    taus_inf,
    ground_temp(taus_inf),
    c=colours.durham.purple,
    label="Surface Temperature",
)
plt.semilogx(
    taus_inf,
    [isa.get_temperature(0) for i in taus_inf],
    linestyle="--",
    c=colours.durham.ink,
    label="Earth Surface Temperature",
)
plt.legend()
plt.xlabel(r"Total atmosphere optical thickness: $OD_{\infty}$")
plt.ylabel("Temperature (K)");
../_images/radiative_transfer_2_0.png

As the optical thickness of the atmosphere increases, the surface temperature and the Grey Gas surface temperature approach each other, reducing the need for convection.

This does highlight a problem at low altitudes, where the total optical thickness is low. As a result, there is a large differential between the surface temperature and the atmosphere temperature. This is a poor model for Earth, as the ground is in equilibrium with the air at surface level. Clearly, this radiative equilibrium assumption is not close enough to be valid throughout the whole atmosphere for our purposes.

An assumptions change is required. On the next page, we assume being in radiation and convection-based equilibrium. This removes the boundary condition that each layer must have the same amount of flux in as it lets out. When convection is included, it is assumed that hot pockets of gas rise and adiabatically expand, lowering their temperature to maintain thermal equilibrium with their surroundings. This enables the model block atmosphere to transfer energy between atmosphere altitudes. For the models, these convection currents are not included. So energy equilibrium conservation is assumed but flux conservation is not required.