Difference between revisions of "Scale factors for anthropogenic emissions"

From Geos-chem
Jump to: navigation, search
(Known issues)
 
(66 intermediate revisions by 6 users not shown)
Line 1: Line 1:
This page describes the methodology of the interannual scaling factors used for anthropogenic emissions in GEOS-Chem. In general, the various emissions inventories are for a particular year (or in some cases, for multi-years).  A set of scale factors is required to adjust the emissions totals to what the
+
This page describes the methodology of the interannual scaling factors used for anthropogenic emissions in GEOS-Chem.
  
 
== Overview ==
 
== Overview ==
  
From [http://acmg.seas.harvard.edu/publications/vanDonkelaar_2008_ACP.pdf ''van Donkelaar et al.'', <nowiki>[2008]</nowiki>]:
+
'''''Philippe Le Sager / Aaron van Donkelaar wrote:'''''
 +
 
 +
<blockquote>Global annual scaling factors are available for NOx, CO, and SOx from 1985-2006. The approach, as implemented by van Donkelaar et al. (ACP, 2008) builds upon the work of Bey et al. (2001) and Park et al. (2004).  Scale factors are automatically applied to any inventory, if needed, to get as close as possible to simulated year conditions.  This behavior can be overwritten.  Scale factors have been derived from numerous sources and directly use government statistics, where available.</blockquote>
 +
 
 +
<blockquote>Trend data for the following regions is derived from:</blockquote>
 +
 
 +
<blockquote>Canada – Environment Canada National Pollutant Release Inventory Trends (http://www.ec.gc.ca/inrp-npri/)</blockquote>
 +
 
 +
<blockquote>United States – Industrial emissions are scaled according to reported trends from the EPA Acid Rain Program (http://camddataandmaps.epa.gov/gdm/).  Other emissions are scaled according to National Emissions Inventory (NEI) Air Pollutant Emissions Trends Data (http://www.epa.gov/ttn/chief/trends/)</blockquote>
 +
 
 +
<blockquote>Europe – European Monitoring and Evaluation Program (EMEP; UNECE (2009) Present state of emission data, ECE/EB.AIR/97, http://www.emep.int/)</blockquote>
 +
 
 +
<blockquote>Asia  – Regional Emission inventory in Asia (REAS; Ohara et al., ACP, 2007; http://www.jamstec.go.jp/frcgc/research/p3/emission.htm)</blockquote>
 +
 
 +
<blockquote>Beyond these regions, NOx, SOx and CO are scaled according to CO2 trends available from the Carbon Dioxide Information Analysis Center (http://cdiac.ornl.gov/).  NOx, SOx and CO are scaled according to trends in total, solid and liquid CO2, respectively.</blockquote>
 +
 
 +
From [http://acmg.seas.harvard.edu/publications/2008/vanDonkelaar_2008_ACP.pdf ''van Donkelaar et al.'', <nowiki>[2008]</nowiki>]:
  
 
<blockquote>We scale all regional and global inventories from their respective base year to 2003, the last year of available statistics, unless its base year is after 2003. Our approach follows ''Bey et al.'' (2001) and ''Park et al.'' (2004). Emissions are scaled according to estimates provided by individual countries, where available. These countries/regions include the United States, Canada, Japan and Europe. NOx emissions of remaining countries are scaled proportional to changes in total CO2 emissions. SOx emissions are similarly scaled to solid fuel CO2 emissions and CO emissions to liquid fuel CO2 emissions. CO2 emission data are obtained from the Carbon Dioxide Information Analysis Center (CDIAC).</blockquote>
 
<blockquote>We scale all regional and global inventories from their respective base year to 2003, the last year of available statistics, unless its base year is after 2003. Our approach follows ''Bey et al.'' (2001) and ''Park et al.'' (2004). Emissions are scaled according to estimates provided by individual countries, where available. These countries/regions include the United States, Canada, Japan and Europe. NOx emissions of remaining countries are scaled proportional to changes in total CO2 emissions. SOx emissions are similarly scaled to solid fuel CO2 emissions and CO emissions to liquid fuel CO2 emissions. CO2 emission data are obtained from the Carbon Dioxide Information Analysis Center (CDIAC).</blockquote>
  
'''''Philippe Le Sager wrote:'''''
+
=== Extension of scale factors to 2010 ===
  
<blockquote>Note that the scale factor have been updated since ''van Donkelaar et al'' [2008] to go up to 2005, and are based on REAS data now for South East Asia: Relative changes in the [http://www.jamstec.go.jp/frsgc/research/d4/reas_h_a.html REAS inventory] ([http://www.atmos-chem-phys.net/7/4419/2007/acp-7-4419-2007.html ''Ohara et al.'', ACP, <nowiki>[2007]</nowiki>]) over East Asia have been used.</blockquote>
+
<span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem_v9-02_benchmark_history#v9-02m|v9-02m]] and approved on 30 Jul 2013. This update is included in Adjoint [[GEOS-Chem_Adjoint_v35|v35d]].'''''</span>
  
<blockquote>This should be a good improvement as REAS emissions are gridded, rather than national scale emissions, giving us much better spatial detail. Also, these scalars are now based on actual NOx, SOx and CO emission estimates, not strictly
+
Annual anthropogenic scale factors will be extended to 2010 for U.S.A. and Canada in v9-02.  Note that these factors remain fixed past 2006 for other regions unless overwritten by other emission inventory data (e.g. EMEP and Streets).
an assumed proportionality between total, solid and liquid CO2 emissions.</blockquote>
+
  
== Validation ==
+
==Seasonal Scale Factors==
  
== References ==
+
<blockquote>Seasonal scale factors are also applied to many anthropogenic emissions:</blockquote>
  
#[http://acmg.seas.harvard.edu/geos/word_pdf_docs/emissions_v8_02_03.pdf Anthropogenic emissions inventories in GEOS-Chem v8-02-03 (and higher versions)]
+
<blockquote>Europe –  When using EMEP, European monthly variation of NOx, SOx, CO and NH3 is based upon input into the Open Source Unified EMEP model (http://www.emep.int/OpenSource/), that was produced by the GENEMIS project, coordinated by the Institute of Energy Economics and the Rational Use of Energy (IER) at the University of Stuttgart.  These source-specific, country-level variations are combined with sector-level gridded EMEP emissions (www.emep.int) to estimate an overall species specific variation.</blockquote>
#[http://wiki.seas.harvard.edu/geos-chem/images/Anth_emiss_gc_v8-01-04.pdf Anthropogenic emissions inventories used in GEOS-Chem v8-01-04 thru v8-02-02]
+
#[http://acmg.seas.harvard.edu/geos/word_pdf_docs/INVENTORY_GC_public-pre_v8.pdf Anthropogenic emissions inventories prior to GEOS-Chem v8-01-01]
+
#Ohara, T., H. Akimoto, J. Kurokawa, N. Horii, K. Yamaji, X. Yan, and T. Hayasaka, ''An Asian emission inventory of anthropogenic emission sources for the period 1980–2020'', <u>Atmos. Chem. Phys.,</u>, '''7''', 4419-4444, 2007. [http://www.atmos-chem-phys.net/7/4419/2007/ ACP site]
+
#van Donkelaar, A., R. V. Martin, W. R. Leaitch, A.M. Macdonald, T. W. Walker, D. G. Streets, Q. Zhang, E. J. Dunlea, J. L. Jimenez, J. E. Dibb, L. G. Huey, R. Weber, and M. O. Andreae, ''Analysis of Aircraft and Satellite Measurements from the Intercontinental Chemical Transport Experiment (INTEX-B) to Quantify Long-Range Transport of East Asian Sulfur to Canada'', <u>Atmos. Chem. Phys.</u>, 8, 2999-3014, 2008. [http://acmg.seas.harvard.edu/publications/vanDonkelaar_2008_ACP.pdf PDF]
+
  
--[[User:Bmy|Bob Y.]] 09:41, 23 February 2010 (EST)
+
<blockquote>North America – Regional NA emission monthly variation of NOx is based upon the VISTAS inventory. All other species, with the exception of NH3, are based upon NEI99 seasonality.  For NEI2005, NH3 seasonality is based upon Park et al., JGR, 2004.  When using VISTAS or NEI2005, a year-varying scalar in additionally included to account for recent ozone regulations.  These are based upon reported trends from the EPA Acid Rain Program (http://camddataandmaps.epa.gov/gdm/).</blockquote>
  
== Known issues ==
+
<blockquote>Asia – When using the Streets inventory, monthly variation is present for NOx and CO only. As of v9-01-02, there is also seasonality in Asian NH3 emissions (described [[Anthropogenic_emissions#Streets_2000|here]]).</blockquote>
  
=== Large jump in emissions in Asia from 2005 to 2006 ===
+
<blockquote>Canada – There is currently no seasonal variation in the CAC inventory, other than ammonia emission which has monthly emissions data.</blockquote>
  
'''''[mailto:jlogan@seas.harvard.edu Jennifer Logan] wrote:'''''
+
<blockquote>Mexico – There is currently no seasonal variation in the BRAVO inventory.</blockquote>
  
:I've inadvertently discovered what seems to be a problem with the time dependent emissions, that likely relates to when it switches from one base inventory to another. For CO, there is a huge difference between emissions in Asia for 2005 and 2006 - see Junhua's plots below:
+
<blockquote>Global – The seasonal variation of both EDGAR (Olivier and Berdowski, 2001) and GEIA is derived from the seasonality of GEIA.</blockquote>
  
*[http://acmg.seas.harvard.edu/geos/wiki_docs/emissions/ff_em_geos4_geos5_SAsia.pdf GEOS-4 vs GEOS-5 emissions comparison: SE Asia only]
+
==Comparison of Emission Configurations==
*[http://acmg.seas.harvard.edu/geos/wiki_docs/emissions/ff_em_geos4_geos5.pdf GEOS-4 vs GEOS-5 emissions comparison: Larger region]
+
  
:I am guessing that this is because 2005 is scaled forward from the Streets inventory for 2000 (his latest for China I hope, not the original TRACE-P one), while 2006 is taken from the new Streets inventory for INTEX-B.
+
<blockquote>The overall impact of the anthropogenic scale factors upon total emissions of NOx, SOx and CO can be greatly affected by the choice of base emission inventory.  The following figures below illustrate the possible range of annual equivalent global and North American emissions for 6 different emission setups, based upon GEOS-Chem v8-02-03 run at 4 x 5 degrees.  Not all setups affects all species (see table), and lines of equivalent emissions are coloured to match.  Note the large impact of switching between Streets 2000 and Streets 2006 from 2000 to 2001 upon global NOx and CO (see discussion below).</blockquote>
  
 +
<blockquote>An additional factor causing this jump is the inclusion of biofuel CO from Streets 2006 in the fossil fuel CO emission diagnostics (ND29).  This inventory does not separate fossil fuel from biofuel emissions and therefore cannot by separated by source.  This should not effect CO concentrations, but has impacted the CO emission diagnositic over the Asian region.</blockquote>
  
 +
<blockquote>[http://fizz.phys.dal.ca/~atmos/EmissionPlot/Global-NOx.jpg Global NOx Comparison]</blockquote>
 +
<blockquote>[http://fizz.phys.dal.ca/~atmos/EmissionPlot/Global-SO2.jpg Global SO2 Comparison]</blockquote>
 +
<blockquote>[http://fizz.phys.dal.ca/~atmos/EmissionPlot/Global-CO.jpg Global CO Comparison]</blockquote>
 +
<blockquote>[http://fizz.phys.dal.ca/~atmos/EmissionPlot/NorthAmerica-NOx.jpg NA NOx Comparison]</blockquote>
 +
<blockquote>[http://fizz.phys.dal.ca/~atmos/EmissionPlot/NorthAmerica-SO2.jpg NA SO2 Comparison]</blockquote>
 +
<blockquote>[http://fizz.phys.dal.ca/~atmos/EmissionPlot/NorthAmerica-CO.jpg NA CO Comparison]</blockquote>
  
:Jennifer
+
<blockquote>Emission within these figures are configured as:</blockquote>
 +
 
 +
<blockquote>
 +
{| border=1 cellpadding=5 cellspacing=0
 +
|-bgcolor="#CCCCCC"
 +
!Descriptor
 +
!Global Inventory
 +
!CAC
 +
!NEI
 +
!VISTAS
 +
!ICARTT
 +
!EMEP
 +
!BRAVO
 +
!Streets
 +
|-align="center"
 +
|"GEIA"
 +
|GEIA
 +
|NO           
 +
|NO
 +
|NO
 +
|NO
 +
|NO
 +
|NO
 +
|NO
 +
|-align="center"
 +
|"EDGAR"
 +
|EDGAR
 +
|NO
 +
|NO
 +
|NO
 +
|NO
 +
|NO
 +
|NO
 +
|NO
 +
|-align="center"
 +
|"RGN-NEI05"
 +
|EDGAR
 +
|YES
 +
|2005
 +
|NO
 +
|NO
 +
|YES
 +
|YES
 +
|YES
 +
|-align="center"
 +
|"RGN-NEI99"
 +
|EDGAR
 +
|YES
 +
|1999
 +
|NO
 +
|NO
 +
|YES
 +
|YES
 +
|YES
 +
|-align="center"
 +
|"RGN-VISTAS"
 +
|EDGAR
 +
|YES
 +
|1999
 +
|YES
 +
|NO
 +
|YES
 +
|YES
 +
|YES
 +
|-align="center"
 +
|"RGN-ICARTT"
 +
|EDGAR
 +
|YES
 +
|1999
 +
|NO
 +
|YES
 +
|YES
 +
|YES
 +
|YES
 +
|}
 +
</blockquote>
 +
 
 +
==Diurnal Variation==
 +
 
 +
<blockquote>A diurnal variation is applied to all NOx. It is derived from EDGAR hourly variations ''sc(k,H)'' for each source ''k'' and hour ''H'', spatially weighted by the sources. In other words:</blockquote>
 +
                           
 +
                            Σ sc(k,H) * NOx(I,J,k)
 +
      ScaleFactor(H,I,J) = --------------------------
 +
                                  Σ NOx(I,J,k)
 +
 
 +
'''''Philippe Le Sager wrote:'''''
 +
 
 +
<blockquote>Note that the scale factors in ''van Donkelaar et al'', 2008] have been updated to go up to 2005, and are based on REAS data now for South East Asia: Relative changes in the [http://www.jamstec.go.jp/frsgc/research/d4/reas_h_a.html REAS inventory] ([http://www.atmos-chem-phys.net/7/4419/2007/acp-7-4419-2007.html ''Ohara et al.'', ACP, <nowiki>[2007]</nowiki>]) over East Asia have been used.</blockquote>
 +
 
 +
<blockquote>This should be a good improvement as REAS emissions are gridded, rather than national scale emissions, giving us much better spatial detail. Also, these scalars are now based on actual NOx, SOx and CO emission estimates, not strictly
 +
an assumed proportionality between total, solid and liquid CO2 emissions.</blockquote>
 +
 
 +
== References ==
 +
 
 +
#[[Anthropogenic emissions]] wiki page
 +
#Bey, I., D. J. Jacob, R. M. Yantosca, J. A. Logan, B. Field, A. M. Fiore, Q. Li, H. Liu, L. J. Mickley, and M. Schultz, ''Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation'', <u>J. Geophys. Res.</u>, '''106''', 23,073–23,096, 2001. [http://acmg.seas.harvard.edu/publications/bey2001a.pdf PDF]
 +
#Ohara, T., H. Akimoto, J. Kurokawa, N. Horii, K. Yamaji, X. Yan, and T. Hayasaka, ''An Asian emission inventory of anthropogenic emission sources for the period 1980–2020'', <u>Atmos. Chem. Phys.,</u>, '''7''', 4419-4444, 2007. [http://www.atmos-chem-phys.net/7/4419/2007/ ACP site]
 +
#Park, R. J., D. J. Jacob, B. D. Field, R. M. Yantosca, and M. Chin, ''Natural and transboundary pollution influences on sulfate-nitrate-ammonium aerosols in the United States: implications for policy'', <u>J. Geophys. Res.</u>, '''109''', D15204, 10.1029/2003JD004473, 2004. [http://acmg.seas.harvard.edu/publications/park2004.pdf PDF]
 +
#van Donkelaar, A., R. V. Martin, W. R. Leaitch, A.M. Macdonald, T. W. Walker, D. G. Streets, Q. Zhang, E. J. Dunlea, J. L. Jimenez, J. E. Dibb, L. G. Huey, R. Weber, and M. O. Andreae, ''Analysis of Aircraft and Satellite Measurements from the Intercontinental Chemical Transport Experiment (INTEX-B) to Quantify Long-Range Transport of East Asian Sulfur to Canada'', <u>Atmos. Chem. Phys.</u>, 8, 2999-3014, 2008. [http://acmg.seas.harvard.edu/publications/2008/vanDonkelaar_2008_ACP.pdf PDF]
 +
 
 +
--[[User:Bmy|Bob Y.]] 09:41, 23 February 2010 (EST)

Latest revision as of 18:21, 24 October 2023

This page describes the methodology of the interannual scaling factors used for anthropogenic emissions in GEOS-Chem.

Overview

Philippe Le Sager / Aaron van Donkelaar wrote:

Global annual scaling factors are available for NOx, CO, and SOx from 1985-2006. The approach, as implemented by van Donkelaar et al. (ACP, 2008) builds upon the work of Bey et al. (2001) and Park et al. (2004). Scale factors are automatically applied to any inventory, if needed, to get as close as possible to simulated year conditions. This behavior can be overwritten. Scale factors have been derived from numerous sources and directly use government statistics, where available.
Trend data for the following regions is derived from:
Canada – Environment Canada National Pollutant Release Inventory Trends (http://www.ec.gc.ca/inrp-npri/)
United States – Industrial emissions are scaled according to reported trends from the EPA Acid Rain Program (http://camddataandmaps.epa.gov/gdm/). Other emissions are scaled according to National Emissions Inventory (NEI) Air Pollutant Emissions Trends Data (http://www.epa.gov/ttn/chief/trends/)
Europe – European Monitoring and Evaluation Program (EMEP; UNECE (2009) Present state of emission data, ECE/EB.AIR/97, http://www.emep.int/)
Asia – Regional Emission inventory in Asia (REAS; Ohara et al., ACP, 2007; http://www.jamstec.go.jp/frcgc/research/p3/emission.htm)
Beyond these regions, NOx, SOx and CO are scaled according to CO2 trends available from the Carbon Dioxide Information Analysis Center (http://cdiac.ornl.gov/). NOx, SOx and CO are scaled according to trends in total, solid and liquid CO2, respectively.

From van Donkelaar et al., [2008]:

We scale all regional and global inventories from their respective base year to 2003, the last year of available statistics, unless its base year is after 2003. Our approach follows Bey et al. (2001) and Park et al. (2004). Emissions are scaled according to estimates provided by individual countries, where available. These countries/regions include the United States, Canada, Japan and Europe. NOx emissions of remaining countries are scaled proportional to changes in total CO2 emissions. SOx emissions are similarly scaled to solid fuel CO2 emissions and CO emissions to liquid fuel CO2 emissions. CO2 emission data are obtained from the Carbon Dioxide Information Analysis Center (CDIAC).

Extension of scale factors to 2010

This update was tested in the 1-month benchmark simulation v9-02m and approved on 30 Jul 2013. This update is included in Adjoint v35d.

Annual anthropogenic scale factors will be extended to 2010 for U.S.A. and Canada in v9-02. Note that these factors remain fixed past 2006 for other regions unless overwritten by other emission inventory data (e.g. EMEP and Streets).

Seasonal Scale Factors

Seasonal scale factors are also applied to many anthropogenic emissions:
Europe – When using EMEP, European monthly variation of NOx, SOx, CO and NH3 is based upon input into the Open Source Unified EMEP model (http://www.emep.int/OpenSource/), that was produced by the GENEMIS project, coordinated by the Institute of Energy Economics and the Rational Use of Energy (IER) at the University of Stuttgart. These source-specific, country-level variations are combined with sector-level gridded EMEP emissions (www.emep.int) to estimate an overall species specific variation.
North America – Regional NA emission monthly variation of NOx is based upon the VISTAS inventory. All other species, with the exception of NH3, are based upon NEI99 seasonality. For NEI2005, NH3 seasonality is based upon Park et al., JGR, 2004. When using VISTAS or NEI2005, a year-varying scalar in additionally included to account for recent ozone regulations. These are based upon reported trends from the EPA Acid Rain Program (http://camddataandmaps.epa.gov/gdm/).
Asia – When using the Streets inventory, monthly variation is present for NOx and CO only. As of v9-01-02, there is also seasonality in Asian NH3 emissions (described here).
Canada – There is currently no seasonal variation in the CAC inventory, other than ammonia emission which has monthly emissions data.
Mexico – There is currently no seasonal variation in the BRAVO inventory.
Global – The seasonal variation of both EDGAR (Olivier and Berdowski, 2001) and GEIA is derived from the seasonality of GEIA.

Comparison of Emission Configurations

The overall impact of the anthropogenic scale factors upon total emissions of NOx, SOx and CO can be greatly affected by the choice of base emission inventory. The following figures below illustrate the possible range of annual equivalent global and North American emissions for 6 different emission setups, based upon GEOS-Chem v8-02-03 run at 4 x 5 degrees. Not all setups affects all species (see table), and lines of equivalent emissions are coloured to match. Note the large impact of switching between Streets 2000 and Streets 2006 from 2000 to 2001 upon global NOx and CO (see discussion below).
An additional factor causing this jump is the inclusion of biofuel CO from Streets 2006 in the fossil fuel CO emission diagnostics (ND29). This inventory does not separate fossil fuel from biofuel emissions and therefore cannot by separated by source. This should not effect CO concentrations, but has impacted the CO emission diagnositic over the Asian region.
Global NOx Comparison
Global SO2 Comparison
Global CO Comparison
NA NOx Comparison
NA SO2 Comparison
NA CO Comparison
Emission within these figures are configured as:
Descriptor Global Inventory CAC NEI VISTAS ICARTT EMEP BRAVO Streets
"GEIA" GEIA NO NO NO NO NO NO NO
"EDGAR" EDGAR NO NO NO NO NO NO NO
"RGN-NEI05" EDGAR YES 2005 NO NO YES YES YES
"RGN-NEI99" EDGAR YES 1999 NO NO YES YES YES
"RGN-VISTAS" EDGAR YES 1999 YES NO YES YES YES
"RGN-ICARTT" EDGAR YES 1999 NO YES YES YES YES

Diurnal Variation

A diurnal variation is applied to all NOx. It is derived from EDGAR hourly variations sc(k,H) for each source k and hour H, spatially weighted by the sources. In other words:
                            Σ sc(k,H) * NOx(I,J,k)
     ScaleFactor(H,I,J) = --------------------------
                                 Σ NOx(I,J,k)

Philippe Le Sager wrote:

Note that the scale factors in van Donkelaar et al, 2008] have been updated to go up to 2005, and are based on REAS data now for South East Asia: Relative changes in the REAS inventory (Ohara et al., ACP, [2007]) over East Asia have been used.
This should be a good improvement as REAS emissions are gridded, rather than national scale emissions, giving us much better spatial detail. Also, these scalars are now based on actual NOx, SOx and CO emission estimates, not strictly an assumed proportionality between total, solid and liquid CO2 emissions.

References

  1. Anthropogenic emissions wiki page
  2. Bey, I., D. J. Jacob, R. M. Yantosca, J. A. Logan, B. Field, A. M. Fiore, Q. Li, H. Liu, L. J. Mickley, and M. Schultz, Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation, J. Geophys. Res., 106, 23,073–23,096, 2001. PDF
  3. Ohara, T., H. Akimoto, J. Kurokawa, N. Horii, K. Yamaji, X. Yan, and T. Hayasaka, An Asian emission inventory of anthropogenic emission sources for the period 1980–2020, Atmos. Chem. Phys.,, 7, 4419-4444, 2007. ACP site
  4. Park, R. J., D. J. Jacob, B. D. Field, R. M. Yantosca, and M. Chin, Natural and transboundary pollution influences on sulfate-nitrate-ammonium aerosols in the United States: implications for policy, J. Geophys. Res., 109, D15204, 10.1029/2003JD004473, 2004. PDF
  5. van Donkelaar, A., R. V. Martin, W. R. Leaitch, A.M. Macdonald, T. W. Walker, D. G. Streets, Q. Zhang, E. J. Dunlea, J. L. Jimenez, J. E. Dibb, L. G. Huey, R. Weber, and M. O. Andreae, Analysis of Aircraft and Satellite Measurements from the Intercontinental Chemical Transport Experiment (INTEX-B) to Quantify Long-Range Transport of East Asian Sulfur to Canada, Atmos. Chem. Phys., 8, 2999-3014, 2008. PDF

--Bob Y. 09:41, 23 February 2010 (EST)