Difference between revisions of "Aerosol optical properties"
(→Aerosol optical properties update for GEOS-Chem v8-03-01) |
|||
Line 26: | Line 26: | ||
=== Overview === | === Overview === | ||
− | The updated <tt>jv_spec.dat</tt> file was created to reflect current information about aerosol size distributions for sulfate, organic carbon, black carbon and sea salt. The method follows Martin et al., [2003]. Refractive indices are based on GADS. The Mie code is from ftp://ftp.giss.nasa.gov/pub/crmim/spher.f and described in Mischenko et al. [1999]. The current calculation uses an geometric standard deviation (sigma_g) of 1.6 for sulfate, BC and OC following Wang et al. [2003ab] and Drury et al. [submitted]. The sigma_g was also reduced to 1.5 and 1.8 for fine and coarse sea salt respectively following Jaegle et al. [2010]. Specific geometric mean radii (r_g) for each species and relative humidity are in the headers of jv_spec.dat. Consistent aerosol optical properties at 550nm were also created for the file jv_spec_aod.dat. | + | The updated <tt>jv_spec.dat</tt> file was created to reflect current information about aerosol size distributions for sulfate, organic carbon, black carbon and sea salt. The method follows Martin et al., [2003]. Refractive indices are based on GADS. The Mie code is from ftp://ftp.giss.nasa.gov/pub/crmim/spher.f and described in Mischenko et al. [1999]. The current calculation uses an geometric standard deviation (sigma_g) of 1.6 for sulfate, BC and OC following Wang et al. [2003ab] and Drury et al. [submitted]. The sigma_g was also reduced to 1.5 and 1.8 for fine and coarse sea salt respectively following Jaegle et al. [2010]. Specific geometric mean radii (r_g) for each species and relative humidity are in the headers of jv_spec.dat. Consistent aerosol optical properties at 550nm were also created for the file <tt>jv_spec_aod.dat</tt>. |
Some illustrative differences between the "old" <tt>jv_spec.dat</tt> and the "current" one are given below for 70% relative humidity. The most important change is that organic carbon aerosols are larger in the current version, and that sea salt are substantially smaller. | Some illustrative differences between the "old" <tt>jv_spec.dat</tt> and the "current" one are given below for 70% relative humidity. The most important change is that organic carbon aerosols are larger in the current version, and that sea salt are substantially smaller. | ||
Line 186: | Line 186: | ||
--[[User:Bmy|Bob Y.]] 11:02, 1 March 2010 (EST) | --[[User:Bmy|Bob Y.]] 11:02, 1 March 2010 (EST) | ||
+ | |||
+ | == Known issues == | ||
+ | |||
+ | === Fix to jv_spec_aod.dat === | ||
+ | |||
+ | '''''[mailto:gabriele.curci@aquila.infn.it Gabriele Curci] wrote:''''' | ||
+ | |||
+ | :I'd like to report a bug in the <tt>jv_spec_aod.dat</tt> file in the GC run dirs. It contains aerosol optical properties for calculations of AOD at 550 nm, with data that the user manually replace in the <tt>jv_spec.dat</tt> file. | ||
+ | |||
+ | :The bug is related to dust and consist in missing update of the parameters as done in the <tt>jv_spec.dat</tt> for other wavelenghts. The update was done some time ago: | ||
+ | |||
+ | ::http://wiki.seas.harvard.edu/geos-chem/index.php/Aerosol_optical_properties | ||
+ | |||
+ | :and for dust included updates from Sinyuk et al. (2003): | ||
+ | |||
+ | ::http://www.atmos.colostate.edu/~heald/docs/GEOS_Chem_optics_description.pdf | ||
+ | |||
+ | :Well, the <tt>jv_spec.dat</tt> is ok, but the companion <tt>jv_spec_aod.dat</tt> for 550 nm data was not consistently updated for dust. | ||
+ | |||
+ | The current values for dust: | ||
+ | |||
+ | w(nm) Q r-eff ss-alb pi(0) pi(1) pi(2) pi(3) pi(4) pi(5) pi(6) pi(7) | ||
+ | 15 Mdust 0.15 = mineral dust (R.V.Martin) | ||
+ | 550 1.5356 0.150 0.973 1.000 1.866 1.790 1.223 0.730 0.371 0.187 0.076 | ||
+ | 16 Mdust 0.25 = mineral dust (R.V.Martin) | ||
+ | 550 2.8155 0.250 0.971 1.000 2.029 2.349 2.013 1.596 1.146 0.792 0.505 | ||
+ | 17 Mdust 0.4 = mineral dust (R.V.Martin) | ||
+ | 550 3.0860 0.400 0.956 1.000 1.999 2.546 2.333 2.241 1.897 1.653 1.345 | ||
+ | 18 Mdust 0.8 = mineral dust (R.V.Martin) | ||
+ | 550 2.5841 0.800 0.906 1.000 2.065 2.911 2.965 3.513 3.476 3.796 3.692 | ||
+ | 19 Mdust 1.5 = mineral dust (R.V.Martin) | ||
+ | 550 2.3461 1.500 0.848 1.000 2.298 3.352 3.751 4.627 4.954 5.714 5.983 | ||
+ | 20 Mdust 2.5 = mineral dust (R.V.Martin) | ||
+ | 550 2.2405 2.500 0.788 1.000 2.453 3.638 4.308 5.354 5.998 7.011 7.618 | ||
+ | 21 Mdust 4.0 = mineral dust (R.V.Martin) | ||
+ | 550 2.1744 4.000 0.724 1.000 2.571 3.885 4.822 6.014 6.956 8.157 9.073 | ||
+ | |||
+ | Should be changed to: | ||
+ | |||
+ | w(nm) Q r-eff ss-alb pi(0) pi(1) pi(2) pi(3) pi(4) pi(5) pi(6) pi(7) | ||
+ | 15 Mdust 0.15 = mineral dust (R.V.Martin) | ||
+ | 550 1.5736 0.151 0.993 1.000 1.766 1.248 0.460 0.143 0.027 0.003 0.000 | ||
+ | 16 Mdust 0.25 = mineral dust (R.V.Martin) | ||
+ | 550 3.6240 0.253 0.994 1.000 2.101 2.473 2.097 1.416 0.769 0.355 0.076 | ||
+ | 17 Mdust 0.4 = mineral dust (R.V.Martin) | ||
+ | 550 3.6424 0.402 0.990 1.000 2.005 2.702 2.636 2.714 2.500 2.179 1.787 | ||
+ | 18 Mdust 0.8 = mineral dust (R.V.Martin) | ||
+ | 550 2.6226 0.818 0.971 1.000 2.142 3.213 3.457 4.381 4.388 4.957 4.747 | ||
+ | 19 Mdust 1.5 = mineral dust (R.V.Martin) | ||
+ | 550 2.3682 1.491 0.953 1.000 2.256 3.339 3.617 4.635 4.902 5.843 6.053 | ||
+ | 20 Mdust 2.5 = mineral dust (R.V.Martin) | ||
+ | 550 2.2699 2.417 0.930 1.000 2.357 3.499 3.897 4.971 5.373 6.463 6.862 | ||
+ | 21 Mdust 4.0 = mineral dust (R.V.Martin) | ||
+ | 550 2.1247 3.721 0.897 1.000 2.409 3.562 4.009 5.070 5.556 6.708 7.235 | ||
+ | |||
+ | --[[User:Melissa Payer|Melissa Payer]] 17:11, 8 March 2012 (EST) | ||
== References == | == References == |
Revision as of 22:11, 8 March 2012
This page describes the aerosol optical properties that are used as input to the FAST-J photolysis mechanism in GEOS-Chem.
Contents
Aerosol optical properties update for GEOS-Chem v8-03-01
Starting with GEOS-Chem v8-03-01, several of the aerosol optical properties for the FAST-J photolysis code (specified in the input file jv_spec.dat have been updated to take into account new observations from various datasets. Colette Heald has created an overview document which describes the implementation of the aerosol optics for FAST-J.
NOTE: For this update there are no modifications to the source code itself, just to the jv_spec.dat file.
Authors
Name | Affiliation | Date |
---|---|---|
Randall Martin | Dalhousie University | 19 November 2009 |
Colette Heald | Colorado State University | 29 January 2010 |
Overview
The updated jv_spec.dat file was created to reflect current information about aerosol size distributions for sulfate, organic carbon, black carbon and sea salt. The method follows Martin et al., [2003]. Refractive indices are based on GADS. The Mie code is from ftp://ftp.giss.nasa.gov/pub/crmim/spher.f and described in Mischenko et al. [1999]. The current calculation uses an geometric standard deviation (sigma_g) of 1.6 for sulfate, BC and OC following Wang et al. [2003ab] and Drury et al. [submitted]. The sigma_g was also reduced to 1.5 and 1.8 for fine and coarse sea salt respectively following Jaegle et al. [2010]. Specific geometric mean radii (r_g) for each species and relative humidity are in the headers of jv_spec.dat. Consistent aerosol optical properties at 550nm were also created for the file jv_spec_aod.dat.
Some illustrative differences between the "old" jv_spec.dat and the "current" one are given below for 70% relative humidity. The most important change is that organic carbon aerosols are larger in the current version, and that sea salt are substantially smaller.
Sulfate | Old Value | New Value |
---|---|---|
r_g (um) | 0.07 | 0.11 |
sigma_g | 2.0 | 1.6 |
r_eff (um) | 0.24 | 0.19 |
Qext (at 300nm) | 2.4 | 2.6 |
Qext (at 600nm) | 1.3 | 0.93 |
Organic Carbon | Old Value | New Value |
r_g (um) | 0.03 | 0.09 |
sigma_g | 2.0 | 1.6 |
r_eff (um) | 0.095 | 0.15 |
Qext (at 300nm) | 1.2 | 2.3 |
Qext (at 600nm) | 0.34 | 0.75 |
Black Carbon | Old value | New value |
r_g (um) | 0.01 | 0.02 |
sigma_g | 2.0 | 1.6 |
r_eff (um) | 0.039 | 0.035 |
Qext (at 300nm) | 1.06 | 0.95 |
Qext (at 600nm) | 0.43 | 0.33 |
Sea salt (accumulation mode) | Old value | New value |
r_g (um) | 0.38 | 0.15 |
sigma_g | 2.0 | 1.5 |
r_eff (um) | 1.32 | 0.23 |
Qext (at 300nm) | 2.33 | 3.04 |
Qext (at 600nm) | 2.60 | 1.37 |
Sea salt (coarse mode) | Old value | New value |
r_g (um) | 3.2 | 0.73 |
sigma_g | 2.0 | 1.8 |
r_eff (um) | 10.1 | 1.7 |
Qext (at 300nm) | 2.07 | 2.22 |
Qext (at 600nm) | 2.11 | 2.43 |
The updated jv_spec.dat for v8-02-05 has the following top-of-file header:
jv_spec.dat: FAST-J, FJX (jmao 4/09), Updated AOD's (rvm, clh, 3/10)
Issues
It is likely that organic aerosol at ultraviolet wavelengths is more absorbing than included here as based on GADS [Koepke et al., 1997].
Aerosols are often internally mixed, in contast with their implementation as an external mixture in jv_spec.dat.
--Bob Y. 11:02, 1 March 2010 (EST)
Download location
The updated jv_spec.dat file (and corresponding ratj.d file) are be located in the 1-month benchmark run directory for GEOS-Chem v8-02-05. You may also download them from the following location:
ftp ftp.as.harvard.edu cd pub/geos-chem/data/aerosol_optics/v8-03-01-updates/
These files are now part of the standard GEOS-Chem run directory archive. You may download these with the Git version control software. Please see this wiki post for more information.
--Bob Y. 15:32, 15 November 2010 (EST)
Aerosol optical properties at high spectral resolution
Randall Martin has made available files containing aerosol optical properties at high spectral resolution. These files are not specifically for GEOS-Chem, but they may be useful for some satellite applications or radiative forcing calculations. The files are available from:
ftp ftp.as.harvard.edu cd /pub/geos-chem/data/aerosol_optics/hi_spectral_res/
For more information about the data, please see this README file.
--Bob Y. 11:02, 1 March 2010 (EST)
Known issues
Fix to jv_spec_aod.dat
Gabriele Curci wrote:
- I'd like to report a bug in the jv_spec_aod.dat file in the GC run dirs. It contains aerosol optical properties for calculations of AOD at 550 nm, with data that the user manually replace in the jv_spec.dat file.
- The bug is related to dust and consist in missing update of the parameters as done in the jv_spec.dat for other wavelenghts. The update was done some time ago:
- and for dust included updates from Sinyuk et al. (2003):
- Well, the jv_spec.dat is ok, but the companion jv_spec_aod.dat for 550 nm data was not consistently updated for dust.
The current values for dust:
w(nm) Q r-eff ss-alb pi(0) pi(1) pi(2) pi(3) pi(4) pi(5) pi(6) pi(7) 15 Mdust 0.15 = mineral dust (R.V.Martin) 550 1.5356 0.150 0.973 1.000 1.866 1.790 1.223 0.730 0.371 0.187 0.076 16 Mdust 0.25 = mineral dust (R.V.Martin) 550 2.8155 0.250 0.971 1.000 2.029 2.349 2.013 1.596 1.146 0.792 0.505 17 Mdust 0.4 = mineral dust (R.V.Martin) 550 3.0860 0.400 0.956 1.000 1.999 2.546 2.333 2.241 1.897 1.653 1.345 18 Mdust 0.8 = mineral dust (R.V.Martin) 550 2.5841 0.800 0.906 1.000 2.065 2.911 2.965 3.513 3.476 3.796 3.692 19 Mdust 1.5 = mineral dust (R.V.Martin) 550 2.3461 1.500 0.848 1.000 2.298 3.352 3.751 4.627 4.954 5.714 5.983 20 Mdust 2.5 = mineral dust (R.V.Martin) 550 2.2405 2.500 0.788 1.000 2.453 3.638 4.308 5.354 5.998 7.011 7.618 21 Mdust 4.0 = mineral dust (R.V.Martin) 550 2.1744 4.000 0.724 1.000 2.571 3.885 4.822 6.014 6.956 8.157 9.073
Should be changed to:
w(nm) Q r-eff ss-alb pi(0) pi(1) pi(2) pi(3) pi(4) pi(5) pi(6) pi(7) 15 Mdust 0.15 = mineral dust (R.V.Martin) 550 1.5736 0.151 0.993 1.000 1.766 1.248 0.460 0.143 0.027 0.003 0.000 16 Mdust 0.25 = mineral dust (R.V.Martin) 550 3.6240 0.253 0.994 1.000 2.101 2.473 2.097 1.416 0.769 0.355 0.076 17 Mdust 0.4 = mineral dust (R.V.Martin) 550 3.6424 0.402 0.990 1.000 2.005 2.702 2.636 2.714 2.500 2.179 1.787 18 Mdust 0.8 = mineral dust (R.V.Martin) 550 2.6226 0.818 0.971 1.000 2.142 3.213 3.457 4.381 4.388 4.957 4.747 19 Mdust 1.5 = mineral dust (R.V.Martin) 550 2.3682 1.491 0.953 1.000 2.256 3.339 3.617 4.635 4.902 5.843 6.053 20 Mdust 2.5 = mineral dust (R.V.Martin) 550 2.2699 2.417 0.930 1.000 2.357 3.499 3.897 4.971 5.373 6.463 6.862 21 Mdust 4.0 = mineral dust (R.V.Martin) 550 2.1247 3.721 0.897 1.000 2.409 3.562 4.009 5.070 5.556 6.708 7.235
--Melissa Payer 17:11, 8 March 2012 (EST)
References
- Drury, E., D.J. Jacob, R.J.D. Spurr, J. Wang, Y. Shinozuka, B.E. Anderson, A.D. Clarke, J. Dibb, C. McNaughton, and R. Weber, Synthesis of satellite (MODIS), aircraft (ICARTT), and surface (IMPROVE, EPA-AQS, AERONET) aerosol observations over North America to improve MODIS aerosol retrievals and constrain surface aerosol concentrations and sources , J. Geophys. Res., submitted.
- Jaegle et al. Global sea salt emissions: New constraints from in situ, AERONET, and MODIS observations, in preparation for submission to Atm. Chem. Phys. Discuss., 2010.
- Koepke, P., M. Hess, I. Schult, and E. P. Shettle, Global Aerosol Data Set, Report No. 243, Max-Planck-Institut fur Meteorologie, Hamburg, ISSN 0937-1060, 1997.
- Martin, R.V., D.J. Jacob, R.M. Yantosca, M. Chin, and P. Ginoux, Global and regional decreases in tropospheric oxidants from photochemical effects of aerosols, J. Geophys. Res., 108, 4097, doi:10.1029/2002JD002622, 2003.
- Mishchenko, M.I., J.M. Dlugach, E.G. Yanovitskij, and N.T. Zakharova, Bidirectional reflectance of flat optically thick particulate layers: an efficient radiative transfer solution and applications to snow and soil surfaces, J. Quant. Spectrosc. Radiat. Transfer, 63, 409-432, 1999.
- Wang, J., S.A. Christopher, F. Brechtel, J. Kim, B. Schmid, J. Redemann, P.B. Russell, P. Quinn, and B.N. Holben, Geostationary satellite rtrievals of aerosol optical thickness during ACE-Asia, J. Geophys. Res., 108, 8657, doi:10.1029/2003JD003580, 2003.
- Wang, J., S.A. Christopher, J.S. Reid, H. Maring, D. Savoie, B.H. Holben, J.M. Livingston, P.B. Russell, and S.K. Yang, GOES-8 retrieval of dust aerosol optical thickness over the Atlantic Ocean during PRIDE, J. Geophys. Res., 108, 8595, doi:10.1029/2002JD002494, 2003.
--Bob Y. 11:02, 1 March 2010 (EST)