Difference between revisions of "Aerosol emissions"

From Geos-chem
Jump to: navigation, search
(Bug fix: Allocate the OCCONV array for marine POA simulations)
(Bug fix: Allocate the OCCONV array for marine POA simulations)
Line 88: Line 88:
 
<span style="color:darkorange">'''''This fix is slated for inclusion in [[GEOS-Chem v11-01#v11-01j|v11-01j]].'''''</span>
 
<span style="color:darkorange">'''''This fix is slated for inclusion in [[GEOS-Chem v11-01#v11-01j|v11-01j]].'''''</span>
  
While unit testing GEOS-Chem with the [[GNU Fortran compiler]], we discovered an error in the marine POA simulation.  The array <code>OCCONV</code> should have been allocated in routine <code>INIT_SEASALT</code> (in module <code>GeosCore/seasalt_mod.F</code>) but wasn't.  We have now added the following IF statement in green, below the existing code:
+
While unit testing GEOS-Chem with the [[GNU Fortran compiler]], we discovered an error in the marine POA simulation.  The array <code>OCCONV</code> should have been allocated in routine <code>INIT_SEASALT</code> (in module <code>GeosCore/seasalt_mod.F</code>) but wasn't.  We have now added the following IF statement in <span style="color:green">GREEN</span>, below the existing code:
  
 
         IF ( AS /= 0 ) CALL ALLOC_ERR( 'DMID' )
 
         IF ( AS /= 0 ) CALL ALLOC_ERR( 'DMID' )

Revision as of 18:44, 13 October 2016

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

NOTE: In GEOS-Chem v10-01 and higher versions, all emissions are handled by the HEMCO emissions component.

Anthropogenic, biofuel, and natural source emissions

The various sectors of aerosol emissions are now handled by the 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:

  1. Aerosol emissions available via HEMCO
  2. Anthropogenic and biofuel emissions available via HEMCO
  3. Aircraft and ship emissions available via HEMCO

--Bob Y. (talk) 20:19, 26 October 2015 (UTC)

Biomass burning emissions

The GFED4 inventory can be used to compute biomass burning emisisons for the following species:

  1. SO2
  2. NH3
  3. EC = elemental carbon (aka "black carbon")
  4. 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):

  1. GFED2
  2. GFED3
  3. Duncan et al [2003] seasonal and interannual biomass emissions (cf Bey et al [2001])

--Bob Y. (talk) 20:23, 26 October 2015 (UTC)

DMS emissions

In GEOS-Chem v10-01 and higher versions, the sea-air flux of DMS is now handled by the HEMCO seaflux extension.

From Park et al [2004]:

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.

EC and OC emissions

The following emissions inventories for elemental carbon (EC) and organic carbon (OC) are used in GEOS-Chem.

  1. Bond et al [2007]: EC and OC emissions over North America
  2. Spracklen et al: Oceanic OC emissions (NOTE: to be included post GEOS-Chem v8-03-01)

Mineral dust emissions

The mineral dust emissions algorithm is described on a separate wiki page.

Sea salt emissions

Please see this discussion about recent updates to the sea salt aerosol emissions.

Ship emissions

Please see our 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.

--Bob Y. 13:30, 19 March 2010 (EDT)

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 is slated for inclusion in v11-01j.

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)