Skip to main content

Ocean Algorithm

Introduction   

The history, theory and strategy of the aerosol retrieval over the ocean is described in in Tanré et al., (1997), Levy et al., (2003), Remer et al., (2005) and the C6 updates are described in Levy et al. (2013). The mechanics of the ocean algorithm are shown in the flowchart below. The algorithm is based on a ‘look-up table’ (LUT) approach, i.e., radiative transfer calculations are pre-computed for a set of aerosol and surface parameters and compared with the observed radiation field. The algorithm assumes that one fine and one coarse lognormal aerosol modes can be combined with proper weightings to represent the ambient aerosol properties over the target. Spectral reflectance from the LUT is compared with MODIS-measured spectral reflectance to find the ‘best’ (least-squares) fit. This best fit, or an ‘average’ of a set of the best fits is the solution to the inversion.

Below is a flow chart of the MODIS dark target ocean algorithm for the 10 Km product.  The algorithm is also explained below the flow chart in text which is largely excerpted from the MODIS Algorithm Theoretical Basis Document (ATBD).

 

 

Pixel Masking and Selection

The algorithm has the arduous task of separating 'good' pixels from 'cloudy' pixels.  The standard M?D35 cloud mask includes using the brightness in the visible channels to identify clouds.  This procedure will mistake heavy aerosol as 'cloudy', and miss retrieving important aerosol events over ocean.  On the other hand, relying on IR-tests alone permits low altitude, warm clouds to escape and be misidentified as 'clear', introducing cloud contamination in the aerosol products. Thus, our cloud mask over ocean combines spatial variability tests (e.g. Martins et al., 2002) along with tests of brightness in visible and infrared channels. Underwater sediments have proved to be a problem in shallow water (near coastlines) as the sediments can easily have land-like surface properties. Thus, the sediment mask of Li et al. (2003) is used in addition to the cloud mask.

The algorithm sorts the remaining pixels that have evaded all the cloud masks and the sediment mask according to their r0.86 value, discards the darkest and brightest 25%, and thereby leaves the middle 50% of the data.  The filter is used to eliminate residual cloud contamination, cloud shadows, or other unusual extreme conditions in the box. For the 10 km retrieval, at least 10 of the 400 pixels in the original box must remain after masking and filtering. If the minimum number of pixels is not met, no retrieval is attempted and fill values are given for the aerosol products in the retrieval box.

The ocean algorithm was designed to retrieve only over dark ocean, (i.e. away from glint). Glint angles of < 40˚ are masked. There is a special case where heavy dust can be retreived over glint, and the retrieval will be assigned a poor quality.

 

Ocean Aerosol Models

The ocean algorithm aerosol models were derived mainly from data gleaned from sun/sky photometers and from analysis of errors in the products from previous versions of the MODIS algorithm. In the C6 ocean algorithm there are four fine modes and five coarse modes described in the following tables.

Table 1: Refractive Indices, Number Median, Standard Deviation and Effective Radius for the aerosol modes used in the MODIS Lookup Table for the ocean algorithm.  Models 1-4 are fine modes and Models 5-9 are coarse modes.

F

λ=0.466-->0.857 µm

λ=1.241 µm

 λ =1.628 µm

 λ =2.113 µm

rg

σ

reff

Comments

1

1.45-0.0035i

1.45-0.0035i

1.43-0.01i

1.40-0.005i

0.07

0.40

0.10

Water Soluble

2

1.45-0.0035i

1.45-0.0035i

1.43-0.01i

1.40-0.005i

0.06

0.60

0.15

Water Soluble

3

1.40-0.0020i

1.40-0.0020i

1.39-0.005i

1.36-0.003i

0.08

0.60

0.20

Water Soluble with humidity

4

1.40-0.0020i

1.40-0.0020i

1.39-0.005i

1.36-0.003i

0.10

0.60

0.25

Water Soluble with humidity

C

λ=0.466-->0.857 µm

λ=1.241 µm

  λ =1.628 µm

 λ =2.113 µm

rg

σ

reff

Comments

5

1.35-0.001i

1.35-0.001i

1.35-0.001i

1.35-0.001i

0.40

0.60

0.98

Wet sea salt type

6

1.35-0.001i

1.35-0.001i

1.35-0.001i

1.35-0.001i

0.60

0.60

1.48

Wet sea salt type

7

1.35-0.001i

1.35-0.001i

1.35-0.001i

1.35-0.001i

0.80

0.60

1.98

Wet sea salt type

8

1.53-0.003i (0.47)

1.53-0.001i (0.55)

1.53-0.000i (0.65)

1.53-0.000i (0.86)

1.46-0.000i

1.46-0.001i

1.46-0.000i

0.60

0.60

1.48

Dust-like type

 

9

1.53-0.003i (0.47)

1.53-0.001i (0.55)

1.53-0.000i (0.65)

1.53-0.000i (0.86)

1.46-0.000i

1.46-0.001i

1.46-0.000i

0.50

0.80

2.50

Dust-like type

 

Table 2: Spectral Extinction Coefficients for the aerosol modes used in the MODIS Lookup Table for the ocean algorithm.  Models 1-4 are fine modes and Models 5-9 are coarse modes.

λ (µm) /  Mode

0.466

0.554

0.645

0.857

1.241

1.628

2.113

1

1.43E-10

9.33E-11

6.15E-11

2.66E-11

7.91E-12

4.30E-12

1.48E-12

2

3.03E-10

2.33E-10

1.78E-10

9.95E-11

3.93E-11

1.88E-11

6.99E-12

3

6.78E-10

5.45E-10

4.34E-10

2.63E-10

1.15E-10

5.67E-11

2.28E-11

4

1.33E-09

1.12E-09

9.36E-10

6.15E-10

3.01E-10

1.57E-10

6.68E-11

5

2.69E-08

2.78E-08

2.84E-08

2.85E-08

2.55E-08

2.12E-08

1.63E-08

6

5.57E-08

5.76E-08

5.95E-08

6.29E-08

6.44E-08

6.09E-08

5.33E-08

7

9.50E-08

9.72E-08

9.97E-08

1.06E-07

1.13E-07

1.15E-07

1.09E-07

8

5.57E-08

9.72E-08

5.70E-08

6.05E-08

6.60E-08

6.63E-08

6.26E-08

9

6.42E-08

6.54E-08

6.66E-08

6.92E-08

7.31E-08

7.43E-08

7.36E-08

 

Table 3: Spectral Single Scattering Albedos for the aerosol modes used in the MODIS Lookup Table for the ocean algorithm.  Models 1-4 are fine modes and Models 5-9 are coarse modes.

λ (µm) /  Mode

0.466

0.554

0.645

0.857

1.241

1.628

2.113

1

0.974

0.968

0.961

0.94

0.879

0.541

0.499

2

0.978

0.977

0.976

0.97

0.956

0.817

0.822

3

0.987

0.986

0.986

0.984

0.978

0.921

0.916

4

0.986

0.987

0.987

0.985

0.982

0.94

0.941

5

0.978

0.982

0.985

0.989

0.991

0.992

0.993

6

0.966

0.972

0.976

0.983

0.988

0.991

0.992

7

0.955

0.962

0.9673

0.976

0.984

0.988

0.99

8

0.901

0.967

1

1

1

0.983

1

9

0.867

0.953

1

1

1

0.9835

1

Table 4: Spectral Asymmetry Parameters for the aerosol modes used in the MODIS Lookup Table for the ocean algorithm.  Models 1-4 are fine modes and Models 5-9 are coarse modes.

λ (µm) /  Mode

0.466

0.554

0.645

0.857

1.241

1.628

2.113

1

0.576

0.511

0.447

0.322

0.178

0.105

0.063

2

0.683

0.660

0.635

0.575

0.468

0.369

0.265

3

0.735

0.718

0.699

0.651

0.559

0.472

0.372

4

0.751

0.740

0.726

0.690

0.618

0.546

0.458

5

0.785

0.786

0.789

0.794

0.795

0.787

0.769

6

0.795

0.788

0.786

0.787

0.794

0.796

0.792

7

0.810

0.800

0.793

0.786

0.788

0.794

0.796

8

0.753

0.72

0.697

0.679

0.713

0.72

0.719

9

0.780

0.746

0.723

0.706

0.722

0.722

0.715