search

geoschem / geos-chem

GEOS-Chem "Science Codebase" repository. Contains GEOS-Chem science routines, run directory generation scripts, and interface code. This repository is used as a submodule within the GCClassic and GCHP wrappers, as well as in other modeling contexts (external ESMs).
http://geos-chem.org
Other
162 stars 155 forks source link

Correcting Injection height in biomass burning inventories #2330

Open jinlx opened 1 month ago

jinlx commented 1 month ago

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:

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 commented 1 month ago

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

yantosca commented 3 weeks ago

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?

jinlx commented 2 weeks ago

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!