HEMCO extensions¶
Overview¶
Emission inventories sometimes include dynamic source types and nonlinear scale factors that have functional dependencies on local environmental variables such as wind speed or temperature, which are best calculated online during execution of the model. HEMCO includes a suite of additional modules (extensions) that perform online emission calculations for a variety of sources (see list below). Extensions are enabled in section Extension Switches of the HEMCO configuration file.
List of extensions¶
The full list of available extensions is given below. Extensions can be selected individually in the Extension Switches section of the The HEMCO configuration file, as can the species to be considered.
-
DustAlk
¶
Species: DSTAL1, DSTAL2, DSTAL3, DSTAL4
Reference: Fairlie et al (check)
-
DustDead
¶
Emissions of mineral dust from the DEAD dust mobilization model.
Species: DST1, DST2, DST3, DST4
Reference: Zender et al. [2003]
-
DustGinoux
¶
Emissions of mineral dust from the P. Ginoux dust mobilization model.
Species: DST1, DST2, DST3, DST4
Reference: Ginoux et al. [2001]
-
FINN
¶
Biomass burning emissions from the FINN model.
Species: NO, CO, ALK4, ACET, MEK, ALD2, PRPE, C2H2, C2H4, C3H8, CH2O, C2H6, SO2, NH3, BCPI, BCPO, OCPI, OCPO, GLYC, HAC, SOAP
Reference: Wiedinmyer et al. [2014]
-
GC_Rn-Pb-Be
¶
Emissions of radionuclide species as used in the GEOS-Chem model.
Species: Rn222, Be7, Be7Strat, Be10, Be10Strat
-
ZHANG_Rn222
¶
If
ZHANG_Rn222
ison
, then Rn222 emissions will be computed according to Zhang et al. [2021].If
ZHANG_Rn222
isoff
, then Rn222 emissions will be computed according to Jacob et al. [1997].
-
GFED
¶
Biomass burning emissions from the GFED model.
Version: GFED3 and GFED4 are available.
Species: NO, CO, ALK4, ACET, MEK, ALD2, PRPE, C2H2, C2H4, C3H8, CH2O C2H6, SO2, NH3, BCPO, BCPI, OCPO, OCPI, POG1, POG2, MTPA, BENZ, TOLU, XYLE NAP, EOH, MOH, SOAP,
Reference: van der Werf et al. [2010]
-
Inorg_Iodine
¶
Species: HOI, I2
Reference: TBD
-
LightNOx
¶
Emissions of NOx from lightning.
Species: NO
Species: [Murray et al., 2012]
-
MEGAN
¶
Biogenic VOC emissions.
Version: 2.1
Species: ISOP, ACET, PRPE, C2H4, ALD2, CO, OCPI, MONX, MTPA, MTPO, LIMO, SESQ
Reference: Guenther et al. [2012]
-
SeaFlux
¶
Air-sea exchange.
-
SeaSalt
¶
Sea salt aerosol emission.
-
Volcano
¶
Emissions of volcanic SO2 from AEROCOM.
Species: SO2
Reference:
-
TOMAS_Jeagle
¶
Size-resolved sea salt emissions for TOMAS aerosol microphysics simulations.
Species: SS1, SS2, SS3, SS4, SS5, SS6, SS7, SS8, SS9, SS10, SS11, SS12, SS13, SS14, SS15, SS16, SS17, SS18, SS19, SS20, SS21, SS22, SS23, SS24, SS25, SS26, SS27, SS28, SS29, SS30, SS31, SS32, SS33, SS34, SS35, SS36, SS37, SS38, SS39, SS40
Reference: Jaeglé et al. [2011]
-
TOMAS_DustDead
¶
Size-resolved dust emissions for TOMAS aerosol microphysics simulations.
Species: DUST1, DUST2, DUST3, DUST4, DUST5, DUST6, DUST7, DUST8, DUST9, DUST10, DUST11, DUST12, DUST13, DUST14, DUST15, DUST16, DUST17, DUST18, DUST19, DUST20, DUST21, DUST22, DUST23, DUST24, DUST25, DUST26, DUST27, DUST28, DUST29, DUST30, DUST31, DUST32, DUST33, DUST34, DUST35, DUST36, DUST37, DUST38, DUST39, DUST40
Reference: Zender et al. [2003]
Gridded data¶
HEMCO can host all environmentally independent data sets (e.g. source functions) used by the extensions. The environmental variables are either provided by the atmospheric model or directly read from file through the HEMCO configuration file. Entries in the HEMCO configuration file file are given priority over fields passed down from the atmospheric model, i.e. if the HEMCO configuration file contains an entry for a given environmental variable, this field will be used instead of the field provided by the atmospheric model. The field name provided in the HEMCO configuration file must exactly match the name of the HEMCO environmental parameter.
To use the NCEP reanalysis monthly surface wind fields (http:, , www.esrl.noaa.gov, psd, data, gridded, data.ncep.reanalysis.derived.surface.html) in all HEMCO extensions, add the following two lines to the Base Emissions section of the HEMCO configuration file:
* U10M /path/to/uwnd.mon.mean.nc uwnd 1948-2014/1-12/1/0 C xy m/s * - 1 1
* V10M /path/to/vwnd.mon.mean.nc vwnd 1948-2014/1-12/1/0 C xy m/s * - 1 1
This will use these wind fields for all emission calculations, even if the atmospheric model uses a different set of wind fields.
It is legal to assign scale factors (and masks) to the environmental variables read through the HEMCO configuration file. This is particularly attractive for sensitivity studies. For example, a scale factor of 1.1 can be assigned to the NCEP surface wind fields to study the sensitivity of emissions on a 10% increase in wind speed:
In the Base Emissions section:
* U10M /path/to/uwnd.mon.mean.nc uwnd 1948-2014/1-12/1/0 C xy m/s * 123 1 1
* V10M /path/to/vwnd.mon.mean.nc vwnd 1948-2014/1-12/1/0 C xy m/s * 123 1 1
In the Scale Factors section:
123 SURFWIND_SCALE 1.1 - - - xy 1 1
As for any other entry in the HEMCO configuration file, spatially uniform values can be set directly in the HEMCO configuration file. For example, a spatially uniform, but monthly varying surface albedo can be specified by adding the following entry to the Base Emissions section of the HEMCO configuration file:
* ALBD 0.7/0.65/0.6/0.5/0.5/0.4/0.45/0.5/0.55/0.6/0.6/0.7 - 2000/1-12/1/0 C xy 1 * - 1 1
Environmental fields used by HEMCO¶
The following fields can be passed from the atmospheric model to HEMCO for use by the various extensions:
-
AIR
¶
Air mass.
Dim: xyz
Units: kg
Used by:
GC_Rn-Pb-Be
,PARANOx
-
FRAC_OF_PBL
¶
Fraction of grid box within the planetary boundary layer (PBL).
-
FRCLND
¶
Land fraction
Dim: xy
Units: unitless
Used by:
GC_Rn-Pb-Be
,SeaFlux
-
SNOWHGT
¶
Snow height (mm of H2O equivalent).
Dim: xy
Units: kg H2O/m2
Used by:
DustDead
,TOMAS_DustDead
-
SPHU
¶
Specific humidity
Dim: xyz
Units: kg H2O/kg air
Used by:
DustDead
,PARANOx
,TOMAS_DustDead
-
TK
¶
Temperature.
Dim: xyz
Units: K
Used by:
DustDead
,LightNOx
,TOMAS_DustDead
-
TROPP
¶
Tropopause pressure.
Dim: xy
Units: Pa
Used by:
GC_Rn-Pb-Be
,LightNOx
-
U10M
¶
E/W wind speed @ 10 meters above surface.
Dim: xy
Units: m/s
Used by:
DustAlk
,DustDead
,DustGinoux
,PARANOx
,SeaFlux
,SeaSalt
,SoilNOx
,TOMAS_DustDead
,TOMAS_Jeagle
-
USTAR
¶
Friction velocity.
Dim: xy
Units: m/s
Used by:
DustDead
,TOMAS_DustDead
-
V10M
¶
N/S wind speed @ 10 meters above surface.
Dim: xy
Units: m/s
Used by:
DustAlk
,DustDead
,DustGinoux
,PARANOx
,SeaFlux
,SeaSalt
,SoilNOx
,TOMAS_DustDead
,TOMAS_Jeagle
-
WLI
¶
Water-land-ice flags (
0
= water,1
= land,2
= ice).Dim: xy
Units: unitless
Used by: Almost every extension
-
Z0
¶
Roughness height.
Dim: xy
Units: m
Used by:
DustDead
,TOMAS_DustDead
Restart variables¶
Some extensions rely on restart variables, i.e. variables that are
highly dependent on historical information such as previous-day leaf
area index or soil NOx pulsing factor. During a simulation run, the
extensions continuously archive all necessary information and update
estart variables accordingly. The updated variables become
automatically written into the HEMCO restart file
(HEMCO_restart.YYYYMMDDhhmmss.nc
) at the end of a
simulation. The fields from this file can then be read through the
HEMCO configuration file to resume the simulation at this date (“warm”
restart). For example, the soil NOx restart variables can be made
available to the soil NOx extension by adding the following lines to
the Base Emissions section of the HEMCO
configuration file.
104 PFACTOR ./HEMCO_restart.$YYYY$MM$DD$HH00.nc PFACTOR $YYYY/$MM/$DD/$HH E xy unitless NO - 1 1
104 DRYPERIOD ./HEMCO_restart.$YYYY$MM$DD$HH00.nc DRYPERIOD $YYYY/$MM/$DD/$HH E xy unitless NO - 1 1
104 GWET_PREV ./HEMCO_restart.$YYYY$MM$DD$HH00.nc GWET_PREV $YYYY/$MM/$DD/$HH E xy unitless NO - 1 1
104 DEP_RESERVOIR ./HEMCO_restart.$YYYY$MM$DD$HH00.nc DEP_RESERVOIR $YYYY/$MM/$DD/$HH E xy unitless NO - 1 1
Many restart variables are very time and date-dependent. It is therefore recommended to set the time slice selection flag to E to ensure that only data is read that exactly matches the simulation start date (also see Base emissions. HEMCO will perform a “cold start” if no restart field can be found for a given simulation start date, e.g. default values will be used for those restart variables.
Built-in tools for scaling/masking¶
HEMCO has built-in tools to facilitate the application of both uniform and spatiotemporal scale factors to emissions calculated by the extensions. At this point, not all extensions take advantage of these tools yet. A list of extensions that support the built-in scaling tools are given below.
For extensions that support the built-in scaling tools, you can specify the uniform and/or spatiotemporal scale factors to be applied to the extension species of interest in section Extension switches the HEMCO configuration file.
For example, to uniformly scale GFED CO by a factor of 1.05 and GFED NO emissions by a factor of 1.2, add the following two lines to the HEMCO configuration file (highlighted in GREEN):
111 GFED : on CO/NO/ACET/ALK4
--> GFED3 : false
--> GFED4 : true
--> GFED_daily : false
--> GFED_3hourly : false
--> Scaling_CO : 1.05
--> Scaling_NO : 1.20
Similarly, a spatiotemporal field to be applied to the species of
interest can be defined via setting ScaleField
, e.g.
111 GFED : on CO/NO/ACET/ALK4
--> GFED3 : false
--> GFED4 : true
--> GFED_daily : false
--> GFED_3hourly : false
--> Scaling_CO : 1.05
--> Scaling_NO : 1.20
--> ScaleField_NO : GFED_SCALEFIELD_NO
The corresponding scale field needs be defined in section Base emissions . A simple example would be a monthly varying scale factor for GFED NO emissions:
111 GFED_SCALEFIELD_NO 0.9/1.1/1.3/1.4/1.6/1.7/1.7/1.8/1.5/1.3/0.9/0.8 - 2000/1-12/1/0 C xy unitless * - 1 1
It is legal to apply scale factors and/or masks to the extension scale fields (in the same way as the ‘regular’ base emission fields). A more sophisticated example on how to scale soil NOx emissions is given in HEMCO examples.