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
166 stars 160 forks source link

Biomass burning CO emissions from GEOSChem.Emissions.20100101_0000z.nc4 are ~4X higher than raw data (GFAS_201001.nc) #2234

Open herizonarivelo opened 6 months ago

herizonarivelo commented 6 months ago

Name and Institution (Required)

Name: Herizo Narivelo Institution: LAERO/OMP Université de Toulouse - Paul Sabatier

Confirm you have reviewed the following documentation

Description of your issue or question

Dear all,

I am using GCHP v14.2.2 for my project. I previously used GFED4 inventories for biomass burning emission. However, GFED4 is not available for the years prior to 2010 that I would like to start my simulation.So, as Bob recommends in my previous questions in other issue #2212, I have to use either GFAS or QFED for biomass burning emissions instead GFED4 if I would like to start my simulation prior 2010 or also set the time cycle flag for all of the GFED entries in HEMCO_Config.rc to C.

So I ran one month simulation test without the chemiistry mechanism to test emission by using GFAS biomass burning with a c48 resolution. So, I found a larger difference in the amount of the CO emission [Tg/month] in the raw data located in the Extdata... and the output of my simulation test. The CO from biomass burnning simulated in GEOS-Chem is very high compared with the raw data (/../../geoschem/ExtData/HEMCO/GFAS/v2018-09/2010). For information, I followed the recommendations in this issue #492 #515 but I didn't solve my problem. My requirement is how to remove these scaling factors in HEMCO_Config entry. I think that these scaling factors are the main cause of this large difference between the raw data and the GEOS-Chem emission diagnostics.

I tried to remove these scaling factors (75 and 300) and got errors. Maybe I didn't do it right so I would like to know how to remove these factors correctly without getting errors in the output.

Here I attached the HEMCO_Config.rc and HEMCO_Diag.rc

Thank you very much and I look forward to hearing from you !

Kind regards, Herizo Narivelo

HEMCO_Diagn.rc.txt HEMCO_Config.rc.txt

yantosca commented 5 months ago

Thanks for writing @herizonarivelo. In your HEMCO_Config.rc, you have this block where the scale factor 75 is still activated:

#==============================================================================
# --- GFAS biomass burning ---
#==============================================================================
(((GFAS
0 GFAS_CO    $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc cofire        2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s CO   75       5 3
0 GFAS_SOAP  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc cofire        2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s SOAP 75/281   5 3
0 GFAS_NO    $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc noxfire       2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s NO   75       5 3
0 GFAS_BCPI  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc bcfire        2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s BCPI 70/75    5 3
0 GFAS_BCPO  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc bcfire        2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s BCPO 71/75    5 3
0 GFAS_OCPI  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc ocfire        2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s OCPI 72/75    5 3
0 GFAS_OCPO  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc ocfire        2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s OCPO 73/75    5 3
0 GFAS_POG1  -                                          -             -                     - -             -       POG1 74/76/75 5 3
0 GFAS_POG2  -                                          -             -                     - -             -       POG2 74/77/75 5 3
0 GFAS_CO2   $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc co2fire       2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s CO2  75       5 3
0 GFAS_CH4   $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc ch4fire       2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s CH4  75       5 3
0 GFAS_SO2   $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc so2fire       2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s SO2  75       5 3
0 GFAS_pFe   $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc so2fire       2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s pFe  75/66    5 3
0 GFAS_NH3   $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc nh3fire       2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s NH3  75       5 3
0 GFAS_ACET  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c3h6ofire     2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s ACET 75       5 3
0 GFAS_ALD2  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c2h4ofire     2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s ALD2 75       5 3
0 GFAS_ALK4  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc hialkanesfire 2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s ALK4 75       5 3
0 GFAS_PRPE1 $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc hialkenesfire 2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s PRPE 75       5 3
0 GFAS_PRPE2 $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c3h6fire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s PRPE 75       5 3
0 GFAS_C2H6  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c2h6fire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s C2H6 75       5 3
0 GFAS_C3H8  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c3h8fire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s C3H8 75       5 3
0 GFAS_CH2O  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc ch2ofire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s CH2O 75       5 3
0 GFAS_C2H4  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c2h4fire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s C2H4 75       5 3
0 GFAS_ISOP  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c5h8fire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s ISOP 75       5 3
0 GFAS_DMS   $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c2h6sfire     2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s DMS  75       5 3
0 GFAS_TOLU  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c7h8fire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s TOLU 75       5 3
0 GFAS_BENZ  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c6h6fire      2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s BENZ 75       5 3
0 GFAS_XYLE  $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc c8h10fire     2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s XYLE 75       5 3
)))GFAS
#0 GFAS_CH3OH $ROOT/GFAS/v2018-09/$YYYY/GFAS_$YYYY$MM.nc ch3ohfire     2003-2021/1-12/1-31/0 C xyL=1:scal300 kg/m2/s MOH  75       5 3

For those lines where 75 is the only scale factor, you can replace it with a - character. Otherwise you can just delete it (i.e. change 71/75 to 71, etc.

Also have you uncommented the GFAS entries in your ExtData.rc file? This is where GCHP will get the file paths and time refresh info, but the scale factors come from HEMCO_Config.rc.

herizonarivelo commented 5 months ago

Hi Bob,

Thanks for your reply and for your explanation. Of course, I uncommented the GFAS path in the Extdata and for be sure that I'm using GFAS inventories, I commented out the GFDE4 path. However, after replacing 75 with - and deleting it for certain lines, the result does not realy change but it slightly reduces the factors between the GFAS raw data in (Tg CO) and the simulated using GFAS. I have also commented out this which refer the scale factor for GFAS in HEMCO_Config.

==============================================================================

(((QFED2.or.GFAS.or.BB4MIPS 75 QFED2_TOD 0.1392/0.1392/0.1368/0.1368/0.1368/0.1368/0.1368/0.1368/0.1368/0.48/0.96/1.68/2.4/3.12/3.84/4.08/2.88/1.68/0.96/0.1368/0.1368/0.1368/0.1368/0.1368 - - - xy unitless 1 )))QFED2.or.GFAS.or.BB4MIPS

==============================================================================

I wonder if the L=1:300 plays another role in the scale factor. As far as I understand, if I remove the L=1:scal300, the species mapping will place all emissions at the surface? And how about the this line ?

==============================================================================

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

Here is my Extdata.rc and HEMCO_Config.rc files.

Thanks a lot ! 

Best regards, Herizo

HEMCO_Config.rc_n75.txt ExtData.rc.txt