# Schwarzschild's equation for radiative transfer

Schwarzschild's equation is used to calculate radiative transfer – energy transfer – through a medium in local thermodynamic equilibrium that both absorbs and emits electromagnetic radiation.

The incremental change in spectral intensity, (dIλ, W/sr/m2/μm) at a given wavelength as radiation travels an incremental distance (ds) through a non-scattering medium is given by:

$dI_{\lambda }=n\sigma _{\lambda }B_{\lambda }(T)\,ds-n\sigma _{\lambda }I_{\lambda }\,ds=n\sigma _{\lambda }[B_{\lambda }(T)-I_{\lambda }]\,ds$ where n is the density of absorbing/emitting molecules, σλ is their absorption cross-section at wavelength λ, Bλ(T) is the Planck function for temperature T and wavelength λ, and Iλ is the spectral intensity of the radiation entering the increment ds. This equation and a variety of equivalent expressions are known as Schwarzschild's equation. The second term describes absorption of radiation by the molecules in a short segment of the radiation's path (ds) and the first term describes emission by those same molecules. In a non-homogeneous medium, these parameters can vary with altitude and location along the path, formally making these terms n(s), σλ(s), T(s), and Iλ(s). Additional terms are added when scattering is important. Integrating the change in spectral intensity (W/sr/m2/μm) over all relevant wavelengths gives the change in intensity (W/sr/m2). Integrating over a hemisphere then affords the flux perpendicular to a plane (F, W/m2).

Schwarzschild's equation contains the fundamental physics needed to understand and quantify how increasing greenhouse gases (GHGs) in the atmosphere reduce the flux of thermal infrared radiation to space. If no other fluxes change, the law of conservation of energy demands that the Earth warm (from one steady state to another) until balance is restored between inward and outward fluxes. Schwarzschild's equation alone says nothing about how much warming would be required to restore balance. When meteorologists and climate scientists refer to "radiative transfer calculations" or "radiative transfer equations" (RTE), the phenomena of emission and absorption are handled by numerical integration of Schwarzschild's equation over a path through the atmosphere. Weather forecasting models and climate models use versions of Schwarzschild's equation optimized to minimize computation time. Online programs permit anyone to perform computations using Schwarzschild's equation.

## Background

Radiative transfer refers to energy transfer through an atmosphere or other medium by means of electromagnetic waves or (equivalently) photons.  The simplest form of radiative transfer involves a collinear beam of radiation traveling through a sample to a detector.  That flux can be reduced by absorption, scattering or reflection, resulting in energy transmission over a path of less than 100%.  The concept of radiative transfer extends beyond simple laboratory phenomena to include thermal emission of radiation by the medium - which can result in more photons arriving at the end of a path than entering it.  It also deals with radiation arriving at a detector from a large source - such as the surface of the Earth or the sky.  Since emission can occur in all directions, atmospheric radiative transfer (like Planck's Law) requires units involving a solid angle, such as W/sr/m2.

At the most fundamental level, the absorption and emission of radiation are controlled by the Einstein coefficients for absorption, emission and stimulated emission of a photon (B12, A21 and B21) and the density of molecules in the ground and excited states (n1 and n2). However, in the simplest physical situation – blackbody radiation – radiation and the medium through which it is passing are in thermodynamic equilibrium, and the rate of absorption and emission are equal. The spectral intensity (W/sr/m2/μm) and intensity (W/sr/m2) of blackbody radiation are given by the Planck function Bλ(T) and the Stefan–Boltzmann law. These expressions are independent of Einstein coefficients. Absorption and emission often reach equilibrium inside dense, non-transparent materials, so such materials often emit thermal infrared of nearly blackbody intensity. Some of that radiation is internally reflected or scattered at a surface, producing emissivity less than 1. The same phenomena makes the absorptivity of incoming radiation less than 1 and equal to emissivity (Kirchhoff's law).

When radiation has not passed far enough through a homogenous medium for emission and absorption to reach thermodynamic equilibrium or when the medium changes with distance, Planck's Law and the Stefan-Boltzmann equation do not apply. This is often the case when dealing with atmospheres. If a medium is in Local Thermodynamic Equilibrium (LTE), then Schwarzschild's equation can be used to calculate how radiation changes as it travels through the medium. A medium is in LTE when the fraction of molecules in an excited state is determined by the Boltzmann distribution. LTE exists when collisional excitation and collisional relaxation of any excited state occur much faster than absorption and emission. (LTE does not require the rates of absorption and emission to be equal.) The vibrational and rotational excited states of greenhouse gases that emit thermal infrared radiation are in LTE up to about 60 km. Radiative transfer calculations show negligible change (0.2%) due to absorption and emission above about 50 km. Schwarzschild's equation therefore is appropriate for most problems involving thermal infrared in the Earth's atmosphere. The absorption cross-sections (σλ) used in Schwarzschild's equation arise from Einstein coefficients and processes that broaden absorption lines. In practice, these quantities have been measured in the laboratory; not derived from theory.

When radiation is scattered (the phenomena that makes the sky appear blue) or when the fraction of molecules in an excited state is not determined by the Boltzmann distribution (and LTE doesn't exist), more complicated equations are required. For example, scattering from clear skies reflects about 32 W/m2 (about 13%) of incoming solar radiation back to space. Visible light is also reflected and scattered by aerosol particles and water droplets (clouds). Neither of these phenomena have a significant impact on the flux of thermal infrared through clear skies.

Schwarzschild's equation can not be used without first specifying the temperature, pressure, and composition of the medium through which radiation is traveling. When these parameters are first measured with a radiosonde, the observed spectrum of the downward flux of thermal infrared (DLR) agrees closely with calculations and varies dramatically with location. Where dI is negative, absorption is greater than emission, and net effect is to locally warm the atmosphere. Where dI is positive, the net effect is "radiative cooling". By repeated approximation, Schwarzschild's equation can be used to calculate the equilibrium temperature change caused by an increase in GHGs, but only in the upper atmosphere where heat transport by convection is unimportant.

## Derivation

Schwarzschild's equation can be derived from Kirchhoff's law of thermal radiation, which states that absorptivity must equal emissivity at a given wavelength. (Like Schwarzschild's equation, Kirchhoff's law only applies to media in LTE.) Given a thin slab of atmosphere of incremental thickness ds, by definition its absorptivity is dIa/I, where I is the incident radiation and dIa is radiation absorbed by the slab. According to Beer's Law:

$dI_{a}/I=n\sigma _{\lambda }\,ds$ Also by definition, emissivity is equal to dIe/Bλ(T), where dIe is the radiation emitted by the slab and Bλ(T) is the maximum radiation any object in LTE can emit. Setting absorptivity equal to emissivity affords:

$dI_{e}/B_{\lambda }(T)=n\sigma _{\lambda }\,ds$ $dI_{e}=n\sigma _{\lambda }B_{\lambda }(T)\,ds$ The total change in radiation, dI, passing through the slab is given by:

$dI=dI_{e}-dI_{a}=n\sigma _{\lambda }B_{\lambda }(T)\,ds-n\sigma _{\lambda }I\,ds$ Schwarzschild's equation has also been derived from Einstein coefficients by assuming a Maxwell–Boltzmann distribution of energy between a ground and excited state (LTE). The oscillator strength for any transition between ground and excited state depends on these coefficients. The absorption cross-section (σλ) is empirically determined from this oscillator strength and the broadening of the absorption/emission line by collisions, the Doppler effect and the uncertainty principle.

## Equivalent equations

Schwarzschild's equation has been expressed in different forms and symbols by different authors. The quantity λ is known as the absorption coefficient (βa), a measure of attenuation with units of cm−1. The absorption coefficient is fundamentally the product of a quantity of absorbers per unit volume (formally cm−3) times an efficiency of absorption (area/absorber, formally cm2). Several sources replace λ with kλr, where kλ is the absorption coefficient per unit density and r is the density of the gas. The absorption coefficient for spectral flux (a beam of radiation with a single wavelength, W/m2/μm) differs from the absorption coefficient for spectral intensity (W/sr/m2/μm) used in Schwarzschild's equation.

Integration of an absorption coefficient over a path from s1 and s2 affords the optical thickness ($\tau$ ) of that path, a dimensionless quantity that is used in some variants of the Schwarzschild equation. When emission is ignored, the incoming radiation is reduced by a factor for 1/e when transmitted over a path with an optical thickness of 1.

$\tau =\int _{s_{1}}^{s_{2}}n(s)\sigma _{\lambda }(s)\,ds$ When expressed in terms of optical thickness, Schwarzschild's equation becomes:

$dI_{\lambda }=[B_{\lambda }(T)-I_{\lambda }]\,d\tau$ After integrating between a sensor located at $\tau$ = 0 and an arbitrary starting point in the medium, $\tau$ ', the spectral intensity of the radiation reaching the sensor, Iλ(0), is:

$I_{\lambda }(0)=I_{\lambda }(\tau ')e^{-\tau '}+\int _{0}^{\tau '}B_{\lambda }(T)e^{-\tau }\,d\tau$ where I($\tau$ ') is the spectral intensity of the radiation at the beginning of the path, $e^{-\tau '}$ is the transmittance along the path, and the final term is the sum of all of the emission along the path attenuated by absorption along the path yet to be traveled.

## Relationship to Planck's and Beer's laws

Both Beer's Law and Planck's Law can be derived from Schwarzschild's equation. In a sense, they are corollaries of Schwarzschild's equation.

When the spectral intensity of radiation is not changing as it passes through a medium, dIλ = 0. In that situation, Schwarzschild's equation simplifies to Planck's law:

$0=n\sigma _{\lambda }[B_{\lambda }(T)-I_{\lambda }]\,ds$ $I_{\lambda }=B_{\lambda }(T)$ When Iλ > Bλ(T), dI is negative and when Iλ < Bλ(T), dI is positive. As a consequence, the intensity of radiation traveling through any medium is always approaching the blackbody intensity given by Planck's law and the local temperature. The rate of approach depends on the density of absorbing/emitting molecules (n) and their absorption cross-section (σλ).

When the intensity of the incoming radiation, Iλ, is much greater than the intensity of blackbody radiation, Bλ(T), the emission term can be neglected. This is usually the case when working with a laboratory spectrophotometer, where the sample is near 300 °K and the light source is a filament at several thousand °K.

$dI_{\lambda }=-n\sigma I_{\lambda }\,ds$ If the medium is homogeneous, λ doesn't vary with location. Integration over a path of length s affords the form of Beer's Law used most often in the laboratory experiments:

$dI_{\lambda }=I_{\lambda }(0)e^{-n\sigma _{\lambda }s}$ ## Origin of the greenhouse effect

Schwarzschild's equation provides a simple explanation for the existence of the greenhouse effect and demonstrates that it requires a non-zero lapse rate. Rising air in the atmosphere expands and cools as the pressure on it falls, producing a negative temperature gradient in the Earth's troposphere. When radiation travels upward through falling temperature, the incoming radiation, I, (emitted by the warmer surface or by GHGs at lower altitudes) is more intense than that emitted locally by Bλ(T). [Bλ(T) − I] is generally less than zero throughout the troposphere, and the intensity of outward radiation decreases as it travels upward. According to Schwarzschild's equation, the rate of fall in outward intensity is proportional to the density of GHGs (n) in the atmosphere and their absorption cross-sections (σλ). Any anthropogenic increase in GHGs will slow down the rate of radiative cooling to space, i.e. produce a radiative forcing until a saturation point is reached.

At steady state, incoming and outgoing radiation at the top of the atmosphere (TOA) must be equal. When the presence of GHGs in the atmosphere causes outward radiation to decrease with altitude, then the surface must be warmer than it would be without GHGs - assuming nothing else changed. Some scientists quantify the greenhouse effect as the 150 W/m2 difference between the average outward flux of thermal IR from the surface (390 W/m2) and the average outward flux at the TOA.

If the Earth had an isothermal atmosphere, Schwarzschild's equation predicts that there would be no greenhouse effect or no enhancement of the greenhouse effect by rising GHGs. In fact, the troposphere over the Antarctic plateau is nearly isothermal. Both observations and calculations show a slight "negative greenhouse effect" – more radiation emitted from the TOA than the surface. Although records are limited, the central Antarctic Plateau has seen little or no warming.

## Saturation

In the absence of thermal emission, wavelengths that are strongly absorbed by GHGs can be significantly attenuated within 10 m in the lower atmosphere. Those same wavelengths, however, are the ones where emission is also strongest. In an extreme case, roughly 90% of 667.5 cm−1 photons are absorbed within 1 meter by 400 ppm of CO2 at surface density, but they are replaced by emission of an equal number of 667.5 cm−1 photons. The radiation field thereby maintains the blackbody intensity appropriate for the local temperature. At equilibrium, Iλ = Bλ(T) and therefore dIλ = 0 even when the density of the GHG (n) increases. In other words, Schwarzschild's equation predicts no radiative forcing at wavelengths where absorption is "saturated". The radiative forcing from doubling carbon dioxide occurs mostly on the flanks of the strongest absorption band.

However, as density decreases with altitude, even the strongest absorption bands eventually become semi-transparent. Once that happens, radiation can travel far enough that the local emission, Bλ(T), can differ from the absorption of incoming Iλ. Depending on what altitude the transition to semi-transparency occurs and the temperature gradient, dI for a layer of atmosphere can be negative (ordinary greenhouse effect) or positive (negative greenhouse effect). Since the vast majority of absorption and emission occurs in the troposphere, where a negative temperature gradient exists, the Earth's greenhouse effect makes the planet warmer than it would be otherwise. However, temperature rises with altitude in the lower stratosphere, and increasing CO2 there increases radiative cooling to space and is predicted by some to cause cooling above 14–20 km.

## Application to climate science

Schwarzschild's equation is used to calculate the outward radiative flux from the Earth (measured in W/m2 perpendicular to the surface) at any altitude, especially the "top of the atmosphere" or TOA. This flux originates at the surface (I0) for clear skies or cloud tops. dI increments are calculated for layers thin enough to be effectively homogeneous in composition and flux (I). These increments are numerically integrated from the surface to the TOA to give the flux of thermal infrared to space, commonly referred to as outgoing long-wavelength radiation (OLR). OLR is the only mechanism by which the Earth gets rid of the heat delivered continuously by the sun. The net downward radiative flux of thermal IR (DLR) produced by emission from GHGs in the atmosphere is obtained by integrating dI from the TOA (where I0 is zero) to the surface. DLR adds to the energy from the sun. Emission from each layer adds equally to the upward and downward fluxes. In contrast, different amounts of radiation are absorbed, because the upward flux entering any layer is usually greater than the downward flux.

In "line-by-line" methods, the change in spectral intensity (dIλ, W/sr/m2/μm) is numerically integrated using a wavelength increment small enough (less than 1 nm) to accurately describe the shape of each absorption line. The HITRAN database contains the parameters needed to describe 7.4 million absorption lines for 47 GHGs and 120 isotopologues. A variety of programs or radiative transfer codes can be used to process this data, including an online facility, SpectralCalc. To reduce the computational demand, weather forecast and climate models use broad-band methods that handle many lines as a single "band". MODTRAN is a broad-band method available online with a simple interface that anyone can use.

To convert intensity (W/sr/m2) to flux (W/m2), calculations usually invoke the "two-stream" and "plane parallel" approximations. The radiative flux is decomposed into three components, upward (+z), downward (-z), and parallel to the surface. This third component contributes nothing to heating or cooling the planet. $ds=dz/\cos \theta$ , where $\theta$ is the zenith angle (away from vertical). Then the upward and downward intensities are integrated over a forward hemisphere, a process that can be simplified by using a "diffusivity factor" or "average effective zenith angle" of 53°. Alternatively, one can integrate over all possible paths from the entire surface to a sensor positioned a specified height above surface for OLR, or over all possible paths from the TOA to a sensor on the surface for DLR.