Research article 16 Apr 2020
Research article  16 Apr 2020
A new lidar inversion method using a surface reference target applied to the backscattering coefficient and lidar ratio retrievals of a fogoil plume at short range
 ^{1}ONERA, DOTA, Université de Toulouse, 31055 Toulouse, France
 ^{2}Univ. Lille, UMR 8518  LOA  Laboratoire d'Optique Atmosphérique, 59000 Lille, France
 ^{3}LEUKOS SARL, 37 Rue Henri Giffard, 87280 Limoges, France
 ^{1}ONERA, DOTA, Université de Toulouse, 31055 Toulouse, France
 ^{2}Univ. Lille, UMR 8518  LOA  Laboratoire d'Optique Atmosphérique, 59000 Lille, France
 ^{3}LEUKOS SARL, 37 Rue Henri Giffard, 87280 Limoges, France
Correspondence: Florian Gaudfrin (florian.gaudfrin@onera.fr) and Nicolas Riviere (nicolas.riviere@onera.fr)
Hide author detailsCorrespondence: Florian Gaudfrin (florian.gaudfrin@onera.fr) and Nicolas Riviere (nicolas.riviere@onera.fr)
In this paper, a new elastic lidar inversion equation is presented. It is based on the backscattering signal from a surface reference target (SRT) rather than that from a volumetric layer of reference (Rayleigh molecular scatterer) as is usually done. The method presented can be used when the optical properties of such a layer are not available, e.g., in the case of airborne elastic lidar measurements or when the lidar–target line is horizontal Also, a new algorithm is described to retrieve the lidar ratio and the backscattering coefficient of an aerosol plume without any a priori assumptions about the plume. In addition, our algorithm allows a determination of the instrumental constant. This algorithm is theoretically tested, viz. by means of simulated lidar profiles and then using real measurements. Good agreement with available data in the literature has been found.
Atmospheric aerosols are liquid or solid particles dispersed in the air (Glickman and Zenk, 2000) of natural (volcano, biomass burnings, desert, ocean) or anthropogenic origins. They play an important role not only in cloud formation (DeMott et al., 2010), in radiative forcing (Hansen et al., 1997; Haywood and Boucher, 2000) and more generally for research on the climate change but also in the context of air quality and public health (Baltensperger et al., 2008; FinlaysonPitts and Pitts, 2000; Popovicheva et al., 2019; Zhang et al., 2018). Their size varies from the nanometer to the millimeter scale (Robert, unpublished). However, a large majority of aerosols have a size between 0.01 and 3 µm (Clark and Whitby, 1967) for which scattering is dominant in the optical domain. The Mie theory is often used, at least statistically (i.e., for a large population of random sized aerosols), although aerosols are not always spherical. The optical backscattering and extinction properties of aerosols are mainly related to their shape (Ceolato et al., 2018), size distribution (VargasUbera et al., 2007), concentration and chemical composition, which is based to their nature (dust, maritime or urban). Lidars are active remote sensing instruments suitable for aerosol detection and characterization (Sicard et al., 2002) over kilometric distances during both day and nighttime.
The optical properties of aerosols are obtained by means of inversion methods using the simple scattering lidar equation. In the 1980s, a stable onecomponent formulation adapted to lidar applications was proposed by Klett (1981). It has then been extended to a twocomponent formulation, viz. separating molecular and aerosol contributions, by Fernald (1984) and Klett (1985). The elastic lidar equation is an illposed problem in the sense of Hadamard (1908) since one searches for extinction and backscattering coefficients with only a single observable. Several assumptions are therefore required in order to invert the lidar equation.
 i.
A calibration constant is usually determined from a volumetric layer of the upper atmosphere as a reference target (Vande Hey, 2014). This calibration layer can be very high in altitude; it has recently been moved from the around 32 km to around 36–39 km for the CALIPSO spaceborne lidar in order to reduce uncertainties in the inversion procedure (Kar et al., 2018; Getzewich et al., 2018). This volume is considered to be made only of pure molecular constituents whose optical scattering properties are well known (Rayleigh regime). The molecular backscattering coefficient is generally estimated from the standard model of the atmosphere (Anon, 1976; Bodhaine et al., 1999). However, poor estimates of the reference or low signaltonoise ratios (SNRs) can lead to severe uncertainties in the retrieved extinction and backscattering coefficients. Few sensitivity studies have been performed to evaluate such uncertainties (Matsumoto and Takeuchi, 1994; Rocadenbosch et al., 2012). Spatial averaging around the volume of reference in addition to time averaging is thus recommended to increase SNR.
 ii.
Lidar ratio is constant over the distance range of measurements (Sasano et al., 1985). This is also an important source of errors in the retrieval values. Some studies have proposed a variable lidar ratio under the form of a powerlaw relationship between the extinction and backscattering coefficients, but such a method requires a priori knowledge of the medium under study (Klett, 1985).
 iii.
The molecular contribution along the lidar line is known. It is estimated, as for the backscattering coefficient, by means of temperature and pressure vertical profiles, using either the standard model of the atmosphere or radio soundings (Jäger, 2005).
In the case of elastic lidar inversion, the most critical parameter is the lidar ratio (LR). It depends on the wavelength (in vacuum) and the microphysics, morphology and size of the particles (Hoff et al., 2008). The LR ranges from 20 to 100 sr at 532 nm (Ackermann, 1998; Cattrall et al., 2005; Leblanc et al., 2005) according to the aerosol origins (maritime, urban, dust particles and biomass burning). It is therefore difficult to assume an a priori value for LR inasmuch as this information is to be found rather than given.
Several alternatives have been analyzed to constrain the inversion procedure while relaxing assumption (ii). These alternatives are based on the determination of the optical thickness, the one which consists of coupling lidar and sun photometer measurements being the most widely used. The measured optical thickness is then used to constrain extinction profiles (Fernald et al., 1972; Pedrós et al., 2010). A second alternative consists of combining elastic lidar and Raman measurements in order to get the optical depth as a function of range (Ansmann et al., 1990, 1992, 1997; Mattis et al., 2004). In a third technique, the optical depth is retrieved from elastic lidar measurements with different zenith angles (Sicard et al., 2002). It is worth indicating that coupling lidar and sun photometer measurements is possible only for daytime, while Raman measurements are carried out preferentially at nighttime in order to increase the SNR. A fourth method consists of the determination of the optical thickness and lidar ratio of transparent layers located above opaque clouds (Hu et al., 2007; Young, 1995) that are used as reference for calibration in the inversion procedure (O'Connor et al., 2004). This method is used for downlooking lidar measurements capable of measuring depolarization ratios. However, the method is limited to lidar systems in nonpolarized detection and for lidar measurements for which clouds cannot be used as a reference. A fifth approach consists of the determination of the optical thickness of the atmosphere from the sea surface echo by combining lidar and radar measurements (Josset et al., 2010a, b, 2008). This method has been used to find the lidar ratio and the optical depth of aerosol layers over oceans (Dawson et al., 2015; Josset et al., 2012; Painemal et al., 2019).
Another limitation of groundbased lidar measurements is related to the overlap function, which strongly impacts (and prevents) observation close to the instrument, i.e., in the lowest layers of the troposphere where aerosols are emitted. Different studies have proposed to modify the overlap function analytically (Comeron et al., 2011; Halldórsson and Langerholc, 1978; Kumar and Rocadenbosch, 2013; Stelmaszczyk et al., 2005) or empirically (Vande Hey et al., 2011; Wandinger and Ansmann, 2002). Some lidar devices are also equipped with a second telescope of higher overlap at short range (Ansmann et al., 2001). However, current lidar systems are not adapted enough to the monitoring and characterization of volumetric targets at short range, for instance in the industrial context or more generally, for anthropogenic activities (Ceolato and Gaudfrin, 2018).
To meet new industrial emission control requirements and the very recently emitted anthropogenic aerosols characterization, we have developed a shortrange lidar of high spatial resolution (Gaudfrin et al., 2019, 2018b). The lidar inversion cannot be performed by means of the classical Klett–Fernald equation, because the reference layer used for the inversion is either impossible to access (horizontal lidar measurements and skytoground lidar airborne measurements) or inaccessible because of a finite lidar range. In the present paper, a modification of the conventional lidar equation is proposed in order to perform lidar inversions using a surface reference target (SRT) at relatively short range (r_{max}≈100 m). Precisely, a unified lidar equation for surface and volumetric scattering media is suggested, and it is then used for a new inversion equation, inspired by the Klett–Fernald equation, using a SRT.
Also a new technique to retrieve the lidar ratio without using any sun photometer, Raman or radar measurements is presented and applied to an aerosol plume. This new inversion technique is both assessed theoretically and experimentally using real lidar measurements. A discussion and a conclusion follow and close the present paper.
Currently, lidar inversion methods use a volumetric layer of the upper atmosphere (higher than an altitude of 8 km above ground level) as a reference target. This volume is considered as being free of aerosols and made only of pure molecular constituents whose optical scattering properties are known. In our approach, we propose using a SRT of known bidirectional reflectance distribution function (BRDF) f_{r,λ} (sr^{−1}) (Nicodemus, 1965; Kavaya et al., 1983).
This requires modification of the usual lidar equation in order to make it suitable for both surface and volumetric targets.
For the singlescattering lidar equation, for which light has undergone only one scattering event, the measured backscattered power, at range r, can be written in a general way, viz. by considering both a surface target (Bufton, 1989; Hall and Ageno, 1970) and a volumetric target (Collis and Russell, 1976), as
where 𝒫_{p,λ} (W) is the peak power of the laser source, $c\approx \mathrm{3}\times {\mathrm{10}}^{\mathrm{8}}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$ is the Einstein constant, τ_{λ} (s) is the laser pulse duration (full width at half maximum) and A_{ef} (m^{2}) is the telescope effective receiving area θ_{i} (rad) the angle between the normal eigenvector to the SRT and the incident beam direction. It should be noted that in the particular case of a Lambertian surface, ${f}_{\mathrm{r},\mathit{\lambda}}({r}_{\mathrm{s}},\phantom{\rule{0.125em}{0ex}}{\mathit{\theta}}_{i})$ can be easily expressed by spectral bidirectional reflectance factor ρ_{λ} from ρ_{λ} cos θ_{i}∕π (Josset et al., 2018, 2010b; Haner et al., 1998). However, the general form of BRDF (f_{r,λ}) will be considered later in this work in order to not restrict the approach to specific cases. Also, the SRT is located at range r_{s}, with ξ_{λ} the dimensionless overlap function and η_{λ} the dimensionless optical efficiency of the whole receiver. 𝒫_{p,λ} is a rectangularshaped pulse in the volumic lidar equation (Measures, 1992), viz. the ratio between the pulse energy and τ_{λ}. In the case of lidar measurements on a SRT, the backscattered peak power is not proportional to 𝒫_{p,λ}. A corrective factor F_{cor} depending on the real shape of the laser pulse is thus introduced. In the present case, ${\mathcal{P}}_{\mathrm{p},\mathit{\lambda}}^{\mathrm{G}}={\mathcal{P}}_{\mathrm{p},\mathit{\lambda}}^{\mathrm{s}}\phantom{\rule{0.125em}{0ex}}{F}_{\mathrm{cor}}$, with ${\mathcal{P}}_{\mathrm{p},\mathit{\lambda}}^{\mathrm{G}}$ and ${\mathcal{P}}_{\mathrm{p},\mathit{\lambda}}^{\mathrm{s}}$ the peak powers of a Gaussianshaped and a square laser pulse, respectively. Conservation of the pulse energy between these two kinds of pulses gives ${F}_{\mathrm{cor}}=\mathrm{2}\phantom{\rule{0.125em}{0ex}}(\mathrm{ln}\mathrm{2}/\mathit{\pi}{)}^{\mathrm{1}/\mathrm{2}}$ (Paschotta, 2008). The factor does not apply to the volume part of the lidar equation because, in this last part, the pulse profile is assumed to be constant over a rate duration τ_{λ}. This approximation cannot be made for the backscatter peak of a SRT, because the backscattered energy is not integrated over a volume.
In Eq. (1), ${T}_{\mathit{\lambda}}^{\mathrm{2}}$ the back and forth atmospheric transmission throughout the environment between the lidar source and range r (Swinehart, 1962) is
with α_{λ} (m^{−1}) being the total extinction coefficient at wavelength λ and range r, as ${\mathit{\alpha}}_{\mathit{\lambda}}={\mathit{\alpha}}_{\mathrm{b},\mathit{\lambda}}+{\mathit{\alpha}}_{\mathrm{a},\mathit{\lambda}}$. The subscripts “b” and “a” refer, respectively, to the contribution of the background (molecules and aerosols) and to the contribution of the aerosol volumetric target under investigation. The total backscattering coefficient β (m^{−1} sr^{−1}) is ${\mathit{\beta}}_{\mathit{\lambda}}={\mathit{\beta}}_{\mathrm{b},\mathit{\lambda}}+{\mathit{\beta}}_{\mathrm{a},\mathit{\lambda}}$, with the same meaning as just above for the subscripts. By definition, the corresponding lidar ratios are LR${}_{\mathrm{b},\mathit{\lambda}}\left(r\right)={\mathit{\alpha}}_{\mathrm{b},\mathit{\lambda}}/{\mathit{\beta}}_{\mathrm{b},\mathit{\lambda}}$ and LR${}_{\mathrm{a},\mathit{\lambda}}\left(r\right)={\mathit{\alpha}}_{\mathrm{a},\mathit{\lambda}}/{\mathit{\beta}}_{\mathrm{a},\mathit{\lambda}}$.
The fundamental quantity measured by the lidar instrument is a voltage V (in volts) which is proportional to the backscattered power as follows: ${V}_{\mathit{\lambda}}\left(r\right)={R}_{\mathrm{v},\mathit{\lambda}}\phantom{\rule{0.125em}{0ex}}{\mathcal{P}}_{\mathit{\lambda}}\left(r\right)$, where R_{v,λ} is the detection constant (V W^{−1}) which determines the light–voltage conversion. It can be written using the instrumental constant as follows: ${C}_{\mathrm{ins}}={R}_{\mathrm{v},\mathit{\lambda}}{K}_{s}$ (V m^{3}), where ${K}_{s}={\mathcal{P}}_{\mathrm{p},\mathit{\lambda}}c{\mathit{\tau}}_{\mathit{\lambda}}{A}_{\mathrm{ef}}\mathit{\eta}/\mathrm{2}$. In the literature, C_{ins} is obtained from 𝒫_{λ}, while, herein, it comes from the voltage and therefore takes into account all the emission, collection, detection and acquisition chain.
In the sequel, for better readability, the subscript λ and θ_{i} will not be written thereafter.
The rangecorrected lidar signal V_{λ}(r)r^{2} is as follows:
To remove the α dependence in the exponential term, we will replace α_{a} and α_{b} with LR_{a} and LR_{b}, respectively, and introduce the term
as detailed in Ansmann and Müller (2004). With such modifications, the final lidar equation for surface and volumetric scatterers can thus be written as
with $Y\left(r\right)={\mathrm{LR}}_{\mathrm{a}}\left(r\right)\phantom{\rule{0.125em}{0ex}}\left[{\mathit{\beta}}_{\mathrm{b}}\right(r)+{\mathit{\beta}}_{\mathrm{a}}(r\left)\right]$.
Thereafter, in order to highlight the expression to solve, it is convenient to define backgroundcorrected transmission factor as
and W(r)=S(r)LR_{a}(r)D(r). Finally, Eq. (3) becomes
We will now introduce the lidar framework adapted to the radiative parameter retrieval of a volumetric scattering medium with a known SRT.
3.1 Radiative parameters identification
The current Klett–Fernald inversion method consists of determining C_{ins} using the high atmosphere as a reference and fixing the LR_{a} a priori. In this paper, C_{ins} is determined by means of a SRT located at range r_{s}. So,
It is worth mentioning that LR_{a}(r_{s}) is the lidar ratio just before the SRT and Y(r_{s})=0 (only at the SRT). Also, obviously, for r<r_{s}, f_{r}=0. Inserting Eq. (8) into Eq. (7) gives
This equation applies only before the SRT and can be solved by integrating both sides from r to r_{s} (Vande Hey, 2014). The exponential term is (see Appendix) as follows:
Plugging Eq. (10) into Eq. (9), we obtain in the following:
Using the definitions of Y(r) and W(r) (see above), β_{a}(r) can be written as
This equation can be written as
Then, by definition of the lidar ratio, we deduct α_{a}(r)=LR_{a}(r)β_{a}(r). Equation (13) is similar to the one defined by Klett (1981), except that β_{b} in Eq. (13) also contains the contribution of the aerosol background.
Assuming that the properties of the SRT are well known, the most critical parameter is LR_{a}(r). Giving a value for LR_{a} requires a priori knowledge of the volumetric target under study, whereas the main objective of lidar remote sensing is precisely to characterize the medium investigated. A priori information is always a topic of discussions, and it causes more or less severe flaws in lidar measurements.
Equation (13) can also be applied to the important context of airborne observations. In this case, it is necessary to know the ground BRDF.
3.2 Determination of LR_{a} and β_{a}: methodology
The objective is to first retrieve β_{a}(r) and LR_{a} – and then to deduce α_{a}(r) – without any a priori knowledge about the medium considered. Two lidar measurements are performed as follows: the first one (signal V_{s}) in the absence of the volumetric aerosol medium of interest and a subsequent one (signal V_{sv}) in its presence. The SRT is obviously present for both measurements. The two measurements should be performed close in time in order to avoid the background environment evolving too much. The experimental setup of these lidar measurements is illustrated in Fig. 1.
By definition, the halflogarithmic ratio of S_{s} and S_{sv} corresponds to the total extinction of the volumetric media under study: ${\mathit{\alpha}}_{\mathrm{tot}}=\mathrm{ln}\left[{S}_{s}\right({r}_{\mathrm{s}})/{S}_{\mathrm{sv}}({r}_{\mathrm{s}}\left)\right]/\mathrm{2}$. Using S_{s}, C_{ins} can be determined independently of the volumetric medium of interest as follows:
which is Eq. (8) with α_{a}=0. C_{ins} and α_{tot} are used in objective functions to retrieve LR_{a} and are assumed to be uniform – r independent. The first objective function is
where α_{a} is the retrieved profile of extinction using Eq. (13) and LR_{a}. The medium is assumed to be at a range of full overlap (r>r_{0}), so that α_{tot} must correspond to the integrated extinction. A second objective function,
is introduced in order to minimize the difference between S_{sv} and the simulated signal S_{sim} obtained from the retrieved β_{a} and α_{a} and from C_{ins}.
The methodology is presented in Fig. 2. The molecular background contribution is computed from pressure and temperature data as in Bucholtz (1995), while the aerosol background contribution is estimated by means of radiative transfer codes, e.g., MATISSE (Simoneau et al., 2002; Labarre et al., 2010) or MODTRAN (Berk et al., 2008, 2014). Another solution consists of using a realistic value of the visibility 𝒱 (km^{−1}) and Koschmieder's relation (Horvath, 1971; Elias et al., 2009; Hyslop, 2009) at 550 nm (maximum human eye sensitivity): 𝒱α_{b}≈3.9.
The signals V_{s} and V_{sv} are introduced in the inversion procedure, which is organized around three main steps (Fig. 2):

A Gaussian fit is first applied to the backscattered signal from the SRT, i.e., V_{s}(r_{s}) and V_{sv}(r_{s}), which gives the amplitude of the backscattering, the position of this peak and its width in position. From these Gaussian models, one can obtain α_{tot} (from its definition; see above) and C_{ins} from Eq. (14). Note: when the target is tilted with respect to the lidartarget line, the backscatter peak of surface target will not be symmetrical. Another fit should be used as a lognormal function.

A first lidar inversion is realized using Eq. (13) with LR_{a}=50 sr at the beginning of the inversion procedure. This value has been chosen because it corresponds to the average LR_{a} data from the literature. For that, the Gaussian model V_{sv} obtained at step 1 is used for signal S(r_{s}) in Eq. (13). A first rangeprofile β_{a}(r) is thus obtained at the end of this second step.

The above β_{a}(r) and LR_{a} allow us to determine α_{a}(r) whose r integration is then compared with α_{tot} in the minimization procedure of Eq. (15). At each iteration, the LR_{a} is modified in order to reduce ε_{1}. The new β_{a}(r), LR_{a} and α_{a}(r) values are then used to compute a simulated lidar signal S_{sim} whose comparison with S_{sv} is minimized according to Eq. (16). In this algorithm, the iterative procedure ends when ${\mathit{\epsilon}}_{\mathrm{1}}+{\mathit{\epsilon}}_{\mathrm{2}}\le {\mathrm{10}}^{\mathrm{6}}$ is reached. The use of 19 iterations is generally enough, depending on the first value of LR_{a} introduced initially (step 2). At the end of this step, one thus obtains final β_{a}(r), α_{a}(r) and LR_{a}. The minimization procedure used is the one implemented by Kraft (1988). Equation (15) is the most important since it determines the rapidity of convergence. Equation (16) is helpful but not critical.
4.1 Theoretical lidar signals
The inversion method described above is tested using theoretical lidar signals generated by PERFALIS (PERFormence Assesment for LIdar Systems; Gaudfrin et al., 2018a). As summarized in Table 1, the simulated atmosphere is composed of three layers and a SRT of BRDF ${f}_{\mathrm{r}}=\mathrm{0.20}/\mathit{\pi}$ located at r_{s}=100 m. Pressure and temperature are uniform (1040 hPa and 290 K), and the continental aerosol background is chosen so that it corresponds to 𝒱=47 km (Hess et al., 1998). In addition, ${\mathit{\beta}}_{\mathrm{b}}=\mathrm{9.97}\times {\mathrm{10}}^{\mathrm{6}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{sr}}^{\mathrm{1}}$ and LR_{b}=118.56 sr. The signal V_{s} is generated from the background components and the SRT, while the signal V_{sv} is generated considering an aerosol plume aerosol between 20 and 30 m (second layer). The plume backscatter coefficient is ${\mathit{\beta}}_{\mathrm{a}}=\mathrm{7.14}\times {\mathrm{10}}^{\mathrm{5}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{sr}}^{\mathrm{1}}$ and LR_{a}=70 sr. Multiple scattering is assumed to be negligible. For dense atmosphere and wider field of view, Eq. (1) has to be corrected by an appropriate factor (Bissonnette, 1996) in order to consider higher orders of scattering events.
Inversion methods are generally applied to averaged signals in order to increase the SNR. In lidar remote sensing, the noise can be, approximately, considered as a white Gaussian noise (Li et al., 2012; Mao et al., 2013; Sun, 2018). In order to assess the impact of noise in the inversion method (see Sect. 3), a Gaussian noise of null mean value and a standard deviation of $\mathrm{1.5}\times {\mathrm{10}}^{\mathrm{5}}\phantom{\rule{0.125em}{0ex}}\mathrm{a}.\mathrm{u}.$ is introduced in the theoretical lidar signals. Figure 3 displays the theoretical noised signals V_{s} and V_{sv}. As expected, because of light extinction by the plume, V_{sv}(r_{s}) is lower than V_{s}(r_{s}) by 9 %. Four datasets are then generated, with respectively, an averaging over 20, 50, 100 and 200 signals, from V_{s} and V_{sv}, and, in addition, a fifth signal without noise is considered (Fig. 4).
4.2 Noise impact on β_{a} and LR_{a} retrievals
LR_{a} is retrieved using Eq. (15). In addition to the sixlidar dataset described above, four different conditions of inversion are considered. In condition 1 the exact data of the background components are used as an input in the inversion algorithm. For conditions 2 and 3, β_{b} is over and underestimated by 20 % compared to the data used to generate the theoretical signals. In conditions 1 to 3, the inversion technique is performed over the entire signal range. Condition 4 is the same as condition 1, but the aerosol plume is spatially delimited. Table 2 summarizes the four conditions for the six datasets. It is worth mentioning that noised lidar signals obviously result in noised retrieved β_{a}(r). Thus, to quantify the performance of the inversion technique, we consider the average value $\stackrel{\mathrm{\u203e}}{{\mathit{\beta}}_{\mathrm{a}}}$ of the plume. The retrieved value of LR_{a} can be directly compared to the theoretical value.
Figure 5 displays $\stackrel{\mathrm{\u203e}}{{\mathit{\beta}}_{\mathrm{a}}}$ for the six datasets and the four inversion conditions. It varies from $\mathrm{7.11}\times {\mathrm{10}}^{\mathrm{5}}$ to $\mathrm{7.22}\times {\mathrm{10}}^{\mathrm{5}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{sr}}^{\mathrm{1}}$, which means an error of approximately 1 % in comparison to the theoretical value. Conditions 2 and 3 result in a translation of the corresponding curve of ±0.4 % with respect to the curve associated with condition 1 because of the over and underestimation of 20 % introduced in β_{b}. The performance is better for condition 4 whatever the dataset is, since the maximum error is 0.5 % for noised signals. The spatially bounded aerosol layer is often applied in inversion methods and seems to herein improve the inversion method.
For signal lidar without noise, ${\stackrel{\mathrm{\u203e}}{\mathit{\beta}}}_{\mathrm{a}}$ is not exactly equal to the theoretical value, maybe because of numerical computation errors in the inversion algorithm. Such a numerical error is about 0.12 % (condition 1) and 0.04 % (condition 4).
Figure 6 is similar to Fig. 5 but considers LR_{a}. One obtains values ranging from 66 to 74 sr, with a maximum error of 5 % compared to the theoretical value. In conditions 1, 2 and 3, using averaged noised signals has no consequence on the retrieved value of LR_{a}, contrary to what was obtained for ${\stackrel{\mathrm{\u203e}}{\mathit{\beta}}}_{\mathrm{a}}$.
In condition 1, the maximum error is 2.1 %. The graphs corresponding to conditions 2 and 3 are translated with respect to the graph under to condition 1 by about ±3 % and permuted by the same but for ${\stackrel{\mathrm{\u203e}}{\mathit{\beta}}}_{\mathrm{a}}$. Nevertheless, the errors remain low with a maximum of 5 % (condition 2) if 50 signals are averaged. However, under condition 4, the LR_{a} is much better for averaged signals and remains quite good for the noised signal (not averaged) with an error rate of 0.6 %. Again, it seems that the spatial limitation of the plume increases the accuracy of the retrieval LR_{a}. Condition 1 remains, however, efficient for noised signals since deviation is below 2.1 %. In the case of the lidar signal without noise, the retrieved LR_{a} values are not exactly equal to the theoretical LR_{a} values; numerical computation errors are about 0.13 % (condition 1) and 0.05 % (condition 4). An error of ±20 % in β_{b} introduced initially will result in an under or overestimation LR_{a} by ±3 %. Condition 4 is preferable for retrieving LR_{a}.
Note that the formalism and methodology adopted here to retrieve the lidar ratio are efficient as long as the peak backscattering of the SRT is present on the lidar signal. The method has been evaluated, in this paper, for short range, around 100 m, because our research focus on application at this range. However, the algorithm developed does not present any limit with respect to the range provided that measurements are made below 1 km of range (this value depends of the power of laser sources) with respect to our applications. However, at first sight there is no limit to the application of the method for measurements at longer ranges such as those at more than 1 km.
4.3 Plume optical property retrieval
The above study allowed us to test the new inversion method on noised signals, for different conditions of inversion, as a function of the number of signals averaged. Thereafter, lidar inversion was performed by considering a spatially bounded plume and 100 signals for averaging. This last condition has been chosen because it corresponds to the number of signals available in less than 0.1 s with our lidar system (see Sect. 5). The theoretical results obtained by the inversion method with 100 averaged signals is also quite good (see above). Figure 7 displays the retrieved β_{a} if a theoretical lidar signal is introduced as a first guess. Table 3 lists the retrieved ${\stackrel{\mathrm{\u203e}}{\mathit{\beta}}}_{\mathrm{a}}$ and LR_{a}. Compared to theoretical values, errors are less than 0.7 % for LR_{a} and below 0.1 % for ${\stackrel{\mathrm{\u203e}}{\mathit{\beta}}}_{\mathrm{a}}$, although a peak of 2.2 % is observed at r=28.8 m.
Our new inversion technique is now applied to real lidar measurements. The instrument used is named COLIBRIS (COmpact LIdar for Broadband polaRImetric multiStatic measurements; Gaudfrin et al., 2018b; Ceolato and Gaudfrin, 2018). This lidar is able to perform shortrange measurements (r_{0}<5 m) at high spatial resolution (lower than 0.25 m). A Nd:YAG microchip laser source from HORUS, part of the LEUKOS company, is used with a pulse energy peaking at 532 nm of 7.3 µJ and a repetition rate of 1 kHz. The backscattered light is collected by a Cassegrain telescope. In the detection part, a dichroic filter for the elastic channel is used before a photomultiplier tube. The signal is digitized at a sample frequency of 3 GHz after being amplified.
5.1 Description of the experimental operations
The lidar measurements are performed horizontally as illustrated in Fig. 8. A Lambertian Zenithal SRT (SphereOptic, 2020) with a ${f}_{\mathrm{r}}=\mathrm{0.20}/\mathit{\pi}$ is placed 52 m away from the source. Its spectral bidirectional reflectance has been checked using laboratory bench measurements (Ceolato et al., 2012). The mean direction of the laser beam is parallel to the normal of the surface.
The repetition laser source has repetition frequency of 1 kHz. In order to increase the SNR, we preprocess the measurements from three lidar measurements:

signal 1 – the first measurement is made by occulting the emitted laser beam to get a measure of the background scene (contribution of passive illumination);

signal 2 – the second measurement is made by occulting the telescope to estimate the dark noise of the instrument;

signal 3 – the last measurement is made without any occultation.
For a given acquisition period, these three series of measured signals are averaged. The averaged signals of the background radiation and the dark noise (signals 1 and 2) are then subtracted from the signal 3 as follows: signal 3 − signal 2 − (signal 1 − signal 2).
Figure 9 shows the lidar results on a volume/surface target over a period of 2 s. This corresponds to 2000 signals per series of measurements. During this period, we assume that the environment does not evolve significantly. The curves V_{sv} and V_{s} are the measurements in the presence and in the absence of oil smoke with SRT, respectively. The oil plume signal is visible between 37.5 and 41 m.
The highspeed sampling allows a measurement every centimeter along the line of sight. Combined with a short pulse duration of the laser source (1.7 ns), this makes it possible to highlight local variation concentrations on the order of 25 cm inside the plume with the presence of two maxima at 38 and 39 m from lidar. The peak of backscattering of the SRT is also well sampled. The signal amplitude corresponding to the backscatter of the SRT is lower for V_{sv} than for V_{s} because of the presence of the oil plume.
During measurement, the pressure, temperature and visibility are respectively 1016 hPa, 288 K and 30 km. These data are used to compute β_{b} as described in Sect. 4.
5.2 Optical property retrieval: fogoil plume
The signals V_{s} and V_{sv} are used in the inversion procedure as described in Sects. 3 and 4. The plume is spatially bounded (condition 4).
The retrieved β_{a}(r) is displayed in Fig. 10. In the densest range of the plume ${\mathit{\beta}}_{\mathrm{a}}\approx \mathrm{2}\times {\mathrm{10}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{sr}}^{\mathrm{1}}$. Also, the retrieved LR_{a} is around 98 sr. According to Bohlmann et al. (2018), this value corresponds, as expected, to smoke particles (at 532 nm, the lidar ratio ranges from 80 to 100 sr).The optical properties of the oilfog plume of experimental retrieved with inverse method are summarized in Sect. 4.
The lidar signal reproduced from the retrieved β_{a}(r), LR_{a} and the instrumental constant deduced from the Eq. (14) gives a standard deviation from the exact value of $\mathrm{1.5}\times {\mathrm{10}}^{\mathrm{5}}\phantom{\rule{0.125em}{0ex}}\mathrm{a}.\mathrm{u}.$ This shows the consistency and reliability of the new inversion method proposed in this paper.
In this paper, a new method has been introduced for lidar measurement inversion in a situation for which a volumetric layer (molecular Rayleigh scatterers) of the high troposphere is not available (e.g., airborne lidar observations, horizontal configuration of measurements). This method is based on a new expression of the lidar equation which allows us to use a surface reference target of a known BRDF instead of a volumetric one. This new formalism permits the inversion of shortrange lidar measurements for which conventional inversion techniques cannot directly be applied. Similarly to common inversion techniques, our method requires the introduction of a background component (molecular and particulate contributions) that can either be estimated from radiative models or deduced from measurements of temperature, pressure and visibility conditions.
Also, a new algorithm has been developed to retrieve, without any a priori assumptions relative to the medium to be characterized (aerosol plume), the backscattering coefficient (β_{a}) and lidar ratio (LR_{a}) of an aerosol plume between the lidar and the surface target reference. In other words, our technique does not need to introduce a lidar ratio as an input for our inverse algorithm. For that, two lidar measurements are necessary, with and without the aerosol plume under consideration.
Comparing these two signals, one can retrieve the total extinction coefficient of the medium analyzed and the instrumental constant of the lidar instrument. These two pieces of information are used to constrain the inversion algorithm and finally to identify LR_{a}.
This algorithm has been first investigated using theoretical (simulated) lidar signals. The quality of the retrieval has been assessed by introducing noise in the simulated signals and by considering various conditions of inversion differing, in particular, from one another according to the initial error introduced in the backscattering coefficient of the aerosol background. Thus, the robustness of algorithm has been shown, since in all the cases, the error on the retrieved values (viz. in β_{a} and LR_{a}) is less than 5 %. Also, we have found that inversion is better for spatially bounded aerosol plumes.
The inversion algorithm has then been applied to real lidar shortrange measurements of an oilfog plume. The retrieved β_{a} and LR_{a} of the plume agree with values found in the literature for smokelike particles. Moreover, thanks to the determination of the instrumental constant, the measured signal has been computed from the inverted products, and an absolute error of ${\mathrm{10}}^{\mathrm{5}}\phantom{\rule{0.125em}{0ex}}\mathrm{a}.\mathrm{u}.$ between the measure and the postprocessed simulation has been encountered.
However, it is worth mentioning that the method proposed herein to find LR_{a} has some limitations. Precisely, the sensitivity of the lidar must be sufficient to detect the signal of weakly thick or weakly backscattering plume. Indeed, since measurements are performed in the absence and in the presence of the medium, by means of a hard surface target of reference of known reflectance, the algorithm converges less easily for very weakly diffusing plumes.
The new inversion technique presented in this paper suggests new airborne lidar applications operated at low altitude from aircraft (helicopters, airplanes) but requires a priori knowledge of the reflectance of the SRT.
Even if some models exist for the BRDF of surfaces (Bréon et al., 2002; Lobell and Asner, 2002; Mishchenko et al., 1999), their use seems difficult to implement because of the diversity of encountered surfaces during airborne measurements. Nevertheless, it may be possible to identify the reflectance of the ground surface by means of a spectroradiometer imager (Poutier et al., 2002; Miesch et al., 2005; Josset et al., 2018). The combination of these measurements with the herein proposed inversion method would a priori be complementary to establishing of new methods of calibration for downlooking lidar measurements (spaceborne or airborne lidars). The evaluation of the method proposed in this paper, considering the uncertainty in the target reflectance, has not been performed. It will be the topic of future works.
To solve Eq. (9), the exponential term can be written under another form. The method proposed by Vande Hey (2014) consists of integrating both members of the equation from r to r_{s}. So,
Since
where F is the primitive of Y. It ensues that
Therefore,
Finally, the exponential term becomes
These data are the property of ONERA and are not accessible.
This work is part of the PhD thesis of FG who is the main contributor of all of the aspects (theoretical, modeling and experimental) of this article. This thesis has been supervised by NR, OP, RC and GH. OP contributed to the theoretical part of this work while experiments have been conducted by FG, RC, NR and GH. FG and RC realized the lidar system used in measurements. FG and OP wrote a main part of the article.
The authors declare that they have no conflict of interest.
This research work has been performed within the framework of a CIFRE grant (ANRT) for the doctoral work of Florian Gaudfrin. The lidar systems has been funded by the PROMETE project (ONERA). The laser source used in this paper for the experimental setup has been designed by Benoit Faure and PaulHenri Pioger, who is from the LEUKOS company
This research work has been supported within the framework of a CIFRE no 2016/0571 grant sign between LEUKOS company and French National Agency for Research and Technology (ANRT) for the doctoral work of Florian Gaudfrin. The lidar system has been funded by the PROMETE project (ONERA).
This paper was edited by Andrew Sayer and reviewed by two anonymous referees.
Ackermann, J.: The extinctiontobackscatter ratio of tropospheric aerosol: A numerical study, J. Atmos. Ocean. Tech., 15, 1043–1050, https://doi.org/10.1175/15200426(1998)015<1043:TETBRO>2.0.CO;2, 1998. a
Anon: U.S. Standard Atmosphere, Tech. rep., NASA, Washington, DC, US, 1976. a
Ansmann, A., Riebesell, M., and Weitkamp, C.: Measurement of atmospheric aerosol extinction profiles with a Raman lidar, Opt. Lett., 15, 746–748, https://doi.org/10.1364/OL.15.000746, 1990. a
Ansmann, A. and Müller, D.: Lidar and Atmospheric Aerosol Particles, in: RangeResolved Optical Remote Sensing of the Atmosphere, 105–141, SpringerVerlag, New York, 2004.
Ansmann, A., Wandinger, U., Riebesell, M., Weitkamp, C., and Michaelis, W.: Independent measurement of extinction and backscatter profiles in cirrus clouds by using a combined Raman elasticbackscatter lidar, Appl. Optics, 31, 7113, https://doi.org/10.1364/AO.31.007113, 1992. a
Ansmann, A., Neuber, R., Rairoux, P., and Wandinger, U.: Advances in Atmospheric Remote Sensing with Lidar, vol. 34, Springer Berlin Heidelberg, Berlin, Heidelberg, https://doi.org/10.1007/9783642606120, 1997. a
Ansmann, A., Baldasano, J., Calpini, B., Chaikovsky, Anatoly and Flamant, P., Mitev, V., Papayannis, A., Pelon, J., Resendes, D., Schneider, J., Trickl, T., and Vaughan, G.: EARLINET: A European Aerosol Research Lidar Network, Tech. rep., 2001. a
Baltensperger, U., Dommen, J., Alfarra, M. R., Duplissy, J., Gaeggeler, K., Metzger, A., Facchini, M. C., Decesari, S., Finessi, E., Reinnig, C., Schott, M., Warnke, J., Hoffmann, T., Klatzer, B., Puxbaum, H., Geiser, M., Savi, M., Lang, D., Kalberer, M., and Geiser, T.: Combined Determination of the Chemical Composition and of Health Effects of Secondary Organic Aerosols: The POLYSOA Project, J. Aerosol Med. Pulm. D., 21, 145–154, https://doi.org/10.1089/jamp.2007.0655, 2008. a
Berk, A., Acharya, P., and Bernstein, L.: Band model method for modeling atmospheric propagation at arbitrarily fine spectral resolution, US Pat., available at: https://patents.google.com/patent/US7433806 (last access: 6 July 2019), 2008. a
Berk, A., Conforti, P., Kennett, R., Perkins, T., Hawes, F., and van den Bosch, J.: MODTRAN6: a major upgrade of the MODTRAN radiative transfer code, 9088, 90880H, International Society for Optics and Photonics, https://doi.org/10.1117/12.2050433, 2014. a
Bissonnette, L. R.: Multiplescattering lidar equation, Appl. Optics, 35, 6449, https://doi.org/10.1364/AO.35.006449, 1996. a
Bodhaine, B. A., Wood, N. B., Dutton, E. G., and Slusser, J. R.: On Rayleigh Optical Depth Calculations, J. Atmos. Ocean. Tech., 16, 1854–1861, https://doi.org/10.1175/15200426(1999)016<1854:ORODC>2.0.CO;2, 1999. a
Bohlmann, S., Baars, H., Radenz, M., Engelmann, R., and Macke, A.: Shipborne aerosol profiling with lidar over the Atlantic Ocean: from pure marine conditions to complex dust–smoke mixtures, Atmos. Chem. Phys., 18, 9661–9679, https://doi.org/10.5194/acp1896612018, 2018. a
Bréon, F.M., Maignan, F., Leroy, M., and Grant, I.: Analysis of hot spot directional signatures measured from space, J. Geophys. Res., 107, 4282, https://doi.org/10.1029/2001JD001094, 2002. a
Bucholtz, A.: Rayleighscattering calculations for the terrestrial atmosphere, Appl. Optics, 34, 2765, https://doi.org/10.1364/AO.34.002765, 1995. a
Bufton, J. L.: Laser Altimetry Measurements from Aircraft and Spacecraft, P. IEEE, 77, 463–477, https://doi.org/10.1109/5.24131, 1989. a
Cattrall, C., Reagan, J., Kurt, T., and Oleg, D.: Variability of aerosol and spectral lidar and backscatter and extinction ratios of key aerosol types derived from selected Aerosol Robotic Network locations, J. Geophys. Res.Atmos., 110, 1–13, https://doi.org/10.1029/2004JD005124, 2005. a
Ceolato, R. and Gaudfrin, F.: Shortrange characterization of freshly emitted carbonaceous particles by LiDAR, in: General Meeting, Nafplio, 17–19 April 2018. a, b
Ceolato, R., Riviere, N., and Hespel, L.: Reflectances from a supercontinuum laserbased instrument: hyperspectral, polarimetric and angular measurements, Opt. Express, 20, 29413, https://doi.org/10.1364/oe.20.029413, 2012. a
Ceolato, R., Gaudfrin, F., Pujol, O., Riviere, N., Berg, M. J., and Sorensen, C. M.: Lidar crosssections of soot fractal aggregates: Assessment of equivalentsphere models, J. Quant. Spectrosc. Ra., 212, 39–44, https://doi.org/10.1016/j.jqsrt.2017.12.004, 2018. a
Clark, W. E. and Whitby, K. T.: Concentration and Size Distribution Measurements of Atmospheric Aerosols and a Test of the Theory of SelfPreserving Size Distributions, J. Atmos. Sci., 24, 677–687, https://doi.org/10.1175/15200469(1967)024<0677:CASDMO>2.0.CO;2, 1967. a
Collis, R. T. H. and Russell, P. B.: Lidar measurement of particles and gases by elastic backscattering and differential absorption, in: Laser Monitoring of the Atmosphere, 71–151, https://doi.org/10.1007/354007743X_18, 1976. a
Comeron, A., Sicard, M., Kumar, D., and Rocadenbosch, F.: Use of a field lens for improving the overlap function of a lidar system employing an optical fiber in the receiver assembly, Appl. Optics, 50, 5538, https://doi.org/10.1364/AO.50.005538, 2011. a
Dawson, K. W., Meskhidze, N., Josset, D., and Gassó, S.: Spaceborne observations of the lidar ratio of marine aerosols, Atmos. Chem. Phys., 15, 3241–3255, https://doi.org/10.5194/acp1532412015, 2015. a
DeMott, P. J., Prenni, A. J., Liu, X., Kreidenweis, S. M., Petters, M. D., Twohy, C. H., Richardson, M. S., Eidhammer, T., and Rogers, D. C.: Predicting global atmospheric ice nuclei distributions and their impacts on climate, P. Natl. Acad. Sci. USA, 107, 11217–11222, https://doi.org/10.1073/pnas.0910818107, 2010. a
Elias, T., Haeffelin, M., Drobinski, P., Gomes, L., Rangognio, J., Bergot, T., Chazette, P., Raut, J. C., and Colomb, M.: Particulate contribution to extinction of visible radiation: Pollution, haze, and fog, Atmos. Res., 92, 443–454, https://doi.org/10.1016/j.atmosres.2009.01.006, 2009. a
Fernald, F. G.: Analysis of atmospheric lidar observations: some comments, Appl. Optics, 23, 652–653, https://doi.org/10.1364/ao.23.000652, 1984. a
Fernald, F. G., Herman, B. M., and Reagan, J. A.: Determination of Aerosol Height Distributions by Lidar, J. Appl. Meteorol., 11, 482–489, https://doi.org/10.1175/15200450(1972)011<0482:DOAHDB>2.0.CO;2, 1972. a
FinlaysonPitts, B. J. and Pitts, J. N. J.: The Atmospheric System, Chemistry of the Upper and Lower Atmosphere: Theory, Experiments, and Applications, 15–42, https://doi.org/10.1016/B9780122570605/500046, 2000. a
Gaudfrin, F., Ceolato, R., Pujol, O., Riviere, N., and Hu, Q.: Supercontinuum lidar for hydrometeors and aerosols optical properties determination, in: EGU General Assembly 2018, 20, 14436, 2018a. a
Gaudfrin, F., Ceolato, R., Riviere, N., Pujol, O., and Huss, G.: New technics to probe aerosols radiative properties from visible to infraredTitle, in: The 13th International IR Target and Background Modeling & Simulation Workshop, BanyulssurMer, 2018b. a, b
Gaudfrin, F., Ceolato, R., Pujol, O., and Riviere, N.: A newlidar inversion technique using a surface reference target at short – range, in: The 29th International Laser Radar Conference, 3–6, 24–28 June, Hefei, 2019. a
Getzewich, B. J., Vaughan, M. A., Hunt, W. H., Avery, M. A., Powell, K. A., Tackett, J. L., Winker, D. M., Kar, J., Lee, K.P., and Toth, T. D.: CALIPSO lidar calibration at 532 nm: version 4 daytime algorithm, Atmos. Meas. Tech., 11, 6309–6326, https://doi.org/10.5194/amt1163092018, 2018. a
Glickman, T. and Zenk, W.: Glossary of Meteorology, American Meteorological Society, Boston, 855 pp. 2, available at: http://oceanrep.geomar.de/6557/ (last access: 17 July 2019), 2000. a
Hadamard, J.: Théorie des équations aux dérivées partielles linéaires hyperboliques et du problème de Cauchy, Acta Math., 31, 333–380, https://doi.org/10.1007/BF02415449, 1908. a
Hall, F. F. and Ageno, H. Y.: Absolute calibration of a laser system for atmospheric probing, Appl. Optics, 9, 1820–1824, https://doi.org/10.1364/AO.9.001820, 1970. a
Halldórsson, T. and Langerholc, J.: Geometrical form factors for the lidar function, Appl. Optics, 17, 240–244, https://doi.org/10.1364/AO.17.000240, 1978. a
Haner, D. A., McGuckin, B. T., Menzies, R. T., Bruegge, C. J., and Duval, V.: Directional–hemispherical reflectance for Spectralon by integration of its bidirectional reflectance, Appl. Optics, 37, 3996, https://doi.org/10.1364/AO.37.003996, 1998. a
Hansen, J., Sato, M., and Ruedy, R.: Radiative forcing and climate response, J. Geophys. Res.Atmos., 102, 6831–6864, https://doi.org/10.1029/96JD03436, 1997. a
Haywood, J. and Boucher, O.: Estimates of the direct and indirect radiative forcing due to tropospheric aerosols: A review, Rev. Geophys., 38, 513–543, https://doi.org/10.1029/1999RG000078, 2000. a
Hess, M., Koepke, P., Schult, I., Hess, M., Koepke, P., and Schult, I.: Optical Properties of Aerosols and Clouds: The Software Package OPAC, B. Am. Meteorol. Soc., 79, 831–844, https://doi.org/10.1175/15200477(1998)079<0831:OPOAAC>2.0.CO;2, 1998. a
Hoff, R. M., Bösenberg, J., and Pappalardo, G.: The GAW Aerosol Lidar Observation Network (GALION), in: International Geoscience and Remote Sensing Symposium (IGARSS08), 27–29 March, Boston (USA), 6–11, 2008. a
Horvath, H.: On the applicability of the Koschmieder visibility formula, Atmos. Environ., 5, 177–184, https://doi.org/10.1016/00046981(71)900813, 1971. a
Hu, Y., Vaughan, M., Liu, Z., Powell, K., and Rodier, S.: Retrieving Optical Depths and Lidar Ratios for Transparent Layers Above Opaque Water Clouds From CALIPSO Lidar Measurements, IEEE Geosci. Remote S., 4, 523–526, https://doi.org/10.1109/LGRS.2007.901085, 2007. a
Hyslop, N. P.: Impaired visibility: the air pollution people see, Atmos. Environ., 43, 182–195, https://doi.org/10.1016/j.atmosenv.2008.09.067, 2009. a
Jäger, H.: Longterm record of lidar observations of the stratospheric aerosol layer at GarmischPartenkirchen, J. Geophys. Res.Atmos., 110, 1–9, https://doi.org/10.1029/2004JD005506, 2005. a
Josset, D., Pelon, J., Protat, A., and Flamant, C.: New approach to determine aerosol optical depth from combined CALIPSO and CloudSat ocean surface echoes, Geophys. Res. Lett., 35, https://doi.org/10.1029/2008GL033442, 2008. a
Josset, D., Pelon, J., and Hu, Y.: MultiInstrument Calibration Method Based on a Multiwavelength Ocean Surface Model, IEEE Geosci. Remote S., 7, 195–199, https://doi.org/10.1109/LGRS.2009.2030906, 2010a. a
Josset, D., Zhai, P.W., Hu, Y., Pelon, J., and Lucker, P. L.: Lidar equation for ocean surface and subsurface, Opt. Express, 18, 20862, https://doi.org/10.1364/OE.18.020862, 2010b. a, b
Josset, D., Pelon, J., Garnier, A., Hu, Y., Vaughan, M., Zhai, P.W., Kuehn, R., and Lucker, P.: Cirrus optical depth and lidar ratio retrieval from combined CALIPSOCloudSat observations using ocean surface echo, J. Geophys. Res.Atmos., 117, https://doi.org/10.1029/2011jd016959, 2012. a
Josset, D., Pelon, J., Pascal, N., Hu, Y., and Hou, W.: On the Use of CALIPSO Land Surface Returns to Retrieve Aerosol and Cloud Optical Depths, IEEE T. Geosci. Remote, 56, 3256–3264, https://doi.org/10.1109/TGRS.2018.2796850, 2018. a, b
Kar, J., Vaughan, M. A., Lee, K.P., Tackett, J. L., Avery, M. A., Garnier, A., Getzewich, B. J., Hunt, W. H., Josset, D., Liu, Z., Lucker, P. L., Magill, B., Omar, A. H., Pelon, J., Rogers, R. R., Toth, T. D., Trepte, C. R., Vernier, J.P., Winker, D. M., and Young, S. A.: CALIPSO lidar calibration at 532 nm: version 4 nighttime algorithm, Atmos. Meas. Tech., 11, 1459–1479, https://doi.org/10.5194/amt1114592018, 2018. a
Kavaya, M. J., Menzies, R. T., Haner, D. A., Oppenheim, U. P., and Flamant, P. H.: Target reflectance measurements for calibration of lidar atmospheric backscatter data, Appl. Optics, 22, 2619, https://doi.org/10.1364/AO.22.002619, 1983. a
Klett, J. D.: Stable analytical inversion solution for processing lidar returns, Appl. Optics, 20, 211–220, https://doi.org/10.1364/AO.20.000211, 1981. a, b
Klett, J. D.: Lidar inversion with variable backscatter/extinction ratios, Appl. Optics, 24, 1638, https://doi.org/10.1364/AO.24.001638, 1985. a, b
Kraft, D.: Software Package for Sequential Quadratic Programming, 1988. a
Kumar, D. and Rocadenbosch, F.: Determination of the overlap factor and its enhancement for mediumsize tropospheric lidar systems: a raytracing approach, J. Appl. Remote Sens., 7, 073591, https://doi.org/10.1117/1.JRS.7.073591, 2013. a
Labarre, L., Caillault, K., Fauqueux, S., Malherbe, C., Roblin, A., Rosier, B., and Simoneau, P.: An overview of MATISSEv2.0, in: Optics in Atmospheric Propagation and Adaptive Systems XIII, 7828, edited by: Stein, K. and Gonglewski, J. D., 9–18, https://doi.org/10.1117/12.868183, 2010. a
Leblanc, T., Trickl, T., and Vogelmann, H.: Lidar, vol. 102 of Springer Series in Optical Sciences, SpringerVerlag, New York, https://doi.org/10.1007/b106786, 2005. a
Li, J., Gong, W., and Ma, Y.: Atmospheric Lidar Noise Reduction Based on Ensemble Empirical Mode Decomposition, Int. Arch. Photogramm. Remote Sens. Spatial Inf. Sci., XXXIXB8, 127–129, https://doi.org/10.5194/isprsarchivesXXXIXB81272012, 2012. a
Lobell, D. B. and Asner, G. P.: Moisture Effects on Soil Reflectance, Soil Sci. Soc. Am. J., 66, 722–727, https://doi.org/10.2136/sssaj2002.7220, 2002. a
Mao, F., Gong, W., and Li, C.: Antinoise algorithm of lidar data retrieval by combining the ensemble Kalman filter and the Fernald method, Opt. Express, 21, 8286, https://doi.org/10.1364/OE.21.008286, 2013. a
Matsumoto, M. and Takeuchi, N.: Effects of misestimated farend boundary values on two common lidar inversion solutions, Appl. Optics, 33, 6451, https://doi.org/10.1364/AO.33.006451, 1994. a
Mattis, I., Ansmann, A., Müller, D., Wandinger, U., and Althausen, D.: Multilayer aerosol observations with dualwavelength Raman lidar in the framework of EARLINET, J. Geophys. Res.Atmos., 109, 1–15, https://doi.org/10.1029/2004JD004600, 2004. a
Measures, R. M.: Laser Remote Sensing: Fundamentals and Applications, reprint ed edn., edited by: Wiley, J., Krieger publishing compagny, Malabar, Florida, 1992. a
Miesch, C., Poutier, L., Achard, V., Briottet, X., Lenot, X., and Boucher, Y.: Direct and inverse radiative transfer solutions for visible and nearinfrared hyperspectral imagery, IEEE T. Geosci. Remote, 43, 1552–1562, https://doi.org/10.1109/TGRS.2005.847793, 2005. a
Mishchenko, M. I., Dlugach, J. M., Yanovitskij, E. G., and Zakharova, N. T.: Bidirectional reflectance of flat, optically thick particulate layers: an efficient radiative transfer solution and applications to snow and soil surfaces, J. Quant. Spectrosc. Ra., 63, 409–432, https://doi.org/10.1016/S00224073(99)00028X, 1999. a
Nicodemus, F. E.: Directional Reflectance and Emissivity of an Opaque Surface, Appl. Optics, 4, 767–775, https://doi.org/10.1364/AO.4.000767, 1965. a
O'Connor, E. J., Illingworth, A. J., and Hogan, R. J.: A Technique for Autocalibration of Cloud Lidar, J. Atmos. Ocean. Tech., 21, 777–786, 2004. a
Painemal, D., Clayton, M., Ferrare, R., Burton, S., Josset, D., and Vaughan, M.: Novel aerosol extinction coefficients and lidar ratios over the ocean from CALIPSO–CloudSat: evaluation and global statistics, Atmos. Meas. Tech., 12, 2201–2217, https://doi.org/10.5194/amt1222012019, 2019. a
Paschotta, R.: Field Guide to Laser Pulse Generation, vol. FG14, SPIE., https://doi.org/10.1111/j.17496632.1965.tb20241.x, 2008. a
Pedrós, R., Estellés, V., Sicard, M., GómezAmo, J. L., Utrillas, M. P., MartínezLozano, J. A., Rocadenbosch, F., Pérez, C., and Recio Baldasano, J. M.: Climatology of the aerosol extinctiontobackscatter ratio from sunphotometric measurements, IEEE T. Geosci. Remote, 48, 237–249, https://doi.org/10.1109/TGRS.2009.2027699, 2010. a
Popovicheva, O. B., Engling, G., Ku, I.T., Timofeev, M. A., and Shonija, N. K.: Aerosol Emissions from LongLasting Smoldering of Boreal Peatlands: Chemical Composition, Markers, and Microstructure, Aerosol Air Qual. Res., 2, 213–218, https://doi.org/10.4209/aaqr.2018.08.0302, 2019. a
Poutier, L., Miesch, C., Lenot, X., Achard, V., and Boucher, Y.: COMANCHE and COCHISE: two reciprocal atmospheric codes for hyperspectral remote sensing, Tech. rep., available at: http://wwwe.onecert.fr/pirrene/references_docs/poutier_et_al_AVIRIS_2002.pdf (last access: 4 April 2020), 2002. a
Rocadenbosch, F., Frasier, S., Kumar, D., Lange Vega, D., Gregorio, E., and Sicard, M.: Backscatter error bounds for the elastic lidar twocomponent inversion algorithm, IEEE T. Geosci. Remote, 50, 4791–4803, https://doi.org/10.1109/TGRS.2012.2194501, 2012. a
Sasano, Y., Browell, E. V., and Ismail, S.: Error caused by using a constant extinction/backscattering ratio in the lidar solution, Appl. Optics, 24, 3929, https://doi.org/10.1364/AO.24.003929, 1985. a
Sicard, M., Chazette, P., Pelon, J., Won, J., and Yoon, S.C.: Variational method for the retrieval of the optical thickness and the backscatter coefficient from multiangle lidar profiles, Appl. Optics, 41, 493–502, https://doi.org/10.1364/AO.41.000493, 2002. a, b
Simoneau, P., Berton, R., Caillault, K., Durand, G., Huet, T., Labarre, L., Malherbe, C., Miesch, C., Roblin, A., and Rosier, B.: MATISSE: advanced Earth modeling for imaging and scene simulation, in: Proceedings of SPIE, edited by: Kohnle, A., Gonglewski, J. D., and Schmugge, T. J., 4538, 39, https://doi.org/10.1117/12.454410, 2002. a
SphereOptic: Diffusers, Zenith Polymer, Foil, nearly ideal Lambertian SphereOptics EN, available at: https://sphereoptics.de/en/product/zenithpolymerdiffusers/?c=79, last access: 4 April 2020. a
Stelmaszczyk, K., Dell'Aglio, M., Chudzyński, S., Stacewicz, T., and Wöste, L.: Analytical function for lidar geometrical compression formfactor calculations, Appl. Optics, 44, 1323, https://doi.org/10.1364/AO.44.001323, 2005. a
Sun, X.: Lidar Sensors From Space, Elsevier, https://doi.org/10.1016/B9780124095489.103276, 2018. a
Swinehart, D.: The BeerLambert law, J. Chem. Educ., 39, 335–335, https://doi.org/10.1021/ed039p333, 1962. a
Vande Hey, J., Coupland, J., Foo, M. H., Richards, J., and Sandford, A.: Determination of overlap in lidar systems, Appl. Optics, 50, 5791, https://doi.org/10.1364/AO.50.005791, 2011. a
Vande Hey, J. D.: Theory of Lidar, Springer International Publishing, 23–41, https://doi.org/10.1007/9783319126135_2, 2014. a, b, c
VargasUbera, J., Aguilar, J. F., and Gale, D. M.: Reconstruction of particlesize distributions from lightscattering patterns using three inversion methods, Appl. Optics, 46, 124–132, https://doi.org/10.1364/AO.46.000124, 2007. a
Wandinger, U. and Ansmann, A.: Experimental determination of the lidar overlap profile with Raman lidar, Appl. Optics, 41, 511–514, https://doi.org/10.1364/AO.41.000511, 2002. a
Young, S. A.: Analysis of lidar backscatter profiles in optically thin clouds, Appl. Opt., 34, 7019–7031, https://doi.org/10.1364/AO.34.007019, 1995. a
Zhang, J., Wei, E., Wu, L., Fang, X., Li, F., Yang, Z., Wang, T., and Mao, H.: Elemental Composition and Health Risk Assessment of PM_{10} and PM_{2.5} in the Roadside Microenvironment in Tianjin, China, Aerosol Air Qual. Res., 18, 1817–1827, https://doi.org/10.4209/aaqr.2017.10.0383, 2018. a
 Abstract
 Introduction
 Unified lidar equation for surface and volumetric scattering media
 New lidar inversion technique
 Theoretical behavior of the retrieval procedure
 Case of real measurements
 Conclusions
 Appendix A
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Unified lidar equation for surface and volumetric scattering media
 New lidar inversion technique
 Theoretical behavior of the retrieval procedure
 Case of real measurements
 Conclusions
 Appendix A
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References