Hello,
I am running a CMAQ5.5 simulation over CONUS using CRACMMv2 and I was interested in testing the effects of zeroing out certain emissions. Specifically, I would like to test out the effects of having zero NH3 fertilizer emissions or zero soil NOx emissions in two runs.
To zero out NH3_fert, should I simply set CTM_abflux to N? For context, I am using ag emissions files from the EQUATES Project that only include NH3 from livestock.
For soil NOx I am a bit confused on how to proceed. From what I understand, soil NOx is calculated using the inline BEIS biogenic emissions. To adjust the biogenic emissions, would I need to edit the NO split factor from 1.0 to 0.0 in the gspro_biogenics.txt file (screenshots below)?
Hello -
Do you mind sending me your run script so I can better answer your Ag NH3 question?
For Soil NO, you can zero out those emissions by placing this line in your CMAQ_Chem_Control file like so:
And then to ensure that these emissions are being zero’d out you can add this as a stream in your CMAQ_Control.nml file like so:
Hi,
Thanks for your advice! Here is my cctm run script for more clarification. I am modeling over the 12km CONUS domain for Summer 2005 using the CRACMM2 mechanism and BEIS biogenic emissions.
run_cctm_2005_12US1_CRACMM2.csh.txt (39.6 KB)
Hello Again,
I have an addendum to my original response. Hopefully this will be more helpful! And have a solution to your fertilizer NH3 question as well thanks to @nash.skipper.
-
When I refer to the CMAQ_Control.nml file, for you that would be the CMAQ_Control_DESID.nml file and then when I refer to the CMAQ_Chem_Control file, I’m referring to your CMAQ_Control_DESID_cracmm2.nml file.
-
For zero’ing out NO emissions from soil in BEIS you would insert the scaling rule I screenshotted before in your CMAQ_Control_DESID_cracmm2.nml file, but instead of “a” under the Op column, you would change this to “o” (which stands for override).
-
To zero out NH3 emissions from fertilizer (from both Bidirectional Fluxes and from your offline EQUATES emissions file - which should be the emissions stream labeled “GRIDDED_EMIS” in your runscript), you would indeed set CTM_ABFLUX to N and then in your CMAQ_Control_DESID_cracmm2.nml you would add this line:
Note that the “a” is retained (meaning add) as you are negatively adding NH3 from fertilizer from your total NH3 emissions (or in plain speak: subtracting).
As a sanity check you can output these streams as I mentioned before by modifying your CMAQ_Control_DESID.nml file. Your NO from the BIOG DESID stream should be 0, and your NH3 from your GRIDDED_EMIS should be less than if you ran without setting that scaling factor to -1.
Hope this helps and good luck with your project!
-Sara
Thank you for all the help! I’ll let you know if I run into any issues.