Difference between revisions of "GEOS-Chem 13.3.0"

• See geoschem/geos-chem PR #919

• See geoschem/geos-chem PR #904

• See geoschem/geos-chem PR #904

• See geoschem/geos-chem PR #905

• See geoschem/geos-chem PR #719
• See geoschem/geos-chem PR #923

• See geoschem/geos-chem #907

• Fixed typo in CloudHet
• Fixed typos in ClNO2 hetchem reactions: HNO2 should be a product, not HNO3

• See geoschem/geos-chem #906
• See geoschem/geos-chem #841

• See geoschem/geos-chem PR #897

• See geoschem/geos-chem PR #604

• See also geoschem/HEMCO PR #93

• See geoschem/geos-chem PR #921
• See geoschem/HEMCO PR #110

• See geoschem/geos-chem PR #816
• Online dust emissions are used in benchmark simulations

• See geoschem/geos-chem PR #920

• See geoschem/HEMCO PR #111

• See geoschem/geos-chem PR #922

• See geoschem/geos-chem PR #925
• See geoschem/HEMCO PR #111

• See geoschem/GCHP #146

• GEOS-Chem Classic: 13.3.0-rc.0
• GCHP: TBD

• GEOS-Chem Classic: 4° x 5°
• GCHP: c48

• GEOS-Chem Classic: 4° x 5°
• GCHP: 0.5° x 0.625°

• NOTE: GCHP c48 restart file regridded from GEOS-Chem Classic 4° x 5°

• Advection
• BL Mixing
• Convection
• Meteorology Fields
• Dry Deposition
• Wet Deposition
• Emissions
• Photolysis
• Chemistry
• Other (please specify):

O3, CO, OH, NOx

DST1, DST2, DST3, DST4, NH4, NIT, SO4, BCPI, BCPO, OCPI, OCPO, AERI, BrSALA, BrSALC, ISALA, ISALC, NITs, SALA, SALC, SO4s, Simple_SOA, Complex_SOA. pFe

Bry, BrOx, Br, Br2, BrCl, BrNO2, BrNO3, BrO, CH3Br, CH2Br2, CHBr3, HOBr, HBr

Cly, ClOx, Cl, ClO, Cl2, Cl2O2, ClOO, ClNO2, ClNO3, CCl4, CFCs, CH3Cl, CH2Cl2, CH3CCl3, CHCl3, HOCl, HCl, Halons, HCFCs, OClO

Iy, IxOy, I, I2, IBr, ICl, IO, IONO, IONO2, CH3I, CH2I2, CH2ICl, CH2IBr, HI, HOI, OIO

NOy, NOx, HNO2, HNO3, HNO4, MPAN, NIT, NO, NO2, NO3, N2O5, MPN, PAN, PPN, N2O, NHx, NH3, NH4, MENO3, ETNO3, IPRNO3, NPRNO3

EOH, MOH, ISOP, MTPA, MTPO, LIMO, ALK4, BENZ, CH4, C2H6, C3H8, PRPE, TOLU, XYLE

ACTA, ALD2, CH2O, HPALDs, MACR, IEPOX, ACET, HAC, MEK, MVK, ISOPN, GLYC, GLYX, HCOOH, MAP, RCHO, MP

TSOA0, TSOA1, TSOA2, TSOA3, ASOA1, ASOA2, ASOA3, ASOAN, TSOG0, TSOG1, TSOG2, TSOG3, ASOG1, ASOG2, ASOG3, INDIOL, LVOCOA, SOAIE, SOAGX, SOAP, SOAS

SOx, DMS, OCS, SO2, SO4

H2O2, H, H2, H2O, HO2, O1D, OH, RO2

Aircraft, Anthropogenic, Biomass Burning, Biogenic, Lightning, Natural, Ocean, Seabirds, Decaying Plants, Ship, Soil, Degassing Volcanoes, Erupting Volcanoes

Tables by inventory

Troposphere only

1. HMS chemistry (via sulfate_mod.F90)

This update introduces new advected species hydroxymethanesulfonate (HMS) in GEOS-Chem. This new chemistry has a relatively small impact on other species' concentrations and on global mean OH. Complete results from a 1-month benchmark evaluating the impact of this chemistry update may be viewed at this link.

2. C2H2 and C2H4 chemistry

This update introduces new advected species C2H2, C2H4, ETHN, and ETHP and adds one photolysis reaction. The new chemistry has a large (>10% impact) on GLYC, GLYX, H2O2, HCOOH, HMHP, MO2, MOH, MP, SOAGX (see this table). Global mean OH is reduced by -1.49%. Complete results from a 1-month benchmark evaluating the impact of this chemistry update may be viewed at this link.

3. Fixed typo in CloudHet

This fix has a large impact on many halogen species and changed global mean OH by +6.65%. These large changes are to be expected, because this bug was in the CloudHet routine (which is used to adjust heterogeneous reaction rates based on the amount of cloud in the grid box). See this Github issue for further discussion on the changes from this fix.

4. Revise gamma(N2O5) computation; cap A at 3.2e-8

Prior to this fix, gamma(N2O5) (i.e. reaction probability) was too high over the oceans. This update now caps gamma(N2O5) at ~0.036, which is consistent with lab results. This update has a relatively small (<10%) impact on most species, but ClNO2 shows a large decrease. However, this is to be expected because when gamma(N2O5) decreases, so will the production of ClNO2 via the reaction N2O5 + HCl = ClNO2 + HNO3. See this Github issue for further discussion on the changes from this update.

5. Updates to IONO2 hydrolysis from Wang et al. (2021)

This update replaced heterogeneous reactions that converted gas-phase species HOI, IONO, and IONO2 to aerosol species ISALA and ISALC with the single hydrolysis reaction IONO2 + H2O = HOI + HNO3. We observe a large (>10%) impact on IONO2 and HOI, but a relatively small impact on other species (see this table). Global mean OH changed by -0.11%. See the original pull request on Github for more details of this update and its impact on model output.

6. Aromatic VOC chemistry

This update added new advected species AONITA, AROMP4, AROMP5, BALD, BENZP, BZCO3H, BZPAN, CSL, MCT, NPHEN, PHEN, along with ~45 gas-phase reactions, 1 heterogeneous reaction, and 4 photolysis reactions. The impact on other species is summarized in this table, and global mean OH changed by -0.36%. Complete results from a 1-month benchmark evaluating the impact of this chemistry update may be viewed at this link.

7. Updated aerosol dry deposition

This feature brings aerosol dry deposition in line with Emerson et al 2020; several aerosol parameters within the dry deposition module were changed. Global mean OH increased by 0.024%. Most gas-phase species show small changes in mass (except for Br2 and ClNO2, which decreased by more than 5%). On the other hand, aerosol species such as AERI, ASOA1-4, DST1-4, show large decreases. Also, the mass of MSA increased by 7.7%. See this Github feature request for more information about this update. Complete results from a 1-month benchmark evaluating the impact of this chemistry update may be viewed at this link.

8. Update volcano emissions to May 2020 and add climatology

The volcano emissions inventory from GMAO has been updated from v20180523 to v202005. Plots for a 1-month benchmark evaluating the impact of this fix may be viewed at this link. The impact on SO2 emissions from volcanoes is:

    ###################################################################################
    ### Emissions totals for inventory VOLCANOdegas [Tg]                            ###
    ### Ref = GCC_13.3.0-alpha.7; Dev = GCC_13.3.0-alpha.8                          ###
    ###################################################################################
                                        Ref                 Dev     Dev - Ref    % diff
    VOLCANOdegas SO2   :           1.953161            1.826465     -0.126697    -6.487
  
    ###################################################################################
    ### Emissions totals for inventory VOLCANOerupt [Tg]                            ###
    ### Ref = GCC_13.3.0-alpha.7; Dev = GCC_13.3.0-alpha.8                          ###
    ###################################################################################
                                        Ref                 Dev     Dev - Ref    % diff
    VOLCANOerupt SO2   :           0.000000            0.101063      0.101063       nan
    

9. Avoid double counting of acetone source from monoterpenes when using online MEGAN

This update reduced ACET emissions in the MEGAN extension. Plots for a 1-month benchmark evaluating the impact of this fix may be viewed at this link. The impact on total ACET emissions from MEGAN is:

    ###################################################################################
    ### Emissions totals for inventory MEGAN [Tg]                                   ###
    ### Ref = GCC_13.2.1; Dev = GCC_13.3.0-alpha.0                                  ###
    ###################################################################################
                                        Ref                 Dev     Dev - Ref    % diff
    MEGAN ACET         :           4.967504            4.374804     -0.592700   -11.932

The following updates introduced differences on the order of numerical noise:

• Cleaning up heterogeneous chemistry reaction rate computations
• Fixed typos in ClNO2 hetchem reactions: HNO2 should be a product, not HNO3