Difference between revisions of "Tagged CO simulation"

Latest revision as of 17:33, 24 October 2023

Previous | Next | Guide to GEOS-Chem simulations

1. Simulations using KPP-built mechanisms
2. Aerosol-only simulation
3. CH4 simulation
4. CO2 simulation
5. Hg simulation
6. POPs simulation
7. Tagged CO simulation
8. Tagged O3 simulation
9. TransportTracers simulation

On this page, we describe the GEOS-Chem tagged CO simulation.

Overview

Description

The tagged CO simulation is an offline simulation that calculates CO concentrations only. It uses monthly mean OH concentrations archived from a previous full-chemistry simulation (more on that below). Because the simulation is linear, CO can be “tagged” by its source region/type. The regions and types used can be adapted to address different problems with a few simple code modifications.

Assumptions

1. The tagged CO simulation doesn’t include direct emissions of volatile organic compounds (VOCs), so CO sources are scaled to account for co-emitted VOCs. Fossil fuel and biofuel emissions are scaled by 19% and biomass burning emissions are scaled by 11%. More information is given in Duncan et al. (2007). See note below!
   
   
2. Biogenic VOCs:
   
   1. Isoprene: Yield of CO from isoprene is assumed to be 30% based on Miyoshi et al. (1994). Isoprene yield can also be computed as a function of NOx concentration by setting ALPHA_ISOP_FROM_NOX = .TRUE. in CHEM_TAGGED_CO, but this is not the default behavior.
   2. Methanol: The CO flux from methanol is scaled to the isoprene flux
   3. Monoterpene: Yield of CO from monoterpenes is assumed to be 20% based on Hatakeyama et al. (1991) and Vinckier et al. (1998).
   4. Acetone: Yield of CO from acetone is assumed to be 2/3 and accounts for acetone loss from reaction with OH and photolysis.
      
      
3. OH concentrations are taken from a previously run full chemistry simulation. The default is from a much earlier version of the model, when OH was thought to be more realistic. The standard code uses OH from version 5-07-08, with GEOS3 meteorology.
   
   
4. Methane concentrations are calculated based on measurements from the NOAA Global Monitoring Division network and are assumed constant over four latitudinal bands (30-90S, 0-30S, 0-30N, 30-90N). Yield is assumed to be one molecule CO per molecule CH4.

NOTE: As described in bullet point 1 above, make sure that your HEMCO_Config.rc file contains these scale factors:

     52 COPROD_FOSSIL   1.19  - - - xy 1 1
     53 COPROD_BIOFUEL  1.189 - - - xy 1 1
     54 COPROD_BIOMASS  1.11  - - - xy 1 1

(i.e. fossil fuel/biofuel increased by ~19% and biomass by 11%). The HEMCO_Config.rc file that shipped with GEOS-Chem v11-01 contain incorrect values for these scale factors. This will be fixed for v11-02.

Standard Tracers

In a standard run, there are 17 tracers.

The regional definitions used for the fossil fuel and biomass burning tracers can be changed by modifying the HEMCO mask file:

ExtData/HEMCO/MASKS/v2014-07/tagged_CO_masks.generic.0.5x0.5.nc

The methane and VOC tracers can be removed by commenting lines in CHEM_TAGGED_CO (look for LSPLIT). Note that if you change the tracers you will also need to make the appropriate changes in your input.geos and restart files.

Notes about using HEMCO with the tagged CO simulation

In GEOS-Chem v11-01 and higher versions, the HEMCO handles all emissions for the tagged CO simulation. You should be aware of the following:

1. HEMCO by default uses a binary masking (either 0 or 1). If a grid box straddles the mask boundary, then HEMCO will count the entire box as part of the masked region. You can disable this behavior by setting Mask fractions: true in the SETTINGS section of the HEMCO_Config.rc file. But this may lead to some further discrepancies. It may be best to use the binary masking but use as fine resolution mask files e.g. 0.5 x 0.5) as possible.

2. For each file that is read from disk, we add underneath that file listing an entry to apply a regional mask to the total emissions, as shown in the previous section.

3. The Yevich & Logan biofuel emissions are added into the same CATEGORY & HIERARCHY as the EDGAR anthropogenic emissions. Because other inventories often do not separate biofuels from anthropogenic emisisons, it makes sense to lump them together. This will make it easier for the tagged tracers to sum together.

4. At present, there is no way to apply regional masks to emissions that are computed from the GFED or FINN biomass burning emissions (which are implemented as HEMCO extensions). For this reason, we the default biomass burning emissions is the QFED inventory (which is simply read from disk, and thus can be separated into tagged tracers with regional masks, as shown in the previous section). Perhaps in a future HEMCO version we will be able to apply regional masks to extension-computed emissions.

   • If you would like to use either GFED or FINN biomass emissions with the tagged CO simulation, then we recommend that you use the HEMCO standalone code to archive the total CO emissions for a given set of met fields (e.g. GEOS-FP, MERRA) and years. Then you can following the example of QFED (shown above) to apply the regional masks to the total CO biomass emissions.

5. We have removed the +LinStratChem+ block around the GMI_PROD_CO and GMI_LOSS_CO. For full-chemistry simulations, +LinStratChem+ is automatically toggled when stratospheric chemistry is turned on in input.geos. But the tagged CO always reads the GMI_PROD_CO and GMI_LOSS_CO fields directly from HEMCO and applies them without using the normal stratospheric chemistry module. ithout having to invoke the strat chem module.<p>

Studies that used Tagged CO simulation

1. Palmer, P. I., D. J. Jacob, D. B. A. Jones, C. L. Heald, R. M. Yantosca, J. A. Logan, G. W. Sachse, and D. G. Streets (2003), Inverting for emissions of carbon monoxide from Asia using aircraft observations over the western Pacific, Journal of Geophysical Research, 108(D21), 4180, doi: 10.1029/2003JD003397.
2. Heald, C. L., D. J. Jacob, D. B. A. Jones, P. I. Palmer, J. A. Logan, D. G. Streets, G. W. Sachse, J. C. Gille, R. N. Hoffman, and T. Nehrkorn (2004), Comparative inverse analysis of satellite (MOPITT) and aircraft (TRACE-P) observations to estimate Asian sources of carbon monoxide, Journal of Geophysical Research, 109(D15S04), doi: 10.1029/2004JD005185.
3. Arellano, A. F., P. S. Kasibhatla, L. Giglio, G. R. van der Werf, and J. T. Randerson (2004), Top-down estimates of global CO sources using MOPITT measurements, Geophysical Research Letters, 31(L01104), doi: 10.1029/2003GL018609.
4. Arellano, A. F., P. S. Kasibhatla, L. Giglio, G. R. van der Werf, J. T. Randerson, and G. J. Collatz (2006), Time-dependent inversion estimates of global biomass-burning CO emissions using Measurement of Pollution in the Troposphere (MOPITT) measurements, J. Geophys. Res., 111(D09303), doi: 10.1029/2005JD006613.
5. Duncan, B. N., Logan, J. A., Bey, I., Megretskaia, I. A., Yantosca, R. M., Novelli, P. C., Jones, N. B., and Rinsland, C. P., Global budget of CO, 1988–1997: Source estimates and validation with a global model, J. Geophys. Res., 112, D22301, doi:10.1029/2007JD008459, 2007.
6. Duncan, B. N., J. A. Logan, I. Bey, I. A. Megretskaia, R. M. Yantosca, P. C. Novelli, N. B. Jones, and C. P. Rinsland (2008), Model analysis of the factors regulating the trends and variability of carbon monoxide between 1988 and 1997, Atmos. Chem. Phys, 8, 7389-3403.
7. Kopacz, M., D. J. Jacob, D. K. Henze, C. L. Heald, D. G. Streets, and Q. Zhang (2009), Comparison of adjoint and analytical Bayesian inversion methods for constraining Asian sources of carbon monoxide using satellite (MOPITT) measurements of CO columns, J. Geophys. Res., 114(D04305), doi: 10.1029/2007JD009264.
8. Fisher, J.A., D.J. Jacob, M.T. Purdy, M. Kopacz, P. Le Sager, C. Carouge, C.D. Holmes, R.M. Yantosca, R.L. Batchelor, K. Strong, G.S. Diskin, H.E. Fuelberg, J.S. Holloway, E.J. Hyer, W.W. McMillan, J. Warner, D.G. Streets, Q. Zhang, Y. Wang, S. Wu, Source attribution and interannual variability of Arctic pollution in spring constrained by aircraft (ARCTAS, ARCPAC) and satellite (AIRS) observations of carbon monoxide, Atm. Chem. Phys. Discuss., 9, 19035-19080, 2009.
9. Kopacz, M., D.J. Jacob, J.A. Fisher, J.A. Logan, L. Zhang, I.A. Megretskaia, R.M. Yantosca, K. Singh, D.K. Henze, J.P. Burrows, M. Buchwitz, I. Khlystova, W.W. McMillan, J.C. Gille, D.P. Edwards, A. Eldering, V. Thouret, P. Nedelec, Global estimates of CO sources with high resolution by adjoint inversion of multiple satellite datasets (MOPITT, AIRS, SCIAMACHY and TES), Atm. Chem. Phys. Discuss., 9, 19967-20018, 2009.

References

1. Duncan, B. N., Logan, J. A., Bey, I., Megretskaia, I. A., Yantosca, R. M., Novelli, P. C., Jones, N. B., and Rinsland, C. P., Global budget of CO, 1988–1997: Source estimates and validation with a global model, J. Geophys. Res., 112, D22301, doi:10.1029/2007JD008459, 2007.
2. Hatakeyama, S., Izumi, K., Fukuyama, T., Akimoto, H. Washida, N., Reactions of OH with alpha-pinene and beta-pinene in air: Estimate of global CO production from the atmospheric oxidation of terpenes, J. Geophys. Res., 96(D1), 947-958, 1991.
3. Miyoshi, A., Hatakeyama, S., Washida, N., OH radical-initiated photooxidation of isoprene: An estimate of global CO production, J. Geophs. Res., 99(D9), 18779-18787, 1994.
4. Vinckier, C., Compernolle, F., Saleh, A. M., Van Hoof, N., Van Hees, I., Product yields of the alpha -pinene reaction with hydroxyl radicals and the implication on the global emission of trace compounds in the atmosphere, Fresenius Env. Bull., 7(5-6), 361-368, 1998.