Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

[PULL REQUEST] Revise N2O5 gamma calculation #913

Merged
merged 1 commit into from
Sep 27, 2021
Merged

Conversation

yantosca
Copy link
Contributor

This is the corresponding PR for Feature Request #907.

While this is a "one-liner" fix, we will run "alpha" benchmarks on this update, as it has the potential to significantly change the model results.

This commit addresses the issue in geoschem/geos-chem #907.

In the computation of the gamma (rxn probability) for N2O5, the
alpha factor (here called A) in routine N2O5_Inorg_Org can get
too large.  Viral Shah writes:

   The issue is that A scales as Volume/Surface area and becomes large
   for coarse-mode sea salt aerosols. The parameterization is based on
   the assumption that the N2O5 reaction is volume-limited, but Gaston &
   Thornton (2016) show that the reaction is largely diffusion-limited
   for particles larger than ~100 nm radius.

   I discussed this with Chris Holmes @cdholmes, and we propose to fix
   this by capping A to 3.2e-8, which is the value reported in
   Bertram & Thornton (2009). This sets the maximum gamma value
   to ~0.036, and would be consistent with the lab results.

Signed-off-by: Bob Yantosca <yantosca@seas.harvard.edu>
@yantosca yantosca added category: Feature Request New feature or request topic: Chemical Mechanisms Related to KPP and/or GEOS-Chem chemistry mechanisms labels Sep 24, 2021
@yantosca yantosca self-assigned this Sep 24, 2021
@yantosca yantosca linked an issue Sep 24, 2021 that may be closed by this pull request
@yantosca yantosca added this to the 13.3.0 milestone Sep 24, 2021
@yantosca yantosca changed the title [PULL REQUEST] Cap A (alpha) factor at 3.2e-8 when computing gamma(N2O5) [PULL REQUEST] Revise N2O5 gamma calculation Sep 24, 2021
@yantosca yantosca merged commit 5119768 into dev Sep 27, 2021
@yantosca yantosca deleted the feature/gammaN2O5 branch September 27, 2021 15:44
Copy link
Contributor Author

Now merged and tagged with 13.3.0-alpha.5

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
category: Feature Request New feature or request topic: Chemical Mechanisms Related to KPP and/or GEOS-Chem chemistry mechanisms
Projects
None yet
Development

Successfully merging this pull request may close these issues.

[FEATURE REQUEST] Revise N2O5 gamma calculation
2 participants