2010-06-07: Mie Calculations Verification
Insert non-formatted text hereMie 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.
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.
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
Output from Mie Tab 4.1
[ scattering angle considered to be about 150 degree]
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
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 relative intensity given by MieTab 4.1 for almost similar scattering angles, when everything else ( i.e aerosol model, wavelength etc) remained identical.
The reason for this (from my investigations!!!) is that the relative intensity output by MieTab 4.1 is normalized with 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.