**Influence of Greenhouse Gases to Global Warming on Account of Radiative Forcing**

Akira Tomizuka

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/58995

## **1. Introduction**

Radiative forcing is a measure of the size of a greenhouse gas's contribution to global warming. Radiative forcing values are estimated by a numerical process using radiative transfer schemes for terrestrial radiation and data from general circulation models. However, the estimation is complex and difficult to understand for non-specialists, including researchers in other fields. Understanding the essence of the Earth system is important for correctly discussing global environmental issues. Accordingly, in this chapter, radiative forcing values are calculated from a simple, intuitive radiative transfer model using the absorption spectra of greenhouse gases and the Planck formula for terrestrial radiation. [1]

#### **1.1. Radiative forcing**

The global atmospheric carbon dioxide (CO2) concentration has increased from about 278 ppmv in pre-industrial times (defined as 1750) to 390.5 ppmv in 2011. During the same period, the concentrations of methane (CH4) and nitrous oxide (N2O) have also increased from about 0.722 ppmv to 1.803 ppmv and about 0.270 ppmv to 0.324 ppmv, respectively (Table 1). [2]


**Table 1.** Concentrations in 1750 and 2011, radiative forcing, and global warming potential for each greenhouse gas. [2]

© 2015 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Radiative forcing is often referred to as an index of the size of a greenhouse gas's contribution to global warming. When the Earth system is at radiative equilibrium, the energy flux reaching the top of the Earth's atmosphere exactly balances with the outgoing energy flux from the Earth to outer space. However, increasing the concentration of greenhouse gases decreases the energy flux to outer space and changes the energy flux to the Earth into a surplus. Conse‐ quently, the Earth's surface temperature and atmospheric temperature rise, causing the outgoing flux to increase, and the Earth system shifts to a new equilibrium. Radiative forcing is defined as the imbalance of the energy flux density caused by these perturbations. The Intergovernmental Panel on Climate Change (IPCC) has estimated the following radiative forcing values due to increased greenhouse gas concentrations in 2011 relative to their preindustrial levels: CO2, 1.82 W m–2; CH4, 0.48 W m–2; and N2O, 0.17 W m–2. The unit of radiative forcing is the same as that of energy flux density. The uncertainties in these values are all ±10% under 90% confidence intervals. While the magnitude of the positive radiative forcing of greenhouse gases is well understood, the effects of other atmospheric constituents such as aerosols are subject to considerable uncertainty.[3] IPCC estimated values are the sum of the contributions from direct effect (via emissions of gases) and several indirect effects (via atmospheric chemistry). Radiative forcing to be compared with calculated value in this chapter is the direct contribution: CO2, 1.68 W m–2; CH4, 0.64 W m–2; and N2O, 0.17 W m–2.

Global warming potential (GWP) is also used in comprehensive policies regarding the regulation of greenhouse gases. GWP is a measure of how much a given mass of a greenhouse gas contributes to global warming and is usually defined as the radiative forcing resulting from an instantaneous release of 1 kg of the greenhouse gas into the atmosphere relative to that of CO2.[4] The GWP values for the next 20 and 100 years are given in Table 1.

Radiative forcing is estimated by a numerical process using radiative transfer schemes for terrestrial radiation and data from models referred to as general circulation models.[5] However, since the process is difficult to understand for non-specialists, including many citizens and researchers in other fields, they simply accept the results announced by the specialists. Yet, understanding the essence of the Earth system is important for correctly discussing global environmental issues. It is therefore necessary to create models such that anyone who has acquired basic scientific knowledge can intuitively understand the Earth system as well as the essence of the calculations based on the models.

In this chapter, radiative forcing values are calculated from a simple radiative transfer model using the absorption spectra of greenhouse gases and the Planck formula for terrestrial radiation. Furthermore, the GWPs of specific greenhouse gases are derived. Finally, the increase in the Earth's surface temperature due to radiative forcing is estimated.

#### **1.2. How can radiative forcing be calculated?**

The mean vertical temperature of the atmosphere results from the balance between heating and cooling. The Earth's surface and the troposphere are strongly coupled by convective heat transfer processes. At the surface, solar heating is balanced by convective transport of latent and sensible heat to the troposphere. In the troposphere, radiative cooling (infrared emission by molecules) is balanced by the release of latent heat via condensation and precipitation and by convective transport of sensible heat from the surface. This radiative-convective interaction leads to a roughly constant lapse rate in the troposphere. However, at the top of the troposphere (the tropopause), which is at an altitude of about 11 km above the Earth's surface, the tem‐ perature tends to become invariant with altitude.[6]

Radiative forcing is often referred to as an index of the size of a greenhouse gas's contribution to global warming. When the Earth system is at radiative equilibrium, the energy flux reaching the top of the Earth's atmosphere exactly balances with the outgoing energy flux from the Earth to outer space. However, increasing the concentration of greenhouse gases decreases the energy flux to outer space and changes the energy flux to the Earth into a surplus. Conse‐ quently, the Earth's surface temperature and atmospheric temperature rise, causing the outgoing flux to increase, and the Earth system shifts to a new equilibrium. Radiative forcing is defined as the imbalance of the energy flux density caused by these perturbations. The Intergovernmental Panel on Climate Change (IPCC) has estimated the following radiative forcing values due to increased greenhouse gas concentrations in 2011 relative to their preindustrial levels: CO2, 1.82 W m–2; CH4, 0.48 W m–2; and N2O, 0.17 W m–2. The unit of radiative forcing is the same as that of energy flux density. The uncertainties in these values are all ±10% under 90% confidence intervals. While the magnitude of the positive radiative forcing of greenhouse gases is well understood, the effects of other atmospheric constituents such as aerosols are subject to considerable uncertainty.[3] IPCC estimated values are the sum of the contributions from direct effect (via emissions of gases) and several indirect effects (via atmospheric chemistry). Radiative forcing to be compared with calculated value in this chapter

is the direct contribution: CO2, 1.68 W m–2; CH4, 0.64 W m–2; and N2O, 0.17 W m–2.

that of CO2.[4] The GWP values for the next 20 and 100 years are given in Table 1.

system as well as the essence of the calculations based on the models.

**1.2. How can radiative forcing be calculated?**

4 Global Warming - Causes, Impacts and Remedies

increase in the Earth's surface temperature due to radiative forcing is estimated.

Global warming potential (GWP) is also used in comprehensive policies regarding the regulation of greenhouse gases. GWP is a measure of how much a given mass of a greenhouse gas contributes to global warming and is usually defined as the radiative forcing resulting from an instantaneous release of 1 kg of the greenhouse gas into the atmosphere relative to

Radiative forcing is estimated by a numerical process using radiative transfer schemes for terrestrial radiation and data from models referred to as general circulation models.[5] However, since the process is difficult to understand for non-specialists, including many citizens and researchers in other fields, they simply accept the results announced by the specialists. Yet, understanding the essence of the Earth system is important for correctly discussing global environmental issues. It is therefore necessary to create models such that anyone who has acquired basic scientific knowledge can intuitively understand the Earth

In this chapter, radiative forcing values are calculated from a simple radiative transfer model using the absorption spectra of greenhouse gases and the Planck formula for terrestrial radiation. Furthermore, the GWPs of specific greenhouse gases are derived. Finally, the

The mean vertical temperature of the atmosphere results from the balance between heating and cooling. The Earth's surface and the troposphere are strongly coupled by convective heat transfer processes. At the surface, solar heating is balanced by convective transport of latent and sensible heat to the troposphere. In the troposphere, radiative cooling (infrared emission by molecules) is balanced by the release of latent heat via condensation and precipitation and

**Figure 1.** Vertical concentration profiles of four principal greenhouse gases. Concentrations are shown in units of ppmv. These data in the US Standard Atmosphere are obtained from the SpectralCalc website. [7]

The IPCC has defined radiative forcing as the change in net energy flux density at the tropo‐ pause. Figure 1 shows the altitude dependence of the concentrations of the four principal greenhouse gases. Water vapor is the most abundant and important greenhouse gas in the atmosphere. Nevertheless, it is excluded as an objective of the radiative forcing estimation since humans cannot directly control it. However, because of its strong absorption band, water vapor should be considered when calculating the radiative forcing of the other greenhouse gases.

The first step in building the model is to divide the atmosphere into appropriate layers, where the pressure, temperature, and concentration of each greenhouse gas are homogeneous within each layer. The concentrations of CO2, CH4, and N2O are homogeneous between the Earth's surface and the tropopause, whereas the concentration of water vapor changes considerably. Moreover, the temperature of the troposphere decreases with altitude at a roughly constant lapse rate. For intervals with a thickness of 100 m, the change in the concentration of water vapor between the adjacent layers is around 5% and the change in temperature is 0.65 K.

Therefore, up to an altitude of 11 km, the atmosphere is divided into 110 layers, each with a thickness of 100 m. The bottom layer tangent to the surface is referred to as the zeroth layer for descriptive purposes.

The terrestrial energy flux emitted from the Earth's surface enters the zeroth layer. Some of the incident flux is absorbed by the molecules of greenhouse gases in the layer, and the remainder is transmitted. Subsequently, the molecules emit radiation flux both upward and downward. Some of the combined upward flux (transmitted flux and emitted upward flux) is again absorbed in the first layer, and the layer emits radiation.[8] For simplicity, the molecules emit radiation only once. The repetition of this simple radiative transfer process leads to the outgoing flux from the tropopause to outer space.

In the first step, the outgoing energy flux density *<sup>F</sup>*¯ at the tropopause is calculated under the assumption that the concentrations of greenhouse gases in the atmosphere are equal to those in the pre-industrial era. In the next step, the flux density *Fi* is similarly calculated for an atmosphere in which the concentration of a specific greenhouse gas *i* increases up to its level in 2011, keeping everything else constant including the temperature.[4] The radiative forcing is obtained as *ΔFi* <sup>=</sup> *<sup>F</sup>*¯ <sup>−</sup> *Fi* for the greenhouse gas *i* that changes in concentration.

In these processes, the effects of the near-infrared region of incident solar radiation are ignored. Also, the intensities of absorption and emission depend on only the respective number densities of the greenhouse gases in the layer. Pressure and temperature affect these intensities through only a change in number density. Thus, for the same concentration of a given greenhouse gas, the ratio of the absorption or emission of an arbitrary layer to that of the zeroth layer is equal to the ratio of the respective number densities for that gas. When concentration is dependent on altitude, as in the case of water vapor, the number density ratio is multiplied by the factor of altitude dependence. The relation between number density and altitude is discussed in the following section.

## **2. Mathematical relation of number density of gas molecules**

Assume that the atmosphere within the troposphere is composed of an ideal gas with a density *ρ*(*z*) at an altitude *z* and a mean molecular mass *μ*. The pressure *p*(*z*) and the temperature *T*(*z*) are described using the temperature in the zeroth layer *T*0, the lapse rate in the troposphere Γ, and the universal gas constant *R* as follows:

$$p = \frac{R}{\mu} \rho T,\tag{1}$$

$$T = T\_0 + \Gamma \text{ z.} \tag{2}$$

Using the hydrostatic equation

$$\frac{dp}{dz} + \rho g = 0,\tag{3}$$

where *g* is the gravitational acceleration, we arrive at the following equation for the pressure:

$$\frac{dp}{dz} = -\frac{\mu g}{R(T\_0 + \Gamma z)} p.\tag{4}$$

Equation (4) can be easily integrated from a height of 0 to *z* to obtain:

The terrestrial energy flux emitted from the Earth's surface enters the zeroth layer. Some of the incident flux is absorbed by the molecules of greenhouse gases in the layer, and the remainder is transmitted. Subsequently, the molecules emit radiation flux both upward and downward. Some of the combined upward flux (transmitted flux and emitted upward flux) is again absorbed in the first layer, and the layer emits radiation.[8] For simplicity, the molecules emit radiation only once. The repetition of this simple radiative transfer process

In the first step, the outgoing energy flux density *<sup>F</sup>*¯ at the tropopause is calculated under the assumption that the concentrations of greenhouse gases in the atmosphere are equal to those

atmosphere in which the concentration of a specific greenhouse gas *i* increases up to its level in 2011, keeping everything else constant including the temperature.[4] The radiative forcing

In these processes, the effects of the near-infrared region of incident solar radiation are ignored. Also, the intensities of absorption and emission depend on only the respective number densities of the greenhouse gases in the layer. Pressure and temperature affect these intensities through only a change in number density. Thus, for the same concentration of a given greenhouse gas, the ratio of the absorption or emission of an arbitrary layer to that of the zeroth layer is equal to the ratio of the respective number densities for that gas. When concentration is dependent on altitude, as in the case of water vapor, the number density ratio is multiplied by the factor of altitude dependence. The relation between number density and altitude is

Assume that the atmosphere within the troposphere is composed of an ideal gas with a density *ρ*(*z*) at an altitude *z* and a mean molecular mass *μ*. The pressure *p*(*z*) and the temperature *T*(*z*) are described using the temperature in the zeroth layer *T*0, the lapse rate in the troposphere

<sup>=</sup> (1)

*TT z* = +G <sup>0</sup> . (2)

(3)

, *<sup>R</sup> p T*r

m

0, *dp <sup>g</sup> dz* + = r

**2. Mathematical relation of number density of gas molecules**

for the greenhouse gas *i* that changes in concentration.

is similarly calculated for an

leads to the outgoing flux from the tropopause to outer space.

in the pre-industrial era. In the next step, the flux density *Fi*

is obtained as *ΔFi* <sup>=</sup> *<sup>F</sup>*¯ <sup>−</sup> *Fi*

6 Global Warming - Causes, Impacts and Remedies

discussed in the following section.

Using the hydrostatic equation

Γ, and the universal gas constant *R* as follows:

$$\frac{p}{p\_0} = \exp\left[-\frac{\mu \text{g}}{R \,\Gamma} \ln(1 + \Gamma \, z \,/ \, T\_0)\right].\tag{5}$$

Again, using Eq. (1), the number density ratio of the atmosphere *N*(*z*) is as follows:

$$\frac{N}{N\_0} = \frac{\rho}{\rho\_0} = \frac{p}{p\_0} \frac{T\_0}{T\_0 + \Gamma z} = \frac{T\_0}{T\_0 + \Gamma z} \exp\left[ -\frac{\mu g}{R\Gamma} \ln(1 + \Gamma z / T\_0) \right]. \tag{6}$$

The subscript 0 denotes the value for the zeroth layer. The number density of a well-mixed greenhouse gas *n*(*z*) is related by the number density of the atmosphere:

$$m = \mathbb{C} \, N\_{\prime} \tag{7}$$

where *C* is the concentration (volume mixing ratio) of the greenhouse gas. For the three gases (CO2, CH4, and N2O) that have concentrations roughly independent of altitude, the ratio of number density is as follows:

$$\frac{m}{m\_0} = \frac{T\_0}{T\_0 + \Gamma \, z} \exp\left[-\frac{\mu g}{R\Gamma} \ln(1 + \Gamma \, z \,/\, T\_0)\right]. \tag{8}$$

In the case of water vapor where the concentration changes with the altitude, Eq. (8) is multiplied by the factor of altitude dependence.

#### **3. Simple radiative transfer**

To obtain the change of outgoing energy flux density at the tropopause, it is necessary to calculate the radiative transfer of terrestrial radiation.[9] For a wavelength in the range between *λ* and *λ*+*dλ*, the net change in flux density passing through a layer with a thickness Δ*z* is described as follows:

$$
\Delta I \, d\lambda = -\sigma\_\text{a} \mathfrak{n} \, \Delta z \, I \, d\lambda + \sigma\_\text{a} \mathfrak{n} \, \Delta z \, B \, d\lambda,\tag{9}
$$

where *I* (*λ*, *z*) is the intensity, which is defined as the amount of radiant energy leaving a unit area of a body per unit time per unit spectral interval *dλ*; in other words, intensity is the energy flux density per unit spectral interval.

The first term on the right-hand side of Eq. (9) uses the Beer–Lambert law of absorption; σa(*λ*) is the absorption cross section per unit greenhouse gas molecule and *n*(*z*) is the number density of the greenhouse gas in the layer. The second term is blackbody radiation based on the Kirchhoff law, which states that emittance and absorption have identical values. *B*(*λ*, *z*) is the intensity of the Planck blackbody function in the troposphere:

$$B = \frac{2\pi hc^2}{\mathcal{A}^5 \left( \exp\left[\frac{hc}{\mathcal{A}k(T\_0 + \Gamma z)}\right] - 1\right)} \,\mathrm{}^\prime\,\tag{10}$$

where *h* is the Planck constant, *k* is the Boltzmann constant, and *c* is the speed of light.

Scattering by molecules is ignored since the reference radiation is located in the infrared region. The absorption of the zeroth layer *a*0(*λ*) is defined as follows:

$$
\Delta a\_0 = \sigma\_\text{a} \eta\_0 \Delta \underline{z}.\tag{11}
$$

As a result, we can assume that the absorption of an arbitrary layer is described by *a*0(*λ*) and the ratio of number density is written as

$$
\sigma\_{\mathfrak{a}} \mathfrak{n} \,\Delta \mathfrak{z} = a\_0 \,\frac{\mathfrak{n}}{\mathfrak{n}\_0}.\tag{12}
$$

Therefore, the net change in flux density is expressed as a function of altitude as follows:

$$
\Delta I \, d\lambda = \frac{a\_0 n}{n\_0} (B - I) d\lambda = \frac{a\_0 T\_0}{T\_0 + \Gamma z} \exp\left[ -\frac{\mu g}{R \Gamma} \ln(1 + \Gamma z / T\_0) \right] (B - I) d\lambda. \tag{13}
$$

Taking into account the absorption by the four kinds of molecules, the actual change of the flux density passing through the layer is

$$
\Delta Id \lambda = \left( w\_0 f\_\text{w} + \sum\_{i=1}^3 a\_0^i \right) \frac{T\_0}{T\_0 + \Gamma \ z} \exp\left[ -\frac{\mu g}{R \Gamma} \ln(1 + \Gamma \ z / T\_0) \right] (B - I) d\lambda,\tag{14}
$$

where *a*<sup>0</sup> *i* and *w*0 represent the absorption spectra of anthropogenic greenhouse gas and that of water vapor of the zeroth layer, respectively. *f*w (*z*) is the factor of altitude dependence for water vapor, which is obtained by normalizing the concentration of an arbitrary layer by 7750 ppmv, which is the value for the zeroth layer.

The absorption spectra of greenhouse gases are discussed in the following section. The outgoing flux density at the tropopause with a wavelength from *λ* to *λ*+*dλ* can be obtained by the repeated use of Eqs. (14) and (10) with an initial value of *I*(*λ*, 0)=*B*(*λ*, 0). They are summed up for the reference range of wavelengths to give the total flux density. Of course, the total absorption through any layers must not be greater than one.

#### **3.1. Estimation of radiative forcing**

where *I* (*λ*, *z*) is the intensity, which is defined as the amount of radiant energy leaving a unit area of a body per unit time per unit spectral interval *dλ*; in other words, intensity is the energy

The first term on the right-hand side of Eq. (9) uses the Beer–Lambert law of absorption; σa(*λ*) is the absorption cross section per unit greenhouse gas molecule and *n*(*z*) is the number density of the greenhouse gas in the layer. The second term is blackbody radiation based on the Kirchhoff law, which states that emittance and absorption have identical values. *B*(*λ*, *z*) is the

2

<sup>2</sup> ,

(10)

0

Scattering by molecules is ignored since the reference radiation is located in the infrared region.

As a result, we can assume that the absorption of an arbitrary layer is described by *a*0(*λ*) and

0 . *n*

m

Taking into account the absorption by the four kinds of molecules, the actual change of the

m

exp ln(1 / ) ( ) , *<sup>i</sup>*

+G G è ø ë û <sup>å</sup> (14)

*n*

. (11)

D = (12)

0

+G G ë û (13)

l

 l

exp 1 ( )

p

<sup>=</sup> æ ö é ù ç ÷ ê ú - + G è ø ë û

where *h* is the Planck constant, *k* is the Boltzmann constant, and *c* is the speed of light.

*a nz* 0 a0 = D s

a 0

Therefore, the net change in flux density is expressed as a function of altitude as follows:

( ) exp ln(1 / ) ( ) . *a n a T <sup>g</sup> Id B I d z T B Id*

é ù D= - = ê ú - +G -

0w 0 0

æ ö é ù D= + ç ÷ ê ú - +G -

*<sup>T</sup> <sup>g</sup> Id w f a z T B Id Tz R*

*nz a*

s

l

*hc kT z*

flux density per unit spectral interval.

8 Global Warming - Causes, Impacts and Remedies

intensity of the Planck blackbody function in the troposphere:

5

l

The absorption of the zeroth layer *a*0(*λ*) is defined as follows:

0 0 0

 l

*n Tz R*

0

0 0

3

*i*

=

1 0

the ratio of number density is written as

l

l

flux density passing through the layer is

*hc <sup>B</sup>*

The calculations of radiative forcing require the absorption spectra of the four greenhouse gases for the zeroth layer in 1750.[10]

**Figure 2.** Absorption spectra of atmospheric greenhouse gases in 1750 for thickness of 100 m obtained from the Spec‐ tralCalc website: (a) water vapor, 7750 ppmv; (b) CO2, 278 ppmv; (c) CH4, 0.722 ppmv; and (d) N2O, 0.270 ppmv. The concentration of water vapor is assumed to be the same as the present value. Absorption is the intensity ratio of absor‐ bed radiation to total radiation incident on the zeroth layer.

Figure 2 shows the absorption spectra of the zeroth layer for 1750, as calculated using the SpectralCalc database. For selected gases, the SpectralCalc website provides transmittance spectra at an arbitrary pressure, temperature, and thickness of a gas layer; the concentration of the gas in the layer; and the range of wavelengths.[7] The SpectralCalc uses a line-by-line model called LINEPAK to accurately model molecular absorption line spectra.[11] These are based on the HITRAN database,[12] which is a compilation of spectroscopic parameters widely used to simulate the gases' transmission and emission of radiation into the atmosphere. The concentrations are 7750 ppmv for water vapor, 278 ppmv for CO2, 0.722 ppmv for CH4, and 0.270 ppmv for N2O. The pre-industrial concentration of water vapor is assumed to be the same as the present value.

In the zeroth layer, the parameters are

$$T\_0 = 288\,\text{K}, \quad p\_0 = 1.013 \times 10^5 \,\text{Pa}, \quad \Delta z = 100\,\text{m}, \text{ and } \quad \lambda = 1.85 - 20.0\,\mu\text{m}.\tag{15}$$

The unit spectral interval is set to 1 nm. The numerical values used for the calculations are as follows:

$$\begin{aligned} \text{R = } & \text{8.3144 J/K}^{-1}, \quad \mu = 28.964 \times 10^{-3} \text{kg}, \quad \Gamma = -6.5 \times 10^{-3} \text{K m}^{-1}, \quad g = 9.8 \,\text{m} \,\text{s}^{-2}, \\\ h = 6.6261 \times 10^{-34} \text{Js}, \quad k = 1.3806 \times 10^{-23} \text{J K}^{-1}, \quad c = 2.9979 \times 10^{8} \text{m s}^{-1}. \end{aligned} \tag{16}$$

Figure 3 shows the calculated outgoing spectra in the pre-industrial era at altitudes of 3 km, 6 km, 9 km, and 11 km (the tropopause). It is clear that terrestrial radiation is strongly absorbed at wavelengths of 5–8 μm and 13–17 μm. The former is due to water vapor and the latter is due to CO2. In both these ranges, the absorption is nearly one. However, the outgoing flux in these ranges is not zero because of emission from the greenhouse gases. Moreover, the intensity of the flux decreases with altitude. The red area under the spectrum at 11 km, 191.60 W m–2, corresponds to the total outgoing flux density *<sup>F</sup>*¯.

The growth rates of greenhouse gas concentrations from the pre-industrial era to 2011 are given as follows: CO2, 1.40; CH4, 2.50; and N2O, 1.20. For simplicity, we assume that the absorption spectrum in 2011 may be approximated by multiplying the spectrum in 1750 by the growth rate. The absorption values never exceed one.

A similar calculation in which the absorption spectrum of only 278 ppmv CO2 is replaced with that of 390.5 ppmv CO2 gives *F*CO2=189.26 W m–2. Therefore, the radiative forcing of CO2 for the period between 1750 and 2011 is

$$
\Delta F\_{\rm CO2} = \overline{F} - F\_{\rm CO2} = 2.34 \,\mathrm{W} \,\mathrm{m}^{-2}.\tag{17}
$$

Similarly, the outgoing flux density *Fi* for the spectrum of greenhouse gas *i* when its concen‐ tration in 1750 is replace with its concentration in 2011, is *F*CH4=190.74 W m–2 and *F*N2O=191.41 Influence of Greenhouse Gases to Global Warming on Account of Radiative Forcing http://dx.doi.org/10.5772/58995 11

Figure 2 shows the absorption spectra of the zeroth layer for 1750, as calculated using the SpectralCalc database. For selected gases, the SpectralCalc website provides transmittance spectra at an arbitrary pressure, temperature, and thickness of a gas layer; the concentration of the gas in the layer; and the range of wavelengths.[7] The SpectralCalc uses a line-by-line model called LINEPAK to accurately model molecular absorption line spectra.[11] These are based on the HITRAN database,[12] which is a compilation of spectroscopic parameters widely used to simulate the gases' transmission and emission of radiation into the atmosphere. The concentrations are 7750 ppmv for water vapor, 278 ppmv for CO2, 0.722 ppmv for CH4, and 0.270 ppmv for N2O. The pre-industrial concentration of water vapor is assumed to be the same

as the present value.

follows:

In the zeroth layer, the parameters are

10 Global Warming - Causes, Impacts and Remedies

5

0 0 *Tp z* = = ´ D= 288K, 1.013 10 Pa, 100m, and 1.85 20.0 m.

The unit spectral interval is set to 1 nm. The numerical values used for the calculations are as

34 23 1 8 1 8.3144 JK , 28.964 10 kg, 6.5 10 K m , 9.8m s ,

Figure 3 shows the calculated outgoing spectra in the pre-industrial era at altitudes of 3 km, 6 km, 9 km, and 11 km (the tropopause). It is clear that terrestrial radiation is strongly absorbed at wavelengths of 5–8 μm and 13–17 μm. The former is due to water vapor and the latter is due to CO2. In both these ranges, the absorption is nearly one. However, the outgoing flux in these ranges is not zero because of emission from the greenhouse gases. Moreover, the intensity of the flux decreases with altitude. The red area under the spectrum at 11 km, 191.60 W m–2,

The growth rates of greenhouse gas concentrations from the pre-industrial era to 2011 are given as follows: CO2, 1.40; CH4, 2.50; and N2O, 1.20. For simplicity, we assume that the absorption spectrum in 2011 may be approximated by multiplying the spectrum in 1750 by the growth

A similar calculation in which the absorption spectrum of only 278 ppmv CO2 is replaced with that of 390.5 ppmv CO2 gives *F*CO2=189.26 W m–2. Therefore, the radiative forcing of CO2 for the

Similarly, the outgoing flux density *Fi* for the spectrum of greenhouse gas *i* when its concen‐ tration in 1750 is replace with its concentration in 2011, is *F*CH4=190.74 W m–2 and *F*N2O=191.41

2 CO2 CO2 *F FF* 2.34W m . - D =- = (17)


6.6261 10 Js, 1.3806 10 JK , 2.9979 10 m s .

= = ´ G=- ´ =

*R g*

*hk c*

m

corresponds to the total outgoing flux density *<sup>F</sup>*¯.

rate. The absorption values never exceed one.

period between 1750 and 2011 is

1 3 3 1 2

=´ =´ =´ (16)


l

 m= - (15)

**Figure 3.** Outgoing terrestrial radiation spectrum at 3 km (blue), 6 km (green), 9 km (yellow), and 11 km (tropopause; red) in 1750, as calculated using the spectra shown in Figure 2.

W m–2. The radiative forcing of each greenhouse gas is listed in Table 2. Despite of the simplicity of the radiative transfer model used here, the calculated values of Δ*F* are close to those in the IPCC estimations, but are relatively higher: CO2, 39%; CH4, 33%; and N2O, 12%. The 2011 spectrum used in this calculation was approximated by multiplying the greenhouse gas's growth rate by the spectrum in 1750. Nevertheless, the differences between the approximated radiative transfer values and the precise values based on the real spectrum in 2011 in Spec‐ tralCalc are very small except for CO2: CO2, 12%; CH4, 1.4%; and N2O, 0.5%.


\* These estimates are direct contributions through the emissions of the gases.

**Table 2.** Calculated outgoing flux density in 1750 and 2011, radiative forcing Δ*F*, and IPCC estimates for each greenhouse gas.

Figure 4 shows the distribution of radiative forcing. We can see that the saturated absorption regions, 5–7 μm for water vapor and 16–18 μm for CO2, do not affect the radiative forcing, whereas the unsaturated area does. Nevertheless, the CO2 absorption in the range 9–11 μm is negligibly small (Figure 2(b)), contributing 10% of the radiative forcing because terrestrial radiation is maximal in this range.

The robustness of the proposed model must be confirmed by testing the sensitivity of the results to changes in certain parameters. A change in water vapor concentration does not affect the radiative forcing considerably. The increase in radiative forcing is around 1% even when water vapor concentration decreases by 5%. Moreover, the outgoing flux density decreases as the lapse rate increases. For a change of 10% in the lapse rate, the radiative forcing values change between about 3% and 6%.

**Figure 4.** Radiative forcing spectra calculated with the precise absorption spectra: CO2 (red), CH4 (green), and N2O (blue).

There are some sources of the uncertainties in this calculation. The first is assuming that the absorption spectra of greenhouse gases in each layer are proportional to only their number density. However, the absorption spectrum also changes with pressure and temperature in the layer: absorption line width broadening is caused by the thermal motion of the molecules and the collisions between them, both of which depend on pressure and temperature.[13,14]

This effect is captured by *a*0(*λ*) provided by SpectralCalc. However, the effects of the other layer are approximated as the effects of the zeroth layer for simplicity. These differences in the spectra in each layer will produce uncertainties in the outgoing fluxes.

The second source of uncertainty is ignoring the radiative transfer of incident solar radiation in the near-infrared region. The absorption and emission at 2.6–3.5 μm may not be negligible since the solar intensity in the region is 5–8 W m–2μm–1 (the maximum intensity of the terrestrial radiation is about 25 W m–2 μm–1, as shown in Figure 3). Therefore, the outgoing flux contains some uncertainties.

Figure 5 shows the calculated radiative forcing for each increasing greenhouse gas from the pre-industrial era to the present age. This result implies that each radiative forcing becomes slightly rounded due to saturation over the specific wavelength areas and it is linearly approximated by the gas's concentration: the radiative forcing per ppmv of CO2 is 0.025 W m–2; CH4, 0.7 W m–2, and N2O, 3.5 W m–2.

**Figure 5.** Calculated radiative forcings for increasing concentrations of greenhouse gases: (a) CO2, for every 20 ppmv from 280 ppmv to 380 ppmv; (b) CH4, for every 0.1 ppmv from 0.7 ppmv to 1.7 ppmv; (c) N2O, for every 0.005 ppmv from 0.270 ppmv to 0.320 ppmv.

#### **3.2. Estimation of global warming potential**

Figure 4 shows the distribution of radiative forcing. We can see that the saturated absorption regions, 5–7 μm for water vapor and 16–18 μm for CO2, do not affect the radiative forcing, whereas the unsaturated area does. Nevertheless, the CO2 absorption in the range 9–11 μm is negligibly small (Figure 2(b)), contributing 10% of the radiative forcing because terrestrial

The robustness of the proposed model must be confirmed by testing the sensitivity of the results to changes in certain parameters. A change in water vapor concentration does not affect the radiative forcing considerably. The increase in radiative forcing is around 1% even when water vapor concentration decreases by 5%. Moreover, the outgoing flux density decreases as the lapse rate increases. For a change of 10% in the lapse rate, the radiative forcing values

**Figure 4.** Radiative forcing spectra calculated with the precise absorption spectra: CO2 (red), CH4 (green), and N2O

There are some sources of the uncertainties in this calculation. The first is assuming that the absorption spectra of greenhouse gases in each layer are proportional to only their number density. However, the absorption spectrum also changes with pressure and temperature in the layer: absorption line width broadening is caused by the thermal motion of the molecules and the collisions between them, both of which depend on pressure and temperature.[13,14]

This effect is captured by *a*0(*λ*) provided by SpectralCalc. However, the effects of the other layer are approximated as the effects of the zeroth layer for simplicity. These differences in the

The second source of uncertainty is ignoring the radiative transfer of incident solar radiation in the near-infrared region. The absorption and emission at 2.6–3.5 μm may not be negligible since the solar intensity in the region is 5–8 W m–2μm–1 (the maximum intensity of the terrestrial

spectra in each layer will produce uncertainties in the outgoing fluxes.

radiation is maximal in this range.

12 Global Warming - Causes, Impacts and Remedies

change between about 3% and 6%.

(blue).

GWP is defined as the radiative forcing per unit mass of released greenhouse gas *i* relative to that of CO2. It is necessary to integrate the radiative forcing for a reference time since the number of gas molecules decreases with time. The Kyoto Protocol is based on GWPs from pulse emissions over a 100-year time frame. For simplicity, the GWP of greenhouse gas *i* for its 100-year effect is approximated as

$$\text{GWP}\_{i} = \frac{\langle \Delta F \, / \Delta n \rangle\_{i} \, / \, M\_{i} \times \int\_{0}^{100} e^{-t/\tau\_{i}} dt}{\langle \Delta F \, / \Delta n \rangle\_{\text{CO2}} \, / \, M\_{\text{CO2}} \times \int\_{0}^{100} e^{-t/\tau\_{\text{CO2}}} dt} \, \tag{18}$$

where Δ*F* / Δ*n* is the radiative forcing per unit increased concentration and *M* is the molecular mass of the greenhouse gas. The mean residual life of a greenhouse gas over 100 years is given by

$$\int\_0^{100} e^{-t/\tau} dt \,\tag{19}$$

where *τ* is its lifetime.

The lifetimes of CH4 and N2O are 12.4 and 121 years, respectively, according to the IPCC Fifth Assessment Report. [2] For CO2, the mean residual life is given as follows:

$$\int\_0^{100} (0.217 + 0.259 \, e^{-t/172.9} + 0.338 \, e^{-t/18.51}) dt. \tag{20}$$

The GWPs of CH4 and N2O for the 100-year effect are calculated as 27.24 and 242.2, respectively, using Δ*F* as obtained in Section 3.1. They are listed in Table 3 along with the IPCC estimates. The GWP of CH4 is very close to the IPCC estimate whereas the radiative forcing value is not.


\* For carbon dioxide, the IPCC uses the parameters in Eq. (20) which represent the mean residual life, without defining its lifetime. [15]

**Table 3.** Increased concentration Δ*n* from 1750 to 2011, lifetime *τ*, calculated GWP, and IPCC estimate of GWP over an interval of 100 years for each greenhouse gas.

#### **3.3. Radiative forcing and surface temperature**

Finally, we estimate the increase in the Earth's temperature by using radiative forcing. In the pre-industrial era, the incoming solar radiation flux at the tropopause exactly balanced the outgoing solar radiation flux reflected by the Earth system and the outgoing terrestrial radiation.

Let the outgoing terrestrial radiation be approximated roughly as

$$
\overline{F} = (1 - f)\sigma T\_0^4 \,\tag{21}
$$

where *f* is the mean absorption of the terrestrial radiation at the tropopause, *T*<sup>0</sup> is the surface temperature and σ is the Stephan–Boltzmann constant.

Increasing abundance of a greenhouse gas (e.g., CO2) decreases the mean absorption by Δ*f*.

$$F\_{\rm CO2} = \left(1 - f - \Delta f\right) \sigma T\_0^4. \tag{22}$$

The perturbation Δ*f* is related to radiative forcing by using Eqs. (21) and (22):

where Δ*F* / Δ*n* is the radiative forcing per unit increased concentration and *M* is the molecular mass of the greenhouse gas. The mean residual life of a greenhouse gas over 100 years is given

The lifetimes of CH4 and N2O are 12.4 and 121 years, respectively, according to the IPCC Fifth

The GWPs of CH4 and N2O for the 100-year effect are calculated as 27.24 and 242.2, respectively, using Δ*F* as obtained in Section 3.1. They are listed in Table 3 along with the IPCC estimates. The GWP of CH4 is very close to the IPCC estimate whereas the radiative forcing value is not.

Carbon dioxide 112 \* 1 1 Methane 1.08 12.4 27.2 28 Nitrous oxide 0.054 121 242 265 \* For carbon dioxide, the IPCC uses the parameters in Eq. (20) which represent the mean residual life, without defining

**Table 3.** Increased concentration Δ*n* from 1750 to 2011, lifetime *τ*, calculated GWP, and IPCC estimate of GWP over an

Finally, we estimate the increase in the Earth's temperature by using radiative forcing. In the pre-industrial era, the incoming solar radiation flux at the tropopause exactly balanced the outgoing solar radiation flux reflected by the Earth system and the outgoing terrestrial

> <sup>4</sup> *F fT* = - (1 ) , s

where *f* is the mean absorption of the terrestrial radiation at the tropopause, *T*<sup>0</sup> is the surface

<sup>0</sup> (0.217 0.259 0.338 ) . [15] *t t <sup>e</sup> e dt* - - + + ò (20)

ò (19)

**GWP This work IPCC**

<sup>0</sup> (21)

100 / <sup>0</sup> , *<sup>t</sup> e dt* t

Assessment Report. [2] For CO2, the mean residual life is given as follows:

**Greenhouse gas Δn (ppmv)** *τ (years)*

Let the outgoing terrestrial radiation be approximated roughly as

temperature and σ is the Stephan–Boltzmann constant.

100 /172.9 /18.51

by

where *τ* is its lifetime.

14 Global Warming - Causes, Impacts and Remedies

its lifetime. [15]

radiation.

interval of 100 years for each greenhouse gas.

**3.3. Radiative forcing and surface temperature**

$$
\Delta f = \frac{\Delta \text{FCO2}}{\sigma T\_0^4}.\tag{23}
$$

Assuming Δ*f* is maintained for some time, a new equilibrium state will develop in which the surface temperature has increased by Δ*T*. At the new radiative equilibrium, the outgoing terrestrial flux at the tropopause does not change if the incoming solar flux and the outgoing solar flux reflected by the Earth system are unchanged. Therefore,

$$(1 - f)\sigma T\_0^4 = (1 - f - \Delta f)\sigma (T\_0 + \Delta T)^4. \tag{24}$$

Ignoring the second-order and higher terms in Eq. (24), we obtain

$$
\Delta T = \frac{\Delta f}{4(1-f)} T\_0. \tag{25}
$$

Substituting Eqs. (21) and (23) into Eq. (25) gives the relationship between Δ*T* and Δ*F*CO2:

$$
\Delta T = \frac{\Delta F\_{\text{CO2}}}{4\overline{F}} T\_0. \tag{26}
$$

By substituting the numerical values in Table 2, we can estimate that the Earth's surface temperature has risen by Δ*T*=0.88 K as a result of the CO2 increase since 1750. Then, in the new equilibrium state, the surface temperature is *T*0′=*T*0+Δ*T*=288.9 K and the lapse rate should be –6.6 × 10–3 K m–1 due to the temperature invariance at the tropopause. Next, the outgoing flux density at the tropopause with the surface temperature *T*0′ is calculated. This value is close to *<sup>F</sup>*¯.

The radiative forcing of CO2 is dominant, and the radiative forcings of other greenhouse gases nearly balance with negative radiative forcings (e.g., of aerosols). As a result, the influence of all radiative forcings is almost equivalent to the influence of CO2.

Figure 6 shows the surface temperature increase due to CO2 considering radiative forcing when the atmospheric CO2 concentration increases to twice that in the pre-industrial era. In the future (when the CO2 concentration exceeds 390 ppmv), the surface temperature increase per 20 ppmv CO2 is about 0.1 K.

**Figure 6.** Surface temperature increase due to the radiative forcing of CO2 when the atmospheric CO2 concentration increases to twice that in the pre-industrial era.

## **4. Conclusions**

The one-dimensional radiative equilibrium model can estimate the radiative forcing and global warming potential of anthropogenic greenhouse gases by using their absorption spectra at Earth's surface and number densities. Under some approximations, radiative forcing values are calculated higher than the IPCC estimates; CO2, 39%; CH4, 33%; and N2O, 12%. Calculated values of global warming potential are very similar to IPCC estimates. Furthermore, increases in Earth's surface temperature due to radiative forcing can be easily estimated by the model.

## **Author details**

#### Akira Tomizuka

Graduate School of Fisheries Science and Environmental Studies, Nagasaki University, Japan

## **References**

(when the CO2 concentration exceeds 390 ppmv), the surface temperature increase per 20 ppmv

**Figure 6.** Surface temperature increase due to the radiative forcing of CO2 when the atmospheric CO2 concentration

The one-dimensional radiative equilibrium model can estimate the radiative forcing and global warming potential of anthropogenic greenhouse gases by using their absorption spectra at Earth's surface and number densities. Under some approximations, radiative forcing values are calculated higher than the IPCC estimates; CO2, 39%; CH4, 33%; and N2O, 12%. Calculated values of global warming potential are very similar to IPCC estimates. Furthermore, increases in Earth's surface temperature due to radiative forcing can be easily

Graduate School of Fisheries Science and Environmental Studies, Nagasaki University, Japan

CO2 is about 0.1 K.

16 Global Warming - Causes, Impacts and Remedies

increases to twice that in the pre-industrial era.

**4. Conclusions**

estimated by the model.

**Author details**

Akira Tomizuka


K. B. Averyt, M. Tignor, and H. L. Miller (Cambridge University Press, Cambridge, 2007).

**Impacts/Threats/Consequences of Global Warming**

K. B. Averyt, M. Tignor, and H. L. Miller (Cambridge University Press, Cambridge,

2007).

18 Global Warming - Causes, Impacts and Remedies

**Chapter 2**
