TOC

(1) 
Author(s) 
RP 
Gordon et al (1975) 
0.3244 
Morel and Prieur
(1977) 
0.33 
Kirk (1981) 
0.328 
Sathyendranath and
Platt (1997) 
0.343 
AOMC 
0.329 
In this problem, a simple homogeneous water body with a solely absorbing
( a = 1 m1 ) environment is simulated. The light source is a single
point in the sky with a zenith angle of 60° resulting in a (below water)
refracted angle of 40.5°. The bottom boundary is at a depth of 1 m and
has a perfectly specular reflecting bottom with Rb
= 1.0. This simulated environment prevents the diffusion (due to scattering)
of the ray of light in the aquatic medium. This allows for an analytical
solution of the rate of attenuation of the beam of light in the downward
and upward directions. In the absence of scattering, the following equation
can be used to calculate Ed(z) knowing Ed(z  dz) and the
rate of vertical attenuation of Ed:
(2) 
Since the model computes Ed at different depths, equation (2) can
be solved for k:
(3) 
For this simulation, the only attenuating factor is a therefore it
is expected that k = a for both the upwelling and downwelling
light stream. Results follow:
Attenuation coefficient

a (given)

Calculated from AOMC simulation

kEd

1.0

0.999

kEu

1.0

1.000

3 Simulation of a Spectral Pivot Point
Natural waters contain various constituents with differing inherent optical
properties (IOP). It is this difference that allows each constituent within
a water body to acquire a unique radiometric spectral signature such as
surface irradiance reflectance (R). The spectral signature varies
as a function of concentration of the constituent. For example, Bukata
et al. (1995) show how for a twocomponent water body (water and chlorophyll),
there exists a wavelength at which the bulk subsurface irradiance reflectance
is independent of chlorophyll concentration. This feature is called a pivot
point. Using an analytical form of the radiative transfer equation, they
identify the pivot point feature at around 497 nm. Using the specific absorption
and scattering coefficients from Bukata et al. (1995)
and the Rayleigh and particle VSF for the water and chlorophyll components
respectively, spectral irradiance reflectance signatures were generated
for concentrations of chlorophyll of 0.0, 0.5, 1.0, 10.0 and 20.0 mg m3.
The purpose of this simulation is to evaluate the model's ability to generate
realistic spectral signatures. The results of these
simulations are shown below . This simulation shows a pivot point at
around 495 nm which is in very close agreement to Bukata
et al.'s value of 497 nm.
4 Mobley et al.'s Canonical Problems
The next set of simulations uses four of Mobley
et al.'s (1993) canonical problems. Mobley et al.'s
problems have the following parameters:
· refractive index of water equals 1.340,
· airwater interface is flat,
· water is vertically homogeneous,
· light source is that from a point at a zenith angle of 60 degrees,
· the sun provides a spectral irradiance just above the surface of
magnitude 1 W m2 nm1 (perpendicular to the sun's rays) on the water surface,
· light source is monochromatic.
Mobley et al. compare the outputs of seven different
models (2 analytical and 5 stochastic) for each one of these problems and
report the variability between them in terms of a coefficient of variation.
The results of the AOMC are compared to the mean values from Mobley
et al.'s multiple model runs and presented in terms of a percent relative
error in the table at the bottom of this page.
4.1 Rayleigh SPF (Mobley et al.'s problem 1)
Mobley et al.'s first test problem
uses a Rayleigh SPF. They present the results for three different apparent
optical properties: Ed, Eou, and Lu. These results are listed below.
The simulations are run for two different b/a ratios: 0.25 and 0.9.
These ratios are translated to absorption coefficients of 1.0 m1 (for b/a=0.25),
1.0 m1(for b/a=0.9), and scattering coefficients of 0.25 m1(for
b/a=0.25) and 9.0 m1 (for b/a=0.9). In all cases, the relative
error between the AOMC model and mean values from Mobley
et al.'s multimodel simulation fall within the coefficient of variation
between Mobley et al.'s seven models.
4.2 Petzold's SPF (Mobley et al.'s problem 2)
Mobley et al.'s second canonical problem replaces the Rayleigh type SPF with that typical of coastal ocean waters which they calculate from three of Petzold's San Diego harbor measurements. The AOMC simulation was performed using the average of these three measured VSF for San Diego Harbor (Petzold, 1972). The table below compares the AOMC results to those from Mobley et al. (1993). Once again, the results from the AOMC simulation fall within the range of the other models in Mobley et al.'s simulations. The radiance values between the AOMC model and Mobley et al.'s are also in good agreement as shown in the following graph (AOMC result is shown in a solid red line, the mean result from Mobley et al.'s models is shown in dashed black lines):
4.3 Stratified water (Mobley et al.'s problem 3)
Mobley et al.'s problem 3 simulates a stratified water body where the concentrations of total particulate matter vary as a function of depth. This problem tests the AOMC model's ability to adjust the photon's pathlength as a function of the concentration of the attenuating medium. Mobley et al. use an analytical equation to describe the vertical profile of the IOPs. For this simulation, tabulated values of a and b are calculated from the analytical equations provided by Mobley et al. at 0.5 meter intervals. The agreement between the AOMC model simulation outputs and the mean values from Mobley et al.'s simulations are reasonably good at the 5 m depth, but diverges at increasing depth where a relative error value for Eou as high as 16.9% is found at a depth of 60 m. This divergence in simulated values is most likely due to the interval resolution used in the vertical profile of the IOPs. AOMC results are shown in the table below. The following graph plots the AOMC result (red solid line) and the results presented in Mobley et al.'s paper (black dashed lines).
4.4 Bottom effect (Mobley et al.'s problem 6)
Mobley et al.'s problem 6 simulates an optically homogeneous water body with a finite depth bottom. The bottom is a Lambertian reflector with an irradiance reflectance of 0.5. This problem tests the AOMC model's bottom boundary interaction. The relative errors for Ed and Eou between the AOMC model and Mobley et al.'s models falls within Mobley et al.'s CV with the exception for the near surface values of Eou and Lu. These differences could be a result of stochastic noise resulting from a low number of upward traveling photons making their way to the surface. Another cause for the error is the possible difference in bin sizes used in logging the angular direction of travel of the photon at each recording depth.
Summary of AOMC/Mobley et al. model comparison
Optical depth 
E_{d} 
E_{ou} 
L_{u} 

AOMC 
Mob. 
CV 
RD 
AOMC 
Mob. 
CV 
RD 
AOMC 
Mob. 
CV 
RD 

Mobley et al’s problem 1: scattering albedo = 0.2 

1 
1.411 
1.411 
0.001 
0.000 
01.332 
1.342 
0.003 
0.007 
1.6703 
1.723 
0.044 
0.029 

5 
1.073 
1.073 
0.005 
0.000 
0.974 
1.004 
0.039 
0.030 
1.2905 
1.375 
0.288 
0.058 

Mobley et al’s problem 1: scattering albedo = 0.9 

1 
3.6601 
3.661 
0.002 
0.000 
3.721 
3.721 
0.005 
0.000 
4.692 
4.852 
0.015 
0.033 

5 
4.3302 
4.332 
0.003 
0.000 
4.352 
4.352 
0.007 
0.000 
5.753 
5.593 
0.052 
0.028 

10 
3.1603 
3.163 
0.015 
0.000 
3.153 
3.203 
0.038 
0.016 
4.264 
4.374 
0.091 
0.025 

Mobley et al’s problem 2: scattering albedo = 0.9 

1 
4.141 
4.131 
0.001 
0.002 
9.162 
9.312 
0.021 
0.016 
6.803 
6.993 
0.063 
0.028 

5 
1.881 
1.871 
0.005 
0.005 
4.572 
4.632 
0.017 
0.013 
3.123 
3.263 
0.055 
0.043 

10 
6.942 
6.852 
0.010 
0.013 
1.632 
1.652 
0.014 
0.012 
1.163 
1.213 
0.109 
0.041 

Physical depth 
Mobley et al’s problem 3: scattering albedo = 0.9 

5 m 
2.351 
2.301 
0.006 
0.022 
4.352 
4.342 
0.025 
0.002 
3.013 
3.133 
0.054 
0.040 

25 m 
1.743 
1.623 
0.028 
0.074 
3.044 
2.864 
0.038 
0.063 
2.205 
2.125 
0.061 
0.038 

60 m 
5.405 
5.235 
0.071 
0.032 
0.064 
5.136 
0.036 
0.169 
4.007 
3.577 
0.434 
0.120 

Mobley et al’s problem 6: scattering albedo = 0.9 

1 
1.621 
1.621 
0.000 
0.000 
9.604 
9.814 
0.010 
0.021 
7.255 
6.845 
0.020 
0.059 

5 
2.293 
2.283 
0.003 
0.004 
2.283 
2.283 
0.002 
0.000 
3.604 
3.604 
0.010 
0.000 

Coming soon: downloadable test files for all simulations listed on this page
Bukata, RP, J.H. Jerome, K.Y. Kondratyev, and D.V. Pozdynyakov, 1995. Optical Properties and Remote Sensing of Inland and Coastal Waters. CRC Press, New York, 362 pp.
Gordon, H.R., O.B. Brown, and M.M. Jacobs, 1975. Computed relationships between the inherent and apparent optical properties of a flat homogeneous ocean. Applied Optics, 14, p. 417427.
Kirk, J.T.O., 1981. Monte Carlo study of the nature of the underwater light field in, and the relationships between optical properties of, turbid yellow waters. Australian Journal of Marine and Freshwater Resource. 32, p. 517532.
Mobley, C.D., B. Gentili, H.R. Gordon, Z. Jin, G.W. Kattawar, A. Morel, P. Reinersmam, K. Stamnes, and R.H. Stavn, 1993. Comparison of numerical models for computing underwater light fields. Applied Optics, 32(6), p. 7,4847,504.
Morel, A. and L. Prieur, 1977. Analysis of variations in ocean color. Limnology and Oceanography, 22(4), p. 709722.
Petzold, T. J. 1972. Volume Scattering Functions for selected ocean waters, Scripps Institution of Oceanography, San Diego, CA. SIO Ref. 7278, 82 pp.
Sathyendranath, S. and Platt, T., 1997. Analytical model of ocean color. Appl. Opt. 36 (12), pp. 26202629.