From Geos-chem
Revision as of 14:58, 25 September 2017 by Melissa Payer (Talk | contribs) (Bug fixes for EOH and MGLY)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

This page includes Emily Fischer's descriptions about changes to the PAN simulation.


Emily Fischer has updated the PAN simulation. This work was summarized in (Fischer et al., 2014) in ACP. The original code was based on GEOS-Chem v9-01-01 and these updates were introduced into the standard code in GEOS-Chem v11-02a (approved 12 May 2017).

Updates to Chemistry

Isoprene chemical mechanism

I replaced the isoprene chemical mechanism with the Paulot scheme. This scheme has already been implemented into the standard code.

Nighttime chemistry

I added nighttime chemistry from reactions of organic peroxy radicals with NO3 following Stone et al. (2013). This may not be incorporated into the standard chemistry, and needs to be discussed by the GEOS-Chem Steering Committee.

Nighttime reactions of NO3 increase PAN formation. The difference in monthly mean simulated PAN is largest in regions with the largest PAN production (Eastern U.S. and East Asia). The addition of nighttime chemistry increases monthly mean surface PAN mixing ratios in these regions by ~50 pptv in spring and summer. This does propagate to higher altitudes, but the changes are very modest. Springtime PAN mixing ratios increase broadly over the Arctic by 5-10 pptv (~5%). The impact in summer is more heterogeneous at higher altitudes, but the overall increase is on the order of ~2% over most of the Arctic. These changes also DO impact ozone. This chemistry increases ozone over the Arctic during spring. Surface ozone increases by ~0.5 ppbv, and the impact is larger aloft (0.75 ppbv at 4 km). I implemented the RO2 + NO3 reactions prior to the adoption of the Paulot isoprene scheme in my version of the code. We need to decide how this fits in with the nighttime isoprene chemistry that is currently under development.

I updated the rate coefficients for the reactions of HO2 with the >C2 peroxy radicals to Equation (iv) in Saunders et al. (2003). These changes have been incorporated into the standard code.

Introduction of new NMVOCs

I added several new NMVOCs. The extended mechanism includes ethanol, benzene, toluene and ethylbenzene (lumped), xylenes and trimethyl benzenes(lumped), and monoterpenes (lumped). Hydroxyacetone and methylglyoxal, which are species in past versions, are treated as tracers in this simulation. Hydroxyacetone has a 1-2 day lifetime. Methylglyoxal is treated as a tracer so that it can be emitted from biomass burning plumes and so I could track PAN production via this pathway.

The inclusion and treatment of aromatics was motivated by Liu et al. (2010). I calculated the associated yield of methylglyoxal using recommended values for the individual aromatic species(toluene, o-xylene, m-xylene, p-xylene, 1,2,3-trimethylbenzene, 1,2,4-trimethylbenzene, and 1,3,5-trimethylbenzene) from Nishino et al. (2010) and the observed mean aromatic speciation for Chinese cities from Barletta et al. (2006). Thus the treatment is particular to the limited observations of aromatic speciation in China. It would be good to determine how different these ratios would be for the US, Europe or other developed/developing regions.

We adopted the treatment of monoterpene oxidation from theRACM2 chemical mechanism (Goliff et al., 2013), lumping terpenes with one double bond (alpha-pinene, beta-pinene, sabinene and delta-3-carene) into one proxy. Unlike Ito et al. (2007), hydroxyacetone is not a product of terpene oxidation in the revised RACM2 mechanism used here. So I expect that the mechanism used here will make less PAN than that in Ito et al., 2007. This code was provided by Jingqiu Mao.

PAN dry deposition

I increased the deposition flux of PAN by assuming PAN reactivity with surfaces is more similar to O3 than NO2. (I changed F0 from 0.1 to 1.0 in drydep_mod.F. F0 is a 'reactivity factor for oxidation of biological substances')

  • 0.1 is used for slightly reactive species (NO2, HNO2)
  • 1 is used for highly reactive species (O3, H2O2, and recently OVOCs)

This change increases the global annual PAN deposition from 6.5 Tg/year to 10.7 Tg/year. This is a change of 0.47 Tg N/year. During the springtime PAN maximum, faster deposition decreased PAN at the surface by 15 – 20% (30 – 50 pptv) throughout the northern hemisphere. At higher altitudes, there is also a hemispheric impact, with PAN decreases of ~10% (20 – 40 pptv). In summer the surface decrease in PAN is most notable over the Arctic, but here it is on the order og 25%. At 4 km, the PAN decrease is about 10%, on the order of 5 – 10 pptv. These changes do propagate to ozone. The faster deposition causes a springtime decrease of ozone of almost a ppb at the surface and in the lower free troposphere. In summer, the changes to ozone are most pronounced for the Arctic, but still 0.5 ppbv over the eastern U.S. The limited evidence from the literature suggests faster PAN deposition than previously assumed. In our model this means f0=1 is probably better than f0=0.1, though a different deposition model would probably be ideal in the longer term.

Emission updates

Ethane and propane emissions

The PAN simulation uses the RETRO emission inventory as global default for anthropogenic NMVOC emissions aside from ethane (C2H6) and propane (C3H8). Ethane and propane emissions in RETRO were far too low compared to the GEOS-Chem inventories from Xiao et al. (2008) and available observations. Emissions of both species appeared to be missing from the major natural gas production region in Russia. I didn't catch this when I did the acetone updates because these emissions don't produce biased results downwind of major anthropogenic source regions (China for example) compared to the existing data. I used the offline ethane emission inventories developed as in Xiao et al. (2008). I have not been able to find the full reference for the propane inventory that corresponds to the Xiao et al. (2008) ethane inventory. I have tried unsuccessfully to make contact with Yaping Xiao.

Aromatic emissions

RETRO includes anthropogenic emissions for benzene (BENZ), xylene (XYLE) and toluene (TOLU). Based on the observed CO to benzene ratio for TRACE-P, I increased benzene emissions over China by 25%. The CO/benzene ratios were examined by Chris Miller. I then scaled xylene and toluene emissions to benzene based on measurements from 43 Chinese cities from Barletta et al. (2006). In summary, RETRO emissions of toluene were increased by a factor of 4 over China to create the lumped toluene (TOLU), and RETRO emissions xylene were increased by a factor of 8 over China to create the lumped xylene species (XYLE). I suspect that upward scaling may also be appropriate in other developing cities, but this would need to be carefully examined and I am not sure what data is available.

Ethanol emissions

Emissions of ethanol follow Dylan Millet's code.

Acetaldehyde ocean source

I also added emissions of acetaldehyde from the ocean. This code was also provided by Dylan Millet.

Biomass burning emission updates

I made several updates to the way fire emissions are incorporated into the model.

GFED emission factors

NOTE: In GEOS-Chem v11-02, these emission factor updates are only applied to GFED3, and not GFED4.

I updated the emission factors for NMVOCs and NOx from extratropical forests, savannas and agricultural fires from Akagi et al. (2011). These are the only categories that overlap with GFED3, and this is why these are the only ones to be updated. The updated NOx emission factor for extratropical fires is approximately a factor of three lower, and the emission factors for the NMVOCs are generally higher. Sometimes they are much higher.

Partitioning of NO biomass burning emissions

Following Alvarado et al. (2010) I directly partition 40% and 20% NO emissions from fires directly to PAN and HNO3 respectively. The Alvarado et al. (2010) partitioning is based on an observation of fresh boreal fire plumes over North America in summer, but I have applied it to all fire types. To support PA radical formation on faster timescales I also added GFED3 emissions of several shorter lived hydrocarbons to the suite of species emitted from fires: terpenes (APINE), aromatics (XYLE, TOLU) along with additional oxygenated species (hydroxyacetone (HAC), and methylgloyxal (MGLY)).

NOTE: After the 1-year benchmark for v11-02a-Run0, the GCSC recommended we should remove the partitioning of NO biomass burning emissions to PAN and HNO3.

Prasad Kasibhatla wrote:

Are we really comfortable with the partitioning of fire NOx emissions instantaneously to 40% PAN and 20% HNO3 everywhere just based one study (Alvarado et al, 2010) in boreal regions? Also, while Alvarado et al. find 40% and 20% NOx converted to PAN and HNO3 relatively rapidly in fire plumes, to get this they use pseudoemission factors that apportion 65% of the emission to PAN (Table 5 in the their study).

Emily Fischer wrote:

Prasad is correct that this ratio is only based on Boreal fires. Removing this is very easy, however, we may want to leave in commented out code so this is an easy option for anyone trying to simulate north American fires.

Injection height of biomass burning emissions

v9.01.01 releases all fire emissions in the boundary layer, but injecting a fraction of the emissions above the boundary layer is important for PAN because temperatures above the boundary layer enhance its stability. I distribute 35% of biomass burning emissions by mass in the 10 sigma layers (4 km) above the boundary layer, and this improves the comparison with PAN observations at high latitudes.

NOTE: After the 1-month benchmark for v11-02a, the GCSC decided we should remove the vertically distributed biomass burning emissions.

Emily Fischer wrote:

I am concerned about implementing the biomass burning injection heights globally in the standard model. I don't think that is the right approach. It is a good approach for boreal burning - but not a global solution. We (Juliet and I) have a new implementation of Maria Val Martin's global injection heights based on MISR data. She will show a poster at IGC8. My recommendation is to keep biomass burning injection heights as is (i.e. into the boundary layer), and then we could implement a global injection height parameterization later.

Scaling regional biomass burning emissions

NOTE: In GEOS-Chem v11-02, these scale factor are not applied because GFED4 is now the default and these factors were computed for 2008.

Finally, I increased wild fire emissions by 60% in North Asia (30 – 75°N, 60 – 190°E), 25% in Canada and 50% in Alaska. The increase over Asia was motivated by Kaiser et al. (2012). This region is one region of comparison in this paper. For Canada and Alaska, I deferred to the analysis of fire records by Xu Yue.

Previous issues that are now resolved

Bug fixes for EOH and MGLY

This update was included in v11-02c and approved on 21 Sep 2017.

Advected species EOH and MGLY were introduced in v11-02a as part of the PAN chemistry updates described above. It was discovered during v11-02c development that following the v11-02a updates MGLY and EOH were no longer being properly identified as KPP species, and therefore their concentrations were not getting properly updated by FlexChem. This was fixed by adding the following lines in Headers/species_database_mod.F90:

         CASE( 'EOH' )
            CALL Spc_Create( am_I_Root     = am_I_Root,                    &
                             ThisSpc       = SpcData(N)%Info,              &
                             ModelID       = N,                            &
                             KppSpcId      = KppSpcId(N),                  &
                             KppVarId      = KppVarId(N),                  &
                             KppFixId      = KppFixId(N),                  &
                             Name          = NameAllCaps,                  &
                             FullName      = 'Ethanol',                    &
                             Formula       = 'C2H5OH',                     &
                             MW_g          = 46.07_fp,                     &
                             EmMW_g        = 12.0_fp,                      &
                             MolecRatio    = 2.0_fp,                       &
                             Is_Advected   = Is_Advected,                  &
                             Is_Gas        = T,                            &
                             Is_Drydep     = T,                            &
                             Is_Wetdep     = T,                            &
                             DD_F0         = 0.0_fp,                       &
                             DD_Hstar_old  = 1.9e+2_fp,                    &
                             Henry_K0      = 1.9e+2_f8,                    &
                             Henry_CR      = 6600.0_f8,                    &
                             WD_RetFactor  = 2.0e-2_fp,                    &
                             RC            = RC )

         CASE( 'MGLY' )
            CALL Spc_Create( am_I_Root     = am_I_Root,                    &
                             ThisSpc       = SpcData(N)%Info,              &
                             ModelID       = N,                            &
                             KppSpcId      = KppSpcId(N),                  &
                             KppVarId      = KppVarId(N),                  &
                             KppFixId      = KppFixId(N),                  &
                             Name          = NameAllCaps,                  &
                             FullName      = 'Methylglyoxal',              &
                             Formula       = 'CH3COCHO',                   &
                             MW_g          = 72.0_fp,                      &
                             Is_Advected   = Is_Advected,                  &
                             Is_Gas        = T,                            &
                             Is_Drydep     = T,                            &
                             Is_Wetdep     = T,                            &
                             Is_Photolysis = T,                            &
                             DD_F0         = 1.0_fp,                       &
                             DD_Hstar_old  = 3.7e+3_fp,                    &
                             Henry_K0      = 3.7e+3_f8,                    &
                             Henry_CR      = 7500.0_f8,                    &
                             WD_RetFactor  = 2.0e-2_fp,                    &
                             RC            = RC )

In addition, there was still a line in flexchem_mod.F90 for setting the concentration in KPP to zero.

      C(ind_EOH)   = 0.0_dp

That line is no longer needed because EOH is now an active chemical species, so it has been removed.

--Melissa Sulprizio (talk) 14:30, 17 August 2017 (UTC)


  • Akagi, S. K., Yokelson, R. J., Wiedinmyer, C., Alvarado, M. J., Reid, J. S., Karl, T., Crounse, J. D., and Wennberg, P. O., Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039-4072, 10.5194/acp-11-4039-2011, 2011.
  • Alvarado, M. J., Logan, J. A., Mao, J., Apel, E., Riemer, D., Blake, D., Cohen, R. C., Min, K. E., Perring, A. E., Browne, E. C., Wooldridge, P. J., Diskin, G. S., Sachse, G. W., Fuelberg, H., Sessions, W. R., Harrigan, D. L., Huey, G., Liao, J., Case-Hanks, A., Jimenez, J. L., Cubison, M. J., Vay, S. A., Weinheimer, A. J., Knapp, D. J., Montzka, D. D., Flocke, F. M., Pollack, I. B., Wennberg, P. O., Kurten, A., Crounse, J., Clair, J. M. S., Wisthaler, A., Mikoviny, T., Yantosca, R. M., Carouge, C. C., and Le Sager, P., Nitrogen oxides and PAN in plumes from boreal fires during ARCTAS-B and their impact on ozone: an integrated analysis of aircraft and satellite observations, Atmos. Chem. Phys., 10, 9739-9760, 10.5194/acp-10-9739-2010, 2010.
  • Barletta, B., Meinardi, S., Simpson, I. J., Sherwood Rowland, F., Chan, C.-Y., Wang, X., Zou, S., Chan, L. Y., and Blake, D. R., Ambient halocarbon mixing ratios in 45 Chinese cities, Atmos. Env., 40, 7706-7719, 10.1016/j.atmosenv.2006.08.039, 2006.
  • Fischer, E.V., D.J. Jacob, R.M. Yantosca, M.P. Sulprizio, D.B. Millet, J. Mao, F. Paulot, H.B. Singh, A.-E. Roiger, L. Ries, R.W. Talbot, K. Dzepina, and S. Pandey Deolal, Atmospheric peroxyacetylnitrate (PAN): a global budget and source attribution, Atmos. Chem. Phys., 14, 2679-2698, 2014. (pdf)
  • Goliff, W. S., Stockwell, W. R., and Lawson, C. V., The regional atmospheric chemistry mechanism, version 2, Atmos. Env., 68, 174-185, 2013.
  • Ito, A., Sillman, S., and Penner, J. E., Effects of additional nonmethane volatile organic compounds, organic nitrates, and direct emissions of oxygenated organic species on global tropospheric chemistry, J. Geophys. Res., 112, D06309, 10.1029/2005jd006556, 2007.
  • Kaiser, J. W., Heil, A., Andreae, M. O., Benedetti, A., Chubarova, N., Jones, L., Morcrette, J. J., Razinger, M., Schultz, M. G., Suttie, M., and van der Werf, G. R., Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power, Biogeosci., 9, 527-554, 10.5194/bg-9-527-2012, 2012.
  • Liu, Z., Wang, Y., Gu, D., Zhao, C., Huey, L. G., Stickel, R., Liao, J., Shao, M., Zhu, T., Zeng, L., Liu, S.-C., Chang, C.-C., Amoroso, A., and Costabile, F., Evidence of Reactive Aromatics As a Major Source of Peroxy Acetyl Nitrate over China, Environ. Sci. Tech., 44, 7017-7022, 10.1021/es1007966, 2010.
  • Millet, D. B., Guenther, A., Siegel, D. A., Nelson, N. B., Singh, H. B., de Gouw, J. A., Warneke, C., Williams, J., Eerdekens, G., Sinha, V., Karl, T., Flocke, F., Apel, E., Riemer, D. D., Palmer, P. I., and Barkley, M., Global atmospheric budget of acetaldehyde: 3-D model analysis and constraints from in-situ and satellite observations, Atmos. Chem. Phys., 10, 3405-3425, doi:10.5194/acp-10-3405-2010, 2010.
  • Millet, D.B., E. Apel, D.K. Henze, J. Hill, J.D. Marshall, H.B. Singh, and C.W. Tessum, Natural and anthropogenic ethanol sources in North America and potential atmospheric impacts of ethanol fuel use, Environ. Sci. Technol., 46, 8484-8492, 2012.
  • Nishino, N., Arey, J., and Atkinson, R., Formation Yields of Glyoxal and Methylglyoxal from the Gas-Phase OH Radical-Initiated Reactions of Toluene, Xylenes, and Trimethylbenzenes as a Function of NO2 Concentration, J. Phys. Chem. A, 114, 10140-10147, 10.1021/jp105112h, 2010.
  • Saunders, S. M., Jenkin, M. E., Derwent, R. G., and Pilling, M. J., Protocol for the development of the Master Chemical Mechanism, MCM v3 (Part A): tropospheric degradation of non-aromatic volatile organic compounds, Atmos. Chem. Phys., 3, 161-180, 10.5194/acp-3-161-2003, 2003.
  • Stone, D., Evans, M. J., Walker, H., Ingham, T., Vaughan, S., Ouyang, B., Kennedy, O. J., McLeod, M. W., Jones, R. L., Hopkins, J., Punjabi, S., Lidster, R., Hamilton, J. F., Lee, J. D., Lewis, A. C., Carpenter, L. J., Forster, G., Oram, D. E., Reeves, C. E., Bauguitte, S., Morgan, W., Coe, H., Aruffo, E., Dari-Salisburgo, C., Giammaria, F., Di Carlo, P., and Heard, D. E., Radical chemistry at night: comparisons between observed and modelled HOx, NO3 and N2O5 during the RONOCO project, Atmos. Chem. Phys., 14, 1299-1321, doi:10.5194/acp-14-1299-2014, 2014.
  • Xiao, Y., Logan, J. A., Jacob, D. J., Hudman, R. C., Yantosca, R., and Blake, D. R., Global budget of ethane and regional constraints on U.S. sources, J. Geophys. Res., 113, D21306, 10.1029/2007jd009415, 2008.