Difference between revisions of "Mean OH concentration"
(→Mean OH in GEOS-Chem v11) |
(→Mean OH in GEOS-Chem v11) |
||
Line 63: | Line 63: | ||
|- bgcolor="#CCCCCC" | |- bgcolor="#CCCCCC" | ||
!width="125px"|Benchmark | !width="125px"|Benchmark | ||
− | !width=" | + | !width="80px"|Met Fields & Year |
!width="125px"|Mean OH<br>[10<sup>5</sup> molec cm<sup>-3</sup>] | !width="125px"|Mean OH<br>[10<sup>5</sup> molec cm<sup>-3</sup>] | ||
!width="125px"|Change from<br>prev. version<br>[10<sup>5</sup> molec cm<sup>-3</sup>] | !width="125px"|Change from<br>prev. version<br>[10<sup>5</sup> molec cm<sup>-3</sup>] | ||
Line 70: | Line 70: | ||
|-valign="top" | |-valign="top" | ||
− | |[[GEOS- | + | |[[GEOS-Chem v11-01 benchmark history#v11-01k-Run0|v11-01k-Run0]] |
|GEOS-FP 2013 | |GEOS-FP 2013 | ||
− | |<tt></tt> | + | |<tt>12.004</tt> |
− | |<tt> (%)</tt> | + | |<tt>-0.065 (-.538%)</tt> |
+ | |'''Updates introduced in [[GEOS-Chem v11-01#v11-01h|v11-01h]]:''' | ||
+ | *[[GEOS-Chem_v11-01#Incorrect_incorporation_of_moisture_in_GEOS-Chem_air_quantities|Fix to remove moisture signature in dry mixing ratio]] | ||
+ | '''Updates introduced in [[GEOS-Chem v11-01#v11-01i|v11-01i]]:''' | ||
+ | *[[Wet_deposition#Fix_bug_in_GEOS-FP.2FMERRA.2FMERRA-2_re-evaporation_calculation|Fix bug in GEOS-FP/MERRA/MERRA-2 re-evaporation calculation]] | ||
+ | '''Updates introduced in [[GEOS-Chem v11-01#v11-01j|v11-01j]]:''' | ||
+ | *[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-01j|Update to HEMCO v2.0.003]], includes: | ||
+ | **Bug fix for PARANOX extension | ||
+ | **Fix for ocean overlap in timezones mask file | ||
+ | *[[CAC_anthropogenic_emissions#High_resolution_emissions|0.1° x 0.1° CAC emissions]] | ||
+ | *[[Implementation_of_RETRO_Anthropogenic_Emissions#RETRO_anthropogenic_propane_emissions_are_too_low|Replace RETRO C3H8 emissions with emissions from Xiao 2008]], includes: | ||
+ | **Overwrite Xiao C2H6 and C3H8 with MIX and NEI2011 regional inventories | ||
+ | *[[EDGAR_v4.2_anthropogenic_emissions#Correct_bug_in_units_of_EDGAR_v4.2_SO2_emissions|Correct bug in units of EDGAR v4.2 SO2 emissions]] | ||
+ | *[[EDGAR_v4.2_anthropogenic_emissions#Prevent_biofuel_emissions_from_being_double_counted|Prevent biofuel emissions from being double counted in EDGAR v4.2]] | ||
+ | *[[New_isoprene_scheme#Set_ACTA_and_HCOOH_to_active_in_globchem.dat|Set ACTA and HCOOH to active]] | ||
+ | *[[UCX_chemistry_mechanism#Fix_for_mesospheric_NOx_chemistry|Fix for mesospheric NOx chemistry]] | ||
+ | *[[GEOS-Chem v11-01#Fix_molecular_weight_of_N2O5|Fix molecular weight of N2O5]] | ||
+ | '''Updates introduced in [[GEOS-Chem v11-01#v11-01k|v11-01k]]:''' | ||
+ | *[[Lightning_NOx_emissions#Update_WEST_NS_DIV_to_fix_ozone_overestimate_in_southern_US|Lightning update based on Travis et al. (2016, ACP) and Zhang et al. (2014, ACP)]] | ||
+ | | | ||
+ | *[[GEOS-Chem v11-01]] | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem v11-01 benchmark history#v11-01g-Run0|v11-01g-Run0]] | ||
+ | |GEOS-FP 2013 | ||
+ | |<tt>12.069</tt> | ||
+ | |<tt>0.025 (0.207%)</tt> | ||
+ | |'''Updates introduced in [[GEOS-Chem v11-01#v11-01g|v11-01g]]:''' | ||
+ | *[[FlexChem]] | ||
+ | | | ||
+ | *[[GEOS-Chem v11-01]] | ||
+ | |||
+ | |-valign="top" | ||
+ | |[[GEOS-Chem_v11-01_benchmark_history#v11-01f-geosfp-Run0|v11-01f-geosfp-Run0]] | ||
+ | |GEOS-FP 2013 | ||
+ | |<tt>12.044</tt> | ||
+ | |<tt>0.033 (0.275%)</tt> | ||
+ | ||'''Updates introduced in [[GEOS-Chem v11-01#v11-01e|v11-01e]]:''' | ||
+ | *[[Implementation_of_HEMCO_in_GEOS-Chem#Features_added_in_v11-01e|Update HEMCO from version 1.1.005 to version 1.1.016]] | ||
+ | '''Updates introduced in [[GEOS-Chem v11-01#v11-01f|v11-01f]]:''' | ||
+ | *[[EMEP_European_anthropogenic_emissions#Zero_EMEP_emissions_at_30N|New EMEP mask to prevent zero emissions at 30N]] | ||
+ | *[[Centralized_chemistry_time_step#Optimal_Configuration|Optimal timestep settings]] | ||
+ | *Bug fixes in globchem.dat for: | ||
+ | **[[New_isoprene_scheme#Syntax_error_in_globchem.dat_file_for_MOBA_reaction_rate|MOBA + OH --> MOBAOO]] | ||
+ | **[[New_isoprene_scheme#Bug_in_reaction_ISNOOA_.2B_NO2_in_globchem.dat|ISNOOA + NO2 --> PMN]] | ||
| | | | ||
*[[GEOS-Chem v11-01]] | *[[GEOS-Chem v11-01]] |
Revision as of 17:28, 29 March 2018
On this page we track the evolution of the mean OH concentration in GEOS-Chem. This value is an indicator of how "hot" the chemistry is. Please also see the related topic about the methyl chloroform lifetime in GEOS-Chem.
Computation of the mean OH diagnostic
The mean OH diagnostic in GEOS-Chem is computed in diag_pl_mod.f.
Routine DO_DIAG_OH sums up the tropospheric air mass and mass-weighted OH mass after each call to the chemistry solver. The summing happens in these lines of code:
! Sum air mass term into AIR_MASS array XAIRMASS = AIRDENS(JLOOP) * VOLUME(JLOOP) AIR_MASS(I,J,L) = AIR_MASS(I,J,L) + XAIRMASS ! Sum OH mass term into OH_MASS array XOHMASS = CSPEC(JLOOP,IDOH) * XAIRMASS OH_MASS(I,J,L) = OH_MASS(I,J,L) + XOHMASS
where the arrays contain the following quantities:
JLOOP = 1-D array index for SMVGEAR I,J,L = 3-D array indices AIRDENS(JLOOP) = air density [molec/cm3] VOLUME(JLOOP) = grid box volume [cm3] AIR_MASS(I,J,L) = air mass [molec air] CSPEC(JLOOP,IDOH) = OH concentration [molec OH/cm3] OH_MASS(I,J,L) = mass-weighted OH = OH mass * air mass [molec OH/cm3 * molec air]
After the GEOS-Chem simulation finishes, the mean OH value is printed out in subroutines PRINT_MEAN_OH:
! Total Mass-weighted OH [molec OH/cm3] * [molec air] SUM_OHMASS = SUM( OH_MASS ) ! Atmospheric air mass [molec air] SUM_MASS = SUM( AIR_MASS ) ! Divide OH [molec OH/cm3 * molec air] by [molec air] and report as [1e5 molec/cm3] OHCONC = ( SUM_OHMASS / SUM_MASS ) / 1d5 ! Write value to log file WRITE( 6, '(/,a)' ) REPEAT( '=', 79 ) WRITE( 6, * ) 'ND23: Mass-Weighted OH Concentration' WRITE( 6, * ) 'Mean OH = ', OHCONC, ' [1e5 molec/cm3]' WRITE( 6, '( a)' ) REPEAT( '=', 79 )
Typical values vary between 105 - 115 molec/cm3.
--Bob Y. 15:12, 11 February 2010 (EST)
Evolution of mean OH in GEOS-Chem
Plot of mean OH evolution
Here is a plot of the mean OH concentration from several recent GEOS–Chem 1-year benchmark simulations. The mean OH is an indicator how the chemistry mechanism is performing. Typical values are 10–12 x 105 molecules cm-1.
--Bob Y. 12:41, 17 January 2014 (EST)
Mean OH in GEOS-Chem v11
The following table chronicles the evolution of the mean OH value from several different GEOS-Chem simulations. All benchmark simulations were conducted at 4 x 5 resolution unless noted otherwise.
Mean OH in GEOS-Chem v10
Benchmark | Met Fields & Year | Mean OH [105 molec cm-3] |
Change from previous version [105 molec cm-3] |
Important updates in this benchmark simulation | References |
---|---|---|---|---|---|
v10-01i-Run0 | GEOS-FP 2013 | 11.125 | 0.033 (0.298%) |
| |
v10-01h-Run0 | GEOS-FP 2013 | 11.092 | -0.751 (-6.341%) |
|
|
v10-01e-Run1 | GEOS-FP 2013 | 11.843 | -0.042 (-0.353%) |
| |
v10-01c-Run1 | GEOS-FP 2013 | 11.885 | -0.665 (-5.299%) |
| |
v10-01c-Run0 | GEOS-FP 2013 | 12.550 | -0.139 (-1.095%) |
|
Mean OH in GEOS-Chem v9
Mean OH in GEOS-Chem v8
Benchmark | Met Fields & Year | Mean OH [105 molec cm-3] |
Change from previous version [105 molec cm-3] |
Important updates in this benchmark simulation | References | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
GEOS-Chem version 8 | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-02-04-Run2 | GEOS-5.1.0 2005 | 11.067 | +0.019 (+0.17%) |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-02-04-Run1 | GEOS-5.1.0 2005 | 11.048 | +0.598 (+5.72%) |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-02-04-Run0 | GEOS-5.1.0 2005 | 10.450 | -1.362 (-11.53%) |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-02-01-Run0 | GEOS-5.1.0 2005 | 11.812 | +0.713 (+6.42%) |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-01-04-Run2 | GEOS-5.1.0 2005 | 11.099 | +0.805 (+7.82%) |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-01-04-Run1 | GEOS-5.0.1 2005 | 10.294 | -0.395 (-3.69%) |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-01-04-Run0 | GEOS-5.0.1 2005 | 10.689 | -1.351 (-11.22%) |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-01-01-Run1 | GEOS-5.0.1 2005 | 12.040 | +0.034 (+0.28%) |
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||
v8-01-01-Run0 | GEOS-5.0.1 2005 | 12.006 | -0.141 (-1.16%) |
|
Mean OH in GEOS-Chem v7
--Bob Y. 12:13, 3 October 2013 (EDT) Additional documentationGEOS-Chem version 8
GEOS-Chem version 7
--Bob Y. 09:32, 23 February 2010 (EST) |