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

Correcting Injection height in biomass burning inventories #2330

Open
jinlx opened this issue Jun 14, 2024 · 3 comments · May be fixed by #2358
Open

Correcting Injection height in biomass burning inventories #2330

jinlx opened this issue Jun 14, 2024 · 3 comments · May be fixed by #2358
Labels
category: Feature Request New feature or request topic: Emissions Related to emissions inventories used in GEOS-Chem

Comments

@jinlx
Copy link

jinlx commented Jun 14, 2024

Your name

Lixu Jin

Your affiliation

University of Montana

Provide a clear and concise overview of the new feature requested.

Current HEMCO treats injection height in the same manner as emission flux, leading to an artificial decrease in injection height when averaging at 0.1ºx0.1º resolution. Here we provided a workaround to resolve the potential injection height underestimation.

This involves implementing weighted averaging of injection heights based on flux amount (special thanks to Dr. Jeff Pierce for insights). Scripts in ipynb/py formats are provided to generate standalone injection heights for various times and resolutions. This correction rectifies the artificial reduction of injection heights by accurately considering the emission strengths of fire pixels.

GFAS_injs_ipynb.txt
GFAS_injs_py.txt

The code generates netCDF files containing one gridded variable: injection height (mami_avg_weight). To implement the new injection heights in GEOS-Chem, two files are required:

  • Original emissions files (0.1ºx0.1º)
  • New netCDF file with mami_avg_weight

Modify the following lines in the HEMCO.Config file:

From:

(((GFAS
300 GFAS_EMITL $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc mami 2003-2021/1-12/1-31/0 C xy m 1
)))GFAS

To:

(((GFAS
300 GFAS_EMITL $ROOT/GFAS/v2021-09/$YYYY/GFAS_$YYYY$MM_hgts_0.25x0.3125_NA.nc mami_avg_weight 2003-2021/1-12/1-31/0 C xy m 1
)))GFAS

This update can potentially be applied to other offline emission inventories like QFED. For online-calculated inventories such as GFED4 and FINN, the direct implementation may differ. One possible workaround is to calculate GFED/FINN offline using HEMCO standalone and apply the aforementioned method to incorporate correct injections.

Please feel free to let me know if there are any things I can further help or any issues in the processing scripts.

Reference(s)

Jin, L., Permar, W., Selimovic, V., Ketcherside, D., Yokelson, R. J., Hornbrook, R. S., Apel, E. C., Ku, I.-T., Collett Jr., J. L., Sullivan, A. P., Jaffe, D. A., Pierce, J. R., Fried, A., Coggon, M. M., Gkatzelis, G. I., Warneke, C., Fischer, E. V., and Hu, L.: Constraining emissions of volatile organic compounds from western US wildfires with WE-CAN and FIREX-AQ airborne observations, Atmos. Chem. Phys., 23, 5969–5991, https://doi.org/10.5194/acp-23-5969-2023, 2023.

@jinlx jinlx added the category: Feature Request New feature or request label Jun 14, 2024
@jinlx jinlx changed the title Put a one-line description here Injection height in biomass burning inventories Jun 14, 2024
@jinlx jinlx changed the title Injection height in biomass burning inventories Correcting Injection height in biomass burning inventories Jun 15, 2024
@jinlx
Copy link
Author

jinlx commented Jun 22, 2024

I mistakenly created an issue instead of a pull request (PR) for the updates. Could you please convert it into a PR? @yantosca.

@yantosca
Copy link
Contributor

Thanks @jinlx. There isn't an easy way to convert it to a PR, it'd have to be done manually. Would you be able to open a PR again?

@yantosca yantosca added the topic: Emissions Related to emissions inventories used in GEOS-Chem label Jun 25, 2024
@jinlx
Copy link
Author

jinlx commented Jul 1, 2024

Hi @yantosca , just figured out how to create a PR and please feel free to let me know if there are any issues or something I can further help!

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: Emissions Related to emissions inventories used in GEOS-Chem
Projects
None yet
Development

Successfully merging a pull request may close this issue.

2 participants