Difference between revisions of "PAN"

From Geos-chem
Jump to: navigation, search
Line 1: Line 1:
Emily Fischer has updated the PAN simulation.  The code was based on v9.01.01, but has been merged with v9.02.h at this point, and a publication [http://acmg.seas.harvard.edu/publications/fischer2013.pdf (Fischer et al., 2013)] has been submitted to ACP.
+
[http://www.atmos.colostate.edu/faculty/fischer.php Emily Fischer] has updated the PAN simulation.  The code was based on [[GEOS-Chem v9-01-01]], but has been merged with [[GEOS-Chem_v9-02#What.27s_new_in_this_version|v9-02h]] at this point, and a publication [http://acmg.seas.harvard.edu/publications/fischer2013.pdf (Fischer et al., 2013)] has been published in ACP.
== '''Updates to Chemistry''' ==
+
  
1.  I replaced the isoprene chemical mechanism with the [http://wiki.seas.harvard.edu/geos-chem/index.php/New_isoprene_scheme Paulot scheme]. This scheme has already been implemented into the standard code.  Here are some bullets on the effect of the addition of nighttime chemistry:
+
== Overview ==
  
2. I added nighttime chemistry from reactions of organic peroxy radicals with NO3 following [http://www.atmos-chem-phys-discuss.net/13/9519/2013/acpd-13-9519-2013.html Stone et al. (2013)].  This may not be incorporated into the standard chemistry, and needs to be discussed by the GEOS-Chem Steering Committee.
+
=== Updates to Chemistry ===
  
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.
+
'''''[http://www.atmos.colostate.edu/faculty/fischer.php Emily Fischer] wrote:'''''
  
3.  I updated the rate coefficients for the reactions of [http://wiki.seas.harvard.edu/geos-chem/index.php/New_isoprene_scheme#Update_One_-_RO2.2BHO2_Reaction_Rate HO2 with the >C2 peroxy radicals] to Equation (iv) in Saunders et al. (2003).  These changes have been incorporated into the standard code.
+
<ol>
 +
<li><p>I replaced the isoprene chemical mechanism with the [http://wiki.seas.harvard.edu/geos-chem/index.php/New_isoprene_scheme Paulot scheme]. This scheme has already been implemented into the standard code. Here are some bullets on the effect of the addition of nighttime chemistry:</p></li>
  
4. 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. 
+
<li><p>I added nighttime chemistry from reactions of organic peroxy radicals with NO3 following [http://www.atmos-chem-phys.net/14/1299/2014/acp-14-1299-2014.html Stone et al. (2013)].  This may not be incorporated into the standard chemistry, and needs to be discussed by the GEOS-Chem Steering Committee.</p>
  
- 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 [http://pubs.acs.org/doi/abs/10.1021/jp105112h Nishino et al. (2010)] and the observed mean aromatic speciation for Chinese cities from [http://ps.uci.edu/~rowlandblake/publications/barlettaambient.pdf 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.
+
<p>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.</p></li>
  
- We adopted the treatment of monoterpene oxidation from theRACM2 chemical mechanism [http://www.sciencedirect.com/science/article/pii/S1352231012011065 (Goliff et al., 2013)], lumping terpenes with one double bond (alpha-pinene, beta-pinene, sabinene and delta-3-carene) into one proxy. Unlike [http://onlinelibrary.wiley.com/doi/10.1029/2005JD006556/abstract Ito et al. (2007)], hydroxyacetone is not a product of terpene oxidation in the revised RACM2 mechanism used hereSo 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.
+
<li><p>I updated the rate coefficients for the reactions of [http://wiki.seas.harvard.edu/geos-chem/index.php/New_isoprene_scheme#Update_One_-_RO2.2BHO2_Reaction_Rate HO2 with the >C2 peroxy radicals] to Equation (iv) in Saunders et al. (2003).  These changes have been incorporated into the standard code.</p></li>
  
5.  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’)
+
<li><p>I added several new NMVOCsThe 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.</p></li> 
 +
<ul>
 +
<li><p>The inclusion and treatment of aromatics was motivated by [http://pubs.acs.org/doi/abs/10.1021/es1007966 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 [http://pubs.acs.org/doi/abs/10.1021/jp105112h Nishino et al. (2010)] and the observed mean aromatic speciation for Chinese cities from [http://ps.uci.edu/~rowlandblake/publications/barlettaambient.pdf 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.</p></li>
 +
<li><p>We adopted the treatment of monoterpene oxidation from theRACM2 chemical mechanism [http://www.sciencedirect.com/science/article/pii/S1352231012011065 (Goliff et al., 2013)], lumping terpenes with one double bond (alpha-pinene, beta-pinene, sabinene and delta-3-carene) into one proxy. Unlike [http://onlinelibrary.wiley.com/doi/10.1029/2005JD006556/abstract 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.</p></li>
 +
</ul>
  
- 0.1 is used for slightly reactive species (NO2, HNO2)
+
<li><p>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’)</p>
 +
<ul>
 +
<li>0.1 is used for slightly reactive species (NO2, HNO2)</li>
 +
<li>1 is used for highly reactive species (O3, H2O2, and recently OVOCs)</li>
 +
</ul>
 +
<p>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.</p></li>
 +
</ol>
  
- 1 is used for highly reactive species (O3, H2O2, and recently OVOCs)
+
=== Updates to Emissions ===
  
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.
+
'''''[http://www.atmos.colostate.edu/faculty/fischer.php Emily Fischer] wrote:'''''
== '''Updates to Emissions''' ==
+
  
 +
<ol>
 +
<li><p>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 [http://acmg.seas.harvard.edu/publications/2007/xiao.2007b.pdf 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.</p></li>
  
1) 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.
+
<li><p>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.</p></li>
  
2) 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
+
<li><p>Emissions of ethanol follow [http://www.atmoschem.umn.edu/papers/millet_2012.pdf Dylan Millet's code].</p></li>
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.
+
  
3)  Emissions of ethanol follow [http://www.atmoschem.umn.edu/papers/millet_2012.pdf Dylan Millet's code].  
+
<li><p>I also added emissions of acetaldehyde from the ocean. This code was also provided by [http://www.atmos-chem-phys.net/10/3405/2010/acp-10-3405-2010.html Dylan Millet].</p></li>
 +
</ol>
  
4) I also added emissions of acetaldehyde from the ocean. This code was also provided by [http://www.atmos-chem-phys.net/10/3405/2010/acp-10-3405-2010.html Dylan Millet].
 
  
I made several '''updates to the way fire emissions are incorporated''' into the model.   
+
:I made several '''updates to the way fire emissions are incorporated''' into the model.   
  
5) 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
+
<ol>
generally higher. Sometimes they are much higher.
+
<li><p>I updated the emission factors for NMVOCs and NOx from extratropical forests, savannas and agricultural fires from [http://www.atmos-chem-phys.net/11/4039/2011/acp-11-4039-2011.html 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.</p></li>
  
6) 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)
+
<li><p>Following [http://www.atmos-chem-phys.net/10/9739/2010/acp-10-9739-2010.html 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)).</p></li>
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)).
+
  
7) 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.
+
<li><p>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.</p></li>
  
8)  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.
+
<li><p>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 [http://www.biogeosciences.net/9/527/2012/bg-9-527-2012.html 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.</p></li>
 +
</ol>
 +
 
 +
== References ==
 +
 
 +
*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'', <u>Atmos. Chem. Phys.</u>, '''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'', <u>Atmos. Chem. Phys.</u>, '''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'', <u>Atmos. Env.</u>, '''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'', <u>Atmos. Chem. Phys.</u>, '''14''', 2679-2698, 2014. [http://acmg.seas.harvard.edu/publications/fischer2013.pdf (pdf)]
 +
*Goliff, W. S., Stockwell, W. R., and Lawson, C. V., ''The regional atmospheric chemistry mechanism, version 2'', <u>Atmos. Env.</u>, '''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'', <u>J. Geophys. Res.</u>, '''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'', <u>Biogeosci.</u>, '''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'', <u>Environ. Sci. Tech.</u>, '''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'', <u>Atmos. Chem. Phys.</u>, '''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'', <u>Environ. Sci. Technol.</u>, '''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'', <u>J. Phys. Chem. A</u>, '''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'', <u>Atmos. Chem. Phys.</u>, '''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'', <u>Atmos. Chem. Phys.</u>, '''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'', <u>J. Geophys. Res.</u>, '''113''', D21306, 10.1029/2007jd009415, 2008.

Revision as of 18:01, 20 May 2015

Emily Fischer has updated the PAN simulation. The code was based on GEOS-Chem v9-01-01, but has been merged with v9-02h at this point, and a publication (Fischer et al., 2013) has been published in ACP.

Overview

Updates to Chemistry

Emily Fischer wrote:

  1. I replaced the isoprene chemical mechanism with the Paulot scheme. This scheme has already been implemented into the standard code. Here are some bullets on the effect of the addition of nighttime chemistry:

  2. 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.

  3. 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.

  4. 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.

  5. 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.

Updates to Emissions

Emily Fischer wrote:

  1. 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.

  2. 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.

  3. Emissions of ethanol follow Dylan Millet's code.

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


I made several updates to the way fire emissions are incorporated into the model.
  1. 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.

  2. 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)).

  3. 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.

  4. 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.

References

  • 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.