Difference between revisions of "NonCO2 Concentrations"
(Created page with "*model description ** carbon cycle ** nonCO2 concentrations **[[Radiative Forcing  radiative forcing r...") 
(No difference)

Revision as of 11:20, 17 June 2013
 model description
 carbon cycle
 nonCO2 concentrations
 radiative forcing routines
 climate model
 upwellingdiffusionentrainment implementation
Contents
NonCO_{2} concentrations
This section provides the formulas used to convert emissions to concentrations, while the Radiative Forcing section provides details on the derivation of radiative forcings.
Methane
Natural emissions of methane are inferred by balancing the budget for a userdefined historical period, e.g. from 19801990, so that
<m>\label{eq_natural emissions}E^n_{\o} = \theta (\Delta C_{\o}  C_{\o '}/\tau_{\rm tot})E^f_{\o}  E^b_{\o}</m>where <m>E^n_{\o}</m>, <m>E^f_{\o}</m> and <m>E^b_{\o}</m> are the average natural, fossil and land use related emissions, respectively; <m>\theta</m> is the conversion factor between atmospheric concentrations and mass loadings. <m>C_{\o'}</m> (and <m>\Delta C_{\o}</m>) are the average (annual changes in) concentrations. The net atmospheric lifetime <m>\tau_{\rm tot}</m> in the case of methane consists of the atmospheric chemical lifetime and lifetimes that characterize the soil and other (e.g. stratospheric) sink components according to
<m>\label{eq_methane lifetime} \frac{1}{\tau_{\rm tot}} = \frac{1}{\tau_{\rm tropos}} + \frac{1}{\tau_{\rm soil}} + \frac{1}{\tau_{\rm other}}</m>The feedback of methane on tropospheric OH and its own lifetime follows the results of the OxComp work (tropospheric oxidant model comparison) (see Ehhalt et al., 2001 in particular Table 4.11, which provides simple parameterizations for simulating complex threedimensional atmospheric chemistry models. As default, tropospheric OH abundances are assumed to decrease by 0.32 % for every 1 % increase in CH<m>_4</m>. The change in tropospheric OH abundances is thus modeled as:
<m>\label{eq_troposphericOH} \noindent \Delta {\rm ln}\,({\rm trop} {\rm OH}) = { S_{\rm CH_4}^{\rm OH}}\, \Delta{{\rm ln}\,(\rm CH_4)}+ { S_Template:\rm NO x^{\rm OH}} {E_Template:\rm NO x + S_{\rm CO}^{\rm OH} E_{\rm CO}} + {\rm S_{\rm VOC}^{\rm OH} E_{\rm VOC}}</m>where <m>S_x^{\rm OH}</m> is the sensitivity of tropospheric OH towards CH<m>_4</m>, NOx, CO and VOC, with default values of <m></m>0.32, +0.0042, <m></m>1.05e4 and <m></m>3.15e4, respectively. Increases in tropospheric OH abundances decrease the tropospheric lifetime <m>\tau '</m> of methane (default 9.6 yrs^{1}), which is approximated as a simple exponential relationship
<m>\label{eq_tropos_methane lifetime} \tau_{\rm CH_4,tropos}' = \tau_{\rm CH_4,tropos}^0\, {\rm exp}\,^{\Delta {\rm ln}\,({\rm trop OH})}</m>Approximating the temperature sensitivity of the net effect of tropospheric chemical reaction rates, the tropospheric lifetime of CH<m>_4</m> is adjusted:
<m>\label{eq_tropos_methane lifetime_inclTemp} \tau_{\rm CH_4,tropos} = \frac{\tau_{\rm CH_4,tropos}^0}{\frac{\tau_{\rm CH_4,tropos}^0}{\tau_{\rm CH_4,tropos}'} + S_{\tau_{\rm CH_4}}\Delta T}</m>where <m>S_{\tau_{\rm CH_4}}</m> is the temperature sensitivity coefficient (default <m>S_{\tau_{\rm CH_4}}</m>=3.16e2^{<m>^{\circ}</m>}C^{1}) and <m>\Delta T</m> is the temperature change above a userdefinable year, e.g. 1990.
Nitrous oxide
As for methane, natural nitrous oxide emissions are estimated by a budget (A28). For nitrous oxide however, the average concentrations <m>C_{\o'}</m>=<m>C_{\o3}</m> are taken for a period shifted by 3 years to account for a three year delay of transport of tropospheric N<m>_2</m>O to the main stratospheric sink. The feedback of the atmospheric burden C_{N<m>_2</m>O} of nitrous oxide on its own lifetime is approximated by:
<m>\tau_{\rm N_2O} = \tau_{\rm N_2O}^0 (\frac{{\rm C}_{\rm N_2O}}{{\rm C}_{\rm N_2O}^0})^{S_{\tau_{\rm N_2O}}}\\ \label{eq_nitrous_lifetime}</m>where <m>S_{\tau_{\rm N_2O}}</m> is the sensitivity coefficient (default <m>S_{\tau_{\rm N_2O}}</m>=<m></m>5e2) and the superscript ``^{<m>0</m>}´´ indicates a preindustrial reference state.
Tropospheric aerosols
Due to their short atmospheric residence time, changes in hemispheric abundances of aerosols are approximated by changes in their hemispheric emissions. Historical emissions of tropospheric aerosols are extended into the future either by emissions scenarios (SO<m>_{\rm x}</m>, NO<m>_{\rm x}</m>, CO) or, if scenario data are not available, with proxy emissions, e.g. using CO as a proxy emission for OC and BC. As with many other emissions scenarios, the harmonized IPCC SRES scenarios do not provide black (BC) and organic carbon (OC) emissions. Hence, various adhoc scaling approaches have been applied, often scaling BC and OC synchronously (Takemura et al., 2006), sometimes linearly with CO<m>_2</m> emissions. The MESSAGE emissions scenario modeling group is one of the few explicitly including BC and OC emissions in their multigas emissions scenarios (Rao et al, 2005; Rao and Riahi, 2006). By analyzing MESSAGE scenarios, a scaling factor was derived for this study in relation to carbon monoxide emissions (CO), varying linearly in time to 0.4 by 2100 relative to current BC/CO or OC/CO emission ratios.
Halogenated gases
The derivation of concentrations of halogenated gases controlled under either the Kyoto or Montreal Protocol assumes timevariable lifetimes. The net atmospheric lifetime <m>\tau_i</m> of each halogenated gas is calculated by summing the inverse lifetimes related to stratospheric, OHrelated and other sinks. Stratospheric lifetimes are assumed to decrease 15 % per degree of global mean surface temperature warming, due to an increased BrewerDobson circulation (Butchart and Scaife, 2001). Tropospheric OHrelated losses are scaled by parameterized changes in OHabundances, matching the respective changes in the lifetime of methane. The concentration <m>C_{t,i}</m> for the beginning of each year <m>t</m> is updated, using a central differencing formulation, according to:
<m>C_{t+1,i}{=}\tau_i E_{t,i} \frac{\rho_{\rm atm}}{m_{\rm atm}\mu_i}(1e^{\frac{1}{\tau_i}})+C_{t,i}(1e^{\frac{1}{\tau_i}})\label{eq_halogas_conc}</m>where <m>E_{t,i}</m> is the average emissions of gas <m>i</m> through year <m>t</m>, <m>C_{t,i}</m> the atmospheric concentration of gas <m>i</m> in year <m>t</m>, <m>\rho_{\rm atm}</m> the average density of air, <m>m_{\rm atm}</m> the total mass of the atmosphere (Trenberth and Guillemot, 1994), and <m>\mu_i</m> is the mass per mol of gas <m>i</m>. For hydrogenated halocarbons, the tropospheric OHrelated lifetimes are assumed to vary in proportion to the changes in methane lifetime.