README -- for HEMCO/MERCURY/v2022-10 directory
21 Oct 2022
GEOS-Chem Support Team
geos-chem-support@g.harvard.edu

------------------------------------------------------------------------------

Mercury emissions inventory from AMAP/UNEP (2018) Global Mercury Assessment

Annual emissions for 2015 for Hg0, Hg2, and HgP. The original emissions
inventory is uploaded as text files and can be found here:

   Steenhuisen, F.; Wilson, S.J., 2022, "Geospatially distributed (gridded) 
   global mercury emissions to air from anthropogenic sources in 2015",
   https://doi.org/10.34894/SZ2KOI, DataverseNL, V1.

This inventory corresponds to the "standard speciation factors" gridded
inventory described in the publication:

   F. Steenhuisen and S. J. Wilson, Development and application of an
   updated geospatial distribution model for gridding 2015 global mercury
   emissions, Atmos. Environ., 2019, 211, 138–150,
   https://doi.org/10.1016/j.atmosenv.2019.05.003.

The emissions are divided into four source types, used in the 2018 Global
Mercury Assessment:

   1) ASGM - artisanal and small-scale gold mining;
   2) PG or POWERGEN - Combustion of fuels for energy and in industry and
      domestic/residential uses, also called stationary combustion in some
      papers;
   3) INDS - Industry sectors;
   4) INTW - Intentional uses (other than ASGM).

All emissions are in units kg m-2 s-1 of Hg. 

Helene Angot processed these text files to netcdf for use in GEOS-Chem
simulations. The employed R script for converting ASGM emission text files
to netcdf can be found in emissions_ASGM.R. Ari Feinberg concatenated
different source files and managed the upload. 

Included variables are:

   emi_hg_0_asgm = "Hg0 emissions from ASGM" ;
   emi_hg_2_asgm = "Hg2 emissions from ASGM" ;
   emi_hg_p_asgm = "HgP emissions from ASGM" ;
   emi_hg_0_inds = "Hg0 emissions from Industry" ;
   emi_hg_2_inds = "Hg2 emissions from Industry" ;
   emi_hg_p_inds = "HgP emissions from Industry" ;
   emi_hg_0_intw = "Hg0 emissions from Intentional uses" ;
   emi_hg_2_intw = "Hg2 emissions from Intentional uses" ;
   emi_hg_p_intw = "HgP emissions from Intentional uses" ;
   emi_hg_0_pg   = "Hg0 emissions from Stationary combustion" ;
   emi_hg_2_pg   = "Hg2 emissions from Stationary combustion" ;
   emi_hg_p_pg   = "HgP emissions from Stationary combustion" ;

Total Hg emissions:
Hg0 emissions (Mg/yr): 1826
Hg2 + HgP emissions (Mg/yr): 399 


------------------------------------------------------------------------------
Additional note by Bob Yantosca (22 Oct 2022)

Also note: The following commands were needed to make the netCDF file
GMA_emissions_Hg.0.25x0.25.2015.nc COARDS-compliant:

   ncap2 -O -h -s 'defdim("time", 1); time[time]={0D}; time@long_name="Time"; time@units="hours since 2015-01-01 00:00:00"; time@calendar="standard"; emi_hg_0_asgm[$time,$lat,$lon]=emi_hg_0_asgm; emi_hg_2_asgm[$time,$lat,$lon]=emi_hg_2_asgm; emi_hg_p_asgm[$time,$lat,$lon]=emi_hg_p_asgm; emi_hg_0_inds[$time,$lat,$lon]=emi_hg_0_inds; emi_hg_2_inds[$time,$lat,$lon]=emi_hg_2_inds; emi_hg_p_inds[$time,$lat,$lon]=emi_hg_p_inds; emi_hg_0_intw[$time,$lat,$lon]=emi_hg_0_intw; emi_hg_2_intw[$time,$lat,$lon]=emi_hg_2_intw; emi_hg_p_intw[$time,$lat,$lon]=emi_hg_p_intw; emi_hg_0_pg[$time,$lat,$lon]=emi_hg_0_pg; emi_hg_2_pg[$time,$lat,$lon]=emi_hg_2_pg; emi_hg_p_pg[$time,$lat,$lon]=emi_hg_p_pg' GMA_emissions_Hg.0.25x0.25.2015.nc tmp.nc

   nccopy -d1 tmp.nc GMA_emissions_Hg.0.25x0.25.2015.nc