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] NEI2016 emissions #920

Merged
merged 2 commits into from Sep 29, 2021
Merged

[PULL REQUEST] NEI2016 emissions #920

merged 2 commits into from Sep 29, 2021

Conversation

msulprizio
Copy link
Contributor

@msulprizio msulprizio commented Sep 28, 2021

This PR corresponds to Feature Request #547.

This inventory will be optional (off by default) and CEDS emissions (1750-2019) will remain the default anthropogenic emissions inventory globally.

Barron Henderson, Lyssa Freese, and Patrick Campbell provided updated
EPA NEI2016 emissions for incorporation into GEOS-Chem. These are to
replace the NEI2011 emissions currently in HEMCO_Config.rc. This inventory
(denoted NEI2016_MONMEAN) is off by default in favor of the CEDS v2 global
anthropogenic emissions which extend to 2019.

For this inventory, the text from the EPA_Config.rc file provided by Barron
has been copied into HEMCO_Config.rc. Fixes made to that text include:
(1) adding entries for SOAP and pFe and (2) removing unnecessary
scale factors 256 and 256 for BCPO and OCPO.

NOTE: Currently mask 1007 (CONUS_MASK) is used for the NEI2016 entries, but
there are fields named for "can" and "mex" containing emissions over Canada
and Mexico. We have added a note to HEMCO_Config.rc instructing users to
revert to the old mask files if they would like to use the emissions over
Canada and Mexico.

For more details, see feature request #547

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
Annual scaling factors for NEI2016 emissions were computed from the EPA
Trends Report for Tier 1 CAPS (obtained 21 Sep 2021). Calculations can be
viewed in ExtData/HEMCO/NEI2016/v2021-06/national_tier1_caps+HEMCOscaling.xlsx.
Values were computed from the "Total without wildfires" field. NH3 and PM2.5
only had annual emission totals without wildfires for 2002-2020, while the
remaining species had annual values for 1990-2020. We include only 2002-2020
scale factors here so all species are consistent. Cl2 and HCl were not
included in the Tier 1 CAPS file, so they are not scaled to other years and
are only applied for 2016.

Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
@msulprizio msulprizio added this to the 13.3.0 milestone Sep 28, 2021
@msulprizio msulprizio self-assigned this Sep 28, 2021
@msulprizio msulprizio changed the title NEI2016 emissions [PULL REQUEST] NEI2016 emissions Sep 28, 2021
Copy link
Contributor

@yantosca yantosca left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This all looks straightforward. Will have to test with NEI emissions turned on but it looks like a straight-up swap of NEI2016 for NEI2011.

Question: why does Ordonez_Ionocarb now use a hierarchy of 1 instead of 50? Is it because the EPA2016 does not have the same emissions? I am sure it is OK but just wanted to know.

@msulprizio
Copy link
Contributor Author

msulprizio commented Sep 29, 2021

Question: why does Ordonez_Ionocarb now use a hierarchy of 1 instead of 50? Is it because the EPA2016 does not have the same emissions?

This appeared to be a typo in the original HEMCO_Config.rc and/or isn't needed anymore. This fix came in with the HEMCO_Config.rc file provided by Barron and Lyssa. This only applies to the ORDONEZ_CH3I entry and caused issues saving out the InvORDONEX_CH3I diagnostic in HEMCO_Diagn.rc. The remaining ORDONEZ_IODOCARB entries have always used hierarchy of 1. Changing the hierarchy value to 1 for ORDONEZ_CH3I helps us save out all species for this inventory in HEMCO_Diagn.rc consistently.

I believe the hierarchy of ORDONEZ_CH3I was 50 before to avoid conflicts with the CH3I_SEAWATER entry which also uses Cat=1 and Hier=1. The CH3I_SEAWATER entry is commented out so there should not be conflicts with the ORDONEZ_IODOCARB entry now.

Copy link
Contributor

Thanks @msulprizio for the explanation. Good to merge.

@msulprizio msulprizio merged commit b7c4b24 into dev Sep 29, 2021
@msulprizio msulprizio deleted the feature/NEI2016 branch September 29, 2021 18:59
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants