˙

         

 

     
   

MainDescription of the codesComparisonBenchmarkEffects of polarizationContacts

 

       
   

The comparison was performed separately for molecular, aerosol and mixed (molecular + aerosol) atmospheres. For a molecular atmosphere, TOA (top-of-atmosphere) reflectances produced by the RT codes were compared to Coulson's tabulated values. In two other cases, the 6SV1 simulations were compared to those produced by the Monte Carlo code and then used as benchmarks for other codes.

(please double click on each item with an asterisk to see or hide the contents)

 

   

1. Molecular atmosphere

  • Conditions of the comparison (6SV1, RT3, MODTRAN, & SHARM vs. Coulson's)*

      Table 1. Conditions of the comparison for a pure molecular atmosphere.

      Optical thickness

      t = 0.1 (l = 530 nm)

       t = 0.25 (l = 440 nm)

      t = 0.5 (l = 360 nm)

      Ground reflectance

      r = 0.0

      r = 0.25

      r = 0.0

      r = 0.25

      r = 0.0

      r = 0.25

      SZA

      23.0739
      53.1301
      78.4630

      0.0
      36.8699
      66.4218

      0.0
      36.8699
      66.4218

      23.0739
      53.1301
      78.4630

      23.0739
      53.1301
      78.4630

      0.0
      36.8699
      66.4218

      AZ

      0.0; 90.0; 180

      VZA

      as in Coulson-s tables (see the Excel file)

      Here SZA denotes the solar zenith angle, VZA is the view zenith angle, and AZ is the relative azimuth.

      The relative differences between the simulations of each code and Coulson's values are calculated as |Coulson-Code|*100%/Coulson.

  • Results of the comparison (6SV1, RT3, MODTRAN, & SHARM vs. Coulson's)*

      Figure 1. Some results of the comparison with Coulson's tabulated values for a molecular atmosphere.

       

      Download the Excel file with the results. An Excel file with VPD results is available by request.

    (Note: VPD was dropped from the comparison after the completion of the molecular case.)

  • Conditions of the comparison (6SV1 vs. Monte Carlo)*

      tmol = 0.25, rsurf = {0.0; 0.25}, SZA = {0.0-; 23.0-; 50.0-}, 1010 photons.
  • Results of the comparison (6SV1 vs. Monte Carlo)*

      Figure 2. Results of the comparison of 6SV1 and Monte Carlo for a pure molecular atmosphere.

      Instead of calculating an exact value of TOA reflectance for a given geometric configuration, Monte Carlo calculates the average value of all TOA reflectances confined within a given solid angle. The hemispherical space at the top of atmosphere is divided into a number of solid angles, specified by VZA and AZ values. In this case, the angular sampling for the relative AZ space (0- - 180-) is 22.5-. The sampling for VZA space (0- - 90-) is shown in Figure 2. The boundary VZA values are presented as angular coordinates. The relative difference varies as the radius coordinate from 0 to 0.4%. A special integration method was applied to match the outputs of 6SV1 and Monte Carlo.

      Download the Excel file with the results.

 
   

2. Aerosol atmosphere

  • Conditions of the comparison (RT3, MODTRAN, & SHARM vs. 6SV1)*

      Table 2. Parameters of the aerosol models used in the comparison.

      Model (short name)

      Urban-industrial and mixed
      (UI)

      Biomass burning
      (AF)

      Biomass burning
      (AS)

      Location

      GSFC,
      Greenbelt, MD(1993-2000)

      Amazonian forest,
      Brazil (1993-1994),
      Bolivia (1998-1999)

      African savanna, Zambia (1995-2000)

      Range of optical thickness t; <t>

      0.1 ≤ t(440) ≤ 1.0
      <
      t(440)> = 0.24

      0.1 ≤ t(440) ≤ 3.0
      <
      t(440)> = 0.74

      0.1 ≤ t(440) ≤ 1.5
      <
      t(440)> = 0.38

      Values of t selected for this study

      0.2; 0.8

      0.2; 0.8; 2.0

      0.2; 0.8

      Real and imaginary parts of refractive index

      1.41-0.03t(440); 0.003

      1.47;
      0.0093

      1.51;
      0.021

      SSA (412/440/670)

      0.97/0.98/0.97

      0.94/0.94/0.93

      0.88/0.88/0.84

      rVf  (mm)

      0.12+0.11t(440)

      0.14+0.013t(440)

      0.12+0.025t(440)

      sf (mm)

      0.38

      0.40

      0.40

      rVc  (mm)

      3.03+0.49t(440)

      3.27+0.58t(440)

      3.22+0.71t(440)

      sc (mm)

      0.75

      0.79

      0.73

      CVf  (mm3/mm2)

      0.15t(440)

      0.12t(440)

      0.12t(440)

      CVc  (mm3/mm2)

      0.01+0.04t(440)

      0.05t(440)

      0.09t(440)

      Reference: O. Dubovik, B. Holben, T. F. Eck, A. Smirnov, Y. J. Kaufman, M. D. King, D. Tanr-, & I. Slutsker, Variability of absorption and optical properties of key aerosol types observed in worldwide locations, J. Atmos. Sci. 59, 590-608 (2002).

      Table 3. Sets of angles for the comparison.

      VZA, deg. (view zenith angle)

      SZA, deg.
      (solar zenith angle)

      AZ, deg.
      (relative azimuth)

      t = 0.2

      t = 0.8

      t = 2.0

      0.0000

      11.4783

      16.2602

      23.0739

      32.8599

      43.9455

      50.2082

      58.6677

      66.4218

      71.3371

      73.7398

      78.4630

      80.7931

      84.2608

      86.5602

      88.8540

      5.0

      10.0

      15.0

      20.0

      25.0

      30.0

      35.0

      40.0

      45.0

      50.0

      55.0

      60.0

      65.0

      70.0

      75.0

      80.0

      3.0

      8.0

      13.0

      18.0

      23.0

      28.0

      33.0

      38.0

      43.0

      48.0

      53.0

      58.0

      63.0

      68.0

      73.0

      78.0

      0.0

      10.0

      23.0709

      45.0

      58.6677

      75.0

      0.0

      90.0

      180.0

  • Results of the comparison (RT3, MODTRAN, & SHARM vs. 6SV1)*

    (Note: In addition to the Excel files with the results, we also provide 6SV1 Mie-files with the calculated phase function values (251 Gaussian angles) for each value of optical thickness.)

  • Conditions of the comparison (6SV1 vs. Monte Carlo)*

      UI: l = 412 nm, taer = {0.2; 0.8}, SZA = {0.0-; 23.0-; 50.0-}, 2 x 1010 photons

      AF: l = 670 nm, taer = {0.2; 0.8}, SZA = {0.0-; 23.0-; 50.0-}, 2 x 1010 photons

      AS: l = 440 nm, taer = {0.2; 0.8}, SZA = {0.0-; 23.0-; 50.0-}, 2 x 1010 photons

  • Results of the comparison (6SV1 vs. Monte Carlo)*

      Figure 4. The results of the comparison of 6SV1 and Monte Carlo for the UI aerosol model.

      The MC geometric configuration is explained in captions for Fig. 2 (Double-click on 'Results of the comparison (6SV1 vs. Monte Carlo)' in the 'Molecular atmosphere' section.

      UI: Results

      AF: Results

      AS: Results

       

       

       

 
   

 

3. Mixed atmosphere

  • Conditions of the comparison (RT3, MODTRAN, & SHARM vs. 6SV1)*
       To create a mixed atmosphere, a molecular US62 atmosphere was added to each aerosol model. The molecular optical thickness depended on the wavelength:

      l = 412 nm, tmol = 0.30319;

      l = 440 nm, tmol = 0.2322;

      l = 670 nm, tmol = 0.04172.

      The aerosol profile was exponential with the scale height of 2 m for all codes except MODTRAN. The molecular profile was exponential with the scale height of approximately 8 m. The surface was simulated as Lambertian with the reflectance rsurf = {0.0; 0.05}.

  • Results of the comparison (RT3, MODTRAN, & SHARM vs. 6SV1)*

    (Note:  )

      Figure 5. Some results of the comparison of RT3, MODTRAN, and SHARM for a mixed atmosphere with the AS aerosol constituent, bounded by black surface.

      Files are organized depending on the aerosol constituent:

      UI: Results (rsurf = 0.0); Results (rsurf = 0.05).

      AF: Results (rsurf = 0.0); Results (rsurf = 0.05).

      AS: Results (rsurf = 0.0); Results (rsurf = 0.05).

  • Conditions of the comparison (6SV1 vs. Monte Carlo)*

      To create a mixed atmosphere, a molecular US62 atmosphere was added to each aerosol model. The molecular optical thickness depended on the wavelength:

      l = 412 nm, tmol = 0.30319;

      l = 440 nm, tmol = 0.2322;

      l = 670 nm, tmol = 0.04172.

      The aerosol profile was exponential with the scale height of 2 m for all codes except MODTRAN. The molecular profile was exponential with the scale height of approximately 8 m. Only black surface cases were simulated for this part of the comparison.

  • Results of the comparison (6SV1 vs. Monte Carlo)*

      Figure 6. Results of the comparison of 6SV1 and Monte Carlo for a mixed atmosphere with the AS aerosol constituent.

      The MC geometric configuration is explained in captions for Fig. 2 (Double-click on 'Results of the comparison (6SV1 vs. Monte Carlo)' in the 'Molecular atmosphere' section.

      UI: Results (rsurf = 0.0)

      AF: Results (rsurf = 0.0)

      AS: Results (rsurf = 0.0)

       

 
   

 

        Official code comparison Web site of the MODIS atmospheric correction group