Page updated: May 19, 2021
Author: Curtis Mobley
View PDF

# The Lidar Equation

“Lidar” is an acronym for “light detection and ranging,” but should it be written as LIDAR, LiDAR, Lidar, lidar, or something else? Deering and Stoker (2014) argue that lidar has now achieved the same status as radar (radio detection and ranging) and sonar (sound navigation and ranging) and should be in lower case (except as the ﬁrst word of a sentence), so that is what is used here.

The purpose of a “lidar equation” is to compute the power returned to a receiver for given transmitted laser power, optical properties of the medium through which the lidar beam passes, and target properties. There are, however, many versions of lidar equations, each of which is tailored to a particular application. For example, Measures (1992) develops lidar equations for elastic and inelastic backscattering scattering by the medium, for ﬂuorescent targets, for topographical targets, for long-path absorption, for broad-band lasers, and so on.

The lidar equation developed here applies to the detection of a scattering layer or in-water target as seen by a narrow-beam laser imaging the ocean from an airborne platform. This equation explicitly shows the eﬀects of atmospheric and sea-surface transmission, the water volume scattering function and beam spread function, water-column diﬀuse attenuation, and transmitter and receiver optics.

(Acknowledgment: I learned this derivation from Richard C. Honey, one of the pioneering geniuses of optical oceanography, and of many other ﬁelds including antenna design and laser eye surgery. Dick Honey is unfortunately little known to the general community because he spent much of his career doing classiﬁed work.)

### Preliminaries

Table 1 lists for reference the variables involved with the derivation of the present form of the lidar equation.

 variable deﬁnition units $H$ Height of the airplane above the sea surface m $z$ Depth of the water layer being imaged m $\Delta z$ Thickness of the water layer being imaged m ${P}_{t}$ Power transmitted by the laser W $\Delta {P}_{r}$ Power detected from water layer $\Delta z$ W ${T}_{a}$ Transmittance by the air nondimen ${T}_{s}$ Transmittance by the water surface nondimen ${A}_{r}$ Receiver aperture area ${m}^{2}$ ${\Omega }_{FOV}$ Receiver ﬁeld of view solid angle sr ${A}_{z}$ Area at depth $z$ seen by the receiver ${m}^{2}$ $\Omega$ Solid angle of the receiver aperture as seen from depth$z$ sr ${\beta }_{\pi }$ Water VSF for 180-deg backscatter ${m}^{-1}\phantom{\rule{2.6108pt}{0ex}}s{r}^{-1}$ $BSF$ Water beam spread function ${m}^{-2}$ ${E}_{i}$ Irradiance incident (downward) onto the water layer at $z$ $W\phantom{\rule{2.6108pt}{0ex}}{m}^{-2}$ $\Delta {E}_{r}$ Irradiance reﬂected (upward) by the water layer $\Delta z$ $W\phantom{\rule{2.6108pt}{0ex}}{m}^{-2}$ ${\overline{K}}_{up}$ Depth-averaged (over 0 to $z$) attenuation coeﬃcient for upwelling (returning) irradiance $\Delta {E}_{r}$ ${m}^{-1}$

Table 1: Variables occurring in the lidar equation.

Figure 1 shows the geometry of the lidar system for detection of a scattering layer. Figure 1: Geometry of the lidar system for detection of a scattering layer. The left panel illustrates the physical variables. The right panel illustrates the system optics; this panel ignores refraction at the sea surface.

#### The volume scattering function (VSF)

Recall from Eq. (1) of the Volume Scattering Function page that the volume scattering function (VSF) is operationally deﬁned by

 $\beta \left(\psi \right)=\frac{\Delta {\Phi }_{s}\left(\psi \right)}{{\Phi }_{i}\phantom{\rule{0.3em}{0ex}}\Delta z\phantom{\rule{0.3em}{0ex}}\Omega }\phantom{\rule{0.3em}{0ex}},$

where ${\Phi }_{i}$ is the power incident onto an element of volume deﬁned by a surface area $A$ and thickness $\Delta z$, and $\Delta {\Phi }_{s}\left(\psi \right)$ is the power scattered through angle $\psi$ into solid angle $\Omega$. These quantities are illustrated in Fig. 2. In this derivation, quantities directly proportional to the layer thickness $\Delta z$ will be labeled with a $\Delta$. Thus $\Delta {\Phi }_{s}\left(\psi \right)$ is the power scattered by the water layer of thickness $\Delta z$. The incident power ${\Phi }_{i}$ onto area $A$ gives an incident irradiance ${E}_{i}={\Phi }_{i}∕A$. Figure 2: Quantities associated with the deﬁnition of the volume scattering function.

Now consider exact backscatter, which is a scattering angle of $\psi =\pi$, or $\psi =$ 180 deg. The backscattered power $\Delta {\Phi }_{s}\left(\pi \right)$ exits the scattering volume through the same area $A$, so the backscattered irradiance is $\Delta {E}_{s}\left(\pi \right)=\Delta {\Phi }_{s}\left(\pi \right)∕A$. The VSF for exact backscatter can then be written as

 ${\beta }_{\pi }\equiv \beta \left(\psi =\pi \right)=\frac{\Delta {\Phi }_{s}\left(\pi \right)∕A}{\left({\Phi }_{i}∕A\right)\phantom{\rule{0.3em}{0ex}}\Delta z\phantom{\rule{0.3em}{0ex}}\Omega }=\frac{\Delta {E}_{s}\left(\pi \right)}{{E}_{i}\phantom{\rule{0.3em}{0ex}}\Delta z\phantom{\rule{0.3em}{0ex}}\Omega }\phantom{\rule{0.3em}{0ex}}.$

This gives

 $\Delta {E}_{s}\left(\pi \right)={E}_{i}\phantom{\rule{0.3em}{0ex}}{\beta }_{\pi }\phantom{\rule{0.3em}{0ex}}\Omega \phantom{\rule{0.3em}{0ex}}\Delta z\phantom{\rule{0.3em}{0ex}}.$ (1)

(Do not confuse this backscattered or reﬂected irradiance with ${E}_{u}=R\phantom{\rule{0.3em}{0ex}}{E}_{d}$, which is deﬁned for any incident radiance distribution and for an arbitrarily thick layer of water. Here ${E}_{i}$ is the irradiance for a collimated incident laser beam and a thin layer of water.)

#### The beam spread function (BSF)

Suppose a collimated beam is emitting power $P$ in direction $𝜃=0$. Then scattering and absorption in the medium will give some irradiance $E\left(r,𝜃\right)$ on the surface of a sphere of radius $r$ at an angle of $𝜃$ relative to the direction of the emitted beam, as illustrated by the green arrow in Fig. 3. The beam spread function (BSF) is then deﬁned as

 $BSF\left(r,𝜃\right)\equiv \frac{E\left(r,𝜃\right)}{P}\phantom{\rule{2em}{0ex}}\left({m}^{-2}\right)\phantom{\rule{0.3em}{0ex}}.$ (2) Figure 3: Geometry for deﬁnition of the beam spread function.

The beam spread function and its equivalent, the point spread function, are discussed in more detail on the Beam and Point Spread Functions page.

### Derivation of the Lidar Equation

The derivation of the lidar equation for the stated application now proceeds via the following steps:

1.
As illustrated in Fig. 1, the laser pulse has transmitted power ${P}_{t}$ to start. After transmission through the atmosphere and sea surface, the pulse has power
 ${P}_{w}\left(0\right)={P}_{t}\phantom{\rule{0.3em}{0ex}}{T}_{a}\phantom{\rule{0.3em}{0ex}}{T}_{s}$

just below the water surface.

2.
The laser pulse is still a narrow beam when it enters the water, but it then begins to spread out because of scattering, and it is attenuated by absorption. This process is quantiﬁed by the beam spread function. We are interested in the “on-axis” irradiance incident onto a layer of water at depth $z$. From Eq. (2), this is given by
 ${E}_{i}\left(z\right)={P}_{w}\left(0\right)\phantom{\rule{0.3em}{0ex}}BSF\left(z,𝜃=0\right)\phantom{\rule{0.3em}{0ex}}.$
3.
The irradiance that is reﬂected by a layer of thickness $\Delta z$ at depth $z$ is given by Eq. (1):
 $\Delta {E}_{r}\left(z\right)={E}_{i}\left(z\right)\phantom{\rule{0.3em}{0ex}}{\beta }_{\pi }\left(z\right)\phantom{\rule{0.3em}{0ex}}\Omega \phantom{\rule{0.3em}{0ex}}\Delta z\phantom{\rule{0.3em}{0ex}}.$

The solid angle $\Omega$ is determined by the aperture size of the receiver and the distance from the receiver to the layer. This is discussed in step 7 below.

4.
The irradiance $\Delta {E}_{r}\left(z\right)$ heading upward from the layer $\Delta z$ will be attenuated by some diﬀuse attenuation function ${\overline{K}}_{up}$ as $\Delta {E}_{r}$ propagates back to the sea surface:
 $\Delta {E}_{r}\left(0\right)=\Delta {E}_{r}\left(z\right)exp\left(-{\overline{K}}_{up}z\right)\phantom{\rule{0.3em}{0ex}}.$

The receiver sees only an area ${A}_{z}$ depth depth $z$. It is assumed that ${A}_{z}$ is less than the total area illuminated by the spreading laser beam at depth $z$. Thus we can multiply $\Delta {E}_{r}\left(z\right)$ by ${A}_{z}$ to get the power leaving the illuminated area that is seen by the receiver. The fraction of the total power reaching the surface, which is seen by the detector, is then

 $\Delta {P}_{r}\left(0\right)=\Delta {E}_{r}\left(z\right)\phantom{\rule{0.3em}{0ex}}{A}_{z}\phantom{\rule{0.3em}{0ex}}exp\left(-{\overline{K}}_{up}\phantom{\rule{0.3em}{0ex}}z\right)\phantom{\rule{0.3em}{0ex}}.$

${A}_{z}$ is determined by the receiver ﬁeld of view and the distance from the receiver to the layer.

5.
The power just below the surface is now transmitted through the water surface and through the atmosphere to the receiver. The power detected is then
 $\Delta {P}_{r}=\Delta {P}_{r}\left(0\right)\phantom{\rule{0.3em}{0ex}}{T}_{s}\phantom{\rule{0.3em}{0ex}}{T}_{a}\phantom{\rule{0.3em}{0ex}}.$
6.
Combining the above results gives
 $\Delta {P}_{r}={P}_{t}\phantom{\rule{0.3em}{0ex}}{T}_{a}^{2}\phantom{\rule{0.3em}{0ex}}{T}_{s}^{2}\phantom{\rule{0.3em}{0ex}}{A}_{z}\phantom{\rule{0.3em}{0ex}}\Omega \phantom{\rule{0.3em}{0ex}}BSF\left(z,0\right)\phantom{\rule{0.3em}{0ex}}{\beta }_{\pi }\phantom{\rule{0.3em}{0ex}}exp\left(-{\overline{K}}_{up}z\right)\phantom{\rule{0.3em}{0ex}}\Delta z\phantom{\rule{0.3em}{0ex}}.$
7.
${A}_{z}$ and $\Omega$ can be written in terms of known parameters. We are interested only in power that is reﬂected by the layer $\Delta z$ into the solid angle $\Omega$ that will take the light into the receiver. As shown in the right panel of Fig. 1, the receiver aperture ${A}_{r}$ and the range $H+z$ determine $\Omega ={A}_{r}∕{\left(H+z\right)}^{2}$. Only power heading into $\Omega$ from area ${A}_{z}$, which is seen by the receiver, reaches the receiver. ${A}_{z}$ is determined by the receiver ﬁeld of view solid angle ${\Omega }_{FOV}$ and the range: ${A}_{z}={\Omega }_{FOV}{\left(H+z\right)}^{2}$. Thus the ${A}_{z}\Omega$ factor in the preceding equation can be rewritten as
 ${A}_{z}\Omega ={\Omega }_{FOV}{\left(H+z\right)}^{2}\frac{{A}_{r}}{{\left(H+z\right)}^{2}}={A}_{r}{\Omega }_{FOV}\phantom{\rule{0.3em}{0ex}}.$

This equation shows how the receiver optics aﬀects the detected power. (This is an application of the “$A\Omega$” theorem of optical engineering. $A\Omega$ is also called the throughput or étendue of the system. Strictly speaking, the in-air solid angle ${\Omega }_{FOV}$ decreases by a factor of $1∕{n}^{2}$ upon entering the water, where $n$ is the water real index of refraction. However, the in-water $\Omega$ increases by a factor of ${n}^{2}$ when entering the air. For a round trip, air to water to air, these ${n}^{2}$ factors cancel. They are thus ignored here and in Fig. 1.)

8.
Collecting the above results gives the desired lidar equation:
 $\Delta {P}_{r}={P}_{t}\phantom{\rule{0.3em}{0ex}}{T}_{a}^{2}\phantom{\rule{0.3em}{0ex}}{T}_{s}^{2}\phantom{\rule{0.3em}{0ex}}{A}_{r}\phantom{\rule{0.3em}{0ex}}{\Omega }_{FOV}\phantom{\rule{0.3em}{0ex}}BSF\left(z,0\right)\phantom{\rule{0.3em}{0ex}}{\beta }_{\pi }\phantom{\rule{0.3em}{0ex}}exp\left(-{\overline{K}}_{up}z\right)\phantom{\rule{0.3em}{0ex}}\Delta z\phantom{\rule{0.3em}{0ex}}.$ (3)

Equation (3) nicely shows the eﬀects of the transmitted power (${P}_{t}$), atmospheric and surface transmission (${T}_{a}^{2}\phantom{\rule{0.3em}{0ex}}{T}_{s}^{2}$), receiver-optics (${A}_{r}\phantom{\rule{0.3em}{0ex}}{\Omega }_{FOV}$), water-column ($BSF\left(z,0\right)\phantom{\rule{0.3em}{0ex}}{\beta }_{\pi }\phantom{\rule{0.3em}{0ex}}exp\left(-{\overline{K}}_{up}z\right)$), and layer thickness ($\Delta z$). The take-home message from this equation is that in order to understand lidar data, the water inherent optical properties you need to know are the beam spread function and ${\beta }_{\pi }$. This observation was in part the incentive for the work of Mertens and Replogle (1977), Voss and Chapin (1990), Voss (1991), McLean and Voss (1991), Maﬃone and Honey (1992), Gordon (1994b), McLean et al. (1998), Sanchez and McCormick (2002), Dolin (2013), Xu and Yue (2015), and others. These papers present several models for beam/point spread functions in terms of the water absorption and scattering properties. Several of those models are reviewed in Hou et al. (2008).

There are various arguments about what to use for ${\overline{K}}_{up}$, which depends on both the water optical properties and on the imaging system details. It is intuitively expected that ${\overline{K}}_{up}$, which is deﬁned for a horizontally small patch of upwelling irradiance, will be greater than the diﬀuse attenuation coeﬃcient for upwelling irradiance, ${K}_{u}$, which is deﬁned for a horizontally inﬁnite light ﬁeld. Likewise, we expect that ${\overline{K}}_{up}$ will be less than $c$, the beam attenuation coeﬃcient. Thus ${K}_{u}<{\overline{K}}_{up}. Because ${\stackrel{̄}{K}}_{up}$ is an attenuation function for a ﬁnite patch of reﬂected irradiance, computing its value is inherently a three-dimensional radiative transfer problem. To pin down the value of ${\overline{K}}_{up}$ more accurately thus requires either actual measurements for a particular system and water body, or three-dimensional raditative transfer simulations (usually Monte Carlo simulations) tailored to a given system and water properties.

#### Example calculation

To develop some intuition about Eq. (3), consider the following example application. Suppose a 532 nm laser is being used to look for objects in the water that have an area of $1\phantom{\rule{2.6108pt}{0ex}}{m}^{2}$. The receiver FOV must be small enough that the object can be distinguished from its surroundings. For $H=100\phantom{\rule{2.6108pt}{0ex}}m$ and $z=10\phantom{\rule{2.6108pt}{0ex}}m$, this requires that

 ${\Omega }_{FOV}\le \frac{1\phantom{\rule{2.6108pt}{0ex}}{m}^{2}}{\left(100+10\right)\phantom{\rule{2.6108pt}{0ex}}{m}^{2}}\approx 8×1{0}^{-5}\phantom{\rule{1em}{0ex}}sr\phantom{\rule{0.3em}{0ex}}.$

For a 15 cm radius receiving telescope, ${A}_{r}=\pi {\left(0.15\phantom{\rule{2.6108pt}{0ex}}m\right)}^{2}\approx 0.07\phantom{\rule{2.6108pt}{0ex}}{m}^{2}$. For normal incidence at the sea surface, ${T}_{s}\approx 0.97$, and for a clear atmosphere, ${T}_{a}\approx 0.98$. Suppose the water is Jerlov type 1 coastal water for which ${K}_{d}\left(532\right)\approx 0.15\phantom{\rule{2.6108pt}{0ex}}{m}^{-1}$ (Light and Water, page 130), and assume that ${\overline{K}}_{up}={K}_{d}$. Further assume that $BSF\approx exp\left(-0.2z\right)$, since the lidar beam attenuation will be more “beam like” than diﬀuse attenuation, and $c>{K}_{d}$. Finally, take ${\beta }_{\pi }=1{0}^{-3}\phantom{\rule{2.6108pt}{0ex}}{m}^{-1}\phantom{\rule{2.6108pt}{0ex}}s{r}^{-1}$ (Light and Water, Table 3.10 for “coastal ocean” water). Then for $\Delta z=1\phantom{\rule{2.6108pt}{0ex}}m$ at a depth of 10 m, the water returns the fraction

$\begin{array}{llll}\hfill \frac{\Delta {P}_{r}}{{P}_{t}}=& {\left(0.98\right)}^{2}\phantom{\rule{0.3em}{0ex}}{\left(0.97\right)}^{2}\phantom{\rule{0.3em}{0ex}}\left(0.07\right)\phantom{\rule{0.3em}{0ex}}\left(8×1{0}^{-5}\right)\phantom{\rule{0.3em}{0ex}}{e}^{-\left(0.2\right)10}\phantom{\rule{0.3em}{0ex}}\left(1{0}^{-3}\right)\phantom{\rule{0.3em}{0ex}}{e}^{-\left(0.15\right)10}\phantom{\rule{0.3em}{0ex}}\left(1\right)\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill \approx & 1.5×1{0}^{-10}\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$

of the transmitted power.

Now suppose that the laser beam hits an object at $z$ = 10 m whose surface is a Lambertian reﬂector of 2% (irradiance) reﬂectance. Then 0.02 of ${E}_{i}$ is reﬂected into $2\pi \phantom{\rule{2.6108pt}{0ex}}sr$. The layer backscatter ${\beta }_{\pi }\Delta z=0.001\phantom{\rule{2.6108pt}{0ex}}s{r}^{-1}$ is then replaced by

 $\frac{0.02}{2\pi \phantom{\rule{2.6108pt}{0ex}}sr}\approx 0.003\phantom{\rule{2.6108pt}{0ex}}s{r}^{-1}\phantom{\rule{0.3em}{0ex}}.$

If all other terms remain the same, the object would return about three times the signal as the water itself.