Dear Sir or Madam,
I am running the WRF–CMAQ modeling system and would like to include biomass burning emissions in CMAQ.
I found through literature that the GFED dataset is commonly used, but it is only available up to 2017, while my simulation year is 2023.
I would like to ask:
What other approaches or datasets can be used to incorporate biomass burning emissions into CMAQ?
Is there any official guidance or tutorial on processing GFED emissions for CMAQ?
I sincerely appreciate your assistance and guidance.
Thank you very much.
Hi,
GFED data is available through 2024 ( Data - Global Fire Emissions Database (GFED) ). The GFED inventory, like many global biomass burning inventories, is distributed as a pre-gridded netCDF. There is currently no official guidance or tutorial on preparing GFED for CMAQ. At a minimum you will need to:
- Interpolate from the 0.25 degree native grid to your model projection and grid
- Temporalize the emissions from daily to hourly
- Map and speciate the inventory species to your chemical mechanism
- Vertically allocate the emissions
This will all need to be done in such a way that the resulting emissions can be used by CMAQ. One potential option is to process the gridded inventory through SMOKE with methods similar to those used to process global HTAP emissions ( C. Processing HTAPv3 gridded global emission inventories · CEMPD/SMOKE Wiki · GitHub ).
The Fire INventory from NCAR is often used in EPA modeling platforms for global and hemispheric biomass burning emissions. EPA also develops biomass burning inventories over the US. Could you please further describe what region you are planning to model?
Thank you very much for your reply. The region I am simulating is actually mainland China, not the United States. I have another question: when using SMOKE to process the GFED dataset, how can I determine its temporal and species allocation coefficients? Thank you in advance for any guidance you can provide.
Diurnal temporal allocation factors can be extracted from the GFED 3-hourly files, averaged, and interpolated. You may also consider using an average profile that peaks mid-afternoon such as profile hourly profile 210 distributed in the EPA EMP ancillary files. The previously linked HTAP processing instructions provides temporal profile assignment examples.
To map the GFED species you will need to first determine corresponding individual and lumped species in your target chemical mechanism, then distribute the species mass using fractions and other factors. This can either be done by adding species in the netCDF inventory or through the SMOKE inventory table and speciation profiles. For example, formaldehyde in GFED can be directly mapped to FORM in the CB6 chemical mechanism by creating an entry in the inventory table that directly maps the inventory species to the model species FORM.