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

Revision as of 12:21, 17 June 2013
 model description
 carbon cycle
 nonCO2 concentrations
 radiative forcing routines
 climate model
 upwellingdiffusionentrainment implementation
Contents
Radiative forcing
The following section highlights the key parameterizations used for estimating the radiative forcing due to humaninduced changes in greenhouse gas concentrations, tropospheric ozone and aerosols. The radiative forcing applied in MAGICC is in general the forcing at tropopause level after stratospheric temperature adjustment. Efficacies of the forcings, as discussed by Hansen et al. (2005) and Meehl et al. (2007) can be applied.
Carbon dioxide
Taking into account the "saturation" effect of CO_{2} forcing, i.e., the decreasing forcing efficiency for a unit increases of CO_{2} concentrations with higher background concentrations, the first IPCC Assessment (Shine et al., 1990) presented the simplified expression of the form:
<m>\label{eq_CO2_forcing} \Delta Q_{\rm CO_2}=\alpha_{\rm CO_2}{\rm ln} (\rm C/C_0)</m>where <m>\Delta Q_{\rm CO_2}</m> is the adjusted radiative forcing by CO_{2} (Wm^{2}) for a CO_{2} concentration <m>C</m> (ppm) above the preindustrial concentration <m>C</m>_{0} (278 ppm). This expression proved to be a good approximation, although the scaling parameter <m>\alpha_{\rm CO_2}</m> has since been updated to a bestestimate of 5.35 Wm^{2} (=<m>\frac{3.71}{ln(2)} </m>Wm^{2})Myhre et al. (1998), used as default in MAGICC. When applying AOGCMspecific CO_{2} forcing, <m>\alpha_{\rm CO_2}</m> is set to:
<m>\label{eq_CO2_forcing_dQ2x} \alpha_{\rm CO_2} = \frac{\Delta Q_{\rm 2\times}}{\rm ln(2)}</m>Methane and nitrous oxide
Methane and nitrous oxide have overlapping absorption bands so that higher concentrations of one gas will reduce the effective absorption by the other and vice versa. This is reflected in the standard simplified expression for methane and nitrous oxide forcing, <m>\Delta Q_{\rm CH4}</m> and <m>\Delta Q_{\rm N2O}</m>, respectively (see Ramaswamy et al., 2001, Myhre et al. 1998):
<m>\label{eq_methane_forcing}\Delta Q_{\rm CH_4} = \alpha_{\rm CH_4}(\sqrt{\rm C_{\rm CH_4}}\sqrt{\rm C_{\rm CH_4}^0}f(\rm C_{\rm CH_4},{\rm C}_{\rm N_2O}^0)f({\rm C}_{\rm CH_4}^0,{\rm C}_{\rm N_2O}^0)</m>
where the overlap is captured by the function
<m>\label{eq_methane_forcing_overlap}f(\rm M,N)=0.47{\rm ln}\,(1+0.6356(\frac{\rm MN}{10^6})^{0.75}+0.007\frac{\rm M}{10^3}(\frac{\rm MN}{10^6})^{1.52})</m>with M and N being CH_{4} and N_{2}O concentrations in ppb. For methane, an additional forcing factor due to methaneinduced enhancement of stratospheric water vapor content is included. This enhancement is assumed to be proportional to (default <m>\beta</m>=15 %) the ``pure´´ methane radiative forcing, i.e., without subtraction of N_{2}O absorption band overlaps:
<m>\Delta Q_{\rm CH_4}^{\rm stratoH2O} =\beta\alpha_{\rm CH_4}(\sqrt{\rm C_{\rm CH_4}}\sqrt{\rm C_{\rm CH_4}^0}). </m>Tropospheric ozone
From the tropospheric ozone precursor emissions and following the updated parameterizations of OxComp as given in footnote a of Table 4.11 in Ehhalt et al. (2001), the change in hemispheric tropospheric ozone concentrations (in DU) is parameterized as:
<m>\label{eq_tropospheric_ozone} \Delta (\rm trop O_3)={\rm S}_{\rm CH_4}^{O_3}\Delta ln(\rm CH_4)+ {\rm S}_{\rm NOx}^{\rm O_3} E_{\rm NOx} + {\rm S}_{\rm CO}^{\rm O_3} E_{\rm CO} + {\rm S}_{\rm VOC}^{\rm O_3} E_{\rm VOC}</m>where <m>S_{\rm x}^{\rm O_3}</m> are the respective sensitivity coefficients of tropospheric ozone to methane concentrations and precursor emissions. The radiative forcing is then approximated by a linear abundance to forcing relationship so that <m>\Delta Q_{\rm trop O_3}=\alpha_{\rm trop O3} \Delta (\rm trop O_3)</m> with <m>\alpha_{\rm trop O3}</m> being the radiative efficiency factor (default 0.042).
Halogenated gases
The globalmean radiative forcing <m>\Delta Q_{t,i}</m> of halogenated gases is simply derived from their atmospheric concentrations C (see NonCO2 concentrations) and radiative efficiencies <m>\varrho_i</m> Ehhalt et al. (2001) table 4.11.
<m>\Delta {Q_{t,i}} = \varrho_i (C_{t,i}C_{0,i})\label{eq_halogas_RF}</m>The landocean forcing contrast in each hemisphere for halogenated gases is assumed to follow the one Hansen et al. (2005) estimated for CFC11. The hemispheric forcing contrast is dependent on the lifetime of the gas. For shortlived gases (<m>{<}</m>1\,yr) the hemispheric forcing contrast is assumed to equal the timevariable hemispheric emission ratio. For longer lived gases (default <m>{>}</m>8\,yrs), the hemispheric forcing contrast is assumed to equal the one from CFC11 with linear scaling in between these two approaches for gases with a medium lifetime.
Stratospheric ozone
Depletion of the stratospheric ozone layer causes a negative globalmean radiative forcing <m>\Delta Q_{t}</m>. The depletion and hence radiative forcing is assumed to be dependent on the equivalent effective stratospheric chlorine (EESC) concentrations as follows:
<m>\Delta Q_{t} = \eta_1 (\eta_2 \times \Delta {\rm EESC}_{t})^{\eta_3} \label{eq_stratospheric_ozone} </m>where <m>\eta_1</m> is a sensitivity scaling factor (default <m></m>4.49e4 Wm^{2}), <m>\Delta {\rm EESC}_{t}</m> the EESC concentrations above 1980 levels (in ppb), the factor <m>\eta_2</m> equals <m>\frac{1}{100}</m> (ppb<m>^{1}</m>) and <m>\eta_3</m> is the sensitivity exponent (default 1.7).
EESC concentrations are derived from the modeled concentrations of 16 ozone depleting substances controlled under the Montreal Protocol, their respective chlorine and bromine atoms, fractional release factors and a bromine versus chlorine ozone depletion efficiency (default 45) (Daniel et al., 1999).
Tropospheric aerosols
\label{section_TropAerosolParameterization}
The direct effect of aerosols is approximated by simple linear forcingabundance relationships for sulfate, nitrate, black carbon and organic carbon. Timevariable hemispheric abundances of these shortlived aerosols are in turn approximated by their hemispheric emissions, justifiable because of their very short lifetimes. The ratio of direct forcing over land and ocean areas in each hemisphere is taken from Hansen et al. (2005) (available at http://data.giss.nasa.gov/efficacy/). Specifying the direct radiative forcing patterns for one particular year, and knowing the hemispheric emissions in that year, allows us to define the future forcing as a function of future emissions.
The indirect radiative forcing, formerly modeled as dependent on SO<m>_{\rm x}</m> abundances only Wigley, 1991a, is now estimated by taking into account timeseries of sulfate, nitrate, black carbon and organic carbon optical thickness:
<m>\label{eq_forcing_indirectAer} \Delta Q_{\rm Alb,i} = r\times {\rm P}_{\rm Alb,i}\times \log(\frac{\sum_g w_g {\rm N}_{g,i}}{\sum_g w_g {\rm N}_{g,i}^0})</m>where <m>\Delta Q_{\rm Alb,i}</m> is the first indirect aerosol forcing in the four atmospheric boxes <m>i</m>, representing land and ocean areas in each hemisphere; P<m>_{\rm Alb}</m> is the fourelement pattern of aerosol indirect effects related to albedo (HTwomey, 1977) in a reference year. The second indirect effect on cloud cover changes (Albrecht, 1989) is modeled equivalently  using a reference year pattern P<m>_{\rm Cvr,i}</m>. The respective default patterns are derived from data displayed in Fig. 13 of Hansen et al. (2005). The scaling factor <m>r</m> allows one to specify a globalmean first or second indirect forcing for a specific reference year. The timevariable number concentrations of soluble aerosols N<m>_{g,i}</m> relative to their preindustrial level in each hemisphere N<m>_{g,i}^0</m> are normed to unity in that reference year. This is done separately for sulfates, nitrates, black carbon and organic carbon. For the latter, the differential solubility from industrial (fossil fuel) and biomass burning sources is taken into account (default solubility ratio 0.6/0.8) (Hansen et al., 2005). The default contribution shares w<m>_g</m> of the individual aerosol types g to the indirect aerosol effect were assigned to reflect the preliminary results by Hansen et al. (2005), namely 36 % for sulfates, 36 % for organic carbon, 23 % for nitrates and 5 % for black carbon. Note, however, that these estimates of the importance of nonSOx aerosol contributions are very uncertain, not least because the solubility, e.g. for organic carbon and nitrates have large uncertainties. The number concentrations <m>N_{g,i}</m> are here approximated by historical optical thickness estimates (as provided on http://data.giss.nasa.gov/efficacy/ see as well Supplement) and extrapolated into the future by scaling with hemispheric emissions. The general logarithmic relation between number concentrations and forcing is based on the findings by Wigley and Raper (1992); Wigley (1991a); Gultepe and Isaac(1999) and used in Hansen et al. (2005)}.