|
|
| (55 intermediate revisions by 4 users not shown) |
| Line 1: |
Line 1: |
| − | On this page, we provide an overview of the chemistry mechanisms used in GEOS-Chem. | + | __FORCETOC__ |
| | + | '''''[[Guide to GEOS-Chem simulations|Previous]] | [[Aerosol-only simulation|Next]] | [[Guide to GEOS-Chem simulations]]''''' |
| | + | #<span style="color:blue">'''Simulations using KPP-built mechanisms'''</span> |
| | + | #[[Aerosol-only simulation]] |
| | + | #[[CH4 simulation]] |
| | + | #[[CO2 simulation]] |
| | + | #[[Mercury|Hg simulation]] |
| | + | #[[POPs simulation]] |
| | + | #[[Tagged CO simulation]] |
| | + | #[[Tagged O3 simulation]] |
| | + | #[[TransportTracers simulation]] |
| | + | |
| | + | On this page, we provide information about GEOS-Chem simulations that use chemistry mechanism solver code built by the [https://kpp.readthedocs.io Kinetic PreProcessor (KPP)]. |
| | | | |
| | == Overview == | | == Overview == |
| | | | |
| − | The following table provides links to information about the available chemistry mechanisms in GEOS-Chem. Please contact the relevant [http://acmg.seas.harvard.edu/geos/geos_working_groups.html GEOS-Chem Working Group] for more information. | + | The following table provides links to information about the available full-chemistry mechanisms in GEOS-Chem. |
| − | | + | |
| | {| border=1 cellspacing=0 cellpadding=5 | | {| border=1 cellspacing=0 cellpadding=5 |
| | |-bgcolor="#CCCCCC" | | |-bgcolor="#CCCCCC" |
| − | !width="200px"|Category | + | !width="100px"|Mechanism |
| − | !width="175px"|Simulation(s) | + | !width="300px"|Description |
| − | !wigth="350px"|Mechanism file | + | !width="250px"|Mechanism file |
| − | !width="225px"|Contact | + | !width="175px"|Extra options |
| | | | |
| | |-valign="top" | | |-valign="top" |
| − | |Full-chemistry<br>([[Tropospheric_chemistry_mechanism#Recent_chemistry_updates|troposphere]] + [[UCX_chemistry_mechanism|stratosphere]]) | + | |fullchem |
| − | |Standard, Benchmark<sup>1</sup>
| + | |NOx + Ox + Br + Cl + I + aerosols chemistry in the [[Tropospheric_chemistry_mechanism|troposphere]] and [[UCX_chemistry_mechanism|stratosphere]] |
| − | |[https://github.com/geoschem/geos-chem/blob/master/KPP/Standard/Standard.eqn <tt>KPP/Standard/Standard.eqn</tt>] | + | |[https://github.com/geoschem/geos-chem/blob/main/KPP/fullchem/fullchem.eqn <tt>KPP/fullchem/fullchem.eqn</tt>] |
| − | |[[Chemistry Working Group]] | + | | |
| | + | *[[GEOS-Chem_benchmarking|benchmark]] |
| | + | *complex SOA |
| | + | *complex SOA + SVPOA |
| | + | *[[#TOMAS|TOMAS]] |
| | + | *[[#APM|APM]] |
| | + | *[[Coupling_GEOS-Chem_with_RRTMG|RRTMG]] |
| | + | *Aciduptake |
| | + | *Marine POA |
| | | | |
| | |-valign="top" | | |-valign="top" |
| − | |Full-chemistry<br>([[Tropospheric_chemistry_mechanism#Recent_chemistry_updates|troposphere]] only) | + | |Hg |
| − | |Tropchem<br>[[Secondary_organic_aerosols#Simple_SOA_scheme|SOA]]<br>[[Secondary_organic_aerosols#Complex_SOA_scheme|complexSOA]]<br>[[#TOMAS|TOMAS]]<br>[[#APM|APM]]<br>[[Coupling_GEOS-Chem_with_RRTMG|RRTMG]] | + | |[[Mercury|Mercury chemistry]] |
| − | |[https://github.com/geoschem/geos-chem/blob/master/KPP/Tropchem/Tropchem.eqn <tt>KPP/Tropchem/Tropchem.eqn</tt>] | + | *Introduced in [[GEOS-Chem 13.4.0|13.4.0]] as a KPP mechanism |
| − | |[[Chemistry Working Group]] | + | |[https://github.com/geoschem/geos-chem/blob/main/KPP/fullchem/fullchem.eqn <tt>KPP/Hg/Hg.eqn</tt>] |
| | + | | |
| | | | |
| | |-valign="top" | | |-valign="top" |
| − | |SOA+semivoatile POA | + | |carboncycle |
| − | |[[Secondary_organic_aerosols#SOA_simulation_with_semi-volatile_POA|complexSOA_SVPOA]]
| + | |[https://github.com/geoschem/geos-chem/blob/main/KPP/carboncycle/carboncycle.eqn <tt>KPP/carboncycle/carboncycle.eqn</tt>] |
| − | |[https://github.com/geoschem/geos-chem/blob/master/KPP/SOA_SVPOA/SOA_SVPOA.eqn <tt>KPP/SOA_SVPOA/SOA_SVPOA.eqn</tt>] | + | *Will debut in [[GEOS-Chem 14.1.0|14.1.0]] as a KPP mechanism |
| − | |[[Aerosols Working Group]]
| + | | |
| − | | + | | |
| − | |-valign="top"
| + | |
| − | |Carbon Gases
| + | |
| − | |[[CH4 simulation|CH4]]<br>[[CH4 simulation|tagCH4]]
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/global_ch4_mod.F <tt>GeosCore/global_ch4_mod.F</tt>]
| + | |
| − | |[[Carbon Cycle Working Group]] | + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Carbon Gases
| + | |
| − | |[[Tagged CO simulation|tagCO]]
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/tagged_co_mod.F <tt>GeosCore/tagged_co_mod.F</tt>]
| + | |
| − | |[[Carbon Cycle Working Group]]
| + | |
| − | | + | |
| − | |-valign="top" | + | |
| − | |Carbon Gases | + | |
| − | |[[CO2 simulation|CO2]]
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/co2_mod.F <tt>GeosCore/co2_mod.F</tt>]
| + | |
| − | |[[Carbon Cycle Working Group]]
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Mercury
| + | |
| − | |[[Mercury|Hg]]<br>tagHg<br>Hg+[[Global Terrestrial Mercury Model|GTMM]]
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/mercury_mod.F <tt>GeosCore/mercury_mod.F</tt>]
| + | |
| − | |[[Hg and POPs Working Group]]
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Persistent Organic Pollutants
| + | |
| − | |[[POPs_simulation|POPs]]
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/pops_mod.F <tt>GeosCore/pops_mod.F</tt>]
| + | |
| − | |[[Hg and POPs Working Group]]
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Ozone
| + | |
| − | |[[Tagged O3 simulation|tagO3]]
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/tagged_o3_mod.F <tt>GeosCore/tagged_o3_mod.F</tt>]
| + | |
| − | |[[Chemistry Working Group]]
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Radionuclides
| + | |
| − | |[[Rn-Pb-Be simulation|Rn-Pb-Be]]<br>[[Transport_Working_Group#Transport_Tracers_simulation|TransportTracers]]
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/RnPbBe_mod.F <tt>GeosCore/RnPbBe_mod.F</tt>]
| + | |
| − | |[[Transport Working Group]]
| + | |
| − | | + | |
| − | |-
| + | |
| − | !colspan="5" bgcolor="#CCCCCC"|<span style="color:red">'''The following mechanisms are obsolete and have been removed:'''</span>
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Carbon Gases
| + | |
| − | |C2H6
| + | |
| − | |[https://github.com/geoschem/geos-chem/blob/master/GeosCore/c2h6_mod.F <tt>GeosCore/c2h6_mod.F</tt>]
| + | |
| − | |[[Carbon Cycle Working Group]]
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Carbon Gases
| + | |
| − | |CH3I
| + | |
| − | |<tt>GeosCore/ch3i_mod.F</tt><br>in [[GEOS-Chem v9-02]] and earlier
| + | |
| − | |[[Carbon Cycle Working Group]]
| + | |
| − | | + | |
| − | |-valign="top"
| + | |
| − | |Radionuclides
| + | |
| − | |H2-HD
| + | |
| − | |<tt>GeosCore/h2_h2_mod.F</tt><br>in [[GEOS-Chem v9-02]] and earlier
| + | |
| − | |[[Transport Working Group]]
| + | |
| | | | |
| | |} | | |} |
| | | | |
| − | <sup>1</sup>The benchmark simulation is used for [[GEOS-Chem_benchmarking|1-month and 1-year benchmarks]]. It uses the '''Standard''' chemistry mechanism, but includes both the [[Secondary_organic_aerosols#Simple_SOA_scheme|simple SOA]] and [[Secondary_organic_aerosols#Complex_SOA_scheme|complex SOA]] species.
| + | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 14:22, 20 September 2022 (UTC) |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 17:02, 22 February 2019 (UTC)
| + | |
| − | | + | |
| − | == Chemistry updates ==
| + | |
| − | | + | |
| − | === Updated isoprene and monoterpene chemistry ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02c|v11-02c]] and approved on 21 Sep 2017.'''''</span>
| + | |
| − | | + | |
| − | '''Developers:'''
| + | |
| − | * Katie Travis (MIT, formerly Harvard)
| + | |
| − | * Jenny Fisher (U. Wollongong)
| + | |
| − | * Christopher Chan Miller (Smithsonian Astrophysical Observatory, formerly Harvard)
| + | |
| − | * Eloise Marais (U. Birminghan, formerly Harvard)
| + | |
| − | | + | |
| − | [[Media:NewChemistry 030917.pdf|'''This document''']] compiled by Katie Travis and Josh Cox describes the updated isoprene and monoterpene chemistry to be included in [[GEOS-Chem v11-02#v11-02c|GEOS-Chem v11-02c]] (also see the [[#Modifications to the original updates|list of modifications below]]). These updates include the [[monoterpene nitrate scheme]] and [[Secondary_organic_aerosols#SOA_formation_from_aqueous_isoprene_uptake|aqueous isoprene uptake]] and were originally implemented for simulation of the [https://www.nasa.gov/seac4rs/ SEAC4RS] data.
| + | |
| − | | + | |
| − | '''References'''
| + | |
| − | *Chan Miller, C., D.J.Jacob, E.A. Marais, K. Yu, K.R. Travis, P.S. Kim, J.A. Fisher, L. Zhu, G.M. Wolfe, F.N. Keutsch, J. Kaiser, K.-E. Min, S.S. Brown, R.A. Washenfelder, G. Gonzalez Abad, and K. Chance, Glyoxal yield from isoprene oxidation and relation to formaldehyde: chemical mechanism, constraints from SENEX aircraft observations, and interpretation of OMI satellite data, Atmos. Chem. Phys., 17, 8725-8738, https://doi.org/10.5194/acp-17-8725-2017, 2017. [http://acmg.seas.harvard.edu/publications/2016/miller2016_seac4rs.pdf PDF]
| + | |
| − | *Fisher, J.A., D.J. Jacob, K.R. Travis, P.S. Kim, E.A. Marais, C. Chan Miller, K. Yu, L. Zhu, R.M. Yantosca, M.P. Sulprizio, J. Mao, P.O. Wennberg, J.D. Crounse, A.P. Teng, T.B. Nguyen, J.M. St. Clair, R.C. Cohen, P. Romer, B.A. Nault, P.J. Wooldridge, J.L. Jimenez, P. Campuzano-Jost, D.A. Day, P.B. Shepson, F. Xiong, D.R. Blake, A.H. Goldstein, P.K. Misztal, T.F. Hanisco, G.M. Wolfe, T.B. Ryerson, A. Wisthaler, and T. Mikoviny. Organic nitrate chemistry and its implications for nitrogen budgets in an isoprene- and monoterpene-rich atmosphere: constraints from aircraft (SEAC4RS) and ground-based (SOAS) observations in the Southeast US. Atmos. Chem. Phys., 16, 2961-2990, 2016. [http://acmg.seas.harvard.edu/publications/2016/fisher2016.pdf PDF]
| + | |
| − | *Marais, E. A., D. J. Jacob, J. L. Jimenez, P. Campuzano-Jost, D. A. Day, W. Hu, J. Krechmer, L. Zhu, P. S. Kim, C. C. Miller, J. A. Fisher, K. Travis, K. Yu, T. F. Hanisco, G. M. Wolfe, H. L. Arkinson, H. O. T. Pye, K. D. Froyd, J. Liao, V. F. McNeill, Aqueous-phase mechanism for secondary organic aerosol formation from isoprene: application to the southeast United States and co-benefit of SO2 emission controls, Atmos. Chem. Phys., 16, 1603-1618, 2016. [http://acmg.seas.harvard.edu/publications/2016/Marais_SEUS_isopSOA_ACP_2016.pdf PDF]
| + | |
| − | *Travis, K. R., D. J. Jacob, J. A. Fisher, P. S. Kim, E. A. Marais, L. Zhu, K. Yu, C. C. Miller, R. M. Yantosca, M. P. Sulprizio, A. M. Thompson, P. O. Wennberg, J. D. Crounse, J. M. St. Clair, R. C. Cohen, J. L. Laughner, J. E. Dibb, S. R. Hall, K. Ullmann, G. M. Wolfe, J. A. Neuman, and X. Zhou, Why do models overestimate surface ozone in the Southeast United States, Atmos. Chem. Phys., 16, 13561-13577, doi:10.5194/acp-16-13561-2016, 2016. [http://acmg.seas.harvard.edu/publications/2016/Travis_ACPD_2016.pdf PDF], [http://acmg.seas.harvard.edu/publications/2016/Travis_ACPD_2016_sup.pdf Supplement]
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 18:06, 12 July 2017 (UTC)
| + | |
| − | | + | |
| − | ==== Modifications to the original updates ====
| + | |
| − | | + | |
| − | The following modifications were made to the original updates listed in the [[Media:NewChemistry 030917.pdf|above document]] following conversations with the developers. These modifications were included in [[GEOS-Chem v11-02#v11-02c|v11-02c]].
| + | |
| − | | + | |
| − | (1) Restore H2O2 Henry's law constant for wet deposition. Daniel Jacob wrote:
| + | |
| − | :''For wetdep of H2O2 we should restore the old Henry’s law constant of 8.3E4exp[7400(1/T – 1/298)] because as Dylan points out that’s the physical value. For drydep of H2O2 we can keep the value of 5E7 as parameterized by Nguyen to fit his drydep data.''
| + | |
| − | | + | |
| − | (2) HC187 is advected
| + | |
| − | | + | |
| − | (3) The following species have different names from the original document:
| + | |
| − | :*API is now MTPA (for consistency with [[Secondary_organic_aerosols#Complex_SOA_scheme|existing SOA scheme]])
| + | |
| − | :*APIO2 is now PIO2 (for consistency with [[PAN|PAN updates]] added in [[GEOS-Chem v11-02#v11-02a|v11-02a]])
| + | |
| − | :*LIM is now LIMO (for consistency with [[Secondary_organic_aerosols#Complex_SOA_scheme|existing SOA scheme]])
| + | |
| − | :*PMN is now NPMN and IPMN (PMN from non-isoprene and isoprene sources; from [[Secondary_organic_aerosols#SOA_formation_from_aqueous_isoprene_uptake|aqueous isoprene uptake updates]])
| + | |
| − | :*ONITAam is now IONITA (Jenny Fisher recommended we change the names - they were originally daytime/nighttime species, but changed to isop/monot)
| + | |
| − | :*ONITApm is now MONITA (Jenny Fisher recommended we change the names - they were originally daytime/nighttime species, but changed to isop/monot)
| + | |
| − | | + | |
| − | (4) Fix typos in the original document
| + | |
| − | | + | |
| − | Orig: ISNOHOO + MO2 = 0.660PROPNN + 0.700GLYX + 1.200HO2 + 0.750CH2O + 0.040ISN1OG
| + | |
| − | <span style="color:red">Rate = 2.00e-13</span>
| + | |
| − | v11-02c: ISNOHOO + MO2 = 0.660PROPNN + 0.700GLYX + 1.200HO2 + 0.750CH2O + <span style="color:green">0.250MOH</span> + 0.040ISN1OG
| + | |
| − | <span style="color:green">Rate = 2.06e-13</span>
| + | |
| − |
| + | |
| − | Orig: ISOPNB + OH = <span style="color:red">ISOPNBO2</span> + 0.100IEPOX + 0.100NO2
| + | |
| − | v11-02c: ISOPNB + OH = <span style="color:green">0.900ISOPNBO2</span> + 0.100IEPOX + 0.100NO2
| + | |
| − |
| + | |
| − | Orig: HONIT + OH = NO3 + <span style="color:red">HKET</span>
| + | |
| − | v11-02c: HONIT + OH = NO3 + <span style="color:green">HAC</span>
| + | |
| − |
| + | |
| − | Orig: HONIT + hv = <span style="color:red">HKET</span> + NO2
| + | |
| − | v11-02c: HONIT + hv = <span style="color:green">HAC</span> + NO2
| + | |
| − | | + | |
| − | (5) Completely replace RIP with RIPA, RIPB, RIPD and IEPOX with IEPOXA, IEPOXB, IEPOXD
| + | |
| − | | + | |
| − | Orig: <span style="color:red">RIP + hv = 0.985OH + 0.985HO2 + 0.710CH2O + 0.425MVK + 0.285MACR + 0.275HC5 + 0.005LVOC</span>
| + | |
| − | v11-02c: <span style="color:green">RIPA + hv = 0.985OH + 0.985HO2 + 0.710CH2O + 0.425MVK + 0.285MACR + 0.275HC5 + 0.005LVOC
| + | |
| − | RIPB + hv = 0.985OH + 0.985HO2 + 0.710CH2O + 0.425MVK + 0.285MACR + 0.275HC5 + 0.005LVOC
| + | |
| − | RIPD + hv = 0.985OH + 0.985HO2 + 0.710CH2O + 0.425MVK + 0.285MACR + 0.275HC5 + 0.005LVOC</span>
| + | |
| − |
| + | |
| − | Orig: ISOPND + OH = <span style="color:red">0.100IEPOX</span> + 0.900ISOPNDO2 +0.100NO2
| + | |
| − | v11-02c: ISOPND + OH = <span style="color:green">0.100IEPOXD</span> + 0.900ISOPNDO2 +0.100NO2
| + | |
| − |
| + | |
| − | Orig: ISOPNB + OH = 0.900ISOPNBO2 + <span style="color:red">0.100IEPOX</span> + 0.100NO2
| + | |
| − | v11-02c: ISOPNB + OH = 0.900ISOPNBO2 + <span style="color:green">0.067IEPOXA + 0.033IEPOXB</span> + 0.100NO2
| + | |
| − |
| + | |
| − | Orig: <span style="color:red">IEPOX = SOAIE : HET(ind_IEPOX,1);</span>
| + | |
| − | v11-02c: <span style="color:green">IEPOXA = SOAIE : HET(ind_IEPOXA,1);
| + | |
| − | IEPOXB = SOAIE : HET(ind_IEPOXB,1);
| + | |
| − | IEPOXD = SOAIE : HET(ind_IEPOXD,1);</span>
| + | |
| − | | + | |
| − | (6) Add LVOC to RIP channels
| + | |
| − | | + | |
| − | Orig: RIPA + OH = 0.750 RIO2 + <span style="color:red">0.250 HC5</span> + 0.125 (OH + H2O)
| + | |
| − | v11-02c: RIPA + OH = 0.750 RIO2 + <span style="color:green">0.245 HC5</span> + 0.125 (OH + H2O) + <span style="color:green">0.005 LVOC</span>
| + | |
| − |
| + | |
| − | Orig: RIPA + OH = 0.850 OH + 0.578 IEPOXA + 0.272 IEPOXB + <span style="color:red">0.150 HC5OO</span>
| + | |
| − | v11-02c: RIPA + OH = 0.850 OH + 0.578 IEPOXA + 0.272 IEPOXB + <span style="color:green">0.145 HC5OO + 0.005 LVOC</span>
| + | |
| − |
| + | |
| − | Orig: RIPB + OH = 0.480 RIO2 + <span style="color:red">0.520 HC5</span> + 0.26 (OH + H2O)
| + | |
| − | v11-02c: RIPB + OH = 0.480 RIO2 + <span style="color:green">0.515 HC5</span> + 0.26 (OH + H2O) + <span style="color:green">0.005 LVOC</span>
| + | |
| − |
| + | |
| − | Orig: RIPD + OH = 0.250 RIO2 + <span style="color:red">0.750 HC5</span> + 0.375 (OH + H2O)
| + | |
| − | v11-02c: RIPD + OH = 0.250 RIO2 + <span style="color:green">0.745 HC5</span> + 0.375 (OH + H2O) + <span style="color:green">0.005 LVOC</span>
| + | |
| − |
| + | |
| − | Orig: RIPD + OH = 0.500 OH + 0.500 IEPOXD + <span style="color:red">0.500 HC5OO</span>
| + | |
| − | v11-02c: RIPD + OH = 0.500 OH + 0.500 IEPOXD + <span style="color:green">0.495 HC5OO + 0.005 LVOC</span>
| + | |
| − |
| + | |
| − | The only reaction that wont have LVOC as a product is RIPB + OH = OH + IEPOXA + IEPOXB.
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:26, 7 September 2017 (UTC)
| + | |
| − | | + | |
| − | === Stratospheric chemistry ===
| + | |
| − | | + | |
| − | GEOS-Chem was historically developed as a model of tropospheric chemistry and composition. The above-mentioned [[#Mechanisms in GEOS-Chem v9-01-03 and prior versions|chemistry mechamisms in GEOS-Chem v9-01-03]] and [[#Mechanisms in GEOS-Chem v9-02 and later versions|in GEOS-Chem v9-02]] only solve the chemical reaction matrix within the troposphere. In order to prevent tropospheric species from accumulating in the stratosphere and being transported back into the troposphere, we have implemented the following simple stratospheric chemistry schemes:
| + | |
| − | | + | |
| − | # [[Linoz stratospheric ozone chemistry]]
| + | |
| − | # [[Stratospheric chemistry|Application of monthly-mean prod/loss rates archived from the GMI model]]
| + | |
| − | | + | |
| − | Linoz only applied to ozone. The simple linearized stratospheric chemistry, which uses production and loss rates archived from the GMI model, is applied to all other species. (NOTE: The user has the option to disable Linoz and use the archived GMI prod/loss rates for ozone, but this is typically not done.)
| + | |
| − | | + | |
| − | In [[GEOS-Chem v10-01]] we added the [[UCX_chemistry_mechanism|'''Unified tropospheric-stratospheric Chemistry eXtension''' (UCX) mechanism]] into GEOS-Chem. UCX was developed by Seb Eastham and Steven Barrett at the [http://lae.mit.edu MIT Laboratory for Aviation and the Environment]. This mechanism combines the existing GEOS-Chem [[NOx-Ox-HC-aerosol|"NOx-Ox-HC-aerosol" mechanism]] with several new stratospheric species and reactions.
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Y.]] 12:11, 1 October 2013 (EDT)<br>--[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 17:18, 26 May 2015 (UTC) | + | |
| − | | + | |
| − | === Correcting ozone from the height of the lowest model level to 10m ===
| + | |
| − | | + | |
| − | <span style="color:darkorange">'''''This update is slated for inclusion in [[GEOS-Chem v11-02#v11-02e|GEOS-Chem v11-02e]].'''''</span>
| + | |
| − | | + | |
| − | Katie Travis created a diagnostic to correct daytime ozone values from the lowest model layer, ~60m, to 10m.
| + | |
| − | | + | |
| − | ''C''(''z<sub>C</sub>'') = (1-''R<sub>a</sub>''(''z<sub>1</sub>'',''z<sub>C</sub>'')''v<sub>d</sub>(''z<sub>1</sub>))''C''(''z<sub>1</sub>'') Eq. 1
| + | |
| − | | + | |
| − | where <tt>''R<sub>a</sub>''(''z<sub>1</sub>'',''z<sub>C</sub>'')</tt> is the aerodynamic resistance between <tt>''z<sub>1</sub>''</tt> and <tt>''z<sub>C</sub>''</tt>, and <tt>''v<sub>d</sub>''(''z<sub>1</sub>'')</tt> is the ozone deposition velocity at <tt>''z<sub>1</sub>''</tt>, and <tt>''C''(''z<sub>1</sub>'')</tt> is the ozone concentration at <tt>''z<sub>1</sub>''</tt>.
| + | |
| − | | + | |
| − | <tt>''R<sub>a</sub>''(''z<sub>1</sub>'',''z<sub>C</sub>'')</tt> is calculated to the lowest model level in drydep_mod.F. We recalculate <tt>''R<sub>a</sub>''</tt> using <tt>''z<sub>1</sub>''</tt> = 10 m, which is the height of the CASTNET measurement for ozone. The new <tt>''R<sub>a</sub>''</tt> is added to the diagnostic array AD_RA and passed to diag49.F for use in Equation 1.
| + | |
| − | | + | |
| − | This new diagnostic is called <tt>O3@10m-$</tt>, and can be called with tracer 539 in ND49 in input.geos.
| + | |
| − | | + | |
| − | '''References'''
| + | |
| − | *Travis, K.R., D.J. Jacob, C.A. Keller, S. Kuang, J. Lin, M.J. Newchurch, A.M. Thompson, ''Resolving ozone vertical gradients in air quality models'', <u>Atmos. Chem. Phys. Disc.</u>,2017.
| + | |
| − | *Zhang, L., D.J. Jacob, E.M. Knipping, N. Kumar, J.W. Munger, C.C. Carouge, A. van Donkelaar, Y. Wang, and D. Chen, ''Nitrogen deposition to the United States: distribution, sources, and processes'', <u>Atmos. Chem. Phys.</u>, '''12''', 4,539-4,4554, 2012.
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 22:26, 17 November 2017 (UTC)
| + | |
| − | | + | |
| − | == Analytical tools ==
| + | |
| − | | + | |
| − | === Process analysis diagnostics ===
| + | |
| − | | + | |
| − | [mailto:barronh@ufl.edu Barron Henderson] (U. Florida) has created a [[Process Analysis Diagnostics|software package for process analysis diagnostics]]. He writes:
| + | |
| − | | + | |
| − | <blockquote>Process-based Analysis examines the change in each species due to each process and reaction. Models predict atmospheric state, which in a time-series can be used to create net-change of each species. What this cannot tell us, is which processes led to that change. To supplement state (or concentration), GEOS-Chem has long archived emissions and employed advanced diagnostics to predict gross chemical production or loss. Process Analysis goes a step further archiving grid-cell budgets for each species, and decomposing gross production/loss into individual reaction contributions. Process Analysis extensions are currently available in CAMx, WRF-Chem, CMAQ, and now GEOS-Chem. This allows for direct comparisons of models at a fundamental, process level.</blockquote>
| + | |
| − | | + | |
| − | To obtain this software, [[Process_Analysis_Diagnostics#Installation_and_Application_-_How_can_I_use_it.3F|please contact Barron Henderson directly]].
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Y.]] 12:26, 1 October 2013 (EDT)
| + | |
| − | | + | |
| − | === Linking GEOS-Chem to CMAQ ===
| + | |
| − | | + | |
| − | [mailto:barronh@ufl.edu Barron Henderson] has created Python software that will let you translate GEOS-Chem output to the proper speciation for input to CMAQ. Please see [[Linking_GEOS-Chem_to_CMAQ|our ''Linking GEOS-Chem to CMAQ'' wiki page]] for more information.
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Y.]] ([[User talk:Bmy|talk]]) 16:46, 26 October 2015 (UTC)
| + | |
| − | | + | |
| − | == Previous issues that have now been resolved ==
| + | |
| − | | + | |
| − | === Fixes for carbon creating reactions ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02c|v11-02c]] and approved on 21 Sep 2017.'''''</span>
| + | |
| − | | + | |
| − | '''''Sarah Safieddine wrote:'''''
| + | |
| − | | + | |
| − | <blockquote>Colette, Barron, Mat and myself modified 13 previous "carbon creating" reactions to preserve carbon. The [table below] lists all the corrections for the reactions in globchem.dat V902 that we corrected, with all the details.</blockquote>
| + | |
| − | | + | |
| − | {| border=1 cellspacing=0 cellpadding=5
| + | |
| − | |-bgcolor="#CCCCCC"
| + | |
| − | !Reaction # in globchem.dat v9-02
| + | |
| − | !Unbalanced Reaction
| + | |
| − | !Rate constant
| + | |
| − | !Issue<br>(R=Reactants, P=Products)
| + | |
| − | !Fix and corrected reaction (in <span style="color:green">green</span>)
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |453
| + | |
| − | |R4O2 + NO → NO2 + 0.32ACET + 0.19MEK + 0.18MO2 + 0.27HO2 + 0.32ALD2 + 0.13RCHO + 0.05A3O2 + 0.18B3O2 + 0.32ETO2
| + | |
| − | |K* (1-YN) where YN is returned from fyrno3.f; K=2.7E-12 exp(350/T) (Xcarbn=4.5)
| + | |
| − | |Creates carbon: R=4C P=4.26C
| + | |
| − | |Replace 0.18B3O2 by 0.093B3O2 to achieve carbon closure (as suggested by Barron).<br>
| + | |
| − | <span style="color:green">R4O2 + NO → NO2 + 0.32ACET + 0.19MEK + 0.18MO2 + 0.27HO2 + 0.32ALD2 + 0.13RCHO + 0.05A3O2 + '''0.093B3O2''' + 0.32ETO2</span><br>
| + | |
| − | <span style="color:darkorange">'''Use Barron's [[Chemistry_Issues#Fixes_to_correct_ALK4_lumping_issue|fix to correct ALK4 lumping issue]] instead.'''</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |453
| + | |
| − | |R4N1 + NO → 2NO2 + 0.39CH2O + 0.75ALD2 + 0.57RCHO + 0.3R4O2
| + | |
| − | |2.7E-12 exp(350/T)
| + | |
| − | |Creates carbon: R=4C, P=4.8C
| + | |
| − | |Fix, as suggested by Matt:<br>
| + | |
| − | <span style="color:green">R4N1 + NO → 2NO2 + 0.570RCHO + '''0.86ALD2''' + '''0.57CH2O'''</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |453
| + | |
| − | |ATO2 + NO → 0.96NO2 + 0.96CH2O + 0.96MCO3 + 0.04R4N2
| + | |
| − | |2.8E-12 exp(300/T)
| + | |
| − | |Creates carbon: R=3C, P=3.04
| + | |
| − | |Fix as suggested by Mat: ditch the R4N2 channel<br>
| + | |
| − | <span style="color:green">ATO2 + NO → NO2 + CH2O + MCO3
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |803
| + | |
| − | |RIO2 → 2HO2 + CH2O + 0.5MGLY + 0.5GLYC + 0.5GLYX + 0.5GLYX + 0.5HAC + OH
| + | |
| − | |4.07E+08 exp(-7694/T)
| + | |
| − | |Creates carbon: R=5C, P=7C<br>There was a [[Caltech_isoprene_scheme#Remove_duplicate_GLYX_product_from_RIO2_reaction|fix proposed on the isoprene scheme wiki page]] but still not enough
| + | |
| − | |Fix as suggested by Sarah: remove CH2O<br>
| + | |
| − | <span style="color:green">RIO2 → 2HO2 + 0.5MGLY + 0.5GLYC + 0.5GLYX + 0.5HAC + OH</span><br>
| + | |
| − | <span style="color:darkorange">This reaction was replaced with '''RIO2 → 0.5HPALD + 0.5DHPCARP''' in the [[GEOS-Chem_chemistry_mechanisms#Updated_isoprene_and_monoterpene_chemistry|isoprene chemistry updates]] added in v11-02c</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |453
| + | |
| − | |ISNOOB + NO3 → R4N2 + GLYX + 2NO2
| + | |
| − | |2.3E-12
| + | |
| − | |Creates carbon: R=5C, P=6C
| + | |
| − | |Fix as suggested by Barron: Replace R4N2 by PROPNN<br>
| + | |
| − | <span style="color:green">ISNOOB + NO3 → '''PROPNN''' + GLYX + 2NO2</span><br>
| + | |
| − | <span style="color:darkorange">This reaction was replaced with '''ISNOOB + NO3 → 0.94PROPNN + GLYX + 2NO2 + 0.04ISN1OG''' in the [[GEOS-Chem_chemistry_mechanisms#Updated_isoprene_and_monoterpene_chemistry|isoprene chemistry updates]] added in v11-02c</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |453
| + | |
| − | |ISNOOB+NO → 0.94R4N2 +0.94GLYX +1.88NO2
| + | |
| − | |2.6E-12 exp(380/T)
| + | |
| − | |Creates carbon: R=5C, P=5.64C
| + | |
| − | |Same as above<br>
| + | |
| − | <span style="color:green">ISNOOB + NO → '''0.06R4N2''' + ''''0.94PROPNN'''' + 0.94GLYX + 1.88NO2</span><br>
| + | |
| − | <span style="color:darkorange">This reaction was replaced with '''ISNOOB + NO → 0.9PROPNN + 0.94GLYX + 1.88NO2 + 0.04ISN1OG''' in the [[GEOS-Chem_chemistry_mechanisms#Updated_isoprene_and_monoterpene_chemistry|isoprene chemistry updates]] added in v11-02c</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |453
| + | |
| − | |ISNOHOO + NO → 0.934R4N2 + 0.934HO2 + 0.919GLYX
| + | |
| − | |2.6E-12 exp(380/T)
| + | |
| − | |Creates carbon: R=5C, P=5.574C
| + | |
| − | |Fix by Barron:<br>
| + | |
| − | <span style="color:green">ISNOHOO + NO → '''0.081R4N2''' + '''0.919PROPNN''' + 0.934HO2 + 0.919GLYX</span><br>
| + | |
| − | <span style="color:darkorange">This reaction was replaced with '''ISNOHOO + NO3 → 0.894PROPNN + 0.934HO2 + 0.919GLYX + 0.4ISN1OG''' in the [[GEOS-Chem_chemistry_mechanisms#Updated_isoprene_and_monoterpene_chemistry|isoprene chemistry updates]] added in v11-02c</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |472
| + | |
| − | |MAN2 + HO2 → 0.075PROPNN + 0.075CO + 0.075HO2 + 0.075MGLY + 0.075CH2O + 0.075NO2 + 0.15OH + 0.85ISNP
| + | |
| − | |2.91E-13*exp(1300/T)[1-exp(-0.245*n)],n=4
| + | |
| − | |Creates carbon: R=4C, P=4.85C
| + | |
| − | |Fix by both Mat and Barron: Replace ISNP with 0.85MAOP + 0.85NO2<br>
| + | |
| − | <span style="color:green">MAN2 + HO2 → 0.075PROPNN + 0.075CO + 0.075HO2 + 0.075MGLY + 0.075CH2O + 0.075NO2 + 0.15OH + '''0.85MAOP + 0.85NO2'''</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |719
| + | |
| − | |ATO2 + MCO3 → MEK + ACTA
| + | |
| − | |1.87E-13 exp(500/T)
| + | |
| − | |Creates carbon: R=5C, P=6C
| + | |
| − | |From the WIKI: replace MEK with MGLY<br>
| + | |
| − | <span style="color:green">ATO2 + MCO3 → '''MGLY''' + ACTA</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |817
| + | |
| − | |Br + ALD2 → HBr + MCO3 + CO
| + | |
| − | |1.3E-11 exp(-360/T)
| + | |
| − | |Creates carbon: R=2C, P=3C
| + | |
| − | |Remove CO Following Parrella et al., Table 2a, reactions R7 to R10 (also for the 3 reactions below)<br>
| + | |
| − | <span style="color:green">Br + ALD2 → HBr + MCO3</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |818
| + | |
| − | |Br + ACET → HBr + ATO2 + CO
| + | |
| − | |1.66E-10exp(-7000/T)
| + | |
| − | |Creates carbon: R=3C, P=4C
| + | |
| − | |Remove CO, same as above<br>
| + | |
| − | <span style="color:green">Br + ACET → HBr + ATO2</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |819
| + | |
| − | |Br + C2H6 → HBr + ETO2 + CO
| + | |
| − | |2.36E-10 exp(-6411/T)
| + | |
| − | |Creates carbon: R=2C, P=3C
| + | |
| − | |Remove CO, same as above<br>
| + | |
| − | <span style="color:green">Br + C2H6 → HBr + ETO2</span>
| + | |
| − | | + | |
| − | |-valign="top
| + | |
| − | |820
| + | |
| − | |Br + C3H8 → HBr + A3O2 + CO
| + | |
| − | |8.77E-11 exp(-4330/T)
| + | |
| − | |Creates carbon: R=3C, P=4C
| + | |
| − | |Remove CO, same as above<br>
| + | |
| − | <span style="color:green">Br + C3H8 → HBr + A3O2</span>
| + | |
| − | | + | |
| − | |}
| + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 19:59, 27 July 2017 (UTC)
| + | |
| − | | + | |
| − | === Identification of carbon leaking reactions ===
| + | |
| − | | + | |
| − | '''''Sarah Safieddine wrote:'''''
| + | |
| − | | + | |
| − | <blockquote>76 other reactions leaked carbon, we enforced carbon conservation by tracking the lost carbon as CO2 (labeled as <tt>{CO2}</tt> in the document [http://onlinelibrary.wiley.com/store/10.1002/2017GL072602/asset/supinfo/grl55781-sup-0001-Supplementary.docx?v=1&s=21c96c26c411290d72b64fc4f9ecdc5e806af2c5 ROC_SI.docx], Table 2). This is the supplementary material for Safieddine, Heald and Henderson, 2017. It contains the corrections for both the carbon leaking and carbon creating reactions and all other information. The paper for reference can be found here: http://onlinelibrary.wiley.com/doi/10.1002/2017GL072602/abstract.</blockquote>
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 21:06, 22 May 2017 (UTC)
| + | |
| − | | + | |
| − | === Fixes to correct ALK4 lumping issue ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''These fixes were included in [[GEOS-Chem v11-02#v11-02a|v11-02a]] and approved on 12 May 2017.'''''</span>
| + | |
| − | | + | |
| − | '''''[[User:Barronh|Barron Henderson]] wrote:'''''
| + | |
| − | | + | |
| − | <blockquote>I have a lumping-related issue that I know some of you are already aware of. There is a chemical carbon source (and secondary ETO2 source).
| + | |
| − | | + | |
| − | Right now, ALK4 (via R4O2) produces 4.26 moles carbon products per reaction. The ALK4 representation can be traced back to a paper by Frederick Lurmann. That paper refers to a report that I have been unable to obtain. In fact, Frederick Lurmann no longer has a copy. When we spoke, however, he confirmed my suspicion that ALK4 is based on a 70% butane and 30% pentane mixture. Our 4.26 carbon product appears to be based on two differences (typos?) from the paper that alter the yields.
| + | |
| − | | + | |
| − | If ALK4 emissions are introduced using a 4C assumption, then ALK4 chemistry is acting as a 7% carbon source. From a ozone reactivity standpoint, this is not a major issue. First, the speciation of VOC is highly uncertain and most of the atmosphere is NOx-limited. Even so, it represents another reason to revisit our lumped species.
| + | |
| − | | + | |
| − | I have [https://www.evernote.com/shard/s315/sh/f2ec9589-d827-4ee1-afcb-96ee5a2d2914/d84318450f729cd414e4a6653c03a296 extensive notes] on what I interpret as happening. To the best of my knowledge, we need to make three modifications to R4O2 + NO. The first two are to make R4O2 correctly linked to Lurmann and the third is to correctly connect the mass emissions with the molar conservation.
| + | |
| − | | + | |
| − | #Increase MO2 stoichiometry from 0.18 to 0.19
| + | |
| − | #Increase RCHO stoichiometry from 0.13 to 0.14 (or A3O2 from 0.05 to 0.06 -- it is not clear to me when this was introduced).
| + | |
| − | #Modify the carbon count for ALK4 (i.e. the <tt>MolecRatio</tt> field in the [[GEOS-Chem species database]]) from 4 to 4.3.
| + | |
| − | | + | |
| − | Fixes 1 and 2—which can be applied to the KPP <tt>globchem.eqn</tt> file—will make the carbon conservation consistent with Lurmann's. Right now, it looks like there were a couple changes that could have been inadvertent (i.e., 0.18 instead of 0.19). If there was a reason for these changes, I have been unable to find it.</blockquote>
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:21, 31 January 2017 (UTC)
| + | |
| − | | + | |
| − | === JPL Released 18th Rate Coefficient Evaluation ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02a|v11-02a]] and approved on 12 May 2017.'''''</span>
| + | |
| − | | + | |
| − | JPL has released its 18th evaluation of chemical rate coefficients for atmospheric studies (Burkholder et al., 2015)." A new page ([[Updates in JPL Publication 15-10]]) is being created to compare rates between GEOS-Chem v10 and JPL Publication 15-10.
| + | |
| − | | + | |
| − | :J. B. Burkholder, S. P. Sander, J. Abbatt, J. R. Barker, R. E. Huie, C. E. Kolb, M. J. Kurylo, V. L. Orkin, D. M. Wilmouth, and P. H. Wine "Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation No. 18," JPL Publication 15-10, Jet Propulsion Laboratory, Pasadena, 2015 http://jpldataeval.jpl.nasa.gov.
| + | |
| − | | + | |
| − | ---[[User:barronh|B. Henderson]] 2016-05-03 15:25 (EDT)
| + | |
| − | | + | |
| − | | + | |
| − | === Centralizing chemistry time step===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-02_benchmark_history#v9-01-02q|v9-01-02q]] and approved on 18 Oct 2011.'''''</span>
| + | |
| − | | + | |
| − | Please see the full discussion on the [[Centralized chemistry time step]] wiki page.
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Y.]] 16:01, 4 November 2011 (EDT)
| + | |
| − | | + | |
| − | === Acetone photolysis ===
| + | |
| − | | + | |
| − | [[FAST-J_photolysis_mechanism#v9-02_post-release_patch_to_fix_bug_in_acetone_photolysis_pressure_dependency|This discussion has been moved to our ''FAST-J photolysis mechanism'' wiki page]].
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Y.]] 15:20, 20 May 2014 (EDT)
| + | |
| − | | + | |
| − | == Issues that have been since rendered obsolete by newer code updates ==
| + | |
| − | | + | |
| − | Most of the issues described below pertained to the SMVGEAR chemical solver (which was replaced by FlexChem in [[GEOS-Chem v11-01|v11-01]]) and/or the FAST-J photolysis mechanism (which was replaced by FAST-JX in [[GEOS-Chem v10-01|v10-01]]).
| + | |
| − | | + | |
| − | === NIT should be converted to molec/cm3 in calcrate.F ===
| + | |
| − | | + | |
| − | [[Image:Obsolete.jpg]]
| + | |
| − | | + | |
| − | <span style="color:red">'''''SMVGEAR was removed from [[GEOS-Chem v11-01]] and higher versions. The code in <tt>calcrate.F</tt> will be replaced by the KPP master equation file.'''''</span>
| + | |
| − | | + | |
| − | In <tt>calcrate.F</tt>, we have:
| + | |
| − | | + | |
| − | ! Nitrate effect; reduce the gamma on nitrate by a
| + | |
| − | ! factor of 10 (lzh, 10/25/2011)
| + | |
| − | IF ( N == 8 ) THEN
| + | |
| − | TMP1 = State_Chm%Tracers(IX,IY,IZ,IDTSO4) +
| + | |
| − | & State_Chm%Tracers(IX,IY,IZ,IDTNIT)
| + | |
| − | TMP2 = State_Chm%tracers(IX,IY,IZ,IDTNIT)
| + | |
| − | IF ( TMP1 .GT. 0.0 ) THEN
| + | |
| − | XSTKCF = XSTKCF * ( 1.0e+0_fp - 0.9e+0_fp
| + | |
| − | & *TMP2/TMP1 )
| + | |
| − | ENDIF
| + | |
| − | ENDIF
| + | |
| − | | + | |
| − | Here NIT is added to SO4 but NIT is in different units than SO4. This unit difference can be traced to the definition of IDTRMB, which is only nonzero for species that are in the SMVGEAR mechanism. Since NIT is not a SMVGEAR species, IDTRMB = 0 for NIT and it is therefore skipped in the unit conversion from kg --> molec/cm3 in <tt>partition.F</tt>.
| + | |
| − | | + | |
| − | This issue was discovered during the implementation of [[FlexChem]]. In [[GEOS-Chem v11-01#v11-01g|GEOS-Chem v11-01g]] and later versions, units of NIT are properly accounted for in routine <tt>HETN2O5</tt> (found in <tt>gckpp_HetRates.F90</tt>).
| + | |
| − | | + | |
| − | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 20:25, 12 September 2016 (UTC)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:27, 31 January 2017 (UTC)
| + | |
| − | | + | |
| − | === rate of HNO4 ===
| + | |
| − | | + | |
| − | [[Image:Obsolete.jpg]]
| + | |
| − | | + | |
| − | <span style="color:red">'''''SMVGEAR was removed from [[GEOS-Chem v11-01]] and higher versions. The <tt>globchem.dat</tt> file is now replaced by the KPP master equation file.'''''</span>
| + | |
| − | | + | |
| − | [mailto:ecbrow@berkeley.edu Ellie Browne] found a typo in the globchem.dat ([[GEOS-Chem v8-02-01]] and beyond)
| + | |
| − | <pre>
| + | |
| − | A 73 9.52E-05 3.2E+00 -10900 1 P 0.60 0. 0.
| + | |
| − | 1.38E+15 1.4E+00 -10900 0 0.00 0. 0.
| + | |
| − | HNO4 + M
| + | |
| − | =1.000HO2 +1.000NO2 + +
| + | |
| − | </pre>
| + | |
| − | This should be corrected as:
| + | |
| − | <pre>
| + | |
| − | A 73 9.52E-05 3.4E+00 -10900 1 P 0.60 0. 0.
| + | |
| − | 1.38E+15 1.1E+00 -10900 0 0.00 0. 0.
| + | |
| − | HNO4 + M
| + | |
| − | =1.000HO2 +1.000NO2 + +
| + | |
| − | </pre>
| + | |
| − | The difference is within 2%.
| + | |
| − | | + | |
| − | --[[User:Jmao|J Mao.]] 19:04, 30 Aug 2010 (EDT)<br>
| + | |
| − | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:29, 31 January 2017 (UTC)
| + | |
| − | | + | |
| − | === near-IR photolysis of HNO4 ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was added to [[GEOS-Chem v8-02-04]].'''''</span>
| + | |
| − | | + | |
| − | [[Image:Obsolete.jpg]]
| + | |
| − | | + | |
| − | <span style="color:red">'''''SMVGEAR was removed from [[GEOS-Chem v11-01]] and higher versions. The <tt>globchem.dat</tt> file is now replaced by the KPP master equation file. Also, FAST-JX has now replaced FAST-J photolysis.'''''</span>
| + | |
| − | | + | |
| − | 1. Since FastJX already takes this into account with cross section data at 574nm, we do not need to redo this in <tt>calcrate.f</tt>. We can therefore comment out this entire IF block:
| + | |
| − |
| + | |
| − | !---------------------------------------------------------------------
| + | |
| − | ! Prior to 10/27/09:
| + | |
| − | ! FastJX has taken near-IR photolysis into account with
| + | |
| − | ! cross section at 574nm, so we don't need to add 1e-5 anymore.
| + | |
| − | ! According to Jimenez et al., "Quantum yields of OH, HO2 and
| + | |
| − | ! NO3 in the UV photolysis of HO2NO2", PCCP, 2005, we also
| + | |
| − | ! changed the branch ratio from 0.67(HO2)/0.33(OH) to 0.95/0.05
| + | |
| − | ! This will put most weight of near-IR photolysis on HO2 channel.
| + | |
| − | ! (jmao, bmy, 10/27/09)
| + | |
| − | !
| + | |
| − | !!==============================================================
| + | |
| − | !! HARDWIRE addition of 1e-5 s-1 photolysis rate to
| + | |
| − | !! HNO4 -> HO2+NO2 to account for HNO4 photolysis in near-IR --
| + | |
| − | !! see Roehl et al. 'Photodissociation of peroxynitric acid in
| + | |
| − | !! the near-IR', 2002. (amf, bmy, 1/7/02)
| + | |
| − | !!
| + | |
| − | !! Add NCS index to NKHNO4 for SMVGEAR II (gcc, bmy, 4/1/03)
| + | |
| − | !!==============================================================
| + | |
| − | !IF ( NKHNO4(NCS) > 0 ) THEN
| + | |
| − | !
| + | |
| − | ! ! Put J(HNO4) in correct spot for SMVGEAR II
| + | |
| − | ! PHOTVAL = NKHNO4(NCS) - NRATES(NCS)
| + | |
| − | ! NKN = NKNPHOTRT(PHOTVAL,NCS)
| + | |
| − | !
| + | |
| − | ! DO KLOOP=1,KTLOOP
| + | |
| − | ! RRATE(KLOOP,NKN)=RRATE(KLOOP,NKN) + 1d-5
| + | |
| − | ! ENDDO
| + | |
| − | !ENDIF
| + | |
| − | !---------------------------------------------------------------------
| + | |
| − | | + | |
| − | | + | |
| − | 2. We need to change the branch ratio of HNO4 photolysis in <tt>ratj.d</tt>. Change these lines from:
| + | |
| − | | + | |
| − | 13 HNO4 PHOTON OH NO3 0.00E+00 0.00 33.3 HO2NO2
| + | |
| − | 14 HNO4 PHOTON HO2 NO2 0.00E+00 0.00 66.7 HO2NO2
| + | |
| − | | + | |
| − | to:
| + | |
| − | | + | |
| − | 13 HNO4 PHOTON OH NO3 0.00E+00 0.00 5.0 HO2NO2
| + | |
| − | 14 HNO4 PHOTON HO2 NO2 0.00E+00 0.00 95.0 HO2NO2
| + | |
| − | | + | |
| − | This is based on Jimenez et al. (Quantum yields of OH, HO2 and NO3 in the UV photolysis of HO2NO2, PHYSICAL CHEMISTRY CHEMICAL PHYSICS, 2005) shows that HO2 yield should be 0.95 and OH yield should be 0.05 for wavelength above 290nm.
| + | |
| − | | + | |
| − | This way all the near-IR photolysis will have most weight on HO2 channel(Stark et al., Overtone dissociation of peroxynitric acid (HO2NO2): Absorption cross sections and photolysis products, JOURNAL OF PHYSICAL CHEMISTRY A, 2008).
| + | |
| − | | + | |
| − | This update has now been added to the [http://acmg.seas.harvard.edu/geos/wiki_docs/chemistry/chemistry_updates_v6.pdf chemistry mechanism documentation file].
| + | |
| − | | + | |
| − | --[[User:Jmao|J Mao.]] 11:00, 26 Oct 2009 (EDT)<br>
| + | |
| − | --[[User:Bmy|Bob Y.]] 16:08, 4 November 2011 (EDT)
| + | |
| − | | + | |
| − | === yield of isoprene nitrates ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was added to [[GEOS-Chem v8-03-02]] as a post-release patch, and standardized in [[GEOS-Chem v9-01-01]].'''''</span>
| + | |
| − | | + | |
| − | [[Image:Obsolete.jpg]]
| + | |
| − | | + | |
| − | <span style="color:red">'''''SMVGEAR was removed from [[GEOS-Chem v11-01]] and higher versions. The <tt>globchem.dat</tt> file is now replaced by the KPP master equation file.'''''</span>
| + | |
| − | | + | |
| − | [mailto:paulot@caltech.edu Fabien Paulot] found a problem in current chemistry scheme. In [[GEOS-Chem v8-02-01]] and beyond, isoprene nitrates are produced twice: one through channel A and one through 10% loss in channel B. This makes the loss of NOx larger than it should be (18.7% vs. 10%) and also reduces the yield of MVK/MACR/CH2O by about 13%.
| + | |
| − | | + | |
| − | A 453 2.70E-12 0.0E+00 350 1 B 0.00 0. 0.
| + | |
| − | 5.00E+00 0.0E+00 0 0 0.00 0. 0.
| + | |
| − | RIO2 + NO
| + | |
| − | =0.900NO2 +0.900HO2 +0.340IALD +0.340MVK
| + | |
| − | +0.220MACR +0.560CH2O + +
| + | |
| − |
| + | |
| − | A 453 2.70E-12 0.0E+00 350 1 A 0.00 0. 0.
| + | |
| − | 5.00E+00 0.0E+00 0 0 0.00 0. 0.
| + | |
| − | RIO2 + NO
| + | |
| − | =1.000HNO3 + + +
| + | |
| − | | + | |
| − | So it should be corrected as (no channel A):
| + | |
| − | | + | |
| − | A 453 2.70E-12 0.0E+00 350 0 0 0.00 0. 0.
| + | |
| − | RIO2 + NO
| + | |
| − | =0.900NO2 +0.900HO2 +0.340IALD +0.340MVK
| + | |
| − | +0.220MACR +0.560CH2O + +
| + | |
| − |
| + | |
| − | D 453 2.70E-12 0.0E+00 350 1 A 0.00 0. 0.
| + | |
| − | 5.00E+00 0.0E+00 0 0 0.00 0. 0.
| + | |
| − | RIO2 + NO
| + | |
| − | =1.000HNO3 + + +
| + | |
| − | | + | |
| − | --[[User:Jmao|J Mao.]] 18:04, 30 Aug 2010 (EDT)<br>
| + | |
| − | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:31, 31 January 2017 (UTC)
| + | |
| − | | + | |
| − | === Potential issue with reading restart.cspec file ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-02_benchmark_history#v9-01-02c|v9-01-02c]] and approved on 21 Jul 2011.'''''</span>
| + | |
| − | | + | |
| − | [[Image:Obsolete.jpg]]
| + | |
| − | | + | |
| − | <span style="color:red">'''''The binary-punch format <tt>restart.cspec.YYYYMMDDhh</tt> file is slated to be replaced by a netCDF-format restart file, starting in [[GEOS-Chem v11-01]] and higher versions. But during a transition period, you can still request binary-punch format output.'''''</span>
| + | |
| − | | + | |
| − | Jingqiu Mao discovered a mis-indexing problem when using the <tt>restart.cspec.YYYYMMDDhh</tt> file. Please see [[Restart files#Potential issue with reading restart.cspec_file|this wiki post]] for more information.
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Y.]] 16:02, 4 November 2011 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:33, 31 January 2017 (UTC)
| + | |
| − | | + | |
| − | === GLCO3, GLPAN bug in standard mechanism ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-03_benchmark_history#v9-01-03a|v9-01-03a]] and approved on 08 Dec 2011.'''''</span>
| + | |
| − | | + | |
| − | [[Image:Obsolete.jpg]]
| + | |
| − | | + | |
| − | <span style="color:red">'''''SMVGEAR was removed from [[GEOS-Chem v11-01]] and higher versions. The <tt>globchem.dat</tt> file is now replaced by the KPP master equation file.'''''</span>
| + | |
| − | | + | |
| − | '''''[mailto:fabienpaulot@gmail.com Fabien Paulot] wrote:'''''
| + | |
| − | | + | |
| − | :I think there is a relatively serious bug in the standard chemistry. GLPAN and GLCO3 are set to inactive but their production and loss reactions are active. As a result they never reach equilibrium and this results in an artificial loss of NOx.
| + | |
| − | | + | |
| − | :If this is the only cause of the imbalance between sources and sinks of NOx in my simulations, this would account for ~5% of NOy losses. I don't see that problem in a simulation with a different chemistry that among other changes does not feature those reactions. So hopefully that's it.
| + | |
| − | | + | |
| − | :To fix the error, I made the following modifications in <tt>globchem.dat</tt>:
| + | |
| − | | + | |
| − | :#I set GLPAN and GLCO3 rxns from active to dead. These rxns were causing an artificial loss of NOx.
| + | |
| − | :#I have physically removed GLCO3, GLP, GLPAN, GPAN, ISNO3, MNO3, O2CH2OH, MVN2 and their associated reactions.
| + | |
| − | :#I have made GLYX active. I'm not sure why it's not active by default.
| + | |
| − | | + | |
| − | :and to <tt>ratj.d</tt>:
| + | |
| − | | + | |
| − | :# I deleted photolysis reactions for MNO3 and GLP, since these species have also now been deleted in <tt>globchem.dat</tt>
| + | |
| − | | + | |
| − | --[[User:Bmy|Bob Y.]] 14:51, 10 November 2011 (EST)<br>
| + | |
| − | --[[User:Melissa Payer|Melissa Payer]] 10:49, 15 December 2011 (EST)<br>
| + | |
| − | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:35, 31 January 2017 (UTC)
| + | |
| − | | + | |
| − | === Bug in routine ARSL1K ===
| + | |
| − | | + | |
| − | <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem v9-01-03 benchmark history#v9-01-03m|v9-01-03m]] and approved on 06 Jun 2012.'''''</span>
| + | |
| − | | + | |
| − | [[Image:Obsolete.jpg]]
| + | |
| − | | + | |
| − | <span style="color:red">'''''SMVGEAR was removed from [[GEOS-Chem v11-01]] and higher versions. The <tt>ARSL1K</tt> routine was replaced by an equivalent function in KPP's rate law library.'''''</span>
| + | |
| − | | + | |
| − | A bug in routine ARSL1K became problematic in the implementation of Justin Parrella's [[Bromine_chemistry_mechanism|tropospheric bromine chemistry mechanism]] for [[GEOS-Chem v9-01-03]]. In the bromine chemistry mechanism, a sticking coefficient of 0.0 is passed to the routine ARSL1K for non-sulfate, non-sea salt aerosol. The IF statement modified in [[GEOS-Chem_v8-02-04#Div-by-zero_error_encountered_in_arsl1k.f|GEOS-Chem v8-02-04]] resulted in the reaction rate being set to the default value of 1.0d-3. A 1-month benchmark for July 2005 indicated that the simulated BrO was a little more than twice the expected zonal mean. Modifying the default value from 1.0d-3 to 1.0d-30 resulted in reasonable simulated BrO values.
| + | |
| − | | + | |
| − | '''''[mailto:mat.evans@york.ac.uk Mat Evans] wrote:'''''
| + | |
| − | | + | |
| − | :I've re-run two 2 month simulation [using [[GEOS-Chem v9-01-02]]]. One with the error handling value of 1e-3 (standard) and one with it being 1e-30. There are 5127 time and space points where the model traps the problem and invokes the 1e-3 or 1e-30 value. There are 30*24*2*37*72*46 (roughly 200 million) time and space points when the error could have occurred so we are looking at a relatively infrequent event.
| + | |
| − | | + | |
| − | :The simulations show virtually no difference between the two simulations.
| + | |
| − | | + | |
| − | :mean and stddev ratio of all grid boxes with and without the fix are shown below
| + | |
| − | NOx 0.999996 0.000409291
| + | |
| − | Ox 1.00000 1.27233e-05
| + | |
| − | O3 1.00000 1.52284e-05
| + | |
| − | PAN 0.997849 0.0111997
| + | |
| − | CO 1.00000 4.21768e-06
| + | |
| − | ALK4 0.990514 0.0351941
| + | |
| − | ISOP 0.999979 0.0108033
| + | |
| − | H2O2 0.992067 0.0264659
| + | |
| − | DST1 1.00000 0.00000
| + | |
| − | HO2 0.999996 0.00309464
| + | |
| − | OH 1.00003 0.00767954
| + | |
| − | | + | |
| − | :So although there are some differences they are very minor. For completeness we should put this in as a bug fix (make the error value 1d-30 rather than 1d-3). But it is not a major problem.
| + | |
| | | | |
| − | --[[User:Melissa Payer|Melissa Payer]] 17:52, 14 May 2012 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:35, 31 January 2017 (UTC) | + | ---- |
| | + | '''''[[Guide to GEOS-Chem simulations|Previous]] | [[Aerosol-only simulation|Next]] | [[Guide to GEOS-Chem simulations]]''''' |
On this page, we provide information about GEOS-Chem simulations that use chemistry mechanism solver code built by the Kinetic PreProcessor (KPP).
The following table provides links to information about the available full-chemistry mechanisms in GEOS-Chem.
