Benchmark/GEOS-Chem 12.5.0

From Geos-chem
Revision as of 14:10, 16 August 2019 by Melissa Payer (Talk | contribs) (Validation of netCDF diagnostics)

Jump to: navigation, search


This page includes information about all benchmark simulations for GEOS-Chem 12.5.0.

GEOS-Chem Classic 1-month benchmark

Benchmark plots

Enable PDF bookmarks1 in the files linked below to navigate between species. See the Species in GEOS-Chem wiki page for species definitions and details.

Species Category2 Plots Additional Information
Level Map Zonal Mean Emissions
Oxidants

O3, CO, OH, NOx

sfc 500hPa full column stratosphere sfc
Aerosols

NH4, NIT, SO4, DST1, DST2, DST3, DST4, BC, OC, AERI, BrSALA, BrSALC, ISALA, ISALC, NITs, SALA, SALC, SO4s, SOA_Simple, SOA_Complex

sfc 500hPa full column stratosphere sfc AOD
Bromine

Bry, BrOx, Br, Br2, BrCl, BrNO2, BrNO3, BrO, CH4Br, CH2Br2, CH2Br2, CH3Br3, HOBr, HBr

sfc 500hPa full column stratosphere sfc
Chlorine

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

sfc 500hPa full column stratosphere --
Iodine

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

sfc 500hPa full column stratosphere sfc
Nitrogen

NOy, NOx, HNO2, HNO3, HNO4, IMPN, NIT, NO, NO2, NO3, N2O5, MPN, PAN, PPN, N2O, NHx, NH3, NH4

sfc 500hPa full column stratosphere sfc
Primary Organics

ALK4, BENZ, CH4, C2H6, C2H4, C3H8, PRPE, TOLU, XYLE, ISOP, MTPA, MTPO, LIMO, GLYX, HCOOH, MAP, MBO, MOH, RCHO

sfc 500hPa full column stratosphere sfc
Secondary Organics

ACET, MEK, MVK, ALD2, CH2O, HPALD, MP, ISOPN, IEPOX, ACTA, EOH, MOH

sfc 500hPa full column stratosphere sfc
ROy

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

sfc 500hPa full column stratosphere --
Sulfur

SOx, DMS, OCS, SO2, SO4

sfc 500hPa full column stratosphere sfc
Emissions by HEMCO category

Total, Aircraft, Anthro, BioBurn, Biofuel, Biogenic,
Lightning, Monoterp, Natural, Ocean, Seabirds, DecayingPlants,
Ship, Soil, VolcDegas, VolcErupt,

Emissions tables
J-values sfc 500hPa full column stratosphere -- --
Plots still created by GAMAP
Budget -- -- -- -- -- Ox and CO budgets

1All plots created with GCPy have PDF bookmarks. If you do not see them in your browser try using a different browser (not Safari) or installing an Adobe Acrobat plug-in.

2Species category and lumped species definitions are simple to adjust. Updating them requires only editing the json files.

Benchmark assessment form

Description
New features added into GEOS-Chem
Feature Submitted by
Updates that will affect full-chemistry simulations
Retire obsolete emission inventories Eloise Marais (Leicester)
Emissions and Deposition Working Group
Updated offline biogenic emissions generated with GEOS-Chem 12.3.0 Hongjian Weng (PKU)
Updated volcano emissions to include eruptive emissions through April 2018 Christoph Keller (NASA/GMAO)
Fixes for GFED4 emissions, including:
  • Change the HEMCO cycling flag so that the model crashes if it doesn't find the right GFED year
  • Remove the 27% scaling of POG to POA and emit 1:1 instead
Sidhant Pai (MIT)
Scale DICE-Africa emissions to address errors in inventory Eloise Marais (Leicester)
Separate anthropogenic PM2.5 dust source from the DustDead extension Melissa Sulprizio (GCST)
Corrections for CFCs in the chemistry mechanism Xuan Wang (Harvard)
Fix typo in CLD_PARAMS routine in gckpp_HetRates.F90 Stephen Steenrod (NASA GMAO)
Bug fix: make sure stratospheric BrY concentrations are read properly each month Lyatt Jaeglé (U. Washington)
Updates that will NOT affect full-chemistry simulations
Enable accurate cloud pH diagnostic calculation via post-processing Jonathan Moch (Harvard)
Bug fix for proposed Henry's law constants Jenny Fisher (Wollongong)
Updates that will only affect GCHP
Update GCHP to use MAPL v1.0.0 and ESMF v8.0.0 beta snapshot 28 Lizzie Lundgren (GCST)
Version, resolution, met fields used: GEOS-FP (72L), 4x5, July 2016
1-month benchmark finished on: Sat Aug 10 06:32:29 EDT 2019
Performance statistics:
  • Compiled with ifort 17.0.4
  • Ran on 24 CPUs of holyjacob01.rc.fas.harvard.edu (Intel(R) Xeon(R) CPU E5-2680 v3 "Haswell" @ 2.50 GHz)
  • Wall time: 07:32:26
  • CPU time / wall time: 21.1184
  • % of ideal performance: 87.99%
  • Memory: 13.8486 GB
Compared to previous benchmark: 12.4.0
This update will impact:
(select all that apply with boldface)
Advection, BL Mixing, Convection, Met Fields, Dry Dep, Wet Dep, Stratosphere, Emissions, Photolysis, Chemistry, Other (please specify): Diagnostics
Unit test results may be viewed at: http://ftp.as.harvard.edu/gcgrid/geos-chem/1mo_benchmarks/GC_12/12.5.0/GCC/Plots/GC_12.5.0_Aug08.results.html
Metrics
Global mean OH from log file (% change): 12.3699948768189 x 105 molec/cm3 (-0.82%)
Methyl chloroform lifetime (% change): 4.9566 (1.22%)
Methane lifetime (% change): 8.3442 years (1.26%)
List all SURFACE species, J-Values, and COLUMN AOD values that changed by 10% or more:
  • Oxidants: O3 CO OH NOx
  • Aerosols: NH4 NIT SO4 DST1 AERI BrSALA BrSALC ISALA ISALC NITs SO4s SOA_Complex SOA_Simple
  • Nitrogen: NOy NOx HNO2 HNO3 HNO4 IPMN NIT NO NO2 NO3 N2O5 MPN PAN PPN NHx NH3 NH4
  • ROy: H2O2 H HO2 O1D OH RO2
  • Sulfur: SOx DMS SO2 SO4
  • Iodine: Iy IxOy I I2 IBr ICl IO INO2 IONO2 CH2I2 HI HOI OIO
  • Bromine: Bry BrOx Br Br2 BrCl BrNO2 BrNO3 BrO HOBr HBr
  • Chlorine: Cly ClOx Cl ClO Cl2 Cl2O2 ClOO ClNO2 ClNO3 HOCl HCl HCFCs OClO
  • Primary_Organics: ALK4 BENZ C2H6 C3H8 PRPE TOLU XYLE ISOP MTPA MTPO LIMO GLYX HCOOH MAP RCHO
  • Secondary_Organics: ACET MEK MVK ALD2 CH2O HPALD MP ISOPN IEPOX ACTA EOH
  • J-Values: CCl4 CFC11 CFC113 CFC114 CFC115 CFC12 CH2Cl2 CH3CCl3 CH3Cl HCFC123 HCFC141b HCFC142b N2O O2
  • Column AOD: Dust OCPI SO4 SOA_from_Aqueous_ISOP Total
List all 500 hPa species and J-Values that changed by 10% or more :
  • Oxidants: OH NOx
  • Aerosols: NH4 NIT SO4 DST1 AERI BrSALA BrSALC ISALC NITs SO4s SOA_Complex
  • Nitrogen: NOx HNO2 HNO4 IPMN NIT NO NO2 NO3 N2O5 MPN PAN PPN NHx NH3 NH4
  • ROy: OH
  • Sulfur: SOx DMS SO2 SO4
  • Iodine: I2 ICl INO2 IONO2 CH2I2 HI OIO
  • Bromine: BrOx Br Br2 BrCl BrNO2 BrNO3 BrO HOBr HBr
  • Chlorine: ClOx Cl ClO Cl2 Cl2O2 ClOO ClNO2 ClNO3 HOCl HCl OClO
  • Primary_Organics: ALK4 BENZ C2H6 C3H8 PRPE TOLU XYLE ISOP MTPA MTPO LIMO GLYX HCOOH MAP RCHO
  • Secondary_Organics: ACET MEK MVK ALD2 HPALD MP ISOPN IEPOX ACTA EOH
  • J-Values: None
List all ZONAL MEAN species and J-Values that changed by 10% or more:
  • Oxidants: OH NOx
  • Aerosols: NH4 NIT SO4 DST1 BrSALA BrSALC NITs SO4s SOA_Complex
  • Nitrogen: NOy NOx HNO2 HNO3 HNO4 IPMN NIT NO NO2 NO3 N2O5 MPN PAN PPN NH3 NH4
  • ROy: OH RO2
  • Sulfur: SOx SO2 SO4
  • Iodine: I2 IBr ICl INO2 IONO2 CH2I2 HI OIO
  • Bromine: BrOx Br Br2 BrCl BrNO2 BrNO3 BrO HOBr HBr
  • Chlorine: Cly ClOx Cl ClO Cl2 Cl2O2 ClOO ClNO2 ClNO3 HOCl OClO
  • Primary_Organics: ALK4 C2H6 C3H8 PRPE TOLU XYLE ISOP MTPA MTPO LIMO GLYX HCOOH MAP RCHO
  • Secondary_Organics: ACET MEK MVK ALD2 CH2O HPALD ISOPN IEPOX ACTA EOH
  • J-Values: None
List all EMISSIONS that changed by 10% or more:
  • Anthro: ALD2 ALK4 BCPI BCPO C2H6 C3H8 CO HAC MEK NH3 NO pFe PRPE SO2 SO4 SOAP
  • Biogenic: ACET ALD2 EOH ISOP LIMO PRPE SOAP SOAS
  • Monoterp: CO
  • Ship: CO HNO3 NO2 NO O3 SO2 SOAP
  • Total: ACET ALD2 ALK4 BCPI BCPO C2H6 C3H8 CO EOH GLYX ISOP MEK MGLY MTPA NH3 NO pFe PRPE SO2 SO4 SOAP
  • VolcDegas: SO2
  • VolcErupt: SO2
Comments on observed differences: Below we summarize the notable changes caused by specific updates.

Retire obsolete emission inventories

Retiring the EMEP, Bond, RETRO, and Yevich & Logan biofuel inventories results in large fractional differences in anthropogenic emissions and species concentrations over Europe, South America, Central America, Russia, and Australia, where the global anthropogenic, biofuel, and ship emissions from the CEDS inventory are now used. The impact on emission totals are included below.
                                   Ref            Dev            Dev - Ref
  ACET Anthro               :      0.171664       0.069874      -0.101790 Tg C
  ALD2 Anthro               :      0.042310       0.042140      -0.000169 Tg C
  ALK4 Anthro               :      2.407840       2.572299       0.164460 Tg C
  BCPI Anthro               :      0.106725       0.106747       0.000021 Tg
  BCPO Anthro               :      0.426900       0.426986       0.000086 Tg
  C2H6 Anthro               :      0.473419       0.664990       0.191571 Tg C
  C3H8 Anthro               :      0.399813       0.595361       0.195548 Tg C
  CH2O Anthro               :      0.263482       0.235037      -0.028445 Tg
  CO Anthro                 :     42.984464      40.778317      -2.206146 Tg
  HAC Anthro                :      0.212160       0.042432      -0.169728 Tg
  HCOOH Anthro              :      0.577087       0.571744      -0.005343 Tg
  MEK Anthro                :      0.043368       0.064733       0.021364 Tg C
  NH3 Anthro                :      4.430414       4.592343       0.161929 Tg
  NO Anthro                 :      5.543591       5.293092      -0.250499 Tg
  OCPI Anthro               :      0.628814       0.610537      -0.018278 Tg
  OCPO Anthro               :      0.628815       0.610537      -0.018278 Tg
  pFe Anthro                :      0.006943       0.006628      -0.000315 Tg
  PRPE Anthro               :      0.490600       0.448865      -0.041735 Tg C
  SO2 Anthro                :      6.974655       7.092453       0.117797 Tg
  SO4 Anthro                :      0.218581       0.210372      -0.008208 Tg
  SOAP Anthro               :      2.978745       2.826521      -0.152224 Tg
  TOLU Anthro               :      0.584716       0.553756      -0.030960 Tg C
  CO Ship                   :      0.079328       0.059695      -0.019634 Tg
  HNO3 Ship                 :      1.484222       1.449563      -0.034659 Tg
  NO Ship                   :      0.081681       0.082875       0.001195 Tg
  NO2 Ship                  :      0.771094       0.738918      -0.032176 Tg
  SO2 Ship                  :      0.805693       0.773045      -0.032648 Tg
  SOAP Ship                 :      0.005474       0.004119      -0.001355 Tg


Updated offline biogenic emissions generated with GEOS-Chem 12.3.0

The offline biogenic emissions were updated to include recent fixes to the MEGAN extension. The impact on biogenic emission totals are included below.
  ###############################################################################
  ### Inventory totals for OFFLINEBIOVOC                                      ###
  ### Ref = GC_12.4.0; Dev = GC_12.5.0                                        ###
  ###############################################################################
                                   Ref            Dev            Dev - Ref
  OFFLINEBIOVOC ACET        :      2.677313       2.982478       0.305164 Tg C
  OFFLINEBIOVOC ALD2        :      1.032241       1.037183       0.004942 Tg C
  OFFLINEBIOVOC CO          :      5.504973      11.207163       5.702190 Tg C
  OFFLINEBIOVOC EOH         :      0.987091       0.991817       0.004726 Tg C
  OFFLINEBIOVOC ISOP        :     35.759849      33.723910      -2.035939 Tg C
  OFFLINEBIOVOC LIMO        :      1.022397       0.976173      -0.046225 Tg C
  OFFLINEBIOVOC PRPE        :      2.126864       2.030293      -0.096571 Tg C
  OFFLINEBIOVOC SOAP        :      0.621433       0.102510      -0.518922 Tg C
  OFFLINEBIOVOC SOAS        :      0.621433       0.102510      -0.518922 Tg C


Updated volcano emissions to include eruptive emissions through April 2018

This update increased both eruptive and degassing volcanic emissions, as summarized below. As a result, localized increases in SO2 and SO4 can be seen in the species concentration plots.
  ###############################################################################
  ### Inventory totals for VOLCANO                                            ###
  ### Ref = GC_12.4.0; Dev = GC_12.5.0                                        ###
  ###############################################################################
                                   Ref            Dev            Dev - Ref
  VOLCANOerupt SO2          :      0.000000       0.035000       0.035000 Tg
  VOLCANOdegas SO2          :      1.509524       1.951957       0.442434 Tg


Scale DICE-Africa emissions to address errors in inventory

This update contributes to changes in BC and OC concentrations over and just downwind of Africa. Slight differences can also be seen over Africa in the emission maps. The impact on the inventory emission totals are included below.
  ###############################################################################
  ### Inventory totals for DICEAfrica                                         ###
  ### Ref = GC_12.4.0; Dev = GC_12.5.0                                        ###
  ###############################################################################
                                   Ref            Dev            Dev - Ref
  DICEAfrica BCPI           :      0.004131       0.004153       0.000022 Tg
  DICEAfrica BCPO           :      0.016526       0.016612       0.000086 Tg
  DICEAfrica CH2O           :      0.044793       0.016348      -0.028445 Tg
  DICEAfrica CO             :      5.651943       3.648702      -2.003242 Tg
  DICEAfrica HAC            :      0.212160       0.042432      -0.169728 Tg
  DICEAfrica HCOOH          :      0.006679       0.001336      -0.005343 Tg
  DICEAfrica NH3            :      0.294221       0.284477      -0.009744 Tg
  DICEAfrica NO             :      0.153533       0.151820      -0.001713 Tg
  DICEAfrica OCPI           :      0.035185       0.016907      -0.018278 Tg
  DICEAfrica OCPO           :      0.035185       0.016907      -0.018278 Tg
  DICEAfrica PRPE           :      0.057572       0.019730      -0.037842 Tg
  DICEAfrica SOAP           :      0.389982       0.251758      -0.138224 Tg
  DICEAfrica TOLU           :      0.054817       0.020857      -0.033960 Tg


Separate anthropogenic PM2.5 dust source from the DustDead extension

This update reintroduces the anthropogenic PM2.5 fine dust source and increases concentrations of DST1. The impact on DST1 emission totals are included below.
  ###############################################################################
  ### Inventory totals for AFCID                                              ###
  ### Ref = GC_12.5.0; Dev = GC_12.5.0                                        ###
  ###############################################################################
                                   Ref            Dev            Dev - Ref
  AFCID DST1                :      0.000000       1.087876       1.087876 Tg


The following updates produced small changes, primarily over the southern hemisphere oceans and at high altitudes:

Approval
Requires further investigation: TBD
Approved by:
Date of approval:

GCHP 1-month benchmark

The following sections contain validation plots for GCHP 12.5.0.

Comparison of GCHP versions

This table contains plots comparing GCHP 12.4.0 (Ref) and GCHP 12.5.0 (Dev) for a 1-month simulation.

Category Plots Additional Information
Level Map Zonal Mean Emissions
Oxidants sfc 500hPa full column stratosphere full column
Aerosols sfc 500hPa full column stratosphere full column Column AOD
Bromine sfc 500hPa full column stratosphere full column
Chlorine sfc 500hPa full column stratosphere --
Iodine sfc 500hPa full column stratosphere full column
Nitrogen sfc 500hPa full column stratosphere full column
Primary Organics sfc 500hPa full column stratosphere full column
Secondary Organics sfc 500hPa full column stratosphere full column
Sulfur sfc 500hPa full column stratosphere full column
Emissions by HEMCO category

Aircraft, Anthro, Bioburn, Biofuel, Biogenic,
Lightning, Monoterp, Natural, Ocean, Seabirds, DecayingPlants,
Ship, Soil, Total, VolcDegas, VolcErupt,

Emission totals

Comparison of GCHP and GEOS-Chem Classic

This table contains plots comparing GEOS-Chem Classic 12.5.0 (Ref) against GCHP 12.5.0 (Dev) for a 1-month simulation.

Category Plots Additional Information
Level Map Zonal Mean Emissions
Oxidants sfc 500hPa full column stratosphere full column
Aerosols sfc 500hPa full column stratosphere full column Column AOD
Bromine sfc 500hPa full column stratosphere full column
Chlorine sfc 500hPa full column stratosphere --
Iodine sfc 500hPa full column stratosphere full column
Nitrogen sfc 500hPa full column stratosphere full column
Primary Organics sfc 500hPa full column stratosphere full column
Secondary Organics sfc 500hPa full column stratosphere full column
Sulfur sfc 500hPa full column stratosphere full column
Emissions by HEMCO category

Aircraft, Anthro, Bioburn, Biofuel, Biogenic,
Lightning, Monoterp, Natural, Ocean, Seabirds, DecayingPlants,
Ship, Soil, Total, VolcDegas, VolcErupt,

Emission totals

GCHP and GEOS-Chem Classic difference of differences

This table contains plots comparing the difference [GCHP 12.5.0 - GCHP 12.4.0] against the difference [GC Classic 12.5.0 - GC Classic 12.4.0], taken from 1-month benchmark simulations.

Category Plots
Level Map Zonal Mean
Oxidants sfc 500hPa full column stratosphere
Aerosols sfc 500hPa full column stratosphere
Bromine sfc 500hPa full column stratosphere
Chlorine sfc 500hPa full column stratosphere
Iodine sfc 500hPa full column stratosphere
Nitrogen sfc 500hPa full column stratosphere
Primary Organics sfc 500hPa full column stratosphere
Secondary Organics sfc 500hPa full column stratosphere
Sulfur sfc 500hPa full column stratosphere

NOTE: At present, we have only created the diff-of-diffs plots for species concentration. This can be extended to other categories in the future.

Validation of netCDF diagnostics

In this version, we did a complete comparison of the BPCH diagnostics to the netCDF diagnostics in preparation for retiring the old BPCH diagnostics. Comparison plots can be accessed via the links below.

NetCDF collection BPCH equivalent Plots
AerosolMass ND42 sfc full column
CloudConvFlux ND14 sfc
ConcAfterChem ND43 sfc full column
DryDep ND44 sfc
Emissions ND06, ND07, ND08, ND11,
ND13, ND28, ND29, ND32,
ND34, ND36, ND46
sfc
JValues ND22 sfc full column
ProdLoss ND05, ND07, ND65 sfc full column
SpeciesConc ND45 sfc full column
ND66, StateMet ND67, ND68 sfc full column
WetLossConv ND38, ND37 850hPa full column
WetLossLS ND39 850hPa full column