|
|
(46 intermediate revisions by 2 users not shown) |
Line 1: |
Line 1: |
| This page describes the various ship emissions inventories that have been implemented in GEOS-Chem. | | This page describes the various ship emissions inventories that have been implemented in GEOS-Chem. |
|
| |
| <span style="color:green">'''''NOTE: In [[GEOS-Chem v10-01]] and higher versions, all ship emissions (including the PARANOX plume model) are handled by the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
|
| == Overview == | | == Overview == |
|
| |
|
| At present there are several different ship emissions options in GEOS-Chem. All of these are available for use via the [[HEMCO|HEMCO emissions component]]. | | At present there are several different ship emissions options in GEOS-Chem. All of these are available for use via the [https://hemco.readthedocs.io Harmonized Emissions Component (HEMCO)]. |
|
| |
|
| ;ARCTAS: The ARCTAS pre-mission ship SO2 emissions from David Streets. [''Eyring et al'', 2005a and 2005b]
| | {| border=1 cellspacing=0 cellpadding=5 |
| ;Corbett: Ship emissions inventory as described in ''Corbett et al'' [1999].
| | |-valign="top" bgcolor="#CCCCCC" |
| ;EDGAR: SO2 ship emissions from the EDGAR 3.2 FT2000 global inventory.
| | !Inventory |
| ;EMEP: Ship emissions from the EMEP (European Monitoring and Evaluation Programme).
| | !Species included |
| ;ICOADS: Emission inventory based on data from the International Comprehensive Ocean−Atmosphere Data Set (ICOADS), the Automated Mutual-Assistance Vessel Rescue System (AMVER) data sets [''Wang et al'', 2008]. This is the recommended global ship emissions inventory.
| | !Year(s) |
| | !HEMCO Hierarchy |
| | !Notes |
|
| |
|
| Starting with [[GEOS-Chem v8-01-04]], ship-emitted NOx is no longer emitted directly as NOx, but as the combination HNO3 + 10*O3, following ''Chen et al'' [2005].
| | |-valign="top" |
| | !colspan=5 bgcolor="CCFFFF"|Currently-used ship emission inventories |
|
| |
|
| --[[User:Bmy|Bob Y.]] ([[User talk:Bmy|talk]]) 20:27, 26 October 2015 (UTC) | | |-valign="top" |
| | |[[CEDS_anthropogenic_emissions#Include_CEDS_ship_emissions|CEDS]] |
| | |NO, CO, SOAP, SO2, SO4, pFe, NH3, BC, OC, CO2, CH4, EOH, C2H6, C3H8, ALK4, C2H4, PRPE, C2H2, BENZ, TOLU, XYLE, CH2O, ALD2, MEK, HCOOH |
| | |1950-2014 |
| | |5 |
| | | |
| | *CEDS is the default global ship emission inventory in [[GEOS-Chem 12#12.1.0|GC 12.1.0]] and later versions. |
| | *Now supersedes ARCTAS, ICOADS, and/or Corbett. |
|
| |
|
| == ARCTAS == | | |-valign="top" |
| | |[[Anthropogenic_emissions#HTAP|HTAP]] |
| | |CO, SOAP, SO2 |
| | |2008-2010 |
| | |4 |
| | | |
| | *Inventory commented out in HEMCO config file. |
| | *HTAP is a research option but not a standard option. |
|
| |
|
| The ARCTAS pre-mission inventory of SO2 was first introduced in [[GEOS-Chem v8-01-04]].
| | |-valign="top" |
| | !colspan=5 bgcolor="CCFFFF"|Deprecated ship emission inventories |
|
| |
|
| '''''Philippe Le Sager wrote:'''''
| | |-valign="top" |
| | |[[#ARCTAS|ARCTAS]] |
| | |SO2 |
| | |2008 |
| | |1 |
| | | |
| | *Superseded by [[CEDS_anthropogenic_emissions#Include_CEDS_ship_emissions|CEDS]] in [[GEOS-Chem 12#12.1.0|GC 12.1.0]] and later versions. |
| | *Can still be used as a research option. |
|
| |
|
| <blockquote>The ship emission is based on the work by ''Eyring et al.'', [2005a and 2005b], which estimates the total international ship emissions for 1985, 1990, 2001, and 2020 (projection). The ship emission for each individual year is interpreted based on the above years, and the spatial pattern (gridded) is mapped based on the EDGAR gridded ship emission for 2000 (total amount from EDGAR is scaled to Eyring-based number).</blockquote>
| | |-valign="top" |
| | |[[#ICOADS|ICOADS]] |
| | |CO, SOAP, SO2 |
| | |2002 |
| | |2 |
| | | |
| | *SO2 commented out |
| | *Superseded by [[CEDS_anthropogenic_emissions#Include_CEDS_ship_emissions|CEDS]] in [[GEOS-Chem 12#12.1.0|GC 12.1.0]] and later versions. |
| | *Can still be used as a research option. |
|
| |
|
| <blockquote>If you want to reference the work on publication or website, you may either mention ''Diehl et al.'' [manuscript in preparation, 2009] or refer to the [http://www-lscedods.cea.fr/aerocom/AEROCOM_HC/readme_ship AeroCom readme document] (prepared by Diehl).</blockquote>
| | |-valign="top" |
| | |[[#Corbett|Corbett]] |
| | |SO2 |
| | |1985 |
| | |3 |
| | | |
| | *Inventory commented out in HEMCO config file |
| | *Superseded by [[CEDS_anthropogenic_emissions#Include_CEDS_ship_emissions|CEDS]] in [[GEOS-Chem 12#12.1.0|GC 12.1.0]] and later versions. |
| | *Can still be used as a research option. |
|
| |
|
| === In GEOS-Chem v10-01 and newer versions === | | |-valign="top" |
| | |[[EDGAR_v4.3_anthropogenic_emissions|EDGAR]] |
| | |NO, CO, SO2 |
| | |1970-2010 |
| | |N/A |
| | | |
| | *Chi Li wrote: Aviation and shipping emissions are not included [in EDGAR v4.3], since it is hard to spatially assign the total emissions of international aviation and shipping in the time series data. |
|
| |
|
| In [[GEOS-Chem v10-01]] and newer versions, the ARCTAS ship emissions data files are read with the [[HEMCO|HEMCO emissions component]]. We have created [http://wiki.seas.harvard.edu/geos-chem/index.php/Ship_emissions#ARCTAS new ARCTAS 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]].
| | |-valign="top" |
| | |[[EMEP_European_anthropogenic_emissions#EMEP_European_ship_emissions_inventory|EMEP]] |
| | |CO, SOAP, SO2 |
| | |1990-2012 |
| | |10 |
| | | |
| | *Retired in [[GEOS-Chem 12#12.5.0|GEOS-Chem 12.5.0]]. |
|
| |
|
| --[[User:Bmy|Bob Y.]] 13:21, 3 March 2015 (EST) | | |-valign="top" |
| | |[[EPA/NEI11_North_American_emissions|EPA/NEI2011]] |
| | |HNO2, CO, SOAP, NH3, CH2O, RCHO, MACR, ACET, C3H8, MEK, ALK4, PRPE, EOH, MOH, XYLE, TOLU, BENZ, SO2, SO4, pFe, C2H4, C2H6, ALD2, BC, OC |
| | |2006-2013 |
| | |50 |
| | | |
| | *<tt>NEI2011_SHIP_HOURLY</tt> and <tt>NEI2011_SHIP_MONTHLY</tt> are set to false by default |
| | *Katie Travis wrote, "The ship inventory for NEI11 only goes a few km from land, so it isn’t a complete ‘ship inventory’. Maybe we just shouldn’t use it as a default." |
|
| |
|
| === In GEOS-Chem versions prior to v10-01 ===
| | |} |
|
| |
|
| In [[GEOS-Chem v9-02]] and older versions, [ftp://ftp.as.harvard.edu/pub/geos-chem/data/GEOS_1x1/ARCTAS_SHIP_2008/README the ARCTAS ship emissions data files are stored in binary punch format]. They are read from disk by the source code in module <tt>GeosCore/arctas_ship_emiss_mod.f</tt>.
| | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:21, 18 June 2019 (UTC)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 21:52, 16 September 2022 (UTC) |
|
| |
|
| --[[User:Bmy|Bob Y.]] 14:18, 13 February 2015 (EST)
| | == ARCTAS == |
| | <span style="color:red"><big><strong>ARCTAS is now superseded by [[CEDS_anthropogenic_emissions#Include_CEDS_ship_emissions|CEDS]] in [[GEOS-Chem 12#12.1.0|GEOS-Chem 12.1.0]] and later versions.</strong></big></span> |
|
| |
|
| == Corbett ==
| | The ARCTAS pre-mission inventory of SO2 was first introduced in [[GEOS-Chem v8-01-04]]. |
|
| |
|
| For more information about this inventory please see: Corbett, J. J., P. S. Fischbeck, and S. N. Pandis, ''Global nitrogen and sulfur inventories for oceangoing ships'', <u>J. Geophys. Res.</u>, '''104'''(D3), 3457–3470, 1999.
| | '''''Philippe Le Sager wrote:''''' |
|
| |
|
| === In GEOS-Chem v10-01 and newer versions ===
| | <blockquote>The ship emission is based on the work by ''Eyring et al.'', [2005a and 2005b], which estimates the total international ship emissions for 1985, 1990, 2001, and 2020 (projection). The ship emission for each individual year is interpreted based on the above years, and the spatial pattern (gridded) is mapped based on the EDGAR gridded ship emission for 2000 (total amount from EDGAR is scaled to Eyring-based number).</blockquote> |
|
| |
|
| In [[GEOS-Chem v10-01]] and newer versions, the Corbett et al ship emissions data files are read with the [[HEMCO|HEMCO emissions component]]. We have created [http://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/CORBETT_SHIP/v2014-07/README new Corbett 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]].
| | <blockquote>If you want to reference the work on publication or website, you may either mention ''Diehl et al.'' [manuscript in preparation, 2009] or refer to the [http://www-lscedods.cea.fr/aerocom/AEROCOM_HC/readme_ship AeroCom readme document] (prepared by Diehl).</blockquote> |
|
| |
|
| --[[User:Bmy|Bob Y.]] 13:21, 3 March 2015 (EST) | | In [[GEOS-Chem v10-01]] and newer versions, the ARCTAS ship emissions data files are read with the [[HEMCO|HEMCO emissions component]]. We have created [http://wiki.seas.harvard.edu/geos-chem/index.php/Ship_emissions#ARCTAS ARCTAS 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]]. |
|
| |
|
| === In GEOS-Chem versions prior to v10-01 ===
| | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 22:01, 16 September 2022 (UTC) |
|
| |
|
| In [[GEOS-Chem v9-02]] and older versions, the Corbett et al ship emissions data files are stored in binary punch format at the following resolutions:
| | == Corbett == |
| | <span style="color:red"><big><strong>Corbett is now superseded by [[CEDS_anthropogenic_emissions#Include_CEDS_ship_emissions|CEDS]] in [[GEOS-Chem 12#12.1.0|GEOS-Chem 12.1.0]] and later versions.</strong></big></span> |
|
| |
|
| # [ftp://ftp.as.harvard.edu/gcgrid/data/GEOS_0.5x0.666_CH/sulfate_sim_200508/README 0.5° x 0.666° China nested grid]
| | For more information about this inventory please see: Corbett, J. J., P. S. Fischbeck, and S. N. Pandis, ''Global nitrogen and sulfur inventories for oceangoing ships'', <u>J. Geophys. Res.</u>, '''104'''(D3), 3457–3470, 1999. |
| # [ftp://ftp.as.harvard.edu/gcgrid/data/GEOS_0.5x0.666_NA/sulfate_sim_200508/README 0.5° x 0.666° North America nested grid]
| |
| # [ftp://ftp.as.harvard.edu/gcgrid/data/GEOS_2x2.5/sulfate_sim_200508/README 2° x 2.5° global grid]
| |
| # [ftp://ftp.as.harvard.edu/gcgrid/data/GEOS_4x5/sulfate_sim_200508/README 4° x 5° global grid]
| |
| | |
| These data files are read from disk by routine <tt>READ_SHIP_SO2</tt> in <tt>GeosCore/sulfate_mod.f</tt>.
| |
| | |
| --[[User:Bmy|Bob Y.]] 14:13, 13 February 2015 (EST)
| |
| | |
| == EDGAR ==
| |
| | |
| === EDGAR v4.2 ship emissions===
| |
| | |
| Please see [[EDGAR v4.2 anthropogenic emissions#Ship emissions|this wiki post on our EDGAR v4.2 anthropogenic emissions wiki page]] for more information about EDGAR inventories of ship exhaust.
| |
| | |
| --[[User:Melissa Payer|Melissa Sulprizio]] 16:39, 17 February 2015 (EST)
| |
| | |
| === EDGAR v3 ship emissions ===
| |
| | |
| Please see [[EDGAR v3 anthropogenic emissions#Ship emissions|this wiki post on our EDGAR v3 anthropogenic emissions wiki page]] for more information about EDGAR inventories of ship exhaust.
| |
|
| |
|
| --[[User:Bmy|Bob Y.]] 12:29, 20 December 2012 (EST) | | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 22:02, 16 September 2022 (UTC) |
| | |
| == EMEP ==
| |
| | |
| For a complete discussion of EMEP ship emissions, please see the [[EMEP_European_anthropogenic_emissions#EMEP_European_ship_emissions_inventory|EMEP European ship emissions inventory]] post on our [[EMEP European anthropogenic emissions]] wiki page.
| |
| | |
| --[[User:Bmy|Bob Y.]] 12:03, 20 December 2012 (EST)
| |
|
| |
|
| == ICOADS == | | == ICOADS == |
| | <span style="color:red"><big><strong>ICOADS is now superseded by [[CEDS_anthropogenic_emissions#Include_CEDS_ship_emissions|CEDS]] in [[GEOS-Chem 12#12.1.0|GEOS-Chem 12.1.0]] and later versions.</strong></big></span> |
|
| |
|
| The ICOADS ship inventory was first introduced in [[GEOS-Chem v8-02-03]]. It contains the species NOx, CO, and SO2. [mailto:cklee79@gmail.com Chulkyu Lee] worked on implementing ICOADS into GEOS-Chem. | | The ICOADS ship inventory was first introduced in [[GEOS-Chem v8-02-03]]. It contains the species NOx, CO, and SO2. [mailto:cklee79@gmail.com Chulkyu Lee] worked on implementing ICOADS into GEOS-Chem. |
Line 90: |
Line 126: |
| <blockquote>Ship activity patterns depicted by the International Comprehensive Ocean−Atmosphere Data Set (ICOADS), the Automated Mutual-Assistance Vessel Rescue System (AMVER) data set, and their combination demonstrate different spatial and statistical sampling biases. These differences could significantly affect the accuracy of ship emissions inventories and atmospheric modeling. We demonstrate (using ICOADS) a method to improve global-proxy representativeness by trimming over-reporting vessels that mitigates sampling bias, augment the sample data set, and account for ship heterogeneity. Apparent under-reporting to ICOADS and AMVER by ships near coastlines, perhaps engaged in coastwise (short sea) shipping especially in Europe, indicates that bottom-up regional inventories may be more representative locally. Primarily due to the long time series available publicly for ICOADS data, the improved ICOADS data set may be the most appropriate global ship traffic proxy identified to date to be used for a top-down approach. More generally, these three spatial proxies can be used together to perform uncertainty analyses of ship air-emissions impacts on a global scale (http://coast.cms.udel.edu/GlobalShipEmissions/).</blockquote> | | <blockquote>Ship activity patterns depicted by the International Comprehensive Ocean−Atmosphere Data Set (ICOADS), the Automated Mutual-Assistance Vessel Rescue System (AMVER) data set, and their combination demonstrate different spatial and statistical sampling biases. These differences could significantly affect the accuracy of ship emissions inventories and atmospheric modeling. We demonstrate (using ICOADS) a method to improve global-proxy representativeness by trimming over-reporting vessels that mitigates sampling bias, augment the sample data set, and account for ship heterogeneity. Apparent under-reporting to ICOADS and AMVER by ships near coastlines, perhaps engaged in coastwise (short sea) shipping especially in Europe, indicates that bottom-up regional inventories may be more representative locally. Primarily due to the long time series available publicly for ICOADS data, the improved ICOADS data set may be the most appropriate global ship traffic proxy identified to date to be used for a top-down approach. More generally, these three spatial proxies can be used together to perform uncertainty analyses of ship air-emissions impacts on a global scale (http://coast.cms.udel.edu/GlobalShipEmissions/).</blockquote> |
|
| |
|
| === In GEOS-Chem v10-01 and newer versions ===
| | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 22:02, 16 September 2022 (UTC) |
| | |
| In [[GEOS-Chem v10-01]] and newer versions, the ICOADS ship emissions are read with the [[HEMCO|HEMCO emissions component]]. We have created [http://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/ICOADS_SHIP/v2014-07/README new ICOADS 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]].
| |
| | |
| --[[User:Bmy|Bob Y.]] 13:22, 3 March 2015 (EST) | |
| | |
| === In GEOS-Chem versions prior to v10-01 ===
| |
| | |
| In [[GEOS-Chem v9-02]] and older versions, [http://ftp.as.harvard.edu/gcgrid/data/GEOS_NATIVE/ICOADS_200907/README the ICOADS data files are stored in binary punch format]. They are read from disk by the source code in module <tt>GeosCore/icoads_ship_mod.f</tt>.
| |
| | |
| --[[User:Bmy|Bob Y.]] 12:34, 13 February 2015 (EST)
| |
|
| |
|
| == PARANOX ship plume model == | | == PARANOX ship plume model == |
|
| |
| === In GEOS-Chem v10-01 and higher versions ===
| |
|
| |
|
| An updated version of PARANOX was brought into [[GEOS-Chem v10-01]] via the [[HEMCO|HEMCO emissions component]]. | | An updated version of PARANOX was brought into [[GEOS-Chem v10-01]] via the [[HEMCO|HEMCO emissions component]]. |
|
| |
|
| ==== PARANOX updates by Chris Holmes and Geert Vinken (July 2014) ====
| | === PARANOX updates by Chris Holmes and Geert Vinken (July 2014) === |
|
| |
|
| This post describes the original updates to PARANOX made by Chris Holmes and Geert Vinken in July 2014. These updates have since been brought into GEOS-Chem as an [[#Re-implementation as a HEMCO extension|extension to the HEMCO emissions component]]. | | This post describes the original updates to PARANOX made by Chris Holmes and Geert Vinken in July 2014. These updates have since been brought into GEOS-Chem as an [[#Re-implementation as a HEMCO extension|extension to the HEMCO emissions component]]. |
Line 115: |
Line 139: |
|
| |
|
| :Geert and I have updated the ship plume emissions in GEOS-Chem. A paper based on the updated model is now accepted in ACP ([http://www.atmos-chem-phys.net/14/6801/2014/acp-14-6801-2014.html see this link]). | | :Geert and I have updated the ship plume emissions in GEOS-Chem. A paper based on the updated model is now accepted in ACP ([http://www.atmos-chem-phys.net/14/6801/2014/acp-14-6801-2014.html see this link]). |
| :I am writing to pass along the updates to you for inclusion in the GC standard code.
| |
|
| |
| :I have posted a patch file and new look-up tales that accompany this code update at
| |
| ftp://halo.ess.uci.edu/Public/holmes/shipnox_v9-02/
| |
| :The patch will apply to the git commit tagged as "v9-02-Public-Release" on 2014-2-27. After applying the patch, you will need to update the "101 FORMAT" statement to point to the correct data directory on your system.
| |
|
| |
| :'''Update 7/17/14: The above patch file is now superseded by the patch file [[Media:Patch_paranox_ws.v9-02b.txt|posted here]]. The new patch file also applied to the commit tagged as "v9-02-Public-Release" and includes a fix for high values of O3 dry deposition. See [[Ship_emissions#Fix_for_high_values_of_O3_dry_deposition|this wiki post]] for more information.'''
| |
|
| |
| :The updated features are as follows: | | :The updated features are as follows: |
| :#Ship plume chemical aging now depends on wind speed, which directly affects the plume dispersion rates and indirectly affects the plume chemistry. | | :#Ship plume chemical aging now depends on wind speed, which directly affects the plume dispersion rates and indirectly affects the plume chemistry. |
Line 130: |
Line 146: |
| :#CH4 oxidation in ship plumes is now diagnosed and written in the <tt>paranox_ts</tt> files. This is only a diagnostic and does not affect the model chemistry. | | :#CH4 oxidation in ship plumes is now diagnosed and written in the <tt>paranox_ts</tt> files. This is only a diagnostic and does not affect the model chemistry. |
|
| |
|
| :As with previous versions, <tt>paranox_mod.F</tt> uses a look-up table derived from a Gaussian plume model to provide effective emission factors for NOx, O3, and HNO3 from ships. Both the Gaussian plume model and GEOS-Chem with the embedded look-up tables have been evaluated against observations (Holmes et al. 2014; Vinken et al., 2011). The emission factors are based on 8 environmental variables from the CTM: ambient concentrations of NOx and O3, solar zenith angle at emission time and 5 hours later, photolysis rates of NO2 and O3, temperature, and wind speed. | | :As with previous versions, PARANOX uses a look-up table derived from a Gaussian plume model to provide effective emission factors for NOx, O3, and HNO3 from ships. Both the Gaussian plume model and GEOS-Chem with the embedded look-up tables have been evaluated against observations (Holmes et al. 2014; Vinken et al., 2011). The emission factors are based on 8 environmental variables from the CTM: ambient concentrations of NOx and O3, solar zenith angle at emission time and 5 hours later, photolysis rates of NO2 and O3, temperature, and wind speed. |
|
| |
|
| --[[User:Melissa Payer|Melissa Sulprizio]] 10:21, 18 June 2014 (EDT) | | --[[User:Melissa Payer|Melissa Sulprizio]] 10:21, 18 June 2014 (EDT) |
|
| |
|
| ==== Re-implementation as a HEMCO extension ====
| | === Re-implementation as a HEMCO extension === |
|
| |
|
| [[User:Christoph Keller|Christoph Keller]] and the [[GEOS-Chem Support Team]] have re-implemented [[#PARANOX updates by Chris Holmes and Geert Vinken (July 2014)|the updated implementation of PARANOX described above]] as a HEMCO extension. A new HEMCO extension module (<tt>HEMCO/Extensions/hcox_paranox_mod.F90</tt>) has now replaced the original <tt>paranox_mod.F</tt> module. | | [[User:Christoph Keller|Christoph Keller]] and the [[GEOS-Chem Support Team]] have re-implemented [[#PARANOX updates by Chris Holmes and Geert Vinken (July 2014)|the updated implementation of PARANOX described above]] as a HEMCO extension. A new HEMCO extension module (<tt>HEMCO/Extensions/hcox_paranox_mod.F90</tt>) has now replaced the original <tt>paranox_mod.F</tt> module. |
Line 141: |
Line 157: |
|
| |
|
| --[[User:Bmy|Bob Y.]] 13:23, 3 March 2015 (EST) | | --[[User:Bmy|Bob Y.]] 13:23, 3 March 2015 (EST) |
|
| |
| ==== Discussion following 1-month benchmark v10-01h ====
| |
|
| |
| The [[#PARANOX updates by Chris Holmes and Geert Vinken (July 2014)|PARANOX updates described above]] were validated in [[GEOS-Chem_v10-01_benchmark_history#v10-01h|GEOS-Chem v10-01h]]. The following discussion ensued.
| |
|
| |
| '''''[[User:Chris_Holmes|Chris Holmes]] wrote:'''''
| |
|
| |
| :I sent in the updates to ship NOx chemistry (aka PARANOX). In the v10-01h benchmark plots I see that surface O3 concentrations in the North Pacific Ocean have decreased 20 ppb and concentrations are now 10 ppb or lower. The spatial pattern is consistent with ship emissions causing the change, but the magnitude is much larger than I expected based on the changes that I made. In my tests with v9-02, my updates to ship NOx chemistry change O3 by only 2 ppb.
| |
|
| |
| :Are there other model changes that could have affected MBL chemistry? Did the underlying inventory of ship emissions change?
| |
|
| |
| '''''[[User:Christoph Keller|Christoph Keller]] wrote:'''''
| |
|
| |
| :Thanks for looking into this. Do you think that PARANOX now destroys more O3, or produces more?
| |
| :In anticipation of FlexChem, the implementation of dry deposition has changed in v10-01h compared to previous versions: All dry deposition calculations are now done as part of turbulence (in the non-local PBL scheme, or before applying the full PBL mixing if this option is selected). Also, all dry deposition frequencies are now calculated over the same deposition height (level 1 for non-local PBL, entire PBL for full mixing scheme).
| |
| :In other words, the PARANOX loss rates for O3 (and HNO3) are now calculated based on surface concentrations, and then applied immediately afterwards in the non-local PBL scheme. In your implementation in v9-02, the PARANOX dry deposition rates were calculated (and applied) over the entire PBL height, and then applied in chemistry (e.g. passed to the chemical solver via calcrate.F). The new version now also includes HNO3 deposition from PARANOX.
| |
|
| |
| :Do you think these changes may explain the differences? And is there anything that should be changed?
| |
|
| |
| --[[User:Melissa Payer|Melissa Sulprizio]] 13:15, 27 March 2015 (EDT)
| |
|
| |
| ==== Bug fixes for the PARANOX HEMCO extension ====
| |
|
| |
| <span style="color:green">'''''These fixes were 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>
| |
|
| |
| Christoph Keller found two bugs in <tt>HEMCO/Extensions/hcox_paranox_mod.F</tt>:
| |
|
| |
| :'''1. Bug fix in the calculation of the current date SZA'''
| |
| : In routine <tt>PARANOX_LUT</tt> change:
| |
|
| |
| VARS(5) = ExtState%SUNCOSmid%Arr%Val(I,J)
| |
|
| |
| : to:
| |
|
| |
| VARS(5) = ASIND( ExtState%SUNCOSmid%Arr%Val(I,J) )
| |
|
| |
| :'''2. Bug fix in calculation of H2O ambient air concentration'''
| |
| :'''''[[User:Christoph Keller|Christoph Keller]] wrote:'''''
| |
| ::The relative humidity was treated as kg H2O/kg air, whereas in reality it is g H2O/kg air. The units are wrong in gigc_state_met_mod.F90, which is what confused me. I did a one-day simulation with the updated code, and the ozone production efficiency now looks basically the same as in the version that Chris Holmes provided to us.
| |
| ::To add just the H2O fix, you need to replace the calculation of variable H2O as follows:
| |
|
| |
| - ! H2O, molec/cm3. Get from specific humidity, which is in kg/kg.
| |
| - ! H2O = CSPEC(JLOOP,IH2O)
| |
| - H2O = ExtState%SPHU%Arr%Val(I,J,1) * 1.e3_sp &
| |
| - * ExtState%AIR%Arr%Val(I,J,1) / MWH2O &
| |
| - * HcoState%Phys%Avgdr / 1e6_sp
| |
|
| |
| + ! H2O, molec/cm3. Get from specific humidity, which is in g/kg.
| |
| + H2O = ExtState%SPHU%Arr%Val(I,J,1) / 1.0e3_sp * DENS &
| |
| + * HcoState%Phys%AIRMW / MWH2O
| |
|
| |
| :'''3. Fix for setting option in the HEMCO configuration file'''
| |
| :The 1-month benchmark [[GEOS-Chem_v10-01_benchmark_history#v10-01h|v10-01h]] revealed that ship emissions weren't being read into HEMCO. PARANOX was still called, so this resulted in a large decrease in O3 over the oceans. It turns out the order of the logical switches in the HEMCO configuration file was causing this error. When <tt>NEI2011_SHIP</tt> was listed before <tt>SHIP</tt>, HEMCO automatically set <tt>SHIP</tt> to <tt>NEI2011_SHIP</tt>. To resolve this issue, Christoph Keller modified HEMCO so that routine <tt>GetExtOpt</tt> searches for the option name that exactly matches the desired option.
| |
|
| |
| :'''4. Fix SUNCOS values stores in the HEMCO restart file'''
| |
| :The PARANOX_SUNCOS values saved to the HEMCO restart file were one hour off, which resulted in the following error in PARANOX:
| |
|
| |
| !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
| |
| HEMCO ERROR: LUT error: Fracnox should be between 0 and 1!
| |
| ERROR LOCATION: PARANOX_LUT
| |
| !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
| |
|
| |
| :Christoph Keller updated <tt>HEMCO/Extensions/hcox_paranox_mod.F90</tt> to save the correct historic SUNCOS values to the HEMCO restart file.
| |
|
| |
| --[[User:Melissa Payer|Melissa Sulprizio]] 17:34, 20 April 2015 (EDT)
| |
|
| |
| '''''Update:''''' You can still get this error if the start date of your GEOS-Chem simulation does not correspond with the date of the dat ain the HEMCO restart file. Please see [[GEOS-Chem_v10-01#Avoid_PARANOX_SUNCOS_error_by_using_the_proper_timestamps|this wiki post for more information on how to work around this issue]].
| |
|
| |
| --[[User:Bmy|Bob Y. ]] ([[User talk:Bmy|talk]]) 19:45, 19 May 2015 (UTC)
| |
|
| |
| === In versions of GEOS-Chem prior to v10-01 ===
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-03_benchmark_history#v9-01-03f|v9-01-03f]] and approved on 17 Feb 2012. This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35a]].'''''</span>
| |
|
| |
| <span style="color:red">'''''Module <tt>paranox_mod.F</tt> was removed from [[GEOS-Chem v10-01]] and newer versions. The PARANOX ship plume model is now implemented via the [[HEMCO]] emissions component.'''''</span>
| |
|
| |
| '''''[mailto:g.c.m.vinken@tue.nl Geert Vinken] wrote:'''''
| |
|
| |
| :[At the] past GEOS-Chem user's meeting I had a [http://acmg.seas.harvard.edu/geos/meetings/2011/pdf/Day1/MonF_Photo_vinken_geert_1_pc.pdf presentation about accounting for non-linear chemistry] of ship plumes in the GEOS-Chem model. I think my code is ready to be incorporated in the GEOS-Chem model.
| |
|
| |
| :I've attached [[Media:Vinken_Appendix_D.pdf|an Appendix]] from a report I wrote a year ago, which contains some more information about how the method is programmed. I'm aware that some routines might need to be moved to new files (e.g. the interpolation of the Look-Up-Table is now done in <tt>icoads_ship_mod.f</tt>).
| |
|
| |
| :The [http://www.atmos-chem-phys.net/11/11707/2011/acp-11-11707-2011.html reference for my work] has just been published.
| |
|
| |
| --[[User:Bmy|Bob Y.]] 10:27, 1 December 2011 (EST)
| |
|
| |
|
| == References == | | == References == |
Line 244: |
Line 175: |
|
| |
|
| --[[User:Bmy|Bob Y.]] 10:31, 1 December 2011 (EST) | | --[[User:Bmy|Bob Y.]] 10:31, 1 December 2011 (EST) |
|
| |
| == Previous issues that are now resolved ==
| |
|
| |
| === Fix bug in ARCTAS_SHIP scale factors ===
| |
|
| |
| <span style="color:darkorange">'''''This update will be included in [[GEOS-Chem 12#12.0.0|GEOS-Chem 12.0.0]].'''''</span>
| |
|
| |
| '''''[[User:barronh|Barron Henderson]] wrote:'''''
| |
|
| |
| <blockquote>ARCTAS has a base year of 2008, but the implicit AnnualScalar base year is 2000. So scaling from 2008 to another year would require normalizing the AnnualScalar to the base year. Instead of a ScaleID of just <tt>11</tt>, it would be <tt>11/19</tt>.
| |
|
| |
| In the <tt>HEMCO_Config.rc</tt> files for full-chemistry and aerosol-only simulations, we have this entry:</blockquote>
| |
|
| |
| 0 ARCTAS_SHIP_SO2 ... SO2 2008/1/1/0 C xy kg/m2/s SO2 <span style="color:red">11</span> 10 1
| |
|
| |
| <blockquote>So it is proposed to change this to:</blockquote>
| |
|
| |
| 0 ARCTAS_SHIP_SO2 ... SO2 2008/1/1/0 C xy kg/m2/s SO2 <span style="color:green">11/19</span> 10 1
| |
|
| |
| --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 15:35, 3 July 2018 (UTC)
| |
|
| |
| === Fix for high values of O3 dry deposition ===
| |
|
| |
| <span style="color:green">'''''This update was validated with the 1-month benchmark simulation [[GEOS-Chem v10-01 benchmark history#v10-01h|v10-01h]] and approved on Approved 01 May 2015.'''''</span>
| |
|
| |
| '''''[[User:Tomas_Sherwen|Tomás Sherwen]] wrote:'''''
| |
|
| |
| :I am experiencing an issue with the O3 dry deposition fields in V9-2 (GEOS-FP), and note a similar issue has been [[GEOS-Chem_v9-02#O3_completely_removed_at_Shanghai_grid_box_with_GEOS-FP_4x5_simulation|previously reported]] within the wiki with a response from Bob Yantosca.
| |
|
| |
| :The O3 deposition fields are extremely high over Shanghai, northern Europe and near Lake Superior. A log plot of this is attached below. Masking for anomalous values, yields a annual deposition within expected bounds (821 Tg / yr). This anomaly is not present for other depositional fields tested (e.g. SO2, iodine tracers, etc). The issue is present within the “off-the-shelf” version i have of v9-2. Standard compile settings have been used, the flags implicated in response on the wiki were not used.
| |
|
| |
| The plot below shows the O3 dry deposition field (2012-07-01 to 2013-07-01, GEOS-FP, V9-2) / Tg O3 / yr (capped at 1000000000.0 for plotting)
| |
| :[[Image:O3dep_v9-02.png]]
| |
|
| |
| '''''[[User:Chris_Holmes|Chris Holmes]] responded:'''''
| |
|
| |
| :I have a solution. The [[Media:Patch_paranox_ws.v9-02b.txt|attached patch file]] applies to v9-02-Public-Release with the parent commit "Clean up use of State_Met in several modules". This supersedes and replaces my previous patch file that is currently [[Ship_emissions#Updates_to_ship_NOx_chemistry|posted on the wiki here]].
| |
|
| |
| :The wiki description of my changes since v9-02 is still accurate. The substantive difference between this patch and the previous one is that when O3 is lost in sub-grid-scale ship plumes, the loss is calculated by the chemical solver using a 1st order rate coefficient from paranox.
| |
|
| |
| :The problem that Tomás identified occurred because paranox calculates a 1st order loss rate for O3 in ship plumes by diving the flux by the ambient O3 concentration in a grid box. The 1st order loss rate can be quite high when the O3 concentration is low (divide by a small number). Tomás's problem occurs when the O3 concentration then increases after the paranox call and before the 1st order loss is applied, so that a large 1st order loss rate is applied to a large O3 concentration. Because of the timing and order of operations in GC, this can occur when using NLPBL mixing, but not with instant PBL mixing. With the plume O3 loss now calculated inside chemistry, the problem does not occur.
| |
|
| |
| :My patch requires adding the following new reaction to globchem.dat
| |
|
| |
| A 998 0.00E+00 0.0E+00 0 0 S 0.00 0. 0.
| |
| O3 +
| |
| = + + +
| |
| + + + +
| |
| + + + +
| |
| + + + +
| |
|
| |
| :I have not used KPP and haven't updated the KPP files.
| |
|
| |
| --[[User:Melissa Payer|Melissa Sulprizio]] 09:34, 17 July 2014 (EDT)
| |
|
| |
| == Outstanding issues not yet resolved ==
| |
|
| |
| None at this time.
| |
|
| |
| == Obsolete information ==
| |
|
| |
| These sections pertain to code or data that has been removed from the most recent versions of GEOS-Chem. We shall keep this information here for reference.
| |
|
| |
| === Only call PARANOX with the full-chemistry simulations ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:red">'''''Module <tt>emissions_mod.F</tt> has been removed from [[GEOS-Chem v10-01]] and newer versions. PARANOX is now implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem v9-02 benchmark history#1-month benchmark with GEOS-5 meteorology|v9-02r]] and approved on 14 Nov 2013.'''''</span>
| |
|
| |
| In routine <tt>DO_EMISSIONS</tt> (in module <tt>GeosCore/emissions_mod.F</tt>), the subroutine <tt>READ_PARANOX_LUT</tt> is called for fullchem, offline aerosol, tagged CO, and H2/HD simulations. But PARANOX assumes that you have NO and O3 tracers defined. Calling PARANOX from an offline simulation will result in a seg-fault error, caused by unallocated arrays.
| |
|
| |
| The solution is to only call <tt>READ_PARANOX_LUT</tt> for the full-chemistry simulations and not from any of the offline simulations.
| |
|
| |
| --[[User:Bmy|Bob Y.]] 15:23, 14 November 2013 (EST)
| |
|
| |
| === Bug in ND36 diagnostic when ship emissions are turned off ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:red">'''''This code was removed from [[GEOS-Chem v10-01]] and newer versions. The [[HEMCO|HEMCO emissions component]] now computes all emissions diagnostics internally.'''''</span>
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-02_benchmark_history#v9-02o|v9-02o]] and approved on 03 Sep 2013.'''''</span>
| |
|
| |
| ''''' [mailto:Yanko.Davila@colorado.edu Yanko Davila] wrote:'''''
| |
|
| |
| :It seems there may be a small potential bug in the ND36 diagnostic in diag3.F. Near lines 2435, the following:
| |
|
| |
| IF ( AD36_SHIP_COUNT >= 1 ) THEN
| |
| ARRAY(:,:,1) = AD36(:,:,MM) / SECONDS +
| |
| & AD36_SHIP(:,:,MM) / REAL( AD36_SHIP_COUNT )
| |
| ENDIF
| |
|
| |
| :should be
| |
|
| |
| IF ( AD36_SHIP_COUNT >= 1 ) THEN
| |
| ARRAY(:,:,1) = AD36(:,:,MM) / SECONDS +
| |
| & AD36_SHIP(:,:,MM) / REAL( AD36_SHIP_COUNT )
| |
| ELSE
| |
| ARRAY(:,:,1) = AD36(:,:,MM) / SECONDS
| |
| ENDIF
| |
|
| |
| :so that AD36 will have the correct units even if shipping emissions are turned off.
| |
|
| |
| --[[User:Melissa Payer|Melissa Sulprizio]] 11:36, 22 August 2013 (EDT)
| |
|
| |
| === Bugs in ship emissions and non-local PBL mixing ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-02_benchmark_history#v9-02l|v9-02l]] and approved on 26 Jun 2013. This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35a]].'''''</span>
| |
|
| |
| '''''[mailto:cdholmes@uci.edu Chris Holmes] wrote:'''''
| |
|
| |
| :I have some bug fixes for ship NOx emissions which have some significant impacts on [[NOx-Ox-HC-aerosol|HOx-NOx-O3 chemistry]]. In addition, I found an unresolved bug that affects all emissions, at least when using non-local PBL mixing. [mailto:g.c.m.vinken@tue.nl Geert Vinken] and I have talked and I believe he concurs with my recommended changes.
| |
|
| |
| :One proposal that I will make is that we should add global mass balance checks to the monthly benchmarks, in order to find and fix these kinds of bugs more quickly in the future. For example, the integrated emissions minus deposition and chemical loss over a 1-month run should equal the change in mass between the first and last restart files. This would help QA/QC during model development. Many, but not all, of the bugs below would have been caught by this kind of test.
| |
|
| |
| :BUG FIXES
| |
|
| |
| ::*Ship NOx emissions with the ICOADS-AMVER inventory were 3.3 times too low, due to a NO2 vs. N molecular weight error. Emissions should be 5.0Tg(N)/yr. As a result, O3 and HNO3 concentrations were underestimated, especially over the N. Atlantic and N. Pacific. This bug has probably persisted in all model versions with ICOADS ship emissions.
| |
| ::: '''''NOTE: Yanko Davila found that this bug fix may not be backwards compatible with [[GEOS-Chem v9-01-01]] and earlier when ship emissions were handled in emfossil.F. Applying the bug fix in earlier versions may result in unreasonably high values of ozone over the shipping zones.'''''
| |
| ::All remaining bugs afflict only [[GEOS-Chem v9-01-03|v9-01-03]] and later.
| |
| ::*Some portion (probably 25%) of ship emissions were neglected in the model runs with non-local PBL mixing, due to ship emissions being in calcrate.F, not setemis.F. I moved the ship code to setemis.F and reordered the call to Fast-J in chemdr.F to accommodate this. (The Fast-J move alone has no impact.)
| |
| ::*Molecular weights of NOx, HNO3, and O3 were all mixed up in ship plume code. Affects both simulation and diagnostics only [[GEOS-Chem v9-01-03|v9-01-03]], although there was a similar bug in v9-01-01 that we've already fixed in v9-01-02.
| |
| ::*ND32 and ND36 diagnostics gave incorrect totals because ship NOx is calculated every convective time step, when using non-linear PBL mixing, rather than every emission time step. I created additional diagnostic counters in order to calculate the correct averages, regardless of the mixing scheme chosen.
| |
|
| |
| :MINOR SCIENCE UPDATES
| |
|
| |
| :These fixes enable the code to do what it was supposed to do, although they're more than simple bug fixes.
| |
|
| |
| ::*Net O3 loss in ship plumes (which occurs at night) is now treated as dry deposition, rather than set to zero. We calculate a dry deposition frequency that would give an equivalent O3 loss rate.
| |
| ::*Ship plume NO2 should be included in Ox emissions (along with O3 production) when using the non-local PBL mixing. This is unnecessary when doing instant PBL mixing because then emissions occur with chemistry when there is no Ox family.
| |
| ::*ND63 diagnostics are now averages of user-specified time intervals and reported in units of molec/cm2/s or unitless (mol/mol). This improves consistency with other diagnostics.
| |
|
| |
| :REMAINING EMISSIONS BUG
| |
|
| |
| ::*With [[Boundary_layer_mixing#VDIFF|non-local PBL mixing]], there are *no emissions* during the first mixing and convection time step. Even worse, you can't tell this from the diagnostics, which show emissions as normal. This happens because some emission variables are not initialized until after chemistry runs, so setemis within vdiff_mod can't run. Emissions of all trace gases are affected, at least in v9-01-03. I don't know how long this error has persisted. Resolving this bug requires putting chemistry before PBL mixing in main. For fully decoupling emissions and chemistry code. I tested the reordering solution, but it causes 20% or greater changes in concentration of some species, so I leave this issue to the GCST. The emission and PBL mixing wiki pages should have a warning about this.
| |
|
| |
| --[[User:Bmy|Bob Y.]] 11:57, 3 May 2013 (EDT)<br>
| |
| --[[User:Melissa Payer|Melissa Sulprizio]] 11:20, 22 August 2013 (EDT)
| |
|
| |
| === Bug in ship CO emissions ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:red">'''''Module <tt>emfossil.F</tt> has been removed from [[GEOS-Chem v10-01]] and newer versions. EMEP ship emissions are now implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-02_benchmark_history#v9-02l|v9-02l]] and approved on 26 Jun 2013. This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35a]].'''''</span>
| |
|
| |
| '''''[mailto:cdholmes@uci.edu Chris Holmes] wrote:'''''
| |
|
| |
| :When EMEP ship emissions are selected in input.geos the EMEP ship emissions inventory was previously overwriting all ship CO emissions globally. As a result, ship CO emissions were zero outside of the EMEP region (European waters). The solution is to only overwrite within the EMEP region. This bug has persisted at least since [[GEOS-Chem v9-01-01|v9-01-01]].
| |
|
| |
| To fix this bug, change the following code in <tt>emfossil.F</tt> from:
| |
|
| |
| IF ( LEMEPSHIP ) SHIP =
| |
| $ GET_EMEP_ANTHRO( I, J, NN, SHIP=.TRUE.)
| |
|
| |
| to:
| |
|
| |
| IF ( LEMEPSHIP ) THEN
| |
|
| |
| ! If we are over the European region ...
| |
| ! (bug fix, now overwrite only EMEP region, cdh 1/2/2013)
| |
| IF ( GET_EUROPE_MASK( I, J ) > 0d0 ) THEN
| |
| SHIP = GET_EMEP_ANTHRO( I, J, NN, SHIP=.TRUE.)
| |
| ENDIF
| |
|
| |
| ENDIF
| |
|
| |
| --[[User:Melissa Payer|Melissa Payer]] 14:19, 3 January 2013 (EST)
| |
|
| |
| === NOx ship emissions restored to the ND32, ND36 diagnostics ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:red">'''''Module <tt>emfossil.F</tt> has been removed from [[GEOS-Chem v10-01]] and newer versions. The [[HEMCO|HEMCO emissions component]] now internally computes all diagnostics for emission quantities. The information below is obsolete, but we shall leave it here for reference.'''''</span>
| |
|
| |
| In versions [[GEOS-Chem v9-01-01|v9-01-01]] and [[GEOS-Chem v9-01-02|v9-01-02]], you may have noticed that ship NOx emissions seem to be lacking from the ND36 (<tt>ANTHSRCE</tt>) or ND32 (<tt>NOX-AN-$</tt>) diagnostics. This is because the ship NOx emissions in these versions were partitioned into HNO3 + 10*Ox (in <tt>GeosCore/emfossil.F</tt>), and then saved into separate slots of the ND36 diagnostic:
| |
|
| |
| ! Store as HNO3 and O3
| |
| ! Convert molec/cm2/s to molec/box/s (cdh, 10/20/2011)
| |
| EMISRR(I,J,IDEHNO3) = SHIP * AREA_CM2
| |
| EMISRR(I,J,IDEOX) = 10D0 * SHIP * AREA_CM2
| |
|
| |
| With the advent of the [[#PARANOX ship plume model|PARANOX plume model]] in [[GEOS-Chem v9-01-03]] and later versions, this has now been corrected. The resultant ship emissions are now properly archived as NOx, O3, and HNO3 in the ND36 diagnostics. Ship NOx is also saved into the ND32 diagnostic. This is done in routine <tt>GeosCore/calcrate.F</tt>, where PARANOX is invoked:
| |
|
| |
| !===================================================
| |
| ! Update diagnostics
| |
| !===================================================
| |
| !ND36 = Anthro source diagnostic...store as [molec/cm2]
| |
| !and convert to [molec/cm2/s] in DIAG3.F
| |
| IF ( ND36 > 0 ) THEN
| |
|
| |
| AD36(IX,IY,IDEHNO3) = AD36(IX,IY,IDEHNO3) + SHIP *
| |
| & (1D0 - fraction_nox) *
| |
| & XNUMOL(IDTHNO3) / AREA_CM2
| |
|
| |
| AD36(IX,IY,IDEOX) = AD36(IX,IY,IDEOX) + SHIP *
| |
| & int_ope * (1D0 - fraction_nox)*
| |
| & XNUMOL(IDTOX) / AREA_CM2
| |
|
| |
| AD36(IX,IY,IDENOX) = AD36(IX,IY,IDENOX) + SHIP *
| |
| & fraction_nox *
| |
| & XNUMOL(IDTNOX) / AREA_CM2
| |
| ENDIF
| |
|
| |
| ! ND32 = save anthro NOx for levels L=1,NOXEXTENT [molec/cm2/s]
| |
| IF ( ND32 > 0 ) THEN
| |
| AD32_an(IX,IY,IZ) = AD32_an(IX,IY,IZ) +
| |
| & ( SHIP * fraction_nox * XNUMOL(NN) /
| |
| & ( DTSRCE * AREA_CM2 ) )
| |
| ENDIF
| |
|
| |
| --[[User:Bmy|Bob Y.]] 16:12, 7 February 2013 (EST)
| |
|
| |
|
| |
| === Duplicate ship emissions in Asia in Streets 2000 inventory ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:red">'''''Module <tt>emfossil.F</tt> has been removed from [[GEOS-Chem v10-01]] and newer versions. PARANOX is now implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-03_benchmark_history#v9-01-03f|v9-01-03f]] and approved on 17 Feb 2012.'''''</span>
| |
|
| |
| Please see [[David_Streets_regional_emissions_for_China_and_SE_Asia#Duplicate_ship_emissions_in_Asia_in_Streets_2000_inventory|this wiki post on our David Streets regional emissions wiki page]] for a full description of the issue.
| |
|
| |
| --[[User:Bmy|Bob Y.]] 12:24, 20 December 2012 (EST)
| |
|
| |
| === Bug in ship NOx emissions ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:red">'''''Module <tt>emfossil.F</tt> has been removed from [[GEOS-Chem v10-01]] and newer versions. PARANOX is now implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-02_benchmark_history#v9-01-02s|v9-01-02s]] and approved on 5 Nov 2011. (Included in Adjoint [[GEOS-Chem_Adjoint_v34 | v34e]] )'''''</span>
| |
|
| |
| '''''[mailto:cdholmes@uci.edu Chris Holmes] wrote:'''''
| |
|
| |
| :I found a bug in ship NOx emissions in [[GEOS-Chem_v9-01-02_benchmark_history#v9-01-02n|v9-01-02n]]. Ship NOx emissions are supposed to be emitted as (1 mol HNO3 + 10 mol O3) per mol NOx in the inventory. The code previously used molar masses incorrectly so that the ratios were actually 0.73 mol HNO3 and 9.5 mol O3. Thus, this bug fix will slightly increase O3 and HNO3 over the oceans. I have pasted the corrected code for <tt>emfossil.F</tt> below:
| |
|
| |
| !---------------------------------------------------------------
| |
| ! Prior to cdh 10/20/2011, HNO3 and O3 emissions were incorrectly scaled
| |
| ! to ship NOx, due to molar masses being in the wrong place
| |
| ! ! Convert molec/cm2/s to kg/box/timestep to get same
| |
| ! ! units as EMISRN, ie default GEIA emiss (phs, 7/9/09)
| |
| ! SHIP = SHIP * ( DTSRCE * AREA_CM2 ) / XNUMOL(NN)
| |
| !
| |
| ! ! Store as HNO3 and O3
| |
| ! EMISRR(I,J,IDEHNO3) = SHIP * XNUMOL(IDTHNO3) / DTSRCE
| |
| ! EMISRR(I,J,IDEOX) = 10D0 * SHIP *
| |
| ! & XNUMOL(IDTOX) / DTSRCE
| |
| !
| |
| ! ! ND36 = Anthro source diagnostic...store as [molec/cm2]
| |
| ! ! and convert to [molec/cm2/s] in DIAG3.F
| |
| ! IF ( ND36 > 0 ) THEN
| |
| !
| |
| ! AD36(I,J,IDEHNO3) = AD36(I,J,IDEHNO3) + SHIP *
| |
| ! & XNUMOL(IDTHNO3) / AREA_CM2
| |
| !
| |
| ! AD36(I,J,IDEOX) = AD36(I,J,IDEOX) + SHIP *
| |
| ! & 10D0 * XNUMOL(IDTOX) / AREA_CM2
| |
| !
| |
| ! ENDIF
| |
| !---------------------------------------------------------------
| |
| ! Store as HNO3 and O3
| |
| ! Convert molec/cm2/s to molec/box/s (cdh, 10/20/2011)
| |
| EMISRR(I,J,IDEHNO3) = SHIP * AREA_CM2
| |
| EMISRR(I,J,IDEOX) = 10D0 * SHIP * AREA_CM2
| |
|
| |
| ! ND36 = Anthro source diagnostic...store as [molec/cm2]
| |
| ! and convert to [molec/cm2/s] in DIAG3.F
| |
| IF ( ND36 > 0 ) THEN
| |
|
| |
| AD36(I,J,IDEHNO3) = AD36(I,J,IDEHNO3) + SHIP * DTSRCE
| |
|
| |
| AD36(I,J,IDEOX) = AD36(I,J,IDEOX) + 10D0 * SHIP *
| |
| & DTSRCE
| |
|
| |
| ENDIF
| |
|
| |
| --[[User:Melissa Payer|Melissa Payer]] 13:44, 21 October 2011 (EDT)
| |
|
| |
| === Inconsistent mask file for EMEP ship emissions ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| Please see [[EMEP_European_anthropogenic_emissions#Inconsistent_mask_file_for_EMEP_ship_emissions|this wiki post on our EMEP emissions wiki page]] for a full description of the issue.
| |
|
| |
| --[[User:Bmy|Bob Y.]] 12:21, 20 December 2012 (EST)
| |
|
| |
| === Potential bug in EDGAR ship emissions code ===
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| <span style="color:red">'''''This code was removed from [[GEOS-Chem v10-01]] and newer versions. EDGAR emissions are now implemented via the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
| <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-01-02_benchmark_history#v9-01-02s|v9-01-02s]] and approved on 5 Nov 2011. (Included in Adjoint [[GEOS-Chem_Adjoint_v34 | v34e]] ) '''''</span
| |
|
| |
|
| |
| Please see [[EDGAR_v3_anthropogenic_emissions#Potential_bug_in_EDGAR_ship_emissions_code|this wiki post on our EDGAR anthropogenic emissions page]] for a full description of the issue.
| |
|
| |
| --[[User:Bmy|Bob Y.]] 12:32, 20 December 2012 (EST)
| |
|
| |
| === Minor fixes for EMEP ship emissions ===
| |
|
| |
| <span style="color:red">'''''This code was removed from [[GEOS-Chem v10-01]] and newer versions. EMEP ship emissions are now implemented through the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
| Please see these posts about minor bug fixes for EMEP ship emissions in [[EMEP European anthropogenic emissions#Bug fix for EMEP ship emissions|GEOS-Chem v8-02-04]] and in [[EMEP European anthropogenic emissions#Fix for initialization of EMEP ship emissions|GEOS-Chem v8-03-01]].
| |
|
| |
| --[[User:Bmy|Bob Y.]] 10:08, 25 February 2010 (EST)
| |
|
| |
| === Bug for ship emissions in emfossil.f ===
| |
|
| |
|
| |
| <span style="color:red">'''''Module <tt>emfossil.F</tt> was removed from [[GEOS-Chem v10-01]] and newer versions. EMEP ship emissions are now implemented through the [[HEMCO|HEMCO emissions component]].'''''</span>
| |
|
| |
| <span style="color:green">'''''NOTE: This fix has been incorporated into [[GEOS-Chem v9-01-01]].'''''</span>
| |
|
| |
| [[Image:Obsolete.jpg]]
| |
|
| |
| '''''[mailto:g.c.m.vinken@tue.nl Geert Vinken] wrote:'''''
| |
|
| |
| :'''Problem:''' The ship NOx emissions are as HNO3 and O3. These are stored in the EMISRR array when the other NOx emissions are determined. However, when later on the HNO3 and O3 emissions are determined these entries in the EMISRR arrays are overwritten. This way the ship HNO3 and O3 emissions are lost and consequently no ship emissions are present in the model!
| |
|
| |
| :'''Code problem:''' In <tt>emissdr.f</tt> (line258) the fossil fuel emissions are called for all the emitted species.
| |
|
| |
| IF ( LANTHRO ) THEN
| |
| DO N = 1, NEMANTHRO
| |
| NN = IDEMS(N)
| |
|
| |
| IF ( NN /= 0 ) THEN
| |
| CALL EMFOSSIL( I, J, N, NN, IREF, JREF, JSCEN )
| |
| ENDIF
| |
| ENDDO
| |
| ENDIF
| |
|
| |
| :In <tt>emfossil.f</tt> the work array is initialized in line 525:
| |
|
| |
| ! Initialize work variables
| |
| EMX(:) = 0d0
| |
|
| |
| :And then in line 814 the EMISRR(I,J,N) array is overwritten with the work array:
| |
|
| |
| EMISRR(I,J,N) = EMX(1) * XNUMOL(NN) / DTSRCE
| |
|
| |
| :Here the present HNO3 and O3 emissions from the NOx emissions are thus overwritten.
| |
|
| |
| :'''Fix:''' I think the following should fix this, as the EMISRR array is zeroed in every call to <tt>emissdr.f</tt> this should give no problems:
| |
|
| |
| !--------------------------------------------------------------
| |
| ! Store in EMISRR array and archive diagnostics
| |
| !--------------------------------------------------------------
| |
| !--- Prior to (gvinken, 11/16/10). Emissions already present in EMISRR
| |
| ! no longer get overwritten.
| |
| ! EMISRR(I,J,N) = EMX(1) * XNUMOL(NN) / DTSRCE
| |
| EMISRR(I,J,N) = EMISRR(I,J,N) + EMX(1) * XNUMOL(NN) / DTSRCE
| |
|
| |
| --[[User:Bmy|Bob Y.]] 12:00, 16 November 2010 (EST)
| |