Difference between revisions of "EMEP European anthropogenic emissions"

From Geos-chem
Jump to: navigation, search
(In GEOS-Chem v10-01 and newer versions)
 
(59 intermediate revisions by 3 users not shown)
Line 1: Line 1:
 +
----
 +
<span style="color:red"><big><strong>EMEP emissions have been superseded by [[CEDS_anthropogenic_emissions|CEDS global emissions]].</strong></big></span>
 +
----
 +
 +
 
On this page we provide information about ''The Co-operative Programme for Monitoring and Evaluation of the Long-range Transmission of Air Pollutants in Europe (EMEP)'' [[Anthropogenic emissions|anthropogenic emissions inventory]] for Europe.
 
On this page we provide information about ''The Co-operative Programme for Monitoring and Evaluation of the Long-range Transmission of Air Pollutants in Europe (EMEP)'' [[Anthropogenic emissions|anthropogenic emissions inventory]] for Europe.
  
 
== Overview ==
 
== Overview ==
  
The EMEP inventory also includes [[#EMEP European ship emissions|emissions from ship exhaust]].
+
The EMEP inventory includes annual anthropogenic emissions of CO, NH3, NOx, and SOx for 1990-2005. This inventory also includes [[#EMEP European ship emissions|emissions of CO, NH3, NOx, and SOx from ship exhaust]].
  
NOTE: The EMEP emissions data have been extended to the year 2007 in [[GEOS-Chem v8-03-01]] and higher versions.
+
The data files were obtained from the [http://www.ceip.at/ EMEP website] and processed by Aaron van Donkelaar.
  
=== Update EMEP emissions for 2008-2010 ===
+
In [[GEOS-Chem v10-01]] and newer versions, the EMEP anthropogenic and ship emissions are read with the [[HEMCO|HEMCO emissions component]].  We have created [http://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/EMEP/v2015-03/README new EMEP data files] (in [[Preparing_data_files_for_use_with_HEMCO|COARDS-compliant netCDF format]]) for use with HEMCO.  These new data files are contained in the [[HEMCO data directories|HEMCO data directory tree]].  For detailed instructions on how to download these data files to your disk server, please see our [[HEMCO_data_directories#Downloading_the_HEMCO_data_directories|''Downloading the HEMCO data directories'' wiki post]].
  
'''''This update is slated for inclusion in [[GEOS-Chem v10-01]]'''''
+
--[[User:Bmy|Bob Y.]] 13:13, 3 March 2015 (EST)
  
Aaron van Donkelaar generated EMEP emissions for 2008, 2009, and 2010. The data files are available at:
+
=== Discussion following 1-month benchmark v10-01h ===
  
ftp://rain.ucis.dal.ca/ctm/GEOS_1x1/EMEP_200911/
+
The [[#Update_EMEP_emissions_for_2008-2010|2008-2010 updated EMEP emissions]] were validated in [[GEOS-Chem_v10-01_benchmark_history#v10-01h|GEOS-Chem v10-01h]]. The following discussion ensued.
+
ftp://rain.ucis.dal.ca/ctm/GEOS_0.5x0.666_EU/EMEP_200911/
+
  
[mailto:Aaron.van.Donkelaar@dal.ca Aaron van Donkelaar] wrote:
+
'''''Aaron van Donkelaar wrote:'''''
  
:You will have to change the maximum allowable value of EMEP_YEAR in emep_mod.f from 2007 to 2010.  This variable is set twice, on line 467 and 647, in v8-03-01, but these may have moved slightly in later versions. The first setting is for the global simulation and the second for the nested EU simulation. The lines you are looking for look like:
+
:I'm a little surprised at the drop of more than 50% in anthropogenic SO2 and NO over Europe.  
  
      EMEP_YEAR = MIN( SCALEYEAR, 2007 )
+
:I decided to explore the large drop in emissions and have found the drop in the revised EMEP emissions that are presently available to be much lower in than what EMEP had estimated when I downloaded/reformatted the original dataset back in 2011.
  
:Change them to:
+
:It should be fairly straightforward for me to reprocess new EMEP input files based on their latest estimates, which run until 2012, if the old format is still useful post-HEMCO. 
  
      EMEP_YEAR = MIN( SCALEYEAR, 2010 )
+
:For now I would not recommend using the EMEP files originally slated as part of this update.  Rather, I would suggest that I/we update the emission based on the latest EMEP inventory.
  
--[[User:Melissa Payer|Melissa Sulprizio]] 15:18, 23 January 2014 (EST)
+
<span style="color:green">'''Update 3/31/15: Aaron has reprocessed EMEP emissions for 1990-2012. Please see ''[[#Reprocess EMEP emissions for 1990-2012|this wiki post]]'' for more information.'''</span>
  
== EMEP European ship emissions inventory ==
+
--[[User:Melissa Payer|Melissa Sulprizio]] 16:40, 27 March 2015 (EDT)
  
=== In GEOS-Chem v10-01 and newer versions ===
+
=== Reprocess EMEP emissions for 1990-2012 ===
  
In [[GEOS-Chem v10-01]] and newer versions, the EMEP ship emissions are read with the [[HEMCO|HEMCO emissions component]].  We have created [http://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/EMEP/v2015-01/README new EMEP ship data files] (in [[Preparing_data_files_for_use_with_HEMCO|COARDS-compliant netCDF format]]) for use with HEMCO. These new data files are contained in the [[HEMCO data directories|HEMCO data directory tree]].  For detailed instructions on how to download these data files to your disk server, please see our [[HEMCO_data_directories#Downloading_the_HEMCO_data_directories|''Downloading the HEMCO data directories'' wiki post]].
+
<span style="color:green">'''''This update was validated with the 1-month benchmark simulation [[GEOS-Chem v10-01 benchmark history#v10-01i|v10-01i]] and approved on Approved 01 May 2015.'''''</span>
  
--[[User:Bmy|Bob Y.]] 13:13, 3 March 2015 (EST)
+
The original updated EMEP emissions for 2008-2010 (generated in 2012) were found to be too low in [[GEOS-Chem_v10-01_benchmark_history#v10-01h|1-month benchmark v10-01h]], so Aaron van Donkelaar reprocessed EMEP anthropogenic and ship emissions of CO, NH3, SO2, and NO for 1990-2012. The [http://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/EMEP/v2015-03/README reprocessed data file] is in netCDF format and includes data that has been mapped from the original 50 km grid x 50 km grid onto a 0.1&deg; x 0.1&deg; grid. Regridding to the finer 0.1° x 0.1° resolution allowed these input files to maintain a reasonable representation of their original grid characteristics, thereby permitting HEMCO to more accurately regrid onto a variety of GEOS-Chem grids. This file is read with the [[HEMCO|HEMCO emissions component]] and will be validated with [[GEOS-Chem_v10-01_benchmark_history#v10-01h-Run0|1-year benchmark v10-01h-Run0]].
  
=== In GEOS-Chem versions prior to v10-01 ===
+
--[[User:Melissa Payer|Melissa Sulprizio]] 14:08, 31 March 2015 (EDT)
  
The EMEP inventory contains ship emissions for NOx, CO, SO2. 
+
== EMEP European ship emissions inventory ==
  
The source code is located in file <tt>GeosCore/emep_mod.f</tt>.   
+
The EMEP inventory contains ship emissions for CO, NH3, SOx, and NOx.   
  
Please also note the following fixes:
+
Please also note the following fixes that are relevant for [[GEOS-Chem v9-01-01]] or older versions:
  
#[[#Minor fixes for EMEP ship emissions|Minor fixes for EMEP ship emissions]]
+
#[[#Bug with EMEP ship emissions for nested grids|Bug with EMEP ship emissions for nested grids]]
 +
#[[#Bug fix for EMEP ship emissions|Bug fix for EMEP ship emissions]]
 +
#[[#Fix for initialization of EMEP ship emissions|Fix for initialization of EMEP ship emissions]]
 
#[[#Inconsistent mask file for EMEP ship emissions|Inconsistent mask file for EMEP ship emissions]]
 
#[[#Inconsistent mask file for EMEP ship emissions|Inconsistent mask file for EMEP ship emissions]]
  
==== 1980 to 1989 ====
+
=== 1980 to 1989 ===
  
 
EMEP ship data from 1980-1989 emit NOx as pure NOx (and not HNO3 + 10*O3).  It was not possible to separate the NOx into HNO3 and O3 from this data set.
 
EMEP ship data from 1980-1989 emit NOx as pure NOx (and not HNO3 + 10*O3).  It was not possible to separate the NOx into HNO3 and O3 from this data set.
Line 54: Line 59:
 
See ''Auvray and Bey'' [2005] and [http://www.emep.int/publ/reports/2002/mscw_note_1_2002.pdf ''Vestreng and Klein'' <nowiki>[2002]</nowiki>] for more information.
 
See ''Auvray and Bey'' [2005] and [http://www.emep.int/publ/reports/2002/mscw_note_1_2002.pdf ''Vestreng and Klein'' <nowiki>[2002]</nowiki>] for more information.
  
==== 1990 to 2005 ====
+
=== 1990 to 2005 ===
  
EMEP ship data from 1990-2005 are based on [''Vestreng et al'', 2007].
+
EMEP ship data from 1990-2005 are based on [http://emep.int/publ/reports/2007/emep_technical_1_2007.pdf ''Vestreng et al.'' <nowiki>[2007]</nowiki>].
 
+
For more information about the data, please see the README file: [http://ftp.as.harvard.edu/gcgrid/data/GEOS_1x1/EMEP_200911/README GEOS_1x1/EMEP_200911/README]
+
  
 
--[[User:Bmy|Bob Y.]] 12:00, 20 December 2012 (EST)
 
--[[User:Bmy|Bob Y.]] 12:00, 20 December 2012 (EST)
  
== References ==
+
== Known issues ==
  
#Auvray, M., and I. Bey, ''Long-Range Transport to Europe: Seasonal Variations and Implications for the European Ozone Budget'', <u>J. Geophys. Res.</u>, '''110''', D11303, doi: 10.1029/2004JD005503, 2005.
+
=== Zero EMEP emissions at 30N ===
# Vestreng, V., and H. Klein (2002), ''Emission data reported to UNECE/EMEP: Quality assurance and trend analysis & Presentation of WebDab'', <u>MSC-W Status Report 2002</u>, Norwegian Meteorological Institute, Oslo Norway, July 2002. [http://www.emep.int/publ/reports/2002/mscw_note_1_2002.pdf PDF]
+
#Vestreng et al., ''EMEP 2005 expert emissions'', 2007. (need full citation!)
+
  
--[[User:Bmy|Bob Y.]] 12:13, 20 December 2012 (EST)
+
<span style="color:green">'''''This update was validated with [[GEOS-Chem_v11-01_benchmark_history#v11-01f|1-month benchmark simulation v11-01f]] and [[GEOS-Chem_v11-01_benchmark_history#v11-01f-geosfp-Run0|1-year benchmark simulation v11-01f-geosfp-Run0]]. This version was approved on 16 Apr 2016.'''''</span>
  
== Previous issues that are now resolved ==
+
'''''[[User:Emarais|Eloise Marais]] wrote:'''''
  
=== Bug with EMEP ship emissions for nested grids ===
+
<blockquote>
 +
Anthropogenic emissions are zero where the EMEP and EDGAR emissions coincide at 30 degrees North (using [[GEOS-Chem v10-01]]). This is because the EMEP mask extends to gridsquares where EMEP emissions are zero. Anthropogenic emissions are relatively low there, so this issue only shows up when NOx emissions are plotted on a log scale.
  
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>This code was removed from [[GEOS-Chem v10-01]] and newer versions.  EMEP ship emissions are now implemented through the [[HEMCO|HEMCO emissions component]].</p><p>&mdash; Bob Yantosca, 13 February 2015</p>'''</div>
+
My quick fix is to shift the minimum latitude of the EMEP mask up by 1 degree latitude (so starting at 31 deg N instead of 30 deg N), in an updated mask file <tt>EMEP_mask.geos.1x1.20151222.nc</tt>.  (The boxes that were removed show up as <span style="color:#AA0000">'''RED'''</span> in the plot below.)
  
'''''NOTE: This issue was resolved in [[GEOS-Chem v8-02-01]].'''''
+
[[Image:Boxes_removed_from_EMEP_mask.png]]
  
'''''Win Trivitayanurak wrote:'''''
+
</blockquote>
  
:This bug is in subroutine EMFOSSIL in emfossil.f.  The use of function GET_EUROPE_MASK when LEMEP = .false. crashes the run because, i think, the function is not really valid when EMEP emission is not turned on.  When I ran the China nested-grid, all other emission choices need to be turned off, e.g. EMEP, BRAVO. Only STREETS can be on for the China nested grid run.  Therefore, w/o this fix, the run won't go anywhere.
+
As of 06 Jan 2016, the [[GEOS-Chem Support Team]] has placed the updated mask file into the following HEMCO data directory:
  
:I copy the section of the code and paste it here below. See the change in the line with !Prior to 5/1/09 (win, 5/1/09).
+
  HEMCO/MASKS/v2014-07/EMEP_mask.geos.1x1.20151222.nc
  
  !Prior to 5/1/09 (win, 5/1/09)
+
We will validate Eloise's update with the 1-month benchmark simulation [[GEOS-Chem v11-01#v11-01f|GEOS-Chem v11-01f]].  This will involve changing the EMEP mask filename in the proper <tt>HEMCO_Config.rc</tt> files.
  !              ELSE IF ( LEMEPSHIP  .AND.
+
  !    $                ( GET_EUROPE_MASK( I, J ) > 0d0 ) ) THEN
+
  !
+
  !                  ! Get SHIP EMEP emissions for NOx [molec/cm2/s]
+
  !                  SHIP = GET_EMEP_ANTHRO( I, J, NN, SHIP=.TRUE.)
+
   
+
              ELSE IF ( LEMEPSHIP ) THEN
+
               
+
                  IF ( GET_EUROPE_MASK( I, J ) > 0d0 )
+
                  ! Get SHIP EMEP emissions for NOx [molec/cm2/s]
+
    &            SHIP = GET_EMEP_ANTHRO( I, J, NN, SHIP=.TRUE.)
+
  
 +
NOTE: The old mask file (<tt>EMEP_mask.geos.1x1.nc</tt>) shall be preserved for backwards compatibility. 
  
 +
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 15:23, 30 March 2016 (UTC)
  
--[[User:Ccarouge|Ccarouge]] 15:50, 7 May 2009 (EDT)<br>
+
==== Observed differences in CO emissions attributed to new EMEP mask file ====
--[[User:Bmy|Bob Y.]] 11:52, 20 December 2012 (EST)
+
  
=== Bug fix for EMEP ship emissions ===
+
Note that simulations using the updated EMEP mask file may cause some emissions differences w/r/t simulations using the old EMEP mask:
  
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>This code was removed from [[GEOS-Chem v10-01]] and newer versions.  EMEP ship emissions are now implemented through the [[HEMCO|HEMCO emissions component]].</p><p>&mdash; Bob Yantosca, 13 February 2015</p>'''</div>
+
'''''[[User:psk9|Prasad Kasibhatla]] wrote:'''''
  
'''''NOTE: This issue was resolved in [[GEOS-Chem v8-02-04]].'''''
+
<blockquote>
 +
I am surprised by the differences in CO emissions between [[GEOS-Chem v11-01#v11-01d|v11-01d]] and [[GEOS-Chem v11-01#v11-01f|v11-01f]]. I was struck by how large they seem to be. What first caught my eye was the difference in boundary-layer CO mixing ratios in Cairo - and it looks like it is related to a big change in anthro+bifuel emissions in one grid box....Is there a simple explanation for this?
 +
</blockquote>
  
'''''Monika Kopacz wrote:'''''
+
'''''[[User:Melissa Payer|Melissa Sulprizio]] replied:'''''
  
:I noticed that the code crashes when EMEP is set to false but ship emissions from EMEP are true in the input file. Sure, it makes sense, but I would suggest the following addition to the <tt>input_mod.f</tt>, around say around line 1550 to prevent the code from crashing.
+
<blockquote>
 +
Differences in the anthropogenic+biofuel CO emissions in Cairo were caused by the updated EMEP mask file introduced in v11-01f. Eloise Marais found that EDGAR emissions were zero at 30N because the EMEP mask extended too far south, so she provided a new file that shifts the EMEP mask north by 1 degree.
 +
</blockquote>
  
    ! Make sure we're not using EMEP ship emissions, if EMEP is off
+
--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 16:01, 15 April 2016 (UTC)
    IF (.not. LEMEP ) THEN
+
      LEMEPSHIP = .FALSE.
+
    ENDIF
+
  
'''''[[User:Bmy|Bob Yantosca]] wrote:'''''
+
== References ==
  
:Good catch. We've put this code into READ_EMISSIONS_MENU in <tt>input_mod.f:</tt>
+
#EMEP website: http://www.emep.int/index.html
 +
#Auvray, M., and I. Bey, ''Long-Range Transport to Europe: Seasonal Variations and Implications for the European Ozone Budget'', <u>J. Geophys. Res.</u>, '''110''', D11303, doi: 10.1029/2004JD005503, 2005.
 +
# Vestreng, V., and H. Klein, ''Emission data reported to UNECE/EMEP: Quality assurance and trend analysis & Presentation of WebDab'', <u>MSC-W Status Report 2002</u>, Norwegian Meteorological Institute, Oslo Norway, July 2002. [http://www.emep.int/publ/reports/2002/mscw_note_1_2002.pdf (PDF)]
 +
#Vestreng, V., Mareckova, K., Kakareka, S., Malchykhina, A., and Kukharchyk, T., ''Inventory Review 2007; Emission Data Reported to LRTAP Convention and NEC Directive'', <u>MSC-W Technical Report 1/07</u>, The Norwegian Meteorological Institute, Oslo, Norway, 2007. [http://emep.int/publ/reports/2007/emep_technical_1_2007.pdf (PDF)]
  
    !%%% Bug fix!  If LEMEPSHIP is turned on but LEMEP is turned %%%
+
--[[User:Bmy|Bob Y.]] 12:13, 20 December 2012 (EST)
    !%%% off, this will cause an error (because arrays are not  %%%
+
    !%%% allocated, etc.).  For now, just turn off LEMEPSHIP    %%%
+
    !%%% and print a warning message.  Whoever wants to fix this %%%
+
    !%%% in a more robust way is welcome to do so.              %%%
+
    !%%% (mak, bmy, 10/19/09)                                    %%%
+
    IF ( LEMEPSHIP .and. ( .not. LEMEP ) ) THEN
+
        LEMEPSHIP = .FALSE.
+
        WRITE( 6, '(a)' ) REPEAT( '=', 79 )
+
        WRITE( 6, '(a)' ) 'WARNING! EMEP emissions are turned off,'
+
        WRITE( 6, '(a)' ) 'so also turn off EMEP ship emissions'
+
        WRITE( 6, '(a)' ) 'in order to avoid crashing GEOS-Chem!'
+
        WRITE( 6, '(a)' ) REPEAT( '=', 79 )
+
    ENDIF
+
 
+
:As you pointed out the simple solution is to just turn off <tt>LEMEPSHIP</tt> if <tt>LEMEP=F</tt>.  However, it may still be possible to keep <tt>LEMEPSHIP</tt> turned on...one would have to go thru  <tt>emep_mod.f</tt> to make sure that the ship emissions arrays get allocated even if the other arrays don't.  Whoever wants to can feel free to do that...
+
 
+
--[[User:Bmy|Bob Y.]] 10:57, 19 October 2009 (EDT)
+
 
+
=== Fix for initialization of EMEP ship emissions ===
+
 
+
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>This code was removed from [[GEOS-Chem v10-01]] and newer versions.  EMEP ship emissions are now implemented through the [[HEMCO|HEMCO emissions component]].</p><p>&mdash; Bob Yantosca, 13 February 2015</p>'''</div>
+
 
+
'''''NOTE: This issue was resolved in [[GEOS-Chem v8-03-01]].'''''
+
 
+
'''''[mailto:lechlm@leeds.ac.uk Helen Macintyre] wrote:'''''
+
 
+
:I am trying to do some simulations with 1985 emissions, but I find that I run into an SST error.
+
 
+
:Looking through the log file reveals this:
+
 
+
    ===============================================================================
+
    M O N T H L Y  E M E P  E U R O P E A N    E M I S S I O N S
+
    ( INCL. SHIP )
+
+
    Base Year :1985
+
+
    EMEP anthropogenic NOx  for 1985/07:    -Infinity Tg N
+
    EMEP anthropogenic CO  for 1985/07:          NaN Tg
+
    EMEP anthropogenic SO2  for 1985/07:      2.217651 Tg S
+
    EMEP anthropogenic NH3  for 1985/07:      0.771631 Tg
+
 
+
'''''[mailto:ccarouge@seas.harvard.edu Claire Carouge] wrote:'''''
+
 
+
:It looks like a bug in initialisation of EMEP ship emissions for CO and NOx.
+
 
+
:In INIT_EMEP (in emep_mod.f) there are the lines:
+
 
+
      ALLOCATE( EMEP_SO2_SHIP( IIPAR, JJPAR ), STAT=AS )
+
      IF ( AS /= 0 ) CALL ALLOC_ERR( 'EMEP_SO2_SHIP' )
+
      EMEP_SO2_SHIP = 0d0
+
     
+
      ALLOCATE( EMEP_CO_SHIP( IIPAR, JJPAR ), STAT=AS )
+
      IF ( AS /= 0 ) CALL ALLOC_ERR( 'EMEP_SO2_SHIP' )
+
      EMEP_SO2_SHIP = 0d0
+
     
+
      ALLOCATE( EMEP_NOx_SHIP( IIPAR, JJPAR ), STAT=AS )
+
      IF ( AS /= 0 ) CALL ALLOC_ERR( 'EMEP_SO2_SHIP' )
+
      EMEP_SO2_SHIP = 0d0
+
 
+
:So after allocating EMEP_CO_SHIP or EMEP_NOx_SHIP we initialize EMEP_SO2_SHIP to 0d0. A typical problem of copy/paste.
+
 
+
:You should change the lines to initialize the appropriate ship arrays to 0d0:
+
 
+
      ALLOCATE( EMEP_SO2_SHIP( IIPAR, JJPAR ), STAT=AS )
+
      IF ( AS /= 0 ) CALL ALLOC_ERR( 'EMEP_SO2_SHIP' )
+
      EMEP_SO2_SHIP = 0d0
+
     
+
      ALLOCATE( EMEP_CO_SHIP( IIPAR, JJPAR ), STAT=AS )
+
      IF ( AS /= 0 ) CALL ALLOC_ERR( 'EMEP_CO_SHIP' )
+
      EMEP_CO_SHIP = 0d0
+
     
+
      ALLOCATE( EMEP_NOx_SHIP( IIPAR, JJPAR ), STAT=AS )
+
      IF ( AS /= 0 ) CALL ALLOC_ERR( 'EMEP_NOx_SHIP' )
+
      EMEP_NOx_SHIP = 0d0
+
 
+
--[[User:Ccarouge|Ccarouge]] 12:54, 26 January 2010 (EST)<br>
+
--[[User:Bmy|Bob Y.]] 12:11, 20 December 2012 (EST)
+
 
+
=== Inconsistent mask file for EMEP ship emissions ===
+
 
+
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>In [[GEOS-Chem v10-01]] and newer versions, EMEP ship emissions are now implemented through the [[HEMCO|HEMCO emissions component]].  This fix was implemented into HEMCO.</p><p>&mdash; Bob Yantosca, 13 February 2015</p>'''</div>
+
 
+
'''''NOTE: This issue was resolved in [[GEOS-Chem v9-01-01]].'''''
+
 
+
'''''[mailto:g.c.m.vinken@tue.nl Geert Vinken] wrote:'''''
+
 
+
:When I was looking at the ship emissions I noticed that the current mask that determines where to apply EMEP emissions, does not correspond anymore with the actual EMEP emissions. The mask still corresponds with the 200510 emissions data. This mismatch in the mask and the actual emissions results in no ship emissions being released over a region of the North Atlantic and north of Norway.
+
 
+
:As a solution I updated the mask file to correspond with the new border of the EMEP emissions domain. I cut of the mask west of 30W and north of 70N. This fixes the problem.
+
 
+
--[[User:Bmy|Bob Y.]] 12:20, 20 December 2012 (EST)
+
 
+
=== Fixes in emep_mod.f ===
+
 
+
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>This code was removed from [[GEOS-Chem v10-01]] and newer versions.  EMEP ship emissions are now implemented through the [[HEMCO|HEMCO emissions component]].</p><p>&mdash; Bob Yantosca, 13 February 2015</p>'''</div>
+
 
+
'''''NOTE: This issue was ultimately resolved in [[GEOS-Chem_v9-01-03_benchmark_history#v9-01-03g|GEOS-Chem v9-01-03g]]'''''.
+
 
+
In routines <tt>READ_EMEP_UPDATED</tt> and <tt>READ_EMEP_UPDATED_05x0666</tt> of <tt>emep_mod.f</tt>, there were some IF statements that were used to pick the proper filename and tracer number with which to read data from the bpch file format:
+
 
+
        IF ( TRACER .eq. IDTNOx ) THEN
+
            FILENAME = TRIM( DIR ) // 'SeasonalVariation/'
+
    &        // 'NOx-EMEP-SeasonalScalar.geos.1x1'
+
            RATIOID = 71
+
        ELSEIF ( TRACER .eq. IDTCO ) THEN
+
            FILENAME = TRIM( DIR ) // 'SeasonalVariation/'
+
    &        // 'CO-EMEP-SeasonalScalar.geos.1x1'
+
            RATIOID = 72
+
        ELSEIF ( TRACER .eq. IDTSO2 ) THEN
+
            FILENAME = TRIM( DIR ) // 'SeasonalVariation/'
+
    &        // 'SOx-EMEP-SeasonalScalar.geos.1x1'
+
            RATIOID = 73
+
        ELSEIF ( TRACER .eq. IDTNH3 ) THEN
+
            FILENAME = TRIM( DIR ) // 'SeasonalVariation/'
+
    &        // 'NH3-EMEP-SeasonalScalar.geos.1x1'
+
            RATIOID = 74
+
        ENDIF
+
 
+
where TRACER is passed to these routines via subroutine calls such as these:
+
 
+
        CALL READ_EMEP_UPDATED(  1, EMEP_YEAR, EMEP_NOx, 0 )
+
        CALL READ_EMEP_UPDATED(  4, EMEP_YEAR, EMEP_CO,  0 )
+
        CALL READ_EMEP_UPDATED( 26, EMEP_YEAR, EMEP_SO2, 0 )
+
        CALL READ_EMEP_UPDATED( 30, EMEP_YEAR, EMEP_NH3, 1 )
+
 
+
The problem is the code is set up to assume that <tt>IDTNOX</tt> will always be 1, <tt>IDTCO</tt> will always be 4, <tt>IDTSO2</tt> will always be 26, and <tt>IDTNH3</tt> will always be 30.  This is true only if you are running a full-chemistry simulation.  However, for an offline-aerosol simulation NOx and CO are undefined, and SO2 and NH3 will not equal 26 and 30.  Therefore, the IF statement will not return the proper file name and tracer ID.
+
 
+
The solution is to remove all references to <tt>IDTNOX</tt>, <tt>IDTCO</tt>, <tt>IDTSO2</tt>, <tt>IDTNH3</tt> from that section of code.  The new code looks like:
+
 
+
        ! Pick proper filename & tracer #
+
        SELECT CASE( TRACER )
+
+
            ! NOx
+
            CASE( 1 )
+
              RATIOID  = 71
+
              FILENAME = TRIM( DIR )                        //
+
    &                    'SeasonalVariation/'              //
+
    &                    'NOx-EMEP-SeasonalScalar.geos.1x1'
+
            ! CO
+
            CASE( 4 )
+
              RATIOID  = 72
+
              FILENAME = TRIM( DIR )                        //
+
    &                    'SeasonalVariation/'              //
+
    &                    'CO-EMEP-SeasonalScalar.geos.1x1'
+
            ! SO2
+
            CASE( 26 )
+
              RATIOID  = 73
+
              FILENAME = TRIM( DIR )                        //
+
    &                    'SeasonalVariation/'              //
+
    &                    'SOx-EMEP-SeasonalScalar.geos.1x1'
+
            ! NH3
+
            CASE( 30 )
+
              RATIOID  = 74
+
              FILENAME = TRIM( DIR )                        //
+
    &                    'SeasonalVariation/'              //
+
    &                    'NH3-EMEP-SeasonalScalar.geos.1x1'
+
        END SELECT
+
 
+
There are a couple more instances of this error in <tt>emep_mod.f</tt>, but the fix is the same.
+
 
+
'''''Update 8/1/11:''''' Implementing this change causes the 1-month standard benchmark simulation to have divergent results than the previous version.  Therefore, we should understand the source of these differences before we add this update into
+
the standard repository.
+
 
+
--[[User:Bmy|Bob Y.]] 10:43, 1 August 2011 (EDT)
+
 
+
'''''Update 11/28/11:''''' Patrick Kim tested this and found that the numerical discrepancy was being caused by ISORROPIA.  Therefore we will add this fix into [[GEOS-Chem v9-01-03]].
+
 
+
--[[User:Bmy|Bob Y.]] 10:19, 28 November 2011 (EST)
+
 
+
'''''Update 2/27/12:''''' This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-03_benchmark_history#v9-01-03g|v9-01-03g]] and approved on 27 Feb 2012.
+
 
+
--[[User:Melissa Payer|Melissa Payer]] 17:47, 27 February 2012 (EST)<br>
+
--[[User:Bmy|Bob Y.]] 12:43, 20 December 2012 (EST)
+
 
+
=== Bug fixes for regional emission masks ===
+
 
+
<div style="color: #aa0000; background: #eeeeee;border: 3px solid red; padding: 1em; margin: auto; width: 90%; ">'''<p>This code was removed from [[GEOS-Chem v10-01]] and newer versions.  EMEP emissions are now implemented through the [[HEMCO|HEMCO emissions component]].</p><p>&mdash; Bob Yantosca, 13 February 2015</p>'''</div>
+
 
+
'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem v9-01-03 benchmark history#v9-01-03q|v9-01-03q]].'''''
+
 
+
Please see [[Regridding_in_GEOS-Chem#Bug_fixes_for_regional_emission_masks|this wiki post on our ''Regridding in GEOS-Chem'' wiki page]] for a full description of an issue that affected the regridding of the EMEP European mask.
+
 
+
--[[User:Bmy|Bob Y.]] 15:15, 20 December 2012 (EST)
+

Latest revision as of 20:43, 16 September 2022


EMEP emissions have been superseded by CEDS global emissions.



On this page we provide information about The Co-operative Programme for Monitoring and Evaluation of the Long-range Transmission of Air Pollutants in Europe (EMEP) anthropogenic emissions inventory for Europe.

Overview

The EMEP inventory includes annual anthropogenic emissions of CO, NH3, NOx, and SOx for 1990-2005. This inventory also includes emissions of CO, NH3, NOx, and SOx from ship exhaust.

The data files were obtained from the EMEP website and processed by Aaron van Donkelaar.

In GEOS-Chem v10-01 and newer versions, the EMEP anthropogenic and ship emissions are read with the HEMCO emissions component. We have created new EMEP data files (in COARDS-compliant netCDF format) for use with HEMCO. These new data files are contained in the HEMCO data directory tree. For detailed instructions on how to download these data files to your disk server, please see our Downloading the HEMCO data directories wiki post.

--Bob Y. 13:13, 3 March 2015 (EST)

Discussion following 1-month benchmark v10-01h

The 2008-2010 updated EMEP emissions were validated in GEOS-Chem v10-01h. The following discussion ensued.

Aaron van Donkelaar wrote:

I'm a little surprised at the drop of more than 50% in anthropogenic SO2 and NO over Europe.
I decided to explore the large drop in emissions and have found the drop in the revised EMEP emissions that are presently available to be much lower in than what EMEP had estimated when I downloaded/reformatted the original dataset back in 2011.
It should be fairly straightforward for me to reprocess new EMEP input files based on their latest estimates, which run until 2012, if the old format is still useful post-HEMCO.
For now I would not recommend using the EMEP files originally slated as part of this update. Rather, I would suggest that I/we update the emission based on the latest EMEP inventory.

Update 3/31/15: Aaron has reprocessed EMEP emissions for 1990-2012. Please see this wiki post for more information.

--Melissa Sulprizio 16:40, 27 March 2015 (EDT)

Reprocess EMEP emissions for 1990-2012

This update was validated with the 1-month benchmark simulation v10-01i and approved on Approved 01 May 2015.

The original updated EMEP emissions for 2008-2010 (generated in 2012) were found to be too low in 1-month benchmark v10-01h, so Aaron van Donkelaar reprocessed EMEP anthropogenic and ship emissions of CO, NH3, SO2, and NO for 1990-2012. The reprocessed data file is in netCDF format and includes data that has been mapped from the original 50 km grid x 50 km grid onto a 0.1° x 0.1° grid. Regridding to the finer 0.1° x 0.1° resolution allowed these input files to maintain a reasonable representation of their original grid characteristics, thereby permitting HEMCO to more accurately regrid onto a variety of GEOS-Chem grids. This file is read with the HEMCO emissions component and will be validated with 1-year benchmark v10-01h-Run0.

--Melissa Sulprizio 14:08, 31 March 2015 (EDT)

EMEP European ship emissions inventory

The EMEP inventory contains ship emissions for CO, NH3, SOx, and NOx.

Please also note the following fixes that are relevant for GEOS-Chem v9-01-01 or older versions:

  1. Bug with EMEP ship emissions for nested grids
  2. Bug fix for EMEP ship emissions
  3. Fix for initialization of EMEP ship emissions
  4. Inconsistent mask file for EMEP ship emissions

1980 to 1989

EMEP ship data from 1980-1989 emit NOx as pure NOx (and not HNO3 + 10*O3). It was not possible to separate the NOx into HNO3 and O3 from this data set.

See Auvray and Bey [2005] and Vestreng and Klein [2002] for more information.

1990 to 2005

EMEP ship data from 1990-2005 are based on Vestreng et al. [2007].

--Bob Y. 12:00, 20 December 2012 (EST)

Known issues

Zero EMEP emissions at 30N

This update was validated with 1-month benchmark simulation v11-01f and 1-year benchmark simulation v11-01f-geosfp-Run0. This version was approved on 16 Apr 2016.

Eloise Marais wrote:

Anthropogenic emissions are zero where the EMEP and EDGAR emissions coincide at 30 degrees North (using GEOS-Chem v10-01). This is because the EMEP mask extends to gridsquares where EMEP emissions are zero. Anthropogenic emissions are relatively low there, so this issue only shows up when NOx emissions are plotted on a log scale.

My quick fix is to shift the minimum latitude of the EMEP mask up by 1 degree latitude (so starting at 31 deg N instead of 30 deg N), in an updated mask file EMEP_mask.geos.1x1.20151222.nc. (The boxes that were removed show up as RED in the plot below.)

Boxes removed from EMEP mask.png

As of 06 Jan 2016, the GEOS-Chem Support Team has placed the updated mask file into the following HEMCO data directory:

HEMCO/MASKS/v2014-07/EMEP_mask.geos.1x1.20151222.nc

We will validate Eloise's update with the 1-month benchmark simulation GEOS-Chem v11-01f. This will involve changing the EMEP mask filename in the proper HEMCO_Config.rc files.

NOTE: The old mask file (EMEP_mask.geos.1x1.nc) shall be preserved for backwards compatibility.

--Bob Yantosca (talk) 15:23, 30 March 2016 (UTC)

Observed differences in CO emissions attributed to new EMEP mask file

Note that simulations using the updated EMEP mask file may cause some emissions differences w/r/t simulations using the old EMEP mask:

Prasad Kasibhatla wrote:

I am surprised by the differences in CO emissions between v11-01d and v11-01f. I was struck by how large they seem to be. What first caught my eye was the difference in boundary-layer CO mixing ratios in Cairo - and it looks like it is related to a big change in anthro+bifuel emissions in one grid box....Is there a simple explanation for this?

Melissa Sulprizio replied:

Differences in the anthropogenic+biofuel CO emissions in Cairo were caused by the updated EMEP mask file introduced in v11-01f. Eloise Marais found that EDGAR emissions were zero at 30N because the EMEP mask extended too far south, so she provided a new file that shifts the EMEP mask north by 1 degree.

--Bob Yantosca (talk) 16:01, 15 April 2016 (UTC)

References

  1. EMEP website: http://www.emep.int/index.html
  2. Auvray, M., and I. Bey, Long-Range Transport to Europe: Seasonal Variations and Implications for the European Ozone Budget, J. Geophys. Res., 110, D11303, doi: 10.1029/2004JD005503, 2005.
  3. Vestreng, V., and H. Klein, Emission data reported to UNECE/EMEP: Quality assurance and trend analysis & Presentation of WebDab, MSC-W Status Report 2002, Norwegian Meteorological Institute, Oslo Norway, July 2002. (PDF)
  4. Vestreng, V., Mareckova, K., Kakareka, S., Malchykhina, A., and Kukharchyk, T., Inventory Review 2007; Emission Data Reported to LRTAP Convention and NEC Directive, MSC-W Technical Report 1/07, The Norwegian Meteorological Institute, Oslo, Norway, 2007. (PDF)

--Bob Y. 12:13, 20 December 2012 (EST)