## Abstract

Multiple scattering is an inevitable effect in spaceborne oceanic lidar because of the large footprint size and the high optical density of seawater. The effective attenuation coefficient *k _{lidar}* in the oceanic lidar equation, which indicates the influence of the multiple scattering effect on the formation of lidar returns, is an important parameter in the retrieval of inherent optical properties (IOPs) of seawater. In this paper, the relationships between

*k*of the spaceborne lidar signal and the IOPs of seawater are investigated by solving the radiative transfer equation with an improved semianalytic Monte Carlo model. Apart from the geometric loss factors,

_{lidar}*k*is found to decrease exponentially with the increase of depth in homogeneous waters.

_{lidar}*k*is given as an exponential function of depth and IOPs of seawater. The mean percentage errors between

_{lidar}*k*calculated by the exponential function and the simulated ones in three typical stratified waters are within 0.5%, proving the effectiveness and applicability of this

_{lidar}*k*-IOPs function. The results in this paper can help researchers have a better understanding of the multiple scattering effect of spaceborne lidar and improve the retrieval accuracy of the IOPs and the chlorophyll concentration of case 1 water from spaceborne lidar measurements.

_{lidar}© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

The ocean is a complex biological system that contains a vast number of biological species. Phytoplankton play a very important role in the material recycle and radiation budget of the marine ecosystem [1]. Satellite-based ocean color observations have been widely used to understand the distribution and variable characteristics of global phytoplankton [2]. The current ocean color measurements were made by passive sensors, which have several fundamental limitations: retrievals are highly sensitive to the atmospheric correction errors; information about the vertical distribution of ocean layers is missing; observations during nighttime, at high solar zenith angles and through thick cloud or absorbing aerosols are unavailable [3]. The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) has demonstrated advantages in day-night continuous sampling, observations through absorbing aerosols and thin cloud layers [4], capabilities of vertical profiling through the subsurface [5–7], and measurements at high-latitude regions [8]. However, CALIOP has limited spatial coverage compared to most passive sensors, and can only provide an integrated view of optical properties in the upper ocean. Specifically, CALIOP observation is limited by its coarse vertical resolution (i.e., 22.4m in the water) [9,10]. To improve the retrieval accuracy of spaceborne lidars, more information is needed to better understand factors involved in lidar functions.

Lidars obtain depth-resolved return profiles by collecting backscattered signals from molecules and particles. It is essential to establish an accurate model of the radiative transfer process. The return signal of spaceborne oceanic lidar is strongly affected by multiple scattering due to the large footprint size and the high optical density of seawater [11].To characterize the multiple-scattering radiative transfer process, an effective attenuation coefficient *k _{lidar}* is introduced into the oceanic lidar equation [12]. The relationships between

*k*and the optical properties of water and the lidar system parameters were extensively studied. Gordon [12] found that the ratio of the radius of the spot on the sea surface viewed by the lidar receiver optics to the mean free path of photons in water is a significant parameter which is related to

_{lidar}*k*. When this ratio is near zero,

_{lidar}*k*is given by the beam attenuation coefficient. If the ratio is greater than 5,

_{lidar}*k*is given by the diffuse attenuation coefficient

_{lidar}*K*. Phillips et al. [13] proposed a slightly different assumption that the upper and lower limits of

_{d}*k*is given by the beam attenuation coefficient and the absorption coefficient, respectively, depending on field of view (FOV). Lee et al. [14] and Schulien et al. [15] observed strong correlations between the

_{lidar}*in situ K*and

_{d}*k*derived from airborne lidar experiments. Collister et al. [16] exhibited the similarity between

_{lidar}*k*and the absorption coefficient under wide FOV using shipborne lidar measurements. The aforementioned studies are based on airborne or shipborne platform, and the depth characteristics of

_{lidar}*k*are ignored. The values of

_{lidar}*k*at different depths in homogenous water are not constant, due to the multiple scattering effect [17,18], especially near the ocean surface. This phenomenon leads to errors in the retrieval of optical properties, which further affects the assessment accuracy of the phytoplankton biomass and primary productivity. Thus, the depth characteristics of

_{lidar}*k*and the relationship between

_{lidar}*k*and the optical properties of seawater need to be studied in detail.

_{lidar}In this paper, the potential of using an improved semianalytic Monte Carlo radiative transfer simulation model (SALMON) to accurately identify the characteristics of both homogeneous and inhomogeneous waters is investigated. The use of SALMON is proposed by Poole et al. [19] and widely applied in oceanic lidar, therefore, it would be beneficial to develop an improved SALMON that can better simulate the inhomogeneous waters in the real oceans. In Section.2, we present the principle and method of proposed SALMON. In Section 3, we use this improved model to investigate the characteristics of lidar return profiles and the corresponding *k _{lidar}* with analysis of the influence of the lidar system parameters on

*k*and the depth characteristic of

_{lidar}*k*for homogenous water. In addition, an exact relationship between the depth-resolved

_{lidar}*k*and IOPs of seawater is derived based on the statistical properties of case 1 waters. In Section 4, we compare

_{lidar}*k*calculated with our function and

_{lidar}*K*from bio-optical model with the simulation results in three vertical stratified waters to verify the applicability of this function in real ocean environment.

_{d}## 2. Principle and method

#### 2.1Return signals of oceanic lidar

The schematic diagram of the theoretical model of laser backscatter from water is illustrated in Fig. 1. The laser pulse transmits downward towards the flat horizontal water surface at the height of $H$. The photons are scattered and absorbed by the molecules and particles in seawater, and the backscattered photons are received by the detector forming the lidar return signal, which can be described by the lidar equation [3,20]:

*P*(

*z*) is the power received from a depth of $z$,

*P*is the transmitted pulse power,

_{0}*M*is the lidar instrument constant,

*v*is the speed of light in vacuum, $\tau $is the temporal pulse length,$vt/(2n)$ is the depth resolution in water,

*T*is the transmission through the atmosphere,

_{atm}*T*is the Fresnel transmittance of the air-sea surface,

_{surf}*A*is the receiver area,

*n*is the refractive index of sea water and${\beta}_{\pi}$is the volume scattering coefficient at a scattering angle of$\pi $rad.

Rewriting Eq. (1), we have the ${P}_{norm}\left(z\right)$ expressed by

Thus, the effective attenuation coefficient can be written as

*k _{lidar}* is an important parameter for retrieving optical properties of seawater. The IOPs considered here include the scattering coefficient$b(\lambda )$(m

^{−1}), the absorption coefficient$a(\lambda )$ (m

^{−1}), and the beam attenuation coefficient$c(\lambda )$(m

^{−1}). The relation among those three parameters are expressed as

The IOPs of case 1 water can be parameterized in terms of the chlorophyll concentration$[Chl]$and wavelength [21]. The case 1 water is chlorophyll dominated water which includes almost 98% of the world's open ocean and coastal waters. Therefore, lidar remote sensing of the case 1 water plays an important role in estimating global ocean primary productivity.

According to the bio-optical model of case 1 water [21], the absorption coefficient could be expressed as

When [Chl] is expressed in the unit of mg/m^{3} and $\lambda $ is in the unit of nm, the resulting $a(\lambda )$ and $b(\lambda )$ are in the unit of m^{−1}. In this paper, the working wavelength is at 532 nm, which is commonly used in oceanic lidars.

The diffuse attenuation coefficient *K _{d}*, one of the apparent optical properties (AOPs) of case 1 water, is characterized by [Chl] [23]. For the wavelength of 530 nm,

*K*is given by

_{d}In the case of single scattering, *k _{lidar}* is equal to

*c*[12,13], whereas for the multiple scattering, the express of

*k*is more complicated. In this paper, an improved SALMON is used to trace the multiple scattering behavior of photons in water to get the lidar return profiles and the slope method was introduced to calculate the depth-dependent

_{lidar}*k*.

_{lidar}#### 2.2 An improved semianalytic Monte Carlo radiative transfer simulation model

Monte Carlo method is an a numerical method of solving the radiative transfer which is used in optical communication [24], biological optics [25,26] and optical remote sensing [27]. In 1981, Poole et al. [19] proposed a semianalytic Monte Carlo radiative transfer simulation model (SALMON) for atmosphere lidar, which is widely used in the current oceanic lidars. The improved SALMON proposed here includes a Gaussian beam emission model and it is applicable for both homogeneous and inhomogeneous waters.

### A: Gaussian beam emission model

The beam divergence angle of lidar is always neglected in airborne lidar simulations because it is very small compared with the FOV of airborne lidar receiving system [12,13]. However, it is not negligible in the space lidar simulation, because the platform of spaceborne lidars usually has an altitude of several hundred kilometers; as a result, the spot diameter of the laser beam incident on the sea surface could be several tens of meters. A Gaussian distribution sampling model is introduced in this paper to determine the coordinates and directions of photons incident on the water surface [28].

The intensity perpendicular to the transmission direction of the Gaussian beam is expressed as

*I*is the light intensity, x and y are the coordinates on the plane particular to the transmission direction of the light ray, and ${\sigma}_{s}$ is the root mean square (RMS) radius of the laser beam, commonly called the

*spot size*, which is used to describe the radial distribution of the laser beam. Approximately 86.5% of the total energy of the spot is concentrated in the circle of radius${\sigma}_{s}$ for a Gaussian distribution beam. In order to characterize a Gaussian beam shown in Eq. (9) in our Monte Carlo simulations, a sampling model given in Eqs. (10) and (11) is used to describe the distance from the photon to the center of light spot

*r*and the azimuth angle $\phi $ of photons, respectively. where

*R*and

_{1}*R*are random values between 0 and 1

_{2}Therefore, the coordinates $(x,y,z)$ of photons incident on the water surface can be expressed by

Because the altitude of the spaceborne lidar, *H*, is much greater than the Rayleigh length of the Gaussian beam, the wave front is approximately spherical near the sea surface and the normal direction of each point on the sphere is the direction $({u}_{x},{u}_{y},{u}_{z})$of photons, which is given by

### B: Photon behavior in seawater

The SALMON model simulates the random collision and propagation in water of individual photons from the transmitted beam through a series of absorption and scattering until the photons move outside of the field of view of the lidar receiver or reach the boundaries of medium. For each segment of a photon trajectory, the distance travelled before collision and the propagation direction are determined by making random selections from cumulative probability distribution function of these parameters. These two parameters play a critical role in describing the photon behavior in seawater. The details about the definitions of these two parameters and the formation of the lidar return signals are given in the following paragraphs.

The distance travelled before collision is the step size (*s*) in dimensionless unit and defined as the integration of the attenuation coefficient over the photon pathway [26]. In the SALMON model, $s$ is calculated by $-\mathrm{ln}(R)/c$, where *R* is a random value between 0 and 1. The larger the attenuation coefficient *c*, the shorter distance the photon traveled before it is scattered or absorbed. For the stratified inhomogeneous waters, the IOPs are assigned separately for each layer. If the photon passes through more than one layer, *c* will change along the pathway. Therefore, the actual step size in layer $i\pm 1$could be calculated as ${s}_{i\pm 1}=\left({s}_{i}-{d}_{i}\right){c}_{i\pm 1}/{c}_{i}$,where *d* is the geometrical distance along the ray in layer *i* where the scattering took place. The sign indicates the photon propagates up (-) or down ( + ). The photon travels until the initial step size has been consumed (${s}_{i}<{d}_{i}$).

The propagation direction of the photon after each collision can be expressed by the scattering angle $\theta $ and the azimuth angle$\alpha $. The value of $\alpha $is distributed uniformly between 0 and 2$\pi $, and the distribution of angle $\theta $ is modeled by the scattering phase function [24]. Both the pure water molecules and particles contribute to the total scattering phase function$\tilde{\beta}$ [27], i.e.,

the subscripts*w*and

*p*represent the water molecules and the particles, respectively. The probability that the photon is scattered by molecule is calculated through $\eta \text{=}{b}_{w}/({b}_{w}+{b}_{p})={b}_{w}/b$. And the phase function for water molecules is determined by Rayleigh scattering [22], i.e.,

*g*is the asymmetry parameter controlling the relative amounts of forward and backward scattering in ${\tilde{\beta}}_{HG}$ and is equal to the average of the cosine of the scattering angle over the angular distribution. The HG phase function allows a simple computation of the radiation transfer equation (RTE). However, it only imitates the forward-directed peak of the Petzold’s experimental phase function and describes the backward-directed peak inadequately [24]. Therefore, a weighted sum of HG phase function, called the two-term Henyey-Greenstein (TTHG), is used here to solve the issue, which is given by Kattawar [30]:

In the process of tracking the photons, each photon is conceived as a large packet of identical photons traveling in the same direction. The photon packet has an initial weight $w=1$. After each collision, the current weight is multiplied by the single scattering albedo ${\omega}_{0}$, which is given by ${\omega}_{0}=b/c$. A roulette procedure is employed when the weight *w* is below the threshold of ${\text{10}}^{\text{-4}}$. The expected value of the fraction of photons collected by the lidar receiver without further interactions at each collision point in layer *i* is given by

*j*and ${\theta}^{\prime}$ is the angle between the scattering direction and the direction from photon towards the receiver. The expected value of the fraction of photons reaching the detector is added to the signal; then the packet is reduced by the same value. The remaining photon packet proceeds to the next interaction. All the calculations in this paper are executed with more than 10

^{11}photon packets to ensure a low level of standard deviation of the statistical estimates.

## 3. Results

#### 3.1 Influence of the lidar system parameters on k_{lidar}

The FOV of the lidar receiver is an important parameter that influences the order of scattering in the lidar return signals. For the airborne lidars, the FOV ranges from tens of milliradians to several hundred milliradians depending on different detection purpose. For a spaceborne lidar, smaller FOV is preferred in order to reduce the impact of the background noise from sun light. The FOV of the lidar receiving system should be larger than the divergence angle of the transmitting laser beam due to the divergence of the Gaussian beam. Besides, the altitude of satellite platform also affects the multiple scattering signal. In this paper, we analyze four satellite orbits of the actual spaceborne atmospheric lidar systems: 705 km for CALIOP on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite [9], 550 km for the Japanese Experimental Lidar in Space Equipment (ELISE) [33], 415 km for the Cloud-Aerosol Transport System (CATS), which is a lidar remote sensing instrument on the International Space Station (ISS) [34], as well as 400 km for the Earth Clouds, Aerosol and Radiation Explorer (EarthCARE) satellite to be launched in 2021 [35]. In these simulations, the divergence angle of the laser beam is assumed to be 0.1 mrad and the diameter of the telescope is 1.5 m [9]. The depth resolution of the signals is set to be 0.3 m. The ocean is assumed to have a depth of 50 m with a black bottom, considering the effective attenuation optical depth of three, and a small chlorophyll concentration of 0.35 mg/m^{3} that is uniformly distributed from the surface to the bottom. See Table 1 for the values of input parameters of the calculations.

The *k _{lida}*

_{r}for different system parameters, such as FOV of the receiving telescope and orbit altitude, are compared in Fig. 2. The orbit altitude in Fig. 2(a) is 400 km, and the FOV is in the range of 0.1 to 0.5 mrad, then the corresponding footprint of the receiver is in the range of 20 to 100 m and the diameter of laser spot on the ocean surface is 40 m. Figure 2(a) shows that

*k*at different optical depths decreases with the increase of FOV, but remains almost constant when the FOV is larger than a critical value, indicating that the lidar receives the full-multiple scattering signals. For the orbit altitude of 400 km, FOV of 0.4 mrad is wide enough to satisfy the full-multiple scattering statement at optical depth (${\tau}_{o}$,${\tau}_{o}=c\times z$) from 1 to 5. Figure 2(b) presents the

_{lidar}*k*for four orbit altitudes, namely 400 km, 415 km, 550 km and 705 km, with FOV equals to 0.4 mrad. The result indicates that

_{lidar}*k*of the full-multiple scattering return signals is independent of

_{lidar}*H*. According to our simulation, this result is still valid even when the divergence angle of the incident beam increases to 0.2 mrad because the beam divergence in water is mainly due to the multiple scattering. Therefore,

*k*of spaceborne lidar is only related to the IOPs of water, and system parameters for the full-multiple scattering statement are adopted in the following sections.

_{lidar}#### 3.2 Relationships between IOPs and k_{lidar}

### A. Characteristics of lidar return signals for homogeneous case 1 waters

The effective attenuation coefficient *k _{lidar}* in lidar equation is the reflection of the IOPs of seawater. Determining the value of

*k*is beneficial to improve the retrieve accuracy of IOPs from lidar equation. In terms of the orbit altitude of 400 km, the characteristics of lidar return signals from homogeneous case 1 water with various chlorophyll concentrations are studied. The normalized lidar return signal power and the associated effective attenuation coefficients of eight groups of seawaters with [Chl] between 0 to 2 mg/m

_{lidar}^{3}, which covers most open ocean and productive coastal environment, are provided in Fig. 3. The normalized lidar return profiles in Fig. 3(a) shows that the linearity of the lidar return signal decreases with the increase of [Chl]. For seawaters with high [Chl],

*k*is not constant even if the IOPs are uniform in the vertical direction [Fig. 3(b)].

_{lidar}*k*decreases with the increase of depth at first, then it becomes more stable. The full range of the nonlinearity in

_{lidar}*k*depends on the sensitivity of the receiving system. If this nonlinearity in

_{lidar}*k*is ignored in homogeneous waters, it may cause error in the retrieval of lidar signal.

_{lidar}Figure 4 illustrates the statistical estimates of the normalized lidar return signal power with different order of scattering and the percentage for each order of scattering in the total scattering signal for [Chl] = 0.35 mg/m^{3}. The signal of the first order of scattering decreases exponentially with the increase of depth, and the slope of the signal is *c*. The presence of multiple scattering deforms the waveform of the total lidar return signal, and the rapid increase of high-order scattering signal causes the decrease of *k _{lidar}*. It is shown in Fig. 4(b), that the percentage of high-order scattering in the total signal increases with depth. When the depth is larger than about 10 meters, signal with scattering order more than 1 becomes dominant in the total signal power. Due to the absorption of water, the intensity of the high-order scattering signal declines and

*k*shows a characteristic of stability and invariability.

_{lidar}The nonuniformity of *k _{lidar}* caused by multiple scattering is inevitable and it increases the uncertainty of the retrieval. However, in the real ocean environment, the IOPs is not homogenously distributed. Therefore, it is necessary to make a further research of the relationship between the

*k*(

_{lidar}*z*) and IOPs of water.

### B. The relationships between the *k*_{lidar}(z) and IOPs of water

_{lidar}

Based on the results shown in Fig. 3(b), an exponential formula is proposed to describe the depth characteristic of *k _{lidar}*, which is given by

The fitting results of the effective attenuation coefficients for 16 types of case 1 waters with [Chl] ranging from 0.005 to 2 mg/m^{3} are listed in Table 2. *R*^{2} is the criterion for the goodness of fit. All the *R*^{2} values for the fitting results of each group of seawater are greater than 0.99, which confirms the validity of Eq. (19). As shown in Eq. (19), *k _{lidar}* will be given by is an asymptotic parameter

*p*, when

*z*is large enough to make the exponential term close to zero. The fitting results show that the asymptotic value for each

*k*is a little bit smaller than

_{lidar}*K*.

_{d}The slope and the asymptotic value of *k _{lidar}* in Eq. (19) are determined by

*m*,

*n*and

*p*, respectively. Moreover,

*k*involves the IOPs of seawater, indicating these three parameters are related to the IOPs of seawater.

_{lidar}The influence of the IOPs on the qualitative behavior of *k _{lidar}* is analyzed using a series of calculations with one variable at a time (Fig. 5). Figure 5(a) shows that the curves of

*k*for different${\omega}_{0}$ (with

_{lidar}*b*= 0.2 m

^{−1}) are parallel to each other, which means that the absorption coefficient

*a*only affects the parameter

*p*. Figure 5(b) shows that the scattering coefficient

*b*influences both the slope and the asymptotic value of

*k*. Figure 5(c) shows that the phase function mainly affects the slope of the

_{lidar}*k*but has little influence on the asymptotic value of

_{lidar}*k*. The dependence of

_{lidar}*k*on the phase function and

_{lidar}*b*can be summarized by a backscattering coefficient

*b*, which is obtained by integrating the phase function in the range from $\pi /\text{2}$ to $\pi $.

_{b}In summary, the relationships among the parameters *m*, *n*, *p* and IOPs can be derived based on the fitting results in Table 1 (Fig. 6). Both *m* and *n* show a linear relation with *b _{b}*, and

*p*is a linear function of

*a*and

*b*. The linear functions can be expressed by

_{b}*k _{lidar}* for the phytoplankton bloom case 1 water with [Chl] of 10 mg/m

^{3}and 20 mg/m

^{3}are also calculated, and the fitting results are coincident with the fitted curves in Fig. 6, which broadens the application range of Eq. (20). It is easy to obtain the relationship of IOPs and depth resolved

*k*by substituting Eq. (20) into Eq. (19). For convenience, this function is termed and referred to as

_{lidar}*k*-IOPs model in the following sections. This model depends on the IOPs’ bio-optical model and the scattering phase function of particle, especially the backscattering properties of particle. Since the scattering properties are variable in real waters, the goodness of fit of this

_{lidar}*k*-IOPs model will be sensitive to particles with very strong or weak backward scattering.

_{lidar}## 4. Applicability of the relationship between IOPs and depth resolved *k*_{lidar} in real seawaters

_{lidar}

In the real ocean environment, the distribution of optical properties is normally inhomogeneous, so the validity of the *k _{lidar}*-IOPs model needs to be tested in the practical seawaters. We build a model of inhomogeneous water by defining multiple parallel thin layers, which have the same optical properties in the horizontal direction but different in the vertical direction. The thickness of each layer is assumed to be the same as the resolution of the lidar return profiles. According to Uitz et al. [36], the vertical distribution of chlorophyll

*a*concentration in the global open ocean conforms to the Gaussian distribution given by

_{$\zeta $}is the dimensionless depth obtained by dividing the geometrical depth,

*z*, by the euphotic depth,

*Z*. $C\left(\zeta \right)$ is the dimensionless chlorophyll

_{eu}*a*concentration defined as the ratio of the chlorophyll

*a*concentration (Chl

*a*), to the average concentration within the euphotic layer ${\overline{\text{Chl}a}}_{{Z}_{eu}}$. ${C}_{b}$is the background chlorophyll concentration, which decreases linearly with slope

*s*starting from the surface value ${[Chla]}_{surf}$. ${C}_{\mathrm{max}}$represents the maximum concentration, ${\zeta}_{\mathrm{max}}$is the depth of the maximum concentration and$\Delta \zeta $denotes the width of the peak.

Uitz et al. [36] classified the stratified case 1 waters into nine trophic categories, named from S1 to S9, with respect to the Chlorophyll *a* concentration within the surface layer ${[Chla]}_{surf}$. In order to examine the relationship between IOPs and depth resolved *k _{lidar}* derived in Section 3 in actual seawaters, three trophic categories, with low ${[Chla]}_{surf}$(<0.04 mg/m

^{3}), medium ${[Chla]}_{surf}$(0.2-0.3 mg/m

^{3}), and high ${[Chla]}_{surf}$(2.2-4 mg/m

^{3}), namely S1, S5 and S9, are selected based on the experimental and statistical results of stratified waters. These three waters can represent the trophic states of most open ocean waters in the world. The parameters in Eq. (21) of these three typical trophic class of stratified waters are outlined in Table 3. The distributions of IOPs and Chl

*a*of these three waters are shown in Fig. 7.

Comparisons of the *K _{d}* computed by the bio-optical model given as Eq. (8) and

*k*obtained from the calculation model proposed in this paper with simulated effective attenuation coefficients for S1, S5 and S9 waters are shown in Fig. 8. The lidar is expected to penetrate about 50 m in S1 water, 40 m in S5 water and 20 m in S9 water, with the effective detection depth of three optical thicknesses. For all these three waters, the effective attenuation coefficient

_{lidar}*k*of simulating return signals is found to be different from

_{lidar}*K*(Fig. 8). The values of${k}_{lidar-c}$which is the attenuation coefficient calculated using the

_{d}*k*-IOPs model are very close to the simulated ones.

_{lidar}To quantify the accuracy of *k _{lidar}*-IOPs model, a mean percentage error within the effective detection depth range is defined, which is expressed by

*k*is the parameter to be evaluated,

*n*is the number of depth intervals in the effective detection range. The mean percentage errors between

*k*and

_{lidar-c}*k*are 0.22%, 0.11% and 0.34% for S1, S5 and S9 waters, respectively. To make a primary comparation between the

_{lidar}*k*-IOPs model and

_{lidar}*K*from the bio-optical model, mean percentage errors between

_{d}*K*and

_{d}*k*are given which are 3.07%, 4.05% and 6.24% for S1, S5 and S9 waters, respectively. The accuracy of the

_{lidar}*k*is more than one order of magnitude higher than that of

_{lidar-c}*K*. It indicates that the

_{d}*k*-IOPs model proposed in this paper can precisely describe the relationships between the lidar effective attenuation coefficient and the IOPs of seawater in real ocean environment, which is an improvement of the lidar measurement of IOPs.

_{lidar}## 5. Conclusions

In this paper, the multiple scattering effect on the spaceborne oceanic lidar returns is investigated with an improved SALMON, which includes a Gaussian beam emission model and is applicable for both homogeneous and inhomogeneous waters. The effective attenuation coefficient decays exponentially with the increase of depth in water according to the simulations of lidar return signal from case 1 waters, because high order scattering signals increase with depth and predominate the total signal at shallow depth. An exponential function is proposed to describe the depth characteristic of *k _{lidar}* and the parameters in this exponential function is related to IOPs. The exact function of klidar, depth and the IOPs of seawater is given and examined using real ocean waters simulated by a model of inhomogeneous water.

Comparison of the simulated and calculated results of three typical stratified case 1 waters shows that the errors between *k _{lidar}-c* and

*k*are less than 0.5%. The better agreement shows that the exponential

_{lidar}*k*-IOPs function is feasible for oceanic lidar measurement in theory. Using this function, one can estimate the inherent optical properties and retrieve the chlorophyll concentration of case 1 water, which improves the accuracy of spaceborne lidar data retrieval and assessing global phytoplankton biomass and primary productivity.

_{lidar}## Funding

National Key Research and Development Program of China (2016YFC1400902, 2016YFC1400905, 2016YFC0200700); National Natural Science Foundation of China (NSFC) (41775023, 61475141).

## Acknowledgments

The authors thank Prof. Heather Bouman in the Department of Earth Sciences, University of Oxford for providing the knowledge of chlorophyll concentration.

## References

**1. **J. H. Churnside, “LIDAR detection of plankton in the ocean,” in IEEE International Geoscience and Remote Sensing Symposium (2007), 3174–3177. [CrossRef]

**2. **C. R. McClain, “A decade of satellite ocean color observations,” Annu. Rev. Mar. Sci. **1**(1), 19–42 (2009). [CrossRef] [PubMed]

**3. **C. A. Hostetler, M. J. Behrenfeld, Y. Hu, J. W. Hair, and J. A. Schulien, “Spaceborne Lidar in the Study of Marine Systems,” Ann. Rev. Mar. Sci. **10**, 121–147 (2017).

**4. **M. J. Behrenfeld, Y. Hu, C. A. Hostetler, G. Dall’Olmo, S. D. Rodier, J. W. Hair, and C. R. Trepte, “Space‐based lidar measurements of global ocean carbon stocks,” Geophys. Res. Lett. **40**(16), 4355–4360 (2013). [CrossRef]

**5. **X. Lu, Y. Hu, J. Pelon, C. Trepte, K. Liu, S. Rodier, S. Zeng, P. Lucker, R. Verhappen, J. Wilson, C. Audouy, C. Ferrier, S. Haouchine, B. Hunt, and B. Getzewich, “Retrieval of ocean subsurface particulate backscattering coefficient from space-borne CALIOP lidar measurements,” Opt. Express **24**(25), 29001–29008 (2016). [CrossRef] [PubMed]

**6. **X. Lu, Y. Hu, C. Trepte, S. Zeng, and J. H. Churnside, “Ocean subsurface studies with the CALIPSO spaceborne lidar,” J. Geophys. Res. Oceans **119**(7), 4305–4317 (2015). [CrossRef]

**7. **J. H. Churnside, B. J. Mccarty, and X. Lu, “Subsurface Ocean Signals from an Orbiting Polarization Lidar,” Remote Sens. **5**(7), 3457–3475 (2013). [CrossRef]

**8. **M. J. Behrenfeld, Y. Hu, R. T. O’Malley, E. S. Boss, C. A. Hostetler, D. A. Siegel, J. L. Sarmiento, J. Schulien, J. W. Hair, and X. Lu, “Annual boom-bust cycles of polar phytoplankton biomass revealed by space-based lidar,” Nat. Geosci. **10**, 118 (2016).

**9. **D. M. Winker, M. A. Vaughan, A. Omar, Y. Hu, K. A. Powell, Z. Liu, W. H. Hunt, and S. A. Young, “Overview of the CALIPSO Mission and CALIOP Data Processing Algorithms,” J. Atmos. Ocean. Tech. **26**(11), 2310–2323 (2009). [CrossRef]

**10. **L. Dong, L. Qun, B. Jian, and Z. Yupeng, “Data processing algorithms of the space-borne lidar CALIOP: a review,” Infr. Laser Eng. **46**(12), 1202001 (2017). [CrossRef]

**11. **G. M. Krekov, M. M. Krekova, and V. S. Shamanaev, “Laser sensing of a subsurface oceanic layer. I. Effect of the atmosphere and wind-driven sea waves,” Appl. Opt. **37**(9), 1589–1595 (1998). [CrossRef] [PubMed]

**12. **H. R. Gordon, “Interpretation of airborne oceanic lidar: effects of multiple scattering,” Appl. Opt. **21**(16), 2996–3001 (1982). [CrossRef] [PubMed]

**13. **D. Phillips and B. Koerber, “A Theoretical Study of an Airborne Laser Technique for determining Sea Water Turbidity,” Aust. J. Phys. **37**(1), 75 (1984). [CrossRef]

**14. **J. H. Lee, J. H. Churnside, R. D. Marchbanks, P. L. Donaghay, and J. M. Sullivan, “Oceanographic lidar profiles compared with estimates from in situ optical measurements,” Appl. Opt. **52**(4), 786–794 (2013). [CrossRef] [PubMed]

**15. **J. A. Schulien, M. J. Behrenfeld, J. W. Hair, C. A. Hostetler, and M. S. Twardowski, “Vertically- resolved phytoplankton carbon and net primary production from a high spectral resolution lidar,” Opt. Express **25**(12), 13577–13587 (2017). [CrossRef] [PubMed]

**16. **B. L. Collister, R. C. Zimmerman, C. I. Sukenik, V. J. Hill, and W. M. Balch, “Remote sensing of optical characteristics and particle distributions of the upper ocean using shipboard lidar,” Remote Sens. Environ. **215**, 85–96 (2018). [CrossRef]

**17. **R. E. Walker and J. W. McLean, “Lidar equations for turbid media with pulse stretching,” Appl. Opt. **38**(12), 2384–2397 (1999). [CrossRef] [PubMed]

**18. **V. I. Feygels, “Mathematical modeling of input signals for oceanographic lidar systems,” Proc. SPIE **5155**, 30–39 (2003). [CrossRef]

**19. **L. R. Poole, D. D. Venable, and J. W. Campbell, “Semianalytic Monte Carlo radiative transfer model for oceanographic lidar systems,” Appl. Opt. **20**(20), 3653–3656 (1981). [CrossRef] [PubMed]

**20. **J. H. Churnside, “Review of profiling oceanographic lidar,” Opt. Eng. **53**(5), 051405 (2014). [CrossRef]

**21. **A. Morel, “Light and marine photosynthesis: a spectral model with geochemical and climatological implications,” Prog. Oceanogr. **26**(3), 263–306 (1991). [CrossRef]

**22. **C. D. Mobley, *Light and Water: Radiative Transfer in Natural Waters* (Academic Press, 1994).

**23. **A. Morel and S. Maritorena, “Bio‐optical properties of oceanic waters: A reappraisal,” J. Geophys. Res. Oceans **106**(C4), 7163–7180 (2001). [CrossRef]

**24. **C. Gabriel, M. A. Khalighi, P. Léon, S. Bourennane, and V. Rigaud, “Monte-Carlo-Based Channel Characterization for Underwater Optical Communication Systems,” J. Opt. Commun. Netw. **5**(1), 1–12 (2013). [CrossRef]

**25. **B. H. Hokr, J. N. Bixler, G. Elpers, B. Zollars, R. J. Thomas, V. V. Yakovlev, and M. O. Scully, “Modeling focusing Gaussian beams in a turbid medium with Monte Carlo simulations,” Opt. Express **23**(7), 8699–8705 (2015). [CrossRef] [PubMed]

**26. **L. Wang, S. L. Jacques, and L. Zheng, “MCML—Monte Carlo modeling of light transport in multi-layered tissues,” Comput. Meth. Prog. Bio. **47**, 131-146 (1995).

**27. **K. I. Gjerstad, J. J. Stamnes, B. Hamre, J. K. Lotsberg, B. Yan, and K. Stamnes, “Monte Carlo and discrete-ordinate simulations of irradiances in the coupled atmosphere-ocean system,” Appl. Opt. **42**(15), 2609–2622 (2003). [CrossRef] [PubMed]

**28. **M. Xia, K. Yang, Y. Zheng, and J. Rao, “Influence of Wavy Sea Surface on Airborne Lidar Underwater Beam Quality with Monte Carlo Method,” Chin. J. Lasers **35**(2), 178–182 (2008). [CrossRef]

**29. **L. G. Henyey and J. L. Greenstein, “Diffuse radiation in the Galaxy,” Astrophys. J. **93**, 70–83 (1941).

**30. **G. W. Kattawar, “A three-parameter analytic phase function for multiple scattering calculations,” J. Quant. Spectrosc. Ra. **15**(9), 839–849 (1975). [CrossRef]

**31. **C. D. Mobley, L. K. Sundman, and E. Boss, “Phase function effects on oceanic light fields,” Appl. Opt. **41**(6), 1035–1050 (2002). [CrossRef] [PubMed]

**32. **V. I. Haltrin, “One-parameter two-term Henyey-Greenstein phase function for light scattering in seawater,” Appl. Opt. **41**(6), 1022–1028 (2002). [CrossRef] [PubMed]

**33. **Z. Liu, P. Voelger, and N. Sugimoto, “Simulations of the observation of clouds and aerosols with the Experimental Lidar in Space Equipment system,” Appl. Opt. **39**(18), 3120–3137 (2000). [CrossRef] [PubMed]

**34. **“Official website of Cloud-Aerosol Transport System(CATS),NASA”, retrieved https://cats.gsfc.nasa.gov/.

**35. **“Earth Clouds, Aerosol and Radiation Explorer (EarthCARE), ESA”, retrieved https://earth.esa.int/web/guest/missions/esa-future-missions/earthcare.

**36. **J. Uitz, H. Claustre, A. Morel, and S. B. Hooker, “Vertical distribution of phytoplankton communities in open ocean: An assessment based on surface chlorophyll,” Jo. Geophys. Res. Oceans **111**, C08005 (2006).