Difference between revisions of "2010-06-07: Mie Calculations Verification"

From Earth Science Information Partners (ESIP)
m
 
(7 intermediate revisions by the same user not shown)
Line 1: Line 1:
<nowiki>Insert non-formatted text here</nowiki>Mie scattering calculations that were carried out on MODIS reflectance data [ http://wiki.esipfed.org/index.php/2010-05-30:_MODIS_aerosol_reflectance_Mie_simulation] using '''MieTab7.24''' [downloaded from http://sourceforge.net/projects/mietab/ ] are compared here with similar computations done with '''MieTab 4.1''' available from http://capita.wustl.edu/CAPITA/People/RHusar/MIE/ folder.
+
Mie scattering calculations that were carried out on MODIS reflectance data [ http://wiki.esipfed.org/index.php/2010-05-30:_MODIS_aerosol_reflectance_Mie_simulation] using '''MieTab7.24''' [downloaded from http://sourceforge.net/projects/mietab/ ] are compared here with similar computations done with '''MieTab 4.1''' available from http://capita.wustl.edu/CAPITA/People/RHusar/MIE/ folder.
  
 
Before discussing the results of this comparison, I would like to point out some of the basic differences in the computation abilities of MieTab4.1 and MieTab7.24
 
Before discussing the results of this comparison, I would like to point out some of the basic differences in the computation abilities of MieTab4.1 and MieTab7.24
Line 13: Line 13:
 
5) MieTab 7.24 provides Modified Stokes Vector for given input Stokes Vector. Apart from the four elements {I(parallel), I(perpendicular), U, V} of the stokes vector , the out put also include polarization and ellipticity.
 
5) MieTab 7.24 provides Modified Stokes Vector for given input Stokes Vector. Apart from the four elements {I(parallel), I(perpendicular), U, V} of the stokes vector , the out put also include polarization and ellipticity.
  
Almost all of the Mie calculations carried out earlier was re-checked using the MieTab 4.1. Only a sample calculation is provided in this note here.
+
---------------------------------------------------------------
 +
 
 +
In-order to verify the earlier Mie calculations,most of them was re-checked using the MieTab 4.1. However, only a sample calculation is provided in this note here.
  
 
As MieTab 4.1 only provides “cross sections” and MieTab 7.24 only provides “volume cross sections” for PSD models, they could not be compared.
 
As MieTab 4.1 only provides “cross sections” and MieTab 7.24 only provides “volume cross sections” for PSD models, they could not be compared.
Line 30: Line 32:
 
Wavelength= 0.469μm
 
Wavelength= 0.469μm
  
 +
scattering angle of interest is about 150 degree. Therefore the results of Mie computations are provided around that angle.
  
  
 
== Output from Mie Tab 4.1 ==
 
== Output from Mie Tab 4.1 ==
[ scattering angle considered to be about 150 degree]
+
 
 +
 
 +
 
 +
Cross Sections (cm2):
  
 
  Waveln.    Cext      Csca      Cabs      Cbak      M        K   
 
  Waveln.    Cext      Csca      Cabs      Cbak      M        K   
Line 52: Line 58:
 
== Output from Mie Tab 7.24 ==
 
== Output from Mie Tab 7.24 ==
  
Interval   Rmin      Rmax    k(ext)    k(sca)    k(abs)  k(radar)  NRADI CNVG            
+
Number of radius intervals (bins) = 5;
 +
Number of radii used = 45
 +
           
 +
Interval   Rmin      Rmax    k(ext)    k(sca)    k(abs)  k(radar)  NRADI CNVG
 
   1      5.000E-02 1.040E+00 1.932E-06 1.810E-06 1.222E-07 3.069E-06  9  Yes               
 
   1      5.000E-02 1.040E+00 1.932E-06 1.810E-06 1.222E-07 3.069E-06  9  Yes               
 
   2      1.040E+00 2.030E+00 1.735E-06 1.531E-06 2.043E-07 6.325E-07  9  Yes               
 
   2      1.040E+00 2.030E+00 1.735E-06 1.531E-06 2.043E-07 6.325E-07  9  Yes               
Line 58: Line 67:
 
   4      3.020E+00 4.010E+00 2.671E-07 2.086E-07 5.844E-08 8.772E-08  9  Yes               
 
   4      3.020E+00 4.010E+00 2.671E-07 2.086E-07 5.844E-08 8.772E-08  9  Yes               
 
   5      4.010E+00 5.000E+00 1.161E-07 8.649E-08 2.957E-08 2.215E-08  9  Yes
 
   5      4.010E+00 5.000E+00 1.161E-07 8.649E-08 2.957E-08 2.215E-08  9  Yes
 +
 +
Volume Cross Sections (cm-1):
 +
 +
k(ext)=4.725E-06;
 +
k(sca)=4.194E-06;
 +
k(abs)=5.311E-07;
 +
 +
Assymetry g= 7.585E-01
 +
 +
'''F' Matrix'''[a subset of actual output]
 +
 +
Angle  Phase Fn.    M2        M1        S21        D21
 +
148.8  1.331E-02  1.056E-02  1.607E-02 -3.945E-04 -6.492E-03
 +
149.8  1.411E-02  1.132E-02  1.689E-02 -5.330E-04 -7.159E-03
 +
150.8  1.502E-02  1.211E-02  1.793E-02 -6.863E-04 -7.964E-03
 +
151.8  1.608E-02  1.295E-02  1.921E-02 -8.167E-04 -8.928E-03
 +
152.8  1.731E-02  1.387E-02  2.076E-02 -9.141E-04 -1.007E-02
 +
153.9  1.874E-02  1.485E-02  2.262E-02 -9.947E-04 -1.139E-02
 +
 +
 +
It is clear from these outputs that the Phase Function given by MieTab 7.24 is larger (by two orders of magnitude) than the angular relative intensity given by MieTab 4.1 for similar scattering angles,wavelength and aerosol model.
 +
The reason for this (from my investigations!!!) is that the relative intensity output by MieTab 4.1 is normalized to particle number density. The Phase function computed from both the Mie programs matched when particle number density was used as a multiplier to the relative intensity values provided by MieTab 4.1.

Latest revision as of 09:30, June 7, 2010

Mie scattering calculations that were carried out on MODIS reflectance data [ http://wiki.esipfed.org/index.php/2010-05-30:_MODIS_aerosol_reflectance_Mie_simulation] using MieTab7.24 [downloaded from http://sourceforge.net/projects/mietab/ ] are compared here with similar computations done with MieTab 4.1 available from http://capita.wustl.edu/CAPITA/People/RHusar/MIE/ folder.

Before discussing the results of this comparison, I would like to point out some of the basic differences in the computation abilities of MieTab4.1 and MieTab7.24

1) MieTab 7.24 can calculate volume (extinction/scattering/absorption) cross sections in cm-1 for various particle size distribution (PSD) models. MieTab 4.1 only provides cross sections (extinction/scattering/absorption) in square centimeters.

2) MieTab7.24 calculates the four elements (M2, M1, S21 and D21) of F’ Matrix [ notations are similar to that used by Van de Hulst ] at each specified scattering angle and displays them under the generic term "intensities" together with the “phase function”. This facility is not available in MieTab 4.1. Mie Tab 4.1 only provides relative intensities at specified scattering angles.

3) The number of angles (NANG) at which the phase function(s) are calculated can be upto maximum 361 in MieTab 7.24. In MieTab 4.1, NANG is maximum 181.

4) MieTab7.24 has a function called NBIN which can subdivide the overall computational range (RMIN to RMAX) in to number of radius intervals (or “bins”). The maximum number of bins that are allowed in MieTab7.24 is five. MieTab 4.1 doesn’t have this function. NBIN is only used with analytic particle distribution models.

5) MieTab 7.24 provides Modified Stokes Vector for given input Stokes Vector. Apart from the four elements {I(parallel), I(perpendicular), U, V} of the stokes vector , the out put also include polarization and ellipticity.


In-order to verify the earlier Mie calculations,most of them was re-checked using the MieTab 4.1. However, only a sample calculation is provided in this note here.

As MieTab 4.1 only provides “cross sections” and MieTab 7.24 only provides “volume cross sections” for PSD models, they could not be compared.

Therefore it was decided to compare the Phase functions computed by these two Mie programs.

Following are the out puts from these two programs for one of the aerosol models (Model 23) used in the earlier calculations.

Details of the Model:

Refractive Index: 1.55+0.0035i; RMIN= 0.05 μm; RMAX= 5 μm; Rbar= 0.5 μm; Std.Dev (σ )= 0.5; Wavelength= 0.469μm

scattering angle of interest is about 150 degree. Therefore the results of Mie computations are provided around that angle.


Output from Mie Tab 4.1

Cross Sections (cm2):

Waveln.     Cext      Csca      Cabs      Cbak       M         K  
--------- --------- --------- --------- --------- --------- ----------
4.690E-01 5.258E-06 4.686E-06 5.714E-07 4.730E-06 1.550E+00 3.500E-03


Variation of scattered intensity with specified angle [a subset of actual output]:

 Angle(deg) Relative Intensity                                                           
 --------   ------------------
 148.83     1.317E-04
 149.83     1.305E-04
 150.84     1.279E-04
 151.84     1.246E-04
 152.85     1.237E-04
 153.85     1.307E-04

Output from Mie Tab 7.24

Number of radius intervals (bins) = 5; Number of radii used = 45

Interval   Rmin      Rmax     k(ext)    k(sca)    k(abs)  k(radar)  NRADI  CNVG
 1      5.000E-02 1.040E+00 1.932E-06 1.810E-06 1.222E-07 3.069E-06   9   Yes              
 2      1.040E+00 2.030E+00 1.735E-06 1.531E-06 2.043E-07 6.325E-07   9   Yes              
 3      2.030E+00 3.020E+00 6.749E-07 5.583E-07 1.166E-07 3.686E-07   9   Yes              
 4      3.020E+00 4.010E+00 2.671E-07 2.086E-07 5.844E-08 8.772E-08   9   Yes              
 5      4.010E+00 5.000E+00 1.161E-07 8.649E-08 2.957E-08 2.215E-08   9   Yes

Volume Cross Sections (cm-1):

k(ext)=4.725E-06; k(sca)=4.194E-06; k(abs)=5.311E-07;

Assymetry g= 7.585E-01

F' Matrix[a subset of actual output]

Angle   Phase Fn.     M2         M1        S21        D21
148.8   1.331E-02   1.056E-02   1.607E-02 -3.945E-04 -6.492E-03
149.8   1.411E-02   1.132E-02   1.689E-02 -5.330E-04 -7.159E-03
150.8   1.502E-02   1.211E-02   1.793E-02 -6.863E-04 -7.964E-03
151.8   1.608E-02   1.295E-02   1.921E-02 -8.167E-04 -8.928E-03
152.8   1.731E-02   1.387E-02   2.076E-02 -9.141E-04 -1.007E-02
153.9   1.874E-02   1.485E-02   2.262E-02 -9.947E-04 -1.139E-02


It is clear from these outputs that the Phase Function given by MieTab 7.24 is larger (by two orders of magnitude) than the angular relative intensity given by MieTab 4.1 for similar scattering angles,wavelength and aerosol model. The reason for this (from my investigations!!!) is that the relative intensity output by MieTab 4.1 is normalized to particle number density. The Phase function computed from both the Mie programs matched when particle number density was used as a multiplier to the relative intensity values provided by MieTab 4.1.