diff --git a/pyaerocom/aeroval/glob_defaults.py b/pyaerocom/aeroval/glob_defaults.py
index 0ca07b89a..56f1a5833 100644
--- a/pyaerocom/aeroval/glob_defaults.py
+++ b/pyaerocom/aeroval/glob_defaults.py
@@ -120,6 +120,8 @@ class CategoryType(str, Enum):
     temperature = "Temperature"
     particle_ratio = "Particle ratio"
     vertical_column_density = "Vertical column density"
+    surface_emission = "Surface emission"
+    column_burden = "Column burden"
     UNDEFINED = "UNDEFINED"
 
     def __str__(self):
diff --git a/pyaerocom/data/variables.ini b/pyaerocom/data/variables.ini
index 26d2fff90..e9236d1d4 100644
--- a/pyaerocom/data/variables.ini
+++ b/pyaerocom/data/variables.ini
@@ -1305,7 +1305,7 @@ var_name = emiisop
 description = Total isoprene emissions
 standard_name = tendency_of_atmosphere_mass_content_of_isoprene_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1360,7 +1360,7 @@ var_name = emiso2
 description = Total emission of SO2
 standard_name = tendency_of_atmosphere_mass_content_of_sulfur_dioxide_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1371,7 +1371,7 @@ var_name = emidms
 description = Total emission of DMS
 standard_name = tendency_of_atmosphere_mass_content_of_dimethyl_sulfide_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1382,7 +1382,7 @@ var_name = emiso4
 description = Total emission of particulate SO4
 standard_name = tendency_of_atmosphere_mass_content_of_sulfate_dry_aerosol_particles_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1393,7 +1393,7 @@ var_name = emioa
 description = Total emission of OA
 standard_name = tendency_of_atmosphere_mass_content_of_particulate_organic_matter_dry_aerosol_particles_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1426,7 +1426,7 @@ var_name = emibc
 description = Total emission of BC
 standard_name = tendency_of_atmosphere_mass_content_of_elemental_carbon_dry_aerosol_particles_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1437,7 +1437,7 @@ var_name = emiss
 description = Total emission of seasalt
 standard_name = tendency_of_atmosphere_mass_content_of_seasalt_dry_aerosol_particles_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1448,7 +1448,7 @@ var_name = emidust
 description = Total emission of dust
 standard_name = tendency_of_atmosphere_mass_content_of_dust_dry_aerosol_particles_due_to_emission
 var_type = Emission fluxes
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1812,7 +1812,7 @@ var_name = dryso2
 description = Dry deposition of SO2
 standard_name = tendency_of_atmosphere_mass_content_of_sulfur_dioxide_due_to_dry_deposition
 var_type = dry deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1823,7 +1823,7 @@ var_name = dryso4
 description = Dry deposition of SO4
 standard_name = tendency_of_atmosphere_mass_content_of_sulfate_dry_aerosol_particles_due_to_dry_deposition
 var_type = dry deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1845,7 +1845,7 @@ var_name = drydms
 description = Dry deposition of DMS
 standard_name = tendency_of_atmosphere_mass_content_of_dimethyl_sulfide_due_to_dry_deposition
 var_type = dry deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -1856,7 +1856,7 @@ var_name = dryss
 description = Dry deposition of seasalt
 standard_name = tendency_of_atmosphere_mass_content_of_seasalt_dry_aerosol_particles_due_to_dry_deposition
 var_type = dry deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -2030,7 +2030,7 @@ var_name = drybc
 description = Dry deposition of BC
 standard_name = tendency_of_atmosphere_mass_content_of_elemental_carbon_dry_aerosol_particles_due_to_dry_deposition
 var_type = dry deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -2120,7 +2120,7 @@ var_name = dryoa
 description = Dry deposition of OA
 standard_name = tendency_of_atmosphere_mass_content_of_particulate_organic_matter_dry_aerosol_particles_due_to_dry_deposition
 var_type = dry deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -2406,7 +2406,7 @@ var_name = wetbc
 description = Wet deposition of BC
 standard_name = tendency_of_atmosphere_mass_content_of_elemental_carbon_dry_aerosol_particles_due_to_wet_deposition
 var_type = wet deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -2417,7 +2417,7 @@ var_name = wetss
 description = Wet deposition of seasalt
 standard_name = tendency_of_atmosphere_mass_content_of_seasalt_dry_aerosol_particles_due_to_wet_deposition
 var_type = wet deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -2439,7 +2439,7 @@ var_name = wetdust
 description = Wet deposition of dust
 standard_name = tendency_of_atmosphere_mass_content_of_dust_dry_aerosol_particles_due_to_wet_deposition
 var_type = wet deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -2615,7 +2615,7 @@ var_name = wetoa
 description = Wet deposition of OA
 standard_name = tendency_of_atmosphere_mass_content_of_particulate_organic_matter_dry_aerosol_particles_due_to_wet_deposition
 var_type = wet deposition flux
-unit = kg m-2 s-1
+unit = mg m-2 d-1
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -4471,7 +4471,7 @@ var_name = loadso4
 description = Column SO4 mass load
 standard_name = atmosphere_mass_content_of_sulfate_dry_aerosol_particles
 var_type = mass load
-unit = kg m-2
+unit = mg m-2
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -4482,7 +4482,7 @@ var_name = loadbc
 description = Column BC mass load
 standard_name = atmosphere_mass_content_of_elemental_carbon_dry_aerosol_particles
 var_type = mass load
-unit = kg m-2
+unit = mg m-2
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -4504,7 +4504,7 @@ var_name = loadoa
 description = Column OA mass load
 standard_name = atmosphere_mass_content_of_particulate_organic_matter_dry_aerosol_particles
 var_type = mass load
-unit = kg m-2
+unit = mg m-2
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -4526,7 +4526,7 @@ var_name = loadss
 description = Column sea salt mass load
 standard_name = atmosphere_mass_content_of_seasalt_dry_aerosol_particles
 var_type = mass load
-unit = kg m-2
+unit = mg m-2
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -4537,7 +4537,7 @@ var_name = loaddust
 description = Column dust mass load
 standard_name = atmosphere_mass_content_of_dust_dry_aerosol_particles
 var_type = mass load
-unit = kg m-2
+unit = mg m-2
 minimum = 0
 maximum = 10000
 dimensions = time,lat,lon
@@ -4680,7 +4680,7 @@ var_name = loaddms
 description = Column mass of DMS
 standard_name = atmosphere_mass_content_of_dimethyl_sulfide
 var_type = mass load
-unit = kg m-2
+unit = mg m-2
 minimum = 0
 maximum = 10000
 dimensions = tune,lat,lon
@@ -4691,7 +4691,7 @@ var_name = loadso2
 description = Column mass of SO2
 standard_name = atmosphere_mass_content_of_sulfur_dioxide
 var_type = mass load
-unit = kg m-2
+unit = mg m-2
 minimum = 0
 maximum = 10000
 dimensions = tune,lat,lon