Difference between revisions of "Aerosol-only simulation"

From Geos-chem
Jump to: navigation, search
(GEOS-Chem v9-01-03 and earlier versions)
Line 27: Line 27:
  
 
=== Tracer list ===
 
=== Tracer list ===
 
==== GEOS-Chem v9-01-03 and earlier versions ====
 
 
In [[GEOS-Chem v9-01-03]] and earlier versions, there are 31 tracers in a complete aerosol-only simulation (see input.geos below). 
 
 
The offline aerosol is simulation type 10. For a full set of aerosol tracers (based on the 59-tracer full-chemistry SOA simulation), the input.geos should look like this:
 
 
Type of simulation      : 10
 
Number of Tracers      : 31
 
Tracer Entries -------> : TR#  Name  g/mole Tracer Members; () = emitted
 
Tracer #1              : 1    DMS    62.00
 
Tracer #2              : 2    SO2    64.00
 
Tracer #3              : 3    SO4    96.00
 
Tracer #4              : 4    MSA    96.00
 
Tracer #5              : 5    NH3    17.00
 
Tracer #6              : 6    NH4    18.00
 
Tracer #7              : 7    NIT    62.00
 
Tracer #8              : 8    H2O2  34.00
 
Tracer #9              : 9    BCPI  12.00
 
Tracer #10              : 10  OCPI  12.00
 
Tracer #11              : 11  BCPO  12.00
 
Tracer #12              : 12  OCPO  12.00
 
Tracer #13              : 13  ALPH  136.23
 
Tracer #14              : 14  LIMO  136.23
 
Tracer #15              : 15  ALCO  142.00
 
Tracer #16              : 16  SOG1  150.00
 
Tracer #17              : 17  SOG2  160.00
 
Tracer #18              : 18  SOG3  220.00
 
Tracer #19              : 19  SOG4  130.00
 
Tracer #20              : 20  SOG5  150.00
 
Tracer #21              : 21  SOA1  150.00
 
Tracer #22              : 22  SOA2  160.00
 
Tracer #23              : 23  SOA3  220.00
 
Tracer #24              : 24  SOA4  130.00
 
Tracer #25              : 25  SOA5  150.00
 
Tracer #26              : 26  DST1  29.00
 
Tracer #27              : 27  DST2  29.00
 
Tracer #28              : 28  DST3  29.00
 
Tracer #29              : 29  DST4  29.00
 
Tracer #30              : 30  SALA  36.00
 
Tracer #31              : 31  SALC  36.00
 
 
This list can be reduced by turning off any aerosol species (eg. SOA, dust, SS, etc).  You may also run without the secondary organic aerosols (ALPH, LIMO, ALCO, SOG{1,2,3,4,5}, SOA{1,2,3,4,5}) if you wish (but we recommend that you use them).
 
  
 
==== GEOS-Chem v9-02 and later versions ====
 
==== GEOS-Chem v9-02 and later versions ====

Revision as of 15:45, 25 August 2022

Previous | Next | Guide to GEOS-Chem simulations

  1. Simulations overview
  2. Mechanisms for full-chemistry simulations (e.g. Standard, Tropchem, etc.)
  3. Aerosol-only simulation
  4. CH4 simulation
  5. CO2 simulation
  6. Hg simulation
  7. POPs simulation
  8. Tagged CO simulation
  9. Tagged O3 simulation
  10. TransportTracers simulation


On this page we describe the aerosol-only (a.k.a. offline aerosol) simulation.

Description

The aerosol simulation is an offline simulation for aerosol tracers only. It uses archived monthly mean OH, NO3, O3 and total nitrate concentrations archived from a previous full-chemistry simulation (more on that below), as well as production and loss rates for H2O2. This simulation does not provide "tagged" capabilities, but reduces the suite of tracers from full chemistry.

Assumptions

  1. OH concentrations are taken from a previously run full chemistry simulation. The default is from a much earlier version of the model, when OH thought to be more realistic. The standard code uses OH from version 5-07-08, with GEOS3 meteorology.
  2. O3, NO3 and total nitrate (HNO3+NIT) concentrations, as well as production rates and photolysis rates for H2O2 are taken from a previously run full chemistry simulation. The current default fields were produced from a GEOS-4, v7-02-03 simulation. These files are generally located in the sulfate_sim_200508/offline/ directory in your GEOS-Chem data directories.

Practicalities

As with all other Offline chemistry simulations, the GEOS-Chem user community shall be responsible for the scientific content and validation of the aerosol-only simulation.

Tracer list

GEOS-Chem v9-02 and later versions

The SOA simulation was updated in GEOS-Chem v9-02o to include a better tracer lumping scheme, as shown in Figure 1 of Pye et al. 2010.

The offline aerosol is simulation type 10. For a full set of aerosol tracers (based on the 100-tracer full-chemistry SOA simulation), the input.geos should look like this:

Type of simulation      : 10
Number of Tracers       : 47
Tracer Entries -------> : TR#  Name   g/mole   Tracer Members; () = emitted
Tracer #1               :  1   DMS     62.00
Tracer #2               :  2   SO2     64.00
Tracer #3               :  3   SO4     96.00
Tracer #4               :  4   SO4s    96.00
Tracer #5               :  5   MSA     96.00
Tracer #6               :  6   NH3     17.00
Tracer #7               :  7   NH4     18.00
Tracer #8               :  8   NIT     62.00
Tracer #9               :  9   NITs    62.00
Tracer #10              : 10   H2O2    34.00
Tracer #11              : 11   BCPI    12.00
Tracer #12              : 12   OCPI    12.00
Tracer #13              : 13   BCPO    12.00
Tracer #14              : 14   OCPO    12.00
Tracer #15              : 15   DST1    29.00
Tracer #16              : 16   DST2    29.00
Tracer #17              : 17   DST3    29.00
Tracer #18              : 18   DST4    29.00
Tracer #19              : 19   SALA    36.00
Tracer #20              : 20   SALC    36.00
Tracer #21              : 21   MTPA   136.23
Tracer #22              : 22   LIMO   136.23
Tracer #23              : 23   MTPO   136.23
Tracer #24              : 24   TSOG1  150.00
Tracer #25              : 25   TSOG2  150.00
Tracer #26              : 26   TSOG3  150.00
Tracer #27              : 27   TSOG0  150.00
Tracer #28              : 28   TSOA1  150.00
Tracer #29              : 29   TSOA2  150.00
Tracer #30              : 30   TSOA3  150.00
Tracer #31              : 31   TSOA0  150.00
Tracer #32              : 32   ISOG1  150.00
Tracer #33              : 33   ISOG2  150.00
Tracer #34              : 34   ISOG3  150.00
Tracer #34              : 35   ISOA1  150.00
Tracer #36              : 36   ISOA2  150.00
Tracer #37              : 37   ISOA3  150.00
Tracer #38              : 38   BENZ    12.00   (6C)
Tracer #39              : 39   TOLU    12.00   (7C)
Tracer #40              : 40   XYLE    12.00   (8C)
Tracer #41              : 41   ASOG1  150.00
Tracer #42              : 42   ASOG2  150.00
Tracer #43              : 43   ASOG3  150.00
Tracer #44              : 44   ASOAN  150.00
Tracer #45              : 45   ASOA1  150.00
Tracer #46              : 46   ASOA2  150.00
Tracer #47              : 47   ASOA3  150.00

This list can be reduced by turning off any aerosol species (eg. SOA, dust, SS, etc). You may also run without the secondary organic aerosols (MTPA, LIMO, MTPO, TSOG0–3, TSOA0–3, ISOG1–3, ISOA1–3, ASOG1–3, ASOA1–3, ASOAN) if you wish (but we recommend that you use them).

--Melissa Sulprizio 12:51, 5 November 2013 (EST)

Previous issues that have now been resolved

Fixed undefined molecular weight of HNO3 in routine SEASALT_CHEM

This fix was included in v11-01j and approved on 03 Dec 2016

In routine SEASALT_CHEM (in GeosCore/sulfate_mod.F), the MW_HNO3 variable was undefined for aerosol-only simulations, which do not use HNO3 as a defined species. This issue was detected by a unit test using the GNU Fortran compiler.

To fix this issue, we have added the following block of code, which defines MW_HNO3 whether or not HNO3 is a defined species:

     ! Get the molecular weight of HNO3 from the species database
     ! if it is defined, otherwise set it manually (bmy, 10/11/16)
     IF ( id_HNO3 > 0 ) THEN
        MW_HNO3 = State_Chm%SpcData(id_HNO3)%Info%emMW_g
     ELSE
        MW_HNO3 = 63.0_fp
     ENDIF

NOTE: This is just a quick fix. Ideally, we will want to declare HNO3 as a non-advected species in the GEOS-Chem species database for these simulationns.

--Bob Yantosca (talk) 16:47, 11 October 2016 (UTC)

Error in ISOROPIAII when sulfate aerosols are turned off

This update was tested in the 1-month benchmark simulation v9-02r and approved on 14 Nov 2013.

Chris Butenhoff wrote:

I set ONLINE SULFATE AEROSOLS = F and removed the sulfate aerosols from the tracer list. However I now get the error "IDTSO4 is undefined" from isoropia11_mod.F (line 207).

For full-chemistry simulations, we only call routine DO_ISOROPIAII if sulfate aerosols are used (i.e. LSULF = T). However, no such check exists when running the offline aerosol simulation. This was probably an oversight and we should add the following fix to routine DO_CHEMISTRY (in GeosCore/chemistry_mod.F):

           !*** AEROSOL THERMODYNAMIC EQUILIBRIUM ***
           IF ( LSULF ) THEN
              IF ( LSSALT ) THEN

#if   !defined( NO_ISORROPIA )
                 ! ISOROPIA takes Na+, Cl- into account
                 CALL DO_ISOROPIAII
    &               ( am_I_Root, Input_Opt, State_Met, State_Chm, RC )
#endif

              ELSE

                 ! RPMARES does not take Na+, Cl- into account
                 ! (skip for crystalline & aqueous offline run)
                 IF ( .not. LCRYST ) THEN
                    CALL DO_RPMARES( am_I_Root, Input_Opt, 
    &                                State_Met, State_Chm, RC )
                 ENDIF
              ENDIF
           ENDIF

--Bob Y. 15:29, 14 November 2013 (EST)

Known issues

ND05 diagnostic quantities zeroed unexpectedly

The binary punch format diagnostics will soon be made obsolete in a future GEOS-Chem version.

Helen McIntyre wrote:

I'm running the new and old versions of GEOS-Chem, and both seem to not output all the ND05 diagnostics correctly. There are 10 prod/loss diagnostics in this category, but only the 5th, 6th and 7th work. All the rest come out as zero.
I've had a brief look though the code, and it seems that the ones that work are calculated in one part of the routine (sulfate_mod.f), and the zero ones in another (I don't know if this has anything to do with it).
I've just done a 1 day run, with GEOS-Chem v8-01-01 at 4x5 resolution using GEOS-5 met fields. The old version I'm using is v7-02-04 and I get the same result.
The 'ctm.bpch', 'geos.log' and 'input.geos' files from the v8-01-01 run can be found here: http://homepages.see.leeds.ac.uk/~lechlm/files/

Claire Carouge replied:

There are some problems in the calculations. For the 5th element of AD05 (in Fortran notations and not IDL), the diagnostic come from the value of L1 which is calculated line 1529 (v8-01-04):
   L1     = ( SO20 - SO2_cd + PSO2_DMS(I,J,L) ) * RK1/RK
But on line 1508, we have:
   RK1 = 0.d0
with the previous comment:
  ! For online runs, SMVGEAR deals w/ this computation,
  ! so we can simply set RK1 = 0 (rjp, bmy, 3/23/03)
So L1 is always 0. I have no idea what RK1 is for, you may have to go into the physics/chemistry behind the code and we are not qualified to help you with it.
The other values for AD05 are calculated in the routine CHEM_DMS and this routine is only called in an offline aerosol simulation (l. 545 in sulfate_mod.f).
So my guess is that the ND05 diagnostic was designed for offline simulations and if you want it for online chemistry, you need to implement it.



Previous | Next | Guide to GEOS-Chem simulations