Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

New theoretical formulation for the determination of radiance transmittance at the water-air interface: reply

Open Access Open Access

Abstract

Gordon and Voss in their comment challenge the recently published paper [Opt. Express 25, 27086 (2017)] on the unity factor of radiance transmittance of the assumed hypothetical case (i.e., for the albedo equal to 1) and question the dependence of particulate contribution to the refractive index of water. Here, we provide answers to their comments.

© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Transmittance in the hypothetical water

The comment presumes that we have termed the well-established n2 – law for radiance Eq. (1) equation to be incorrect. We object to this supposition and we do not ignore or term it as incorrect, which can be evidenced by the use of the same equation in our formulation. The geometry ascribed in n2 – law for radiance is for the radiance interacting with the interface only once [1,2], and the multiple interactions are completely omitted. Here, we answer to the comment challenging the transmittance of the assumed hypothetical water is unity. The hypothetical situation assumed by Gordon and Voss is analogous to our hypothetical assumption, in which they used R = 1 and we assumed the water medium as non-absorbing and fully upward scattering. Both assumptions state that no photons are retained within the water body in any form, i.e., fully reflecting medium. Now, we prove that Eq. (2) of the comment itself shows the transmittance of hypothetical waters turns to unity, when multiple surface reflections of the photons at the waterside are considered. For that we provide the proof below. The n2 – law for radiance applied to the air-water medium is given by,

Lu+(θa,ϕa)=Lu(θw,  ϕw)nw2tf(θw,φw).

Upwelling radiance before crossing the interface, inclusive of multiple surface interactions at the waterside (same as represented in the comment), is expressed as,

Lu=Lu0+Lu1+Lu2+=Lu0+r¯RLu0+(r¯R)2Lu0+=Lu01r¯R.

Lu0 is the in-water radiance before striking the interface and Lu1, Lu2, ... are the successive surface reflected components of Lu0. Radiance after crossing the interface [Eq. (2) in Eq. (1)],

Lu+=tfnw2Lu01r¯R=tfnw2[11r¯R]Lu0.

Comparing Eq. (3) and Eq. (1), we get an extra term [11r¯R]. As the Lu0 is the in-water radiance before striking the interface, the Lu0 in Eq. (3) and the Lu(θw, ϕw) in Eq. (1) contain the same number of photons within the solid angle. Therefore, [11r¯R] cannot be the component of radiance as stated in the comment [Eq. (2)]. It is clearly a component of the interface and water. The term [11r¯R] varies with the optical properties of water and interface, which additionally contributes to the existing transmittance factor (tf/nw2) and has nothing to do with the radiance. For the hypothetical case of Gordon and Voss (in the comment) R = 1,

Lu+=tfnw2[11r¯]Lu0=0.541×[110.475]Lu0=1.03×Lu0.

Clearly, the above equation results in unity. The value of r¯ given in Gordon et al. [3] and the comment is 0.475 (approx.). Practically, there cannot be 103% of Lu0. We understand the constant r¯~0.475 is an approximate value, and for a 100% of Lu0, the calculation results in the value of r¯=0.459. The correct value of r¯ can be obtained by considering a uniform upwelling light field above and below the interface.

Eu+=ΩdLu+(θa,ϕa)cosθadΩa=πLu+. 
Eu=ΩdLu(θw,ϕw)cosθwdΩw=πLu.

For this assumed hypothetical case, if no photons are retained within the medium, then Eqs. (5) and (6) will be equal and the Eq. (3) reduces to tfnw2[11r¯]=1. This leads the r¯ value to 0.459 by assuming tf/nw2 = 0.541. Now, we consider the present formulation (Eq. (11) in Dev and Shanmugam [4]), where Lu0 is the in-water radiance before striking the interface. For a non-absorbing and a fully upward scattering medium, μu=1 and ω=1,

Lu+=tfnw2[1μuωrf2+μuωnw2tf]×Lu0=tfnw2[11rf2+nw2tf]×Lu0=1.00×Lu0.

Equations (4) and (7) conclude the same and agree with each other that the transmittance in the hypothetical water turns to be unity, after the inclusion of multiple surface reflections of photons with the water-air interface. The only difference between these two formulations is that the comment deals the processes from the standpoint of internal reflection, whereas the present formulation deals from the standpoint of transmittance of photons on every interaction with the interface. If [1μuωrf2+μuωnw2tf] replaces the [11r¯R] in Eq. (2), all the photons will eventually cross the interface satisfying the Eq. (3) to Eq. (7) (in the comment) as described. Therefore, no violation of physical concepts in both the cases.

Here, we provide an example that allows us to examine the correctness of Gordon and Voss (GV) interpretation [Eq. (2)]. According to GV, Lu is the measurable quantity and Lu0 is non-measurable and a theoretical construct.

Now, let us assume radiance at two depths z1 and z2 in a homogeneous and a diffuse medium as Luz1 and Luz2. Now, let us extrapolate Luz1 to ‘just below-water surface radiance denoted by Luz0. Using Eq. (8), Luz1 can be extrapolated to Luz0.

Luz0=Luz1eKLuz1.

where, KLu=1(z2z1)[ln(Luz1)ln(Luz2)]. As we know, Luz1 and Luz2 are not involved in any interactions with the interface, and the obtained Luz0 though the extrapolation is also not involved in any interaction with the interface. This Luz0 and GV’s Lu0 are analogous to each other because both radiances have not yet interacted with the interface. On propagating both the quantities through the interface, Luz0 can be successfully translated to Lw, but GV’s Lu0 requires 1r¯R (according to the comment) which is not possible to translate though the interface. Therefore, two things we can understand from this:

  • 1) GV’s interpretation failed to explain Eq. (8) and the propagation process through the interface.
  • 2) Lu0 is not a theoretical construct as claimed in the comment, since we have simulated it using Eq. (8) (Luz0 is same as their Lu0).

Now, it became clear that 1r¯R cannot be represented along with the radiance, as it is the property of the interface and the water. For this reason, we interpreted GV’s Lu0 as the radiance before interacting with the interface. The only way to avoid the failure of GV’s interpretation is to consider Lu0 as Lu0 as we have shown above. It means that the usage of the term Lu0 is incorrect, and hence, Eq. (2) becomes invalid.

2. Dependence of particulate contribution to the refractive index of water and the transmittance

We initially suspected that there would be significant differences in the Snell cones of the pure water and the turbid water (especially high scattering water) because of large difference in the particle concentrations between these waters. To find the Snell cones of these two water regimes, we performed an experiment in our laboratory to measure the critical angles for the pure and the turbid water by tracing the incident ray. Pure water was filled inside the glass tank (1 ft x 1 ft x 1 ft) and red laser source was used. Tracing the incident ray for the critical angle measurement, it provided the critical angle of 48 degrees. Then the particles were introduced to it slowly until it becomes very turbid. During this measurement, the light became too diffuse, hindering the observable point, however critical angle hovered around 47 degrees. We were only able to observe one 1 degree of difference between the pure and the turbid water which gives the difference in refractive index about (1.3628-1.34 = 0.0228). Our experiment may contain measurement errors, as the experiment was not performed with sophisticated instruments.

Since the particulates contribute significantly less to the refractive index of waters, n(λ) can be replaced with the nw(λ) in the Eq. (8) of the published paper [4]. The corrected equations are given below. The transmittance equation Eq. (10) in [4] becomes,

τw,a=(1ρw,a)nw2[(1μuω)+μuωnw2(1ρw,a)].
Similarly, the water-leaving radiance (Eq. (11) in [4] and translation of below to above-water remote sensing signal (Eq. (18) in [4]) become,
Lw(air,λ)=(1ρw,a)nw2[(1μuω)+μuωnw2(1ρw,a)]Lu(0,λ).
Rrs(0+,λ)=(1ρa,w)(1ρw,a)nw2[(1μuω)+μuωnw2(1ρw,a)]rrs(0,λ).
The term representing the particulate contribution to the refractive index of the water, rf2, has been removed, and it is the only change between this reply and the original manuscript [4]. The results do not vary significantly in the case of clear oceanic waters, but show some variations (approximately up to 15%) in the case of high scattering waters. However, the multiple scattering effect is significantly high in turbid waters as confirmed by this study.

Funding

Department of Science and Technology (DST) through the Networked Project on Big Data Analytics-Hyperspectral Data (OEC1617130DSTXPSHA).

References and links

1. C. L. Wyatt, Radiometric Calibration: Theory and Methods (Academic Press, 1978).

2. C. D. Mobley, “Radiative Transfer: Across the surface,” in Light and Water: Radiative Transfer in the Natural Waters (Academic Press, 1994), pp. 159–161.

3. H. R. Gordon, J. W. Brown, O. B. Brown, R. H. Evans, R. C. Smith, K. S. Baker, and D. K. Clark, “A semianalytic radiance model of ocean color,” J. Geophys. Res. 93(D9), 10909–10924 (1988). [CrossRef]  

4. P. J. Dev and P. Shanmugam, “New theoretical formulation for the determination of radiance transmittance at the water-air interface,” Opt. Express 25(22), 27086–27103 (2017). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Equations (11)

Equations on this page are rendered with MathJax. Learn more.

L u + ( θ a , ϕ a )= L u ( θ w ,   ϕ w ) n w 2 t f ( θ w , φ w ).
L u = L u0 + L u1 + L u2 += L u0 + r ¯ R L u0 + ( r ¯ R ) 2 L u0 += L u0 1 r ¯ R .
L u + = t f n w 2 L u0 1 r ¯ R = t f n w 2 [ 1 1 r ¯ R ] L u0 .
L u + = t f n w 2 [ 1 1 r ¯ ] L u0 =0.541×[ 1 10.475 ] L u0 =1.03× L u0 .
E u + = Ω d L u + ( θ a , ϕ a )cos θ a d Ω a =π L u + . 
E u = Ω d L u ( θ w , ϕ w )cos θ w d Ω w =π L u .
L u + = t f n w 2 [ 1 μ u ω r f 2 + μ u ω n w 2 t f ]× L u0 = t f n w 2 [ 11 r f 2 + n w 2 t f ]× L u0 =1.00× L u0 .
L u z 0 = L u z 1 e K Lu z 1 .
τ w,a = (1 ρ w,a ) n w 2 [ ( 1 μ u ω )+ μ u ω n w 2 (1 ρ w,a ) ].
L w ( air,λ )= ( 1 ρ w,a ) n w 2 [ ( 1 μ u ω )+ μ u ω n w 2 ( 1 ρ w,a ) ] L u ( 0 ,λ).
R rs ( 0 + ,λ )= (1 ρ a,w )( 1 ρ w,a ) n w 2 [ ( 1 μ u ω )+ μ u ω n w 2 ( 1 ρ w,a ) ] r rs ( 0 ,λ).
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.