Aerosol emissions: Difference between revisions

From Geos-chem
Jump to navigation Jump to search
No edit summary
No edit summary
Line 1: Line 1:
This page contains links to the relevant pages which detail the aerosol emissions inventories used by GEOS-Chem.
This page contains links to the relevant pages which detail the aerosol emissions inventories used by GEOS-Chem.


<span style="color:green">'''''NOTE: In [[GEOS-Chem v10-01]] and higher versions, all emissions are handled by the [[HEMCO|HEMCO emissions component]].'''''</span>
== Emissions by sector ==


== Anthropogenic, biofuel, and natural source emissions ==
*[[HEMCO_data_directories#Aerosol_emissions|Aerosol emission inventories]]
 
*[[HEMCO_data_directories#Anthropogenic_and_biofuel_emissions|Anthropogenic and biofuel emission inventories]]
The various sectors of aerosol emissions are now handled by the [[HEMCO|HEMCO emissions component]] in [[GEOS-Chem v10-01]] and higher versions.  For more information about individual inventories that are available, please visit the following links:
*[[HEMCO_data_directories#Anthropogenic_aircraft_and_ship_emissions|Aircraft and ship emission inventories]]
 
*[[HEMCO_data_directories#Biomass_burning_emissions|Biomass burning emission inventories]]
#[[HEMCO_data_directories#Aerosol_emissions|Aerosol emissions available via HEMCO]]
*[[HEMCO_data_directories#Emissions_implemented_as_HEMCO_extensions|Emissions implemented as HEMCO extensions]]
#[[HEMCO_data_directories#Anthropogenic_and_biofuel_emissions|Anthropogenic and biofuel emissions available via HEMCO]]
**[[Mineral_dust_aerosols#Overview|Mineral dust emissions]]
#[[HEMCO_data_directories#Anthropogenic_aircraft_and_ship_emissions|Aircraft and ship emissions available via HEMCO]]
**[[Sea salt aerosols#Updates to sea salt emissions algorithm|Sea salt aerosol emissions]]
 
**[[Volcanic SO2 emissions from Aerocom]]
--[[User:Bmy|Bob Y.]] ([[User talk:Bmy|talk]]) 20:19, 26 October 2015 (UTC)
**[http://acmg.seas.harvard.edu/publications/2004/park2004.pdf DMS emissions]
 
== Biomass burning emissions ==
 
The [[GFED4 biomass burning emissions|GFED4 inventory]] can be used to compute biomass burning emisisons for the following species:
 
#SO2
#NH3
#EC = elemental carbon (aka "black carbon")
#OC = organic carbon
 
for the years 1997 thru 2014.
 
The following inventories are now considered obsolete (but you may still use them for research purposes):
 
#[[GFED2 biomass burning emissions|GFED2]]
#[[GFED3 biomass burning emissions|GFED3]]
#''Duncan et al'' [2003] seasonal and interannual biomass emissions (cf ''Bey et al'' [2001])
 
--[[User:Bmy|Bob Y.]] ([[User talk:Bmy|talk]]) 20:23, 26 October 2015 (UTC)
 
== DMS emissions ==
 
<span style="color:green">'''''In [[GEOS-Chem v10-01]] and higher versions, the sea-air flux of DMS is now handled by the [[HEMCO]] seaflux extension.'''''</span>
 
From [http://acmg.seas.harvard.edu/publications/2004/park2004.pdf ''Park et al'' <nowiki>[2004]</nowiki>]:
 
<blockquote>The oceanic emission of DMS is calculated calculated as the product of local seawater DMS concentration and sea-to-air transfer velocity. The seawater DMS concentrations are gridded monthly averages from Kettle et al. [1999], and the transfer velocity of DMS is computed using an empirical formula from Liss and Merlivat [1986] as a function of the surface (10 m) wind speed.</blockquote>
 
== EC and OC emissions ==
 
The following emissions inventories for elemental carbon (EC) and organic carbon (OC) are used in GEOS-Chem.
 
# [[Carbonaceous aerosols#EC_and_OC_emissions|''Bond et al'' <nowiki>[2007]</nowiki>: EC and OC emissions over North America]]
# [[Carbonaceous_aerosols#Oceanic_OC_emissions|''Spracklen et al'': Oceanic OC emissions]] '''''(NOTE: to be included post [[GEOS-Chem v8-03-01]])'''''
 
== Mineral dust emissions ==
 
The [[Mineral_dust_aerosols#Overview|mineral dust emissions algorithm]] is described on a separate wiki page.
 
== Sea salt emissions==
 
Please see this discussion about [[Sea salt aerosols#Updates to sea salt emissions algorithm|recent updates to the sea salt aerosol emissions]].
 
== Ship emissions ==
 
Please see our [[Ship emissions|ship emissions wiki page]] for a full description of the SO2 ship emission inventories that are available for use in GEOS-Chem.
 
== Volcanic emissions ==
 
Please see our [[Volcanic SO2 emissions from Aerocom]] wiki page.
 
--[[User:Bmy|Bob Y.]] 13:30, 19 March 2010 (EDT)


== Online emission of marine primary organic aerosol (POA) ==
== Online emission of marine primary organic aerosol (POA) ==

Revision as of 17:11, 11 June 2019

This page contains links to the relevant pages which detail the aerosol emissions inventories used by GEOS-Chem.

Emissions by sector

Online emission of marine primary organic aerosol (POA)

This update was validated with 1-month benchmark simulation v11-01b and 1-year benchmark simulation v11-01b-Run0. This version was approved on 19 Aug 2015.

Matthew Johnson wrote:

The online emission parameterization for sub-micron marine primary organic aerosol (POA) was developed for implementation into GEOS-Chem at North Carolina State University (Gantt et al., 2012). The wind speed dependent size-resolved parameterization calculates the organic mass fraction of sea spray aerosol under varying levels of surface ocean chlorophyll-α concentrations ([Chl-α]). The parameterization was designed to be used for all global and regionally nested spatial resolutions and was tested in GEOS-Chem v9-02. The submicron marine POA emissions proved to improve the model-predicted concentrations of total organic aerosol concentrations over coastal and marine regions (Gantt et al., 2015).

References for this work include:

  • Gantt, B., Johnson, M. S., Meskhidze, N., Sciare, J., Ovadnevaite, J., Ceburnis, D., and O'Dowd, C. D.: Model evaluation of marine primary organic aerosol emission schemes, Atmos. Chem. Phys., 12, 8553-8566, doi:10.5194/acp-12-8553-2012, 2012.
  • Gantt, B., Johnson, M. S., Crippa, M., Prévôt, A. S. H., and Meskhidze, N.: Implementing marine organic aerosols into the GEOS-Chem model, Geosci. Model Dev., 8, 619-629, doi:10.5194/gmd-8-619-2015, 2015.

As part of this update, Matthew Johnson has provided MODIS-Aqua chlorophyll-A data. For more information about these data, please visit this post on our Leaf area indices in GEOS-Chem wiki page.

NOTE: A minor bug fix for this update was included in v11-01d to declare variable CHLR as OMP private in routine HCOX_SeaSalt_Run in file HEMCO/Extensions/hcox_seasalt_mod.F90.

--Lizzie Lundgren (talk) 15:21, 22 October 2015 (UTC)

Bug fix: Allocate the OCCONV array for marine POA simulations

This fix was included in v11-01j and approved on 03 Dec 2016

While unit testing GEOS-Chem with the GNU Fortran compiler, we discovered an error in the marine POA simulation. The array OCCONV should have been allocated in routine INIT_SEASALT (in module GeosCore/seasalt_mod.F) but wasn't. We have now added the following IF statement in GREEN, below the existing code: 

       IF ( AS /= 0 ) CALL ALLOC_ERR( 'DMID' )
       DMID = 0e+0_fp

      ! Allocate OCCONV only for marine-POA simulations (bmy, 10/13/16)
      IF ( Input_Opt%LMPOA ) THEN
         ALLOCATE( OCCONV( IIPAR, JJPAR, LLPAR ), STAT=AS )
         IF ( AS /= 0 ) CALL ALLOC_ERR( 'OCCONV' )
         OCCONV = 0e+0_fp
      ENDIF

This will only allocate OCCONV if we are doing a marine POA simulation.

--Bob Yantosca (talk) 18:41, 13 October 2016 (UTC)

Bug fix for MOPO and MOPI definitions in species database

This update was included in v11-02c and approved on 21 Sep 2017.

Katie Travis wrote:

I believe I have come across a bug in the marine organic simulation. In v11, the marine organic aerosol species (MOPI/MOPO) are given their own tracers (I believe they were added into OCPI/OCPO in v9-2?). There is a note on the wiki that these species in the species database are set to gas, when they should be aerosol. But they should also wet deposit, and in fact the parameters for that are on the wiki table. I ran this by one of the developers, Matt Johnson, and he also agrees this appears to be a bug.

Matthew Johnson wrote:

Yes that is correct. Both marine organic tracers should be dry deposited. For wet deposition, MOPO should be treated the same as OCPO (washout only, if this is still how OCPO is treated in the new model) and MOPI is treated the same as OCPI.

To fix this issue, the following lines in Headers/species_database_mod.F90 should be changed from: 

         CASE( 'MOPI' )
            
            ! Halve the Kc (cloud condensate -> precip) rate
            ! for the temperature range 237 K <= T < 258 K.
            KcScale = (/ 1.0_fp, 0.5_fp, 1.0_fp /)

            ! Turn off rainout only when 237 K <= T < 258K.
            RainEff = (/ 1.0_fp, 0.0_fp, 1.0_fp /)

            CALL Spc_Create( am_I_Root     = am_I_Root,                    &
                             ThisSpc       = SpcData(N)%Info,              &
                             ModelID       = N,                            &
                             KppSpcId      = KppSpcId(N),                  &
                             KppVarId      = KppVarId(N),                  &
                             KppFixId      = KppFixId(N),                  &
                             Name          = NameAllCaps,                  &
                             FullName      = 'Hydrophilic marine OC',      &
                             Formula       = ,                           &
                             MW_g          = 12.01_fp,                     &
                             EmMW_g        = 12.0_fp,                      &
                             MolecRatio    = 1.0_fp,                       &
                             Is_Advected   = Is_Advected,                  &
                             Is_Gas        = T,                            &
                             Is_Drydep     = F,                            &
                             Is_Wetdep     = F,                            &
                             DD_DvzAerSnow = 0.03_fp,                      &
                             DD_Hstar_old  = 0.0_fp,                       &
                             WD_AerScavEff = 1.0_fp,                       &
                             WD_KcScaleFac = KcScale,                      &
                             WD_RainoutEff = RainEff,                      &
                             RC            = RC )

         CASE( 'MOPO' )

            ! Turn off rainout because MOPO is hydrophobic
            KcScale = (/ 1.0_fp, 1.0_fp, 1.0_fp /)
            RainEff = (/ 0.0_fp, 0.0_fp, 0.0_fp /)

            CALL Spc_Create( am_I_Root     = am_I_Root,                    &
                             ThisSpc       = SpcData(N)%Info,              &
                             ModelID       = N,                            &
                             KppSpcId      = KppSpcId(N),                  &
                             KppVarId      = KppVarId(N),                  &
                             KppFixId      = KppFixId(N),                  &
                             Name          = NameAllCaps,                  &
                             FullName      = 'Hydrophobic marine OC',      &
                             Formula       = ,                           &
                             MW_g          = 12.01_fp,                     &
                             EmMW_g        = 12.0_fp,                      &
                             MolecRatio    = 1.0_fp,                       &
                             Is_Advected   = Is_Advected,                  &
                             Is_Gas        = T,                            &
                             Is_Drydep     = F,                            &
                             Is_Wetdep     = F,                            &
                             DD_DvzAerSnow = 0.03_fp,                      &
                             DD_Hstar_old  = 0.0_fp,                       &
                             WD_AerScavEff = 0.0_fp,                       &
                             WD_RainOutEff = RainEff,                      &
                             RC            = RC )

to: 

         CASE( 'MOPI', 'MOPO' )
            
            ! MOPO is treated the same as OCPO and
            ! MOPI is treated the same as OCPI (msj, krt, 8/23/17)

            ! These have mostly identical properties
            ! Turn off rainout for hydrophobic OC, for all temperatures.
            SELECT CASE( NameAllCaps )

               CASE( 'MOPI' )
                  FullName = 'Hydrophilic marine organic carbon aerosol'

                  ! Halve the Kc (cloud condensate -> precip) rate
                  ! for the temperature range 237 K <= T < 258 K.
                  KcScale  = (/ 1.0_fp, 0.5_fp, 1.0_fp /)

                  ! Turn off rainout only when 237 K <= T < 258K.
                  RainEff  = (/ 1.0_fp, 0.0_fp, 1.0_fp /)

               CASE( 'MOPO' )
                  Fullname = 'Hydrophobic marine organic carbon aerosol'

                  ! For all temperatures:
                  ! (1) Halve the Kc (cloud condensate -> precip) rate
                  ! (2) Turn off rainout (OCPO is hydrophobic)
                  KcScale  = (/ 0.5_fp, 0.5_fp, 0.5_fp /)
                  RainEff  = (/ 0.0_fp, 0.0_fp, 0.0_fp /)

            END SELECT

            CALL Spc_Create( am_I_Root     = am_I_Root,                    &
                             ThisSpc       = SpcData(N)%Info,              &
                             ModelID       = N,                            &
                             KppSpcId      = KppSpcId(N),                  &
                             KppVarId      = KppVarId(N),                  &
                             KppFixId      = KppFixId(N),                  &
                             Name          = NameAllCaps,                  &
                             FullName      = FullName,                     &
                             Formula       = ,                           &
                             MW_g          = 12.01_fp,                     &
                             EmMW_g        = 12.0_fp,                      &
                             Is_Advected   = Is_Advected,                  &
                             Is_Gas        = F,                            &
                             Is_Drydep     = T,                            &
                             Is_Wetdep     = T,                            &
                             Density       = 1300.0_fp,                    &
                             DD_DvzAerSnow = 0.03_fp,                      &
                             DD_F0         = 0.0_fp,                       &
                             DD_Hstar_Old  = 0.0_fp,                       &
                             WD_AerScavEff = 1.0_fp,                       &
                             WD_KcScaleFac = KcScale,                      &
                             WD_RainoutEff = RainEff,                      &
                             RC            = RC )

--Melissa Sulprizio (talk) 15:59, 23 August 2017 (UTC)

Retrieved from "https://wiki.seas.harvard.edu/geos-chem/index.php?title=Aerosol_emissions&oldid=42070"