Journal cover
Journal topic
**Atmospheric Measurement Techniques**
An interactive open-access journal of the European Geosciences Union

Journal topic

- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Subscribe to alerts
- Peer review
- For authors
- For reviewers
- EGU publications
- Imprint
- Data protection

- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Subscribe to alerts
- Peer review
- For authors
- For reviewers
- EGU publications
- Imprint
- Data protection

**Research article**
12 Mar 2019

**Research article** | 12 Mar 2019

An algorithm to retrieve ice water content profiles

^{1}Laboratoire d'Optique Atmosphérique, Université de Lille, Villeneuve-d'Ascq, France^{2}Cimel Electronique, Paris, France^{a}now at: Laboratoire de Météorologie Dynamique/Institut Pierre-Simon Laplace (LMD/IPSL), Sorbonne Université, Ecole Polytechnique, CNRS, Paris, France

^{1}Laboratoire d'Optique Atmosphérique, Université de Lille, Villeneuve-d'Ascq, France^{2}Cimel Electronique, Paris, France^{a}now at: Laboratoire de Météorologie Dynamique/Institut Pierre-Simon Laplace (LMD/IPSL), Sorbonne Université, Ecole Polytechnique, CNRS, Paris, France

Abstract

Back to toptop
The algorithm presented in this paper was developed to retrieve ice water content (IWC) profiles in cirrus clouds. It is based on optimal estimation theory and combines ground-based visible lidar and thermal infrared (TIR) radiometer measurements in a common retrieval framework in order to retrieve profiles of IWC together with a correction factor for the backscatter intensity of cirrus cloud particles. As a first step, we introduce a method to retrieve extinction and IWC profiles in cirrus clouds from the lidar measurements alone and demonstrate the shortcomings of this approach due to the backscatter-to-extinction ambiguity. As a second step, we show that TIR radiances constrain the backscattering of the ice crystals at the visible lidar wavelength by constraining the ice water path (IWP) and hence the IWC, which is linked to the optical properties of the ice crystals via a realistic bulk ice microphysical model. The scattering phase function obtained from the microphysical model is flat around the backscatter direction (i.e., there is no backscatter peak). We show that using this flat backscattering phase function to define the backscatter-to-extinction ratio of the ice crystals in the retrievals with the lidar-only algorithm results in an overestimation of the IWC, which is inconsistent with the TIR radiometer measurements. Hence, a synergy algorithm was developed that combines the attenuated backscatter profiles measured by the lidar and the measurements of TIR radiances in a common optimal estimation framework to retrieve the IWC profile together with a correction factor for the phase function of the bulk ice crystals in the backscattering direction. We show that this approach yields consistent lidar and TIR results. The resulting lidar ratios for cirrus clouds are found to be consistent with previous independent studies.

Download & links

How to cite

Back to top
top
How to cite.

Hemmer, F., C.-Labonnote, L., Parol, F., Brogniez, G., Damiri, B., and Podvin, T.: An algorithm to retrieve ice water content profiles in cirrus clouds from the synergy of ground-based lidar and thermal infrared radiometer measurements, Atmos. Meas. Tech., 12, 1545-1568, https://doi.org/10.5194/amt-12-1545-2019, 2019.

1 Introduction

Back to toptop
The importance of clouds for the climate system has been extensively discussed during the last decades (Stephens, 2005). Although their essential role in the Earth's radiative budget is unquestionable, they still remain a major source of uncertainty in climate change estimates (Boucher et al., 2013). In particular, the important but complex impact of cirrus clouds has long been recognized (Liou, 1986) but is still not well quantified. This is due to the large range of varying shapes and sizes of the ice crystals observed in cirrus clouds which may interact in different ways with atmospheric radiation by scattering and absorption processes. The net radiative effect of cirrus clouds is generally positive but can be negative as well (Zhang et al., 1999). It is determined on the one hand by the macrophysical cloud properties, e.g., altitude, geometrical thickness, temperature, and the difference between the temperature of the cloud and the surface (e.g., Stephens and Webster, 1981). On the other hand, it depends on the optical properties of the cloud, which are in turn governed by the microphysics, especially the size, shape and number density of particles.

Recent advances in satellite observational systems, particularly the Cloud Profiling Radar (CPR; Stephens et al., 2002, 2008) aboard CloudSat and the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP; Winker et al., 2009, 2010) aboard CALIPSO as part of the international satellite constellation known as A-Train, have shown that the occurrence of cirrus clouds in the atmosphere is much higher than previously presumed (Stubenrauch et al., 2013). Mace et al. (2009) quantified the global occurrence frequencies to 40 %–60 % whereas earlier estimates expected about 20 %–30 % with a higher coverage of 60 %–70 % in the tropics (e.g., Liou, 1986; Wylie et al., 1994). This underlines the importance of studying the characteristics of cirrus clouds to estimate their influence on the radiation budget.

Lidar systems have proven to be powerful tools to study even the most tenuous cloud layers (e.g., Sassen, 1991). The measured backscatter profiles provide information about the cloud-base and cloud-top altitudes, which can be related to temperature by using atmospheric temperature profiles from model reanalysis or radiosounding. Furthermore, these measurements yield the possibility to retrieve profiles of particle extinction and hence the optical depth of the cirrus cloud by making assumptions about the so-called backscatter-to-extinction ratio. However, there are more advanced lidar systems which do not require such assumptions that can introduce large errors in the estimated cloud optical depth. Raman lidars, for example, can provide particle extinction directly since the inelastic Raman backscatter signal is only sensitive to extinction but not to particle backscattering (Ansmann et al., 1990, 1992). High-spectral-resolution lidars are also capable of measuring particle extinction directly with the help of two channels: one that measures the backscatter originating from the entire atmosphere (molecular plus particle) and one that measures only the molecular contribution by removing the central portion of the signal that is associated with aerosol or cloud particles with a filter. From these two simultaneous measurements the particle extinction can be derived from the change in the slope of the molecular signal relative to a clear-sky atmosphere (Turner and Eloranta, 2008). Other lidars, e.g., CALIOP, include polarization measurements from which the cloud phase can be determined since ice crystals tend to depolarize the incident visible radiation whereas for water droplets no such depolarization is observed (Sassen, 1991). However, most ground-based lidars operated at present are simpler systems (Campbell et al., 2015). Thus, the algorithm presented here was developed for the exploitation of data from a simple micro-pulse lidar (combined with thermal infrared, TIR, radiance measurements) although it might be applied to other lidars in future studies. Our method should be applicable to combined TIR and simple backscatter lidar measurements from ground-based as well as space-based observations. Concerning cirrus clouds, there are already a large number of lidar studies dealing with their occurrence frequencies and characteristics, based on satellite data (e.g., Berthier et al., 2008; Campbell et al., 2015), ground-based data (e.g., Ansmann et al., 1993; Keckhut et al., 2006; Giannakaki et al., 2007; Seifert et al., 2007; Liu et al., 2015) or a combination of both (e.g., Pandit et al., 2015; Córdoba-Jabonero et al., 2017). Recently, Campbell et al. (2016) went even further and characterized the daytime radiative forcing of cirrus clouds at the top of the atmosphere from ground-based lidar measurements at a midlatitude site and thereby underlined the importance of ground-based measurements for the estimation of the radiative effect of cirrus.

In spite of their undenied importance, lidar systems are not the only tools to study cirrus clouds. Another important source of information are measurements from passive TIR radiometers, which are also performed from the ground or from space, e.g., the Imaging Infrared Radiometer (IIR) aboard CALIPSO (Winker et al., 2003) or the MODerate resolution Imaging Spectroradiometer (MODIS) aboard Aqua and Terra (King et al., 1992, 2003). These measurements are sensitive to the optical and integrated properties of the cloud, for example the ice water path (IWP). A well-known method using radiances in the TIR wavelength region is the split window technique (Inoue, 1985, 1987; Parol et al., 1991), which allows the retrieval of the cloud-top temperature and the effective emissivity of semitransparent cirrus clouds from two channels centered around 11 and 12 µm, respectively. The method is based on the fact that the brightness temperature difference (BTD) of these channels is always more important for thin cirrus clouds than for thick clouds or under clear-sky conditions. In addition, the BTD is sensitive to the radiative and microphysical properties of the cloud. Dubuisson et al. (2008) showed, by conducting radiative transfer calculations with different ice crystal models, that it is possible to retrieve microphysical properties of cirrus clouds from passive TIR radiometer measurements alone.

However, in recent years synergistic approaches using independent sets of measurements in a common retrieval framework have become more and more popular. Examples that could be cited here are the raDAR/liDAR (DARDAR) algorithm to retrieve ice cloud properties from the synergy of the CPR and CALIOP measurements (Delanoë and Hogan, 2008, 2010) or the multilayer algorithm to retrieve ice and liquid water cloud properties simultaneously from three TIR radiances measured by the IIR and two MODIS reflectances measured at 0.85 and 2.13 µm (Sourdeval et al., 2015, 2016). There are also a few methods combining lidar and TIR radiometer measurements that have been developed in the past. The lidar and infrared radiometric (LIRAD) method introduced by Platt (1973, 1979) was the first method that combined lidar and infrared radiometer data to retrieve optical properties of cirrus clouds. It has been applied and further developed in several studies (e.g., Platt et al., 1987, 2002; Comstock and Sassen, 2001). In this approach, the lidar backscatter coefficient is related theoretically to the infrared volume absorption coefficient. The emissivity of the cloud is then derived in an iterative process by calculating a theoretical cloud radiance, which is compared to the infrared radiometer measurement and adjusting the backscatter-to-extinction ratio until the theoretical and measured radiances converge. Other studies focused on the combination of lidar measurements and the split-window technique to improve the retrievals of cloud properties from passive sensors alone by integrating the information provided by the active lidar measurements in the radiative transfer calculations. Chiriaco et al. (2004) showed the theoretical potential of this approach to improve the particle size retrieval from the instruments aboard CALIPSO, i.e., the IIR and CALIOP, and Garnier et al. (2012, 2013) developed an algorithm based on this idea to retrieve the effective emissivity, optical depth, effective diameter and IWP from CALIPSO measurements. Nevertheless, in their approach the retrieval is based on the split-window technique in which information from the lidar such as scene identification and cloud altitude have been integrated. Saito et al. (2017) recently demonstrated a method to simultaneously infer the IWP, the cloud effective radius, the surface temperature and two morphological parameters, namely the fraction of plates and the surface roughness of ice crystal aggregates, from a synergistic approach based on optimal estimation. They used the layer-integrated total attenuated backscatter and the depolarization ratio at 532 nm from CALIOP as well as the brightness temperatures at 8.65, 10.6 and 12.0 µm from the IIR in a common retrieval framework to obtain the parameters cited above.

The algorithm proposed in this paper also establishes a synergy between lidar and TIR radiometer measurements, although our lidar is a simple micro-pulse lidar and does not possess depolarization channels. In contrast to Saito et al. (2017), we use the whole backscattering profile measured by the lidar together with two TIR radiances in the measurement vector to retrieve profiles of particle extinction and ice water content (IWC). This allows us to include the profile information from the active lidar measurements in the radiative transfer calculations in the TIR, which is an improvement since common retrieval algorithms often assume plane-parallel and homogeneous conditions. As a first step, we developed an algorithm to retrieve extinction and IWC profiles in thin cirrus clouds from ground-based lidar measurements alone. This algorithm is based on the method of Stephens et al. (2001), who used an optimal estimation approach to invert the lidar equation to retrieve profiles of particle extinction from spaceborne data collected during the Lidar in Space Technology Experiment (LITE, McCormick et al., 1993). To overcome the backscatter-to-extinction ambiguity arising from the combination of scattering and absorption processes when regarding the lidar measurements alone, Stephens et al. (2001) introduced an optical depth constraint in the form of an additional measurement. In contrast to this approach, we developed, in a second step, a synergy algorithm that integrates actual measurements of TIR radiances in the optimal estimation framework. We will show that these radiances constrain the backscattering of the ice crystals at the visible lidar wavelength by constraining the IWP and hence the IWC, which is linked to the optical properties of the ice crystals via the microphysical ice cloud model of Baran et al. (2001, 2014a, b) and Vidot et al. (2015).

The paper is organized as follows: Sect. 2 briefly introduces the instruments and data used in this study. Section 3 presents our approach to the lidar retrieval problem and describes the algorithm for the retrieval of extinction and IWC profiles from lidar measurements as well as the underlying microphysical model for cirrus clouds. In this section, we also discuss the abovementioned backscatter-to-extinction ambiguity before the ability of TIR radiances to constrain the backscatter-to-extinction ratio is outlined. Section 4 presents the new algorithm using the synergy of lidar and TIR radiances, which has been developed on the basis of the lidar-only algorithm. Finally, Sect. 5 concludes this study.

2 Instrumentation and data

Back to toptop
The data used in this study originate from the measurement platform of the Laboratoire d'Optique Atmosphérique (LOA) situated on the campus of the University of Lille in northern France. This platform is equipped with, amongst other instruments, an elastic-backscatter micro-pulse lidar and a TIR radiometer.

The lidar is a Cloud and Aerosol Microlidar (CAML) CE370 (Pelon et al., 2008)
developed by the company CIMEL Electronique. It is an eye-safe lidar system
which operates at a single wavelength of 532 nm and does not include
depolarization. The system is automated and has been operated continuously
since 2007; hence a large archive of data is available for the LOA
measurement site. The type of laser integrated in the instrument is a
frequency-doubled Nd:YAG laser. The divergence of the laser beam as well as
the field of view (FOV) of the receiver are both 55 µrad. The
pulse duration is 100 ns and the repetition rate is 4.7 kHz, which results in a
vertical resolution of 15 m defined by $\mathrm{100}\phantom{\rule{0.125em}{0ex}}\mathrm{ns}\cdot c/\mathrm{2}$, where *c*
is the speed of light. The lidar profiles used in this study are averaged
over 1 min and a vertical binomial filter has been applied in order to
smooth the signal. The lidar pointed directly vertical with a zenith
angle of 0^{∘}.

The radiometer is called Conveyable Low-Noise Infrared Radiometer for
Measurements of Atmosphere and Ground Surface Targets (CLIMAT)
(Sicard et al., 1999; Legrand et al., 2000; Brogniez et al., 2003). It was developed to measure
radiances in the TIR wavelength region in three different spectral bands
centered at 8.7, 10.8 and 12.0 µm. The
full width at half maximum (FWHM) of each of these channels is 1 µm. In the following, we will call them C09, C11 and C12, respectively. There
are two different versions of the instrument: one was designed for
ground-based measurements and the other one for airborne measurements.
Although we exploit ground-based observations in this study, the CLIMAT
instrument currently installed on the LOA measurement platform is of the aircraft
type. This instrument measures the radiances of the three different channels
simultaneously and has a FOV of 3.5^{∘}. It consists of two main
parts: the optical head containing the optical elements as well as the
detector and the control unit containing the electronics and the memory. The
main optic consists of two germanium lenses: the objective, which is a
standard convex-plane lens, and the condenser, which is a “best-shaped”
meniscus designed to minimize the geometrical aberrations
(Legrand et al., 2000). The condenser is situated in the focal plane of the
objective. The optical head is constructed respecting the so-called Köhler
design, which means that the detector is located in the conjugate plane of the
objective with respect to the condenser. The radiation is measured by a
thermopile of which the hot junction is heated by the incident radiation and the
temperature of the cold junction is determined by the ambient temperature of
the cavity. In contrast to the lidar system, the CLIMAT instrument is
operated manually depending on the weather conditions.

It should be noted that due to the larger FOV of the TIR radiometer compared to the FOV of the lidar, the two instruments do not see exactly the same cloud area. This difference also depends on the altitude of the cloud. As in almost all remote-sensing algorithms, we assumed a homogeneous cloud in the instrument FOV and did not take into account any uncertainty due to sub-pixel heterogeneity.

3 Lidar-only algorithm

Back to toptop
The algorithm presented here was developed to retrieve profiles of particle extinction and IWC from measured lidar backscattering profiles. We propose a method to simultaneously retrieve a profile of aerosol extinction in the layers close to the ground and a profile of IWC inside cirrus layers. However, our main focus is the characterization of the microphysical properties of cirrus clouds. There are many techniques to invert the lidar equation including the classical Klett–Fernald method (Klett, 1981, 1985; Fernald, 1984). Our algorithm closely follows the method described by Stephens et al. (2001), which is based on optimal estimation theory introduced by Rodgers (1976, 1990, 2000), which is now a common approach for the inversion of remote-sensing data. One advantage of this approach is that it directly provides an estimation of the uncertainties together with the retrieved quantities. Furthermore, it facilitates the introduction of additional information in a common retrieval framework in order to constrain the retrieved parameters. The additional information could be, for example, measurements of polarization or measurements at other lidar wavelengths, as well as TIR radiometer measurements, as will be discussed in Sect. 4 of this paper. As a first step, we focus in this section on the retrieval of extinction and IWC profiles from the lidar measurements alone.

The relationship between the range-resolved backscattered power, *P*(*r*), and
the atmospheric scattering and attenuation properties is described by the
lidar equation and may be expressed as follows:

$$\begin{array}{}\text{(1)}& C\cdot P\left(r\right){r}^{\mathrm{2}}=\left({\mathit{\beta}}_{\mathrm{m}}\right(r)+{\mathit{\beta}}_{\mathrm{p}}(r\left)\right)\cdot \mathrm{exp}[-\mathrm{2}\underset{\mathrm{0}}{\overset{r}{\int}}({\mathit{\sigma}}_{\mathrm{m}}\left({r}^{\prime}\right)+\mathit{\eta}{\mathit{\sigma}}_{\mathrm{p}}\left({r}^{\prime}\right)\left)\mathrm{d}{r}^{\prime}\right],\end{array}$$

where *C* is a calibration constant depending on the lidar system and the
atmospheric profile. *β*(*r*) and *σ*(*r*) represent the
backscattering and extinction coefficients, respectively, and both contain a
contribution arising from purely molecular backscattering (*β*_{m}(*r*)) or
extinction (*σ*_{m}(*r*)) and a contribution arising from cloud or aerosol
particles that may be present in the atmosphere (*β*_{p}(*r*) and
*σ*_{p}(*r*), respectively). The factor *η* accounts for multiple-scattering processes. For the remainder of this paper we assume that *η*=1
for aerosols and *η*=0.75 for cirrus clouds. The value of 0.75 for cirrus
clouds has been chosen based on the PhD thesis of Nohra (2016) in which the
multiple-scattering factor has been evaluated by comparing the optical depth
retrieved from measurements of the micro-pulse lidar in Lille to the optical
depth retrieved from CALIOP and adjusting *η* to find a coherent
retrieval. The multiple-scattering factor used for the CALIOP version 3
retrievals is *η*=0.6 (Garnier et al., 2015). For ground-based lidars the
multiple-scattering effect is less important since they have a much smaller
FOV in combination with a shorter distance to the cloud, although it should
not be neglected because large ice crystals may considerably increase the
forward scattering of the laser beam (Donovan and van Lammeren, 2001). However, since our
knowledge of this parameter is rather poor we assign a large error to it in
our optimal estimation algorithm (see Sect. 3.2).

The retrieval of particle optical properties from elastic lidar measurements
alone is challenging since there is an intrinsic ambiguity between the
effects of backscattering and extinction arising from the combination of
scattering and absorption processes in the atmosphere. In Eq. (1), the backscattering coefficient of aerosol or cloud
particles, *β*_{p}(*r*), can be replaced by

$$\begin{array}{}\text{(2)}& {\mathit{\beta}}_{\mathrm{p}}\left(r\right)=k\left(r\right)\cdot {\mathit{\sigma}}_{\mathrm{p}}\left(r\right),\end{array}$$

where *k*(*r*) represents the range-dependent backscatter-to-extinction ratio.
Since we use a simple micro-pulse lidar, we need to introduce some
assumptions for *k*(*r*) to retrieve profiles of extinction by aerosol and
cloud particles. Unfortunately, this parameter is highly variable and depends
strongly on the type, size and shape of the atmospheric particles.

In this study we are focusing on the retrieval of cirrus cloud properties. Thus, the backscatter-to-extinction coefficient for aerosols is assumed to be constant and is fixed to 64 sr. This value originates from the Optical Properties of Aerosols and Clouds (OPAC) database (Hess et al., 1998b) and corresponds to a water-soluble urban aerosol. Since our measurement site is located in an urban/industrial area, we used the optical properties for this aerosol type (for the visible lidar wavelength as well as for the TIR as will be discussed in Sect. 4.1) from the OPAC database to test our new algorithm. This parameter needs to be refined in future studies depending on the aerosol type that is actually present during the measurement obtained from additional information. However, as discussed above for the multiple-scattering factor, we also assign a large uncertainty to the lidar ratio of aerosols to account for our rather poor knowledge of it. Unfortunately, this reduces the quality of the information returned by our algorithm.

The backscatter-to-extinction ratio for cirrus clouds is calculated using the definition of Mishchenko et al. (1997):

$$\begin{array}{}\text{(3)}& k={\mathit{\varpi}}_{\mathrm{0}}\cdot {P}_{\mathrm{11}}\left(\mathit{\pi}\right),\end{array}$$

where *ϖ*_{0} is the particle single-scattering albedo and *P*_{11}(*π*)
the phase function in the exact backscattering direction.

We obtain the single-scattering properties (scattering coefficient, absorption coefficient and asymmetry parameter) for each cloud layer from the parametrization of Vidot et al. (2015), which is based on the ensemble model for cirrus introduced by Baran and Labonnote (2007). The idea of this model is to represent the variability of ice crystal sizes and shapes inside a cirrus cloud by assuming a distribution of some idealized shapes rather than assuming just a single geometrical form throughout the whole size spectrum. Observed ice crystal shapes in cirrus clouds range from simple pristine particles such as hexagonal ice columns and bullet rosettes (associated with small particles) over aggregates of these particles to aggregate chains, while the complexity of the crystal tends to increase with increasing size. The ensemble model of Baran and Labonnote (2007) attempts to reproduce these observations by using the six members shown in Fig. 1. The smallest ice crystals are represented by the first two members, which are simple hexagonal ice columns and bullet rosettes. The following members represent larger and more complex ice crystals by arbitrarily attaching up to 10 hexagonal elements to create chain-like structures. For the calculation of the bulk optical properties, the particle size distribution (PSD) of Field et al. (2005, 2007) is assumed, which is independent of assumptions about the ice crystal shape and depends only on the in-cloud temperature and the IWC. This parametrization has been constructed based on a large number of in situ measured PSDs and does not include measurements of ice crystal sizes less than 100 µm due to the shattering problem (Strapp et al., 2001; Field et al., 2003). For particles smaller than 100 µm an exponential PSD is assumed. Baran et al. (2011, 2014a, b) used the Field et al. (2005, 2007) parametrization to calculate the single-scattering properties for the ensemble model as functions of IWC and in-cloud temperature for a total of 20 662 in situ measurements from different aircraft-based field campaigns located in the tropics and in the midlatitudes. They created a database of optical ice cloud properties comprising 145 wavelengths between 0.2 and 120 µm. This database was used by Vidot et al. (2015) to develop a new ice cloud parametrization that predicts the single-scattering properties named above as functions of the in-cloud temperature and IWC without the need of a priori information on the shape and the effective diameter of the ice crystals. For the remainder of this paper, we will call this microphysical model BV2015.

Since our algorithm seeks to retrieve the IWC for cirrus cloud layers, the
extinction required in the lidar equation (Eq. 1) is
calculated from the scattering and absorption coefficients obtained from the
BV2015 parametrization as a function of the IWC of each cloud layer. The
necessary temperature information is obtained from the European Centre for
Medium-Range Weather Forecasts (ECMWF) reanalysis by matching atmospheric
temperature profiles to the corresponding cirrus cloud altitude. The single-scattering albedo required in Eq. (3) is also calculated
from the scattering and absorption coefficients, and the scattering phase
function is generated from the asymmetry parameter using the analytic phase
function of Baran et al. (2001) which is a linear piecewise parametrization of
the Henyey–Greenstein phase function depending only on the asymmetry
parameter. It is kept smooth and featureless since atmospheric ice crystals
may be distorted, be roughened or contain inclusions of air bubbles or aerosols.
All these processes would remove or reduce the optical features of the phase
function like the halos at 22 or 46^{∘} and the
backscattering peak (e.g., Macke et al., 1996; C.-Labonnote et al., 2001). Baran et al. (2001) demonstrated that their parametrization
reproduces short-wave multi-angle satellite and aircraft observations, and
Baran and Francis (2004) showed a good agreement with high-resolution infrared
observations between 3 and 18 µm. It has also been shown to be in
good agreement with the backscattering features observed from POLarization
and Directionality of the Earth's Reflectances (POLDER) measurements
(Baran and Labonnote, 2007). The scattering phase function, especially in the exact
backscattering direction, is a crucial parameter in our algorithm since it
defines the lidar backscatter-to-extinction ratio. It will be discussed in
more detail in Sect. 3.3.1.

As mentioned above, we apply an optimal estimation method to invert the lidar
equation following Stephens et al. (2001). Optimal estimation is based on a
Bayesian approach which uses probability density functions to link the
measurement space to the state space accounting for their uncertainties
(Rodgers, 2000). This approach allows us to find the most likely solution
that is consistent with both the measurement and any given prior knowledge
of the state within the range of their uncertainties. In general, the
measurement vector ** y** can be related to the state vector

$$\begin{array}{}\text{(4)}& \mathit{y}=\mathbf{F}\left(\mathit{x}\right)+\mathit{\u03f5},\end{array}$$

where ** ϵ** represents the uncertainties arising from the
measurements and the forward model. The aim of every inversion method is to
invert the connection between the state vector and the measurement vector,
which is given by the forward model, in order to retrieve the elements of the
state vector using the information provided by the measurement vector.

Following Rodgers (2000), the best estimation of the state vector can be obtained by minimizing the following cost function:

$$\begin{array}{}\text{(5)}& \mathrm{\Phi}=[\mathit{y}-\mathbf{F}(\mathit{x}\left){]}^{T}{\mathbf{S}}_{\mathit{\u03f5}}^{-\mathrm{1}}\right[\mathit{y}-\mathbf{F}\left(\mathit{x}\right)]+[\mathit{x}-{\mathit{x}}_{\mathrm{a}}{]}^{T}{\mathbf{S}}_{\mathrm{a}}^{-\mathrm{1}}[\mathit{x}-{\mathit{x}}_{\mathrm{a}}].\end{array}$$

There are two contributions in this cost function: the first term on the
right-hand side of Eq. (5) represents the contribution
arising from the forward model and the measurement, where
**S**_{ϵ} is the sum of the variance–covariance matrices of the
forward model and the measurement, and the second term represents the
contribution from the so-called a priori state vector which contains the prior knowledge
of the state vector before the measurement has been performed. **S**_{a}
is the variance–covariance matrix of the a priori state vector which, in our case, was
chosen to be sufficiently large to reduce the influence of the a priori assumptions on the
final retrieval.

To find the best estimate of the state vector $\widehat{\mathit{x}}$ that minimizes
the cost function Φ, an iterative method was applied following the
approach of Levenberg–Marquardt (Levenberg, 1944; Marquardt, 1963), which is
described in detail by Rodgers (2000). This approach is based on the
Newton–Gauss method to which the parameter *γ* is added that regulates
the size of each iteration step in order to diminish the cost function
compared to the previous step. The equation for this iteration may be
expressed by

$$\begin{array}{ll}{\displaystyle}{\mathit{x}}_{i+\mathrm{1}}& {\displaystyle}={\mathit{x}}_{i}+\left[\right(\mathrm{1}+\mathit{\gamma}){\mathbf{S}}_{\mathrm{a}}^{-\mathrm{1}}+{\mathbf{K}}_{i}^{T}{\mathbf{S}}_{\mathit{\u03f5}}^{-\mathrm{1}}{\mathbf{K}}_{i}{]}^{-\mathrm{1}}\mathit{\{}{\mathbf{K}}_{i}^{T}{\mathbf{S}}_{\mathit{\u03f5}}^{-\mathrm{1}}\\ \text{(6)}& {\displaystyle}& {\displaystyle}[\mathit{y}-\mathbf{F}({\mathit{x}}_{i}\left)\right]-{\mathbf{S}}_{\mathrm{a}}^{-\mathrm{1}}[{\mathit{x}}_{i}-{\mathit{x}}_{\mathrm{a}}]\mathit{\}},\end{array}$$

where **K** is the Jacobian containing the sensitivities of each of
the parameters of the state vector to each individual measurement.
**K** acts as a matrix of weights in Eq. (6) and may
also be referred to as the weighting matrix or kernel. When convergence is
reached, the variance–covariance matrix of the retrieved state vector
$\widehat{\mathit{x}}$ is given by

$$\begin{array}{}\text{(7)}& {\mathbf{S}}_{\widehat{x}}=({\mathbf{S}}_{\mathrm{a}}^{-\mathrm{1}}+{\mathbf{K}}^{T}{\mathbf{S}}_{\mathit{\u03f5}}^{-\mathrm{1}}\mathbf{K}{)}^{-\mathrm{1}},\end{array}$$

where **K** and **S**_{ϵ} correspond to the last
iteration. The matrix ${\mathbf{S}}_{\widehat{x}}$ allows us to identify the error on
each retrieved parameter. Convergence is obtained when the following
convergence test is true

$$\begin{array}{}\text{(8)}& [\mathit{y}-\mathbf{F}(\widehat{\mathit{x}}\left){]}^{T}{\mathbf{S}}_{\mathit{\u03f5}}^{-\mathrm{1}}\right[\mathit{y}-\mathbf{F}\left(\widehat{\mathit{x}}\right)]<N,\end{array}$$

where *N* is the number of elements in the measurement vector.

The application of this theoretical framework to the lidar retrieval problem
requires the definition of all necessary elements described above. The state
vector ** x** contains the desired quantities to be retrieved. These are
in our case a profile of extinction (denoted by

$$\begin{array}{ll}{\displaystyle}\mathit{x}& {\displaystyle}=\left[\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}\right({r}_{\mathrm{1}}),\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}({r}_{\mathrm{2}}),\phantom{\rule{0.125em}{0ex}}\mathrm{\cdots},\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}({r}_{j\mathrm{\_}\mathrm{bot}-\mathrm{1}}),\phantom{\rule{0.125em}{0ex}}\mathrm{IWC}({r}_{j\mathrm{\_}\mathrm{bot}}),\\ \text{(9)}& {\displaystyle}& {\displaystyle}\mathrm{\cdots},\phantom{\rule{0.125em}{0ex}}\mathrm{IWC}\left({r}_{j\mathrm{\_}\mathrm{top}}\right),\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}\left({r}_{j\mathrm{\_}\mathrm{top}+\mathrm{1}}\right),\phantom{\rule{0.125em}{0ex}}\mathrm{\cdots},\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}\left({r}_{N}\right)\phantom{\rule{0.125em}{0ex}}{]}^{T},\end{array}$$

where the subscripts *j*_bot and *j*_top denote the range
index of the bottom cloud layer and the top cloud layer, respectively. The
measurement vector ** y** consists of the logarithm of the calibrated
range-corrected lidar signal,

$$\begin{array}{ll}{\displaystyle}\mathit{y}& {\displaystyle}=\left[\phantom{\rule{0.125em}{0ex}}\mathrm{ln}\right(C\cdot P\left({r}_{\mathrm{1}}\right){r}_{\mathrm{1}}^{\mathrm{2}}),\phantom{\rule{0.125em}{0ex}}\mathrm{ln}(C\cdot P\left({r}_{\mathrm{2}}\right){r}_{\mathrm{2}}^{\mathrm{2}}),\\ \text{(10)}& {\displaystyle}& {\displaystyle}\mathrm{\cdots},\phantom{\rule{0.125em}{0ex}}\mathrm{ln}(C\cdot P({r}_{N}\left){r}_{N}^{\mathrm{2}}\right)\phantom{\rule{0.125em}{0ex}}{]}^{T}.\end{array}$$

As mentioned in Sect. 2, the measurements provide information about the backscattering particles every 15 m and the same constant vertical resolution is used in the state vector.

Following Stephens et al. (2001), the forward model **F** is given by the
lidar equation in its logarithmic and discretized form:

$$\begin{array}{ll}{\displaystyle}\mathrm{F}({x}_{j},{b}_{j})& {\displaystyle}=\mathrm{ln}\left({\mathit{\beta}}_{\mathrm{m}}\right({r}_{j})+k({r}_{j}\left){\mathit{\sigma}}_{\mathrm{p}}\right({r}_{j}\left)\right)\\ \text{(11)}& {\displaystyle}& {\displaystyle}-\mathrm{2}\sum _{l=\mathrm{1}}^{j}[{\overline{\mathit{\sigma}}}_{\mathrm{m},\mathrm{l}}+\mathit{\eta}{\overline{\mathit{\sigma}}}_{\mathrm{p},\mathrm{l}}]\mathrm{\Delta}R,\end{array}$$

defined at each range *r*_{j}, where $j=\mathrm{2},\mathrm{\dots},N$. The overline indicates layer
mean values. As discussed in Sect. 3.1, the multiple-scattering factor *η* is set to unity for aerosols and 0.75 for cirrus
clouds, and Δ*R* is the range resolution of 15 m of the lidar system.
The state vector ** x** defined in Eq. (9) contains the
IWC inside the cloud; hence the extinction

The Jacobian **K** given in Eq. (12) contains the sensitivities of the forward model to
each element of the state vector,
for which the terms F(*x*_{j},*b*_{j}), *σ*_{p}(*r*_{j}) and
IWC(*r*_{j}) have been shortened to F_{j}, *σ*_{p,j} and
IWC_{j}, respectively, to increase the readability. This short
notation will be used for all variables that are a function of range for the
remainder of this article. Inside cirrus cloud layers, the partial
derivatives are expressed by

$$\begin{array}{}\text{(13)}& {\mathrm{K}}_{ij}={\displaystyle \frac{\partial {\mathrm{F}}_{i}}{\partial {\mathrm{IWC}}_{j}}}={\displaystyle \frac{\partial {\mathrm{F}}_{i}}{\partial {\mathit{\sigma}}_{\mathrm{p},j}}}\phantom{\rule{0.33em}{0ex}}{\displaystyle \frac{\partial {\mathit{\sigma}}_{\mathrm{p},j}}{\partial {\mathrm{IWC}}_{j}}},\end{array}$$

and can be obtained from the BV2015 parametrization. The partial derivatives
with respect to extinction are calculated by differentiating Eq. (11) with respect to *σ*_{p}:

$$\begin{array}{}\text{(14)}& {\mathrm{K}}_{ij}=\left\{\begin{array}{ll}\mathrm{0}& \phantom{\rule{1em}{0ex}}\mathrm{for}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}i<j\\ -\mathrm{2}\mathit{\eta}\mathrm{\Delta}R& \phantom{\rule{1em}{0ex}}\mathrm{for}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}i>j\\ \frac{{k}_{i}}{{\mathit{\beta}}_{\mathrm{m},i}+{k}_{i}{\mathit{\sigma}}_{\mathrm{p},i}}-\mathrm{2}\mathit{\eta}\mathrm{\Delta}R& \phantom{\rule{1em}{0ex}}\mathrm{for}\phantom{\rule{0.125em}{0ex}}\phantom{\rule{0.125em}{0ex}}i=j\end{array}\right.\phantom{\rule{1em}{0ex}}.\end{array}$$

It should be noted that in the case of opaque cirrus clouds that completely
attenuate the lidar signal, the size of the measurement vector and
consequently the size of the state vector are reduced. In this case, only the
altitudes until full attenuation of the lidar signal are considered. Thus,
the size of the Jacobian is reduced as well and it contains only
*N*_{att} lines (and columns), where *N*_{att} is the number of
levels until the altitude of full attenuation.

Following Stephens et al. (2001), all variance–covariance matrices are assumed to be diagonal. Hence, the variance–covariance matrix of the a priori state vector can be defined by

$$\begin{array}{}\text{(15)}& {\mathrm{S}}_{a,ii}={\mathit{\sigma}}_{\mathrm{a},i}^{\mathrm{2}},\end{array}$$

where *σ*_{a,i} represents the variances of each of the elements of the a
priori state vector. In this study, we have assigned sufficiently large
variances to the a priori state vector in order to mainly rely on the information
contained in the measurement vector.

As mentioned above, vector ** b** in Eq. (11) represents
the non-retrieved parameters, which are each a function of altitude,

$$\begin{array}{}\text{(16)}& {\displaystyle}{b}_{\mathrm{1}}& {\displaystyle}=\mathit{\{}\phantom{\rule{0.125em}{0ex}}{\mathit{\beta}}_{\mathrm{m},j}\phantom{\rule{0.33em}{0ex}};j=\mathrm{1},\mathrm{\cdots}N\mathit{\}},\text{(17)}& {\displaystyle}{b}_{\mathrm{2}}& {\displaystyle}=\mathit{\{}\phantom{\rule{0.125em}{0ex}}{k}_{j}\phantom{\rule{0.33em}{0ex}};j=\mathrm{1},\mathrm{\cdots}N\mathit{\}},\text{(18)}& {\displaystyle}{b}_{\mathrm{3}}& {\displaystyle}=\mathit{\{}\phantom{\rule{0.125em}{0ex}}{\mathit{\eta}}_{j}\phantom{\rule{0.33em}{0ex}};j=\mathrm{1},\mathrm{\cdots}N\mathit{\}}.\end{array}$$

The molecular extinction *σ*_{m} does not need to be considered a
non-retrieved parameter because it can be obtained from *β*_{m} by
multiplication with the constant molecular backscatter-to-extinction ratio of
3∕8*π* (Fernald, 1984). The variance–covariance matrix of the forward
model and the measurement is then defined by

$$\begin{array}{}\text{(19)}& {S}_{\mathit{\u03f5},ii}={\mathit{\sigma}}_{y,i}^{\mathrm{2}}+{\mathit{\sigma}}_{b\mathrm{1},i}^{\mathrm{2}}+{\mathit{\sigma}}_{b\mathrm{2},i}^{\mathrm{2}}+{\mathit{\sigma}}_{b\mathrm{3},i}^{\mathrm{2}},\end{array}$$

where *σ*_{y,i} represents the measurement error and *σ*_{b1,i},
*σ*_{b2,i} and *σ*_{b3,i} represent the errors on the non-retrieved
parameters in the forward model calculated via

$$\begin{array}{}\text{(20)}& {\displaystyle}{\mathit{\sigma}}_{b\mathrm{1},i}& {\displaystyle}={\displaystyle \frac{{p}_{\mathit{\beta}}\left(\mathit{\%}\right)\cdot {\mathit{\beta}}_{\mathrm{m},i}}{{\mathit{\beta}}_{\mathrm{m},i}+{k}_{i}{\mathit{\sigma}}_{\mathrm{p},i}}},\text{(21)}& {\displaystyle}{\mathit{\sigma}}_{b\mathrm{2},i}& {\displaystyle}={\displaystyle \frac{{p}_{k}\left(\mathit{\%}\right)\cdot {k}_{i}{\mathit{\sigma}}_{\mathrm{p},i}}{{\mathit{\beta}}_{\mathrm{m},i}+{k}_{i}{\mathit{\sigma}}_{\mathrm{p},i}}},\text{(22)}& {\displaystyle}{\mathit{\sigma}}_{b\mathrm{3},i}& {\displaystyle}={p}_{\mathit{\eta}}\left(\mathit{\%}\right)\cdot (-\mathrm{2}\mathit{\eta}\mathrm{\Delta}R),\end{array}$$

where *p*_{β}(*%*), *p*_{k}(*%*) and *p*_{η}(*%*) represent the percentage
errors assumed for the molecular backscattering profile, the
backscatter-to-extinction ratio and the multiple-scattering factor,
respectively. As discussed in Sect. 3.1, we chose large
errors on the multiple-scattering factor for ice clouds and the
backscatter-to-extinction ratio for aerosols and quantified them to
*p*_{η}(*%*)=25 *%* and ${p}_{k,\mathrm{aer}}\left(\mathit{\%}\right)=\mathrm{25}\phantom{\rule{0.125em}{0ex}}\mathit{\%}$, respectively. The same
error has been attributed to the backscatter-to-extinction ratio for ice
clouds since the knowledge of the phase function in the exact backscattering
direction which is used in Eq. (3) is rather poor as well
(${p}_{k,\mathrm{ice}}\left(\mathit{\%}\right)=\mathrm{25}\phantom{\rule{0.125em}{0ex}}\mathit{\%}$). The error on the molecular backscattering
profile, which is obtained from the empirical equation of Flamant et al. (2008)
using the atmospheric temperature and pressure profiles from ECMWF
reanalysis, is set to *p*_{β}(*%*)=2 *%*. The error on the lidar measurement
depends on the altitude since the measurement noise increases with increasing
altitude. It is calculated as the standard deviation around the mean over a
vertically sliding window of 20 gates.

It should be noted that the characterization of the errors related to the BV2015 parametrization is very challenging. Thus, no error for the microphysical model is currently taken into account. However, this issue needs to be addressed in future studies, and an evaluation of the uncertainty arising from the Vidot et al. (2015) parametrization, which has to be integrated in future versions of our algorithm, is planned. This evaluation could be performed by comparing the single-scattering properties calculated directly from the ensemble model of Baran and Labonnote (2007) with the results from the Vidot et al. (2015) parametrization. Unfortunately, this is very costly to realize for all couples of IWC and temperature, particularly since the parametrization has not been developed by us. Furthermore, the uncertainty arising from this parametrization is assumed to be smaller than 5 % (Anthony J. Baran, personal communication, 2018). Consequently, the present study neglects an error related to the microphysical model.

To start the iteration, a first guess is required and in this study we chose to use the a priori state vector as a first guess. In order to reach faster convergence of the algorithm, the elements of the a priori state vector for the layers close to the ground where aerosols are present are calculated from a one-step solution of the lidar equation following the approach of Stephens et al. (2001),

$$\begin{array}{}\text{(23)}& {\mathit{\sigma}}_{\mathrm{p},\mathrm{a},i}=\mathit{\left\{}\mathrm{exp}\right[{y}_{i}+\mathrm{2}\sum _{l=\mathrm{1}}^{i-\mathrm{1}}({\mathit{\sigma}}_{\mathrm{m},\mathrm{l}}+\mathit{\eta}{\mathit{\sigma}}_{\mathrm{p},\mathrm{a},\mathrm{l}})\mathrm{\Delta}R]-{\mathit{\beta}}_{\mathrm{m},i}\mathit{\}}/{k}_{i}.\end{array}$$

For the layers above the boundary layer, the algorithm converges fast enough
when the molecular signal is used in the a priori state vector. For ice cloud layers, we start
the iteration from a small IWC of 0.001 g m^{−3}.

Figure 2 shows an example of a measured lidar profile (represented by the red lines) containing a cirrus cloud in altitudes between 8865 and 10 200 m measured on 30 November 2016 at 18:11 UTC. The black lines show the calculated forward model, in Fig. 2a for the a priori state vector and in Fig. 2b after the last iteration step. Figure 2c shows the relative difference (as a percentage) between the forward model and the measurement after the last iteration. Since the elements of the a priori state vector for the lowest layers have been pre-calculated based on the lidar equation, the forward model of the a priori state vector is already close to the measurement for the layers close to the ground. After the last iteration the forward model and the measured lidar signal overlay each other almost perfectly and the relative difference between the forward model and the measurement is less than 1 % for all altitudes. This indicates that the retrieval was successful and that the cost function has been reduced by reducing the difference between the measurement and the forward model.

Figure 3 presents the corresponding retrieved IWC (Fig. 3a) and extinction (Fig. 3b) profiles. As explained in Sect. 3.2, we retrieve the IWC for the layers containing a cirrus cloud and the particle extinction for the rest of the profile. Hence, the extinction profile inside the cloud in Fig. 3b is not retrieved directly but recalculated from the IWC using the BV2015 parametrization. In the layers close to the ground, an enhanced extinction due to aerosols can be observed. Above these layers in the middle portion of the profile, the particle extinction is close to zero because there were no or very few particles present in this zone. In the ice cloud an important increase in extinction due to the ice crystals can be observed.

The retrieval results for the whole afternoon of 30 November 2016 are shown
in Fig. 4. Figure 4a shows the measured lidar signal. The cloud-base and cloud-top altitudes are defined based on the threshold method described by
Platt et al. (1994). Cloud base is defined as altitude at which the lidar signal
increases above the clear background level and this increase is larger than
*n* times the standard deviation of the background fluctuations. As a second
condition it is required that the signal continues to increase for *m*
following altitude gates to assure that sudden maxima in the signal due to
measurement noise are not misinterpreted as clouds. We chose values of 4 and
5 for *n* and *m*, respectively, which are suitable for our lidar system. The
cloud top is found in the same way but by starting the search from the far
end of the measurement range (about 15 km) and moving downwards.

Figure 4b shows the retrieved extinction profiles and Fig. 4c the retrieved IWC in cirrus cloud layers. The height limit up to which the retrieval is performed depends on the profile. Since the signal above the cloud becomes noisy due to attenuation by the cloud particles, the height limit for the retrieval is defined for each profile as retrieved cloud-top altitude plus 500 m. If the lidar signal is completely attenuated or if the measured power becomes negative due to noise in lower altitudes, the upper limit of the measurement vector (and hence the state vector) is fixed at the uppermost measurable layer. In Fig. 4d the cloud optical thickness (COT) calculated from the retrieved extinction profile is compared to the COT derived from the transmission method introduced by Young (1995) and Chen et al. (2002) in which the COT is derived from the shift of the signal below and above the cloud due to the extinction of the cloud. There are phases when the COT from both methods coincides quite well, e.g., between 16:36 and 18:12 UTC. However, during this period the retrieved COT is very small. When the cloud becomes geometrically and optically thicker, our algorithm tends to overestimate the COT compared to the transmission method, e.g., between 16:06 and 16:36 UTC. At the end of the presented period between 19:24 and 20:00 UTC, the lidar signal increases importantly and for this cloud the convergence of our algorithm is less good. Figure 4e shows the value of the cost function after the iteration (normalized by the size of the measurement vector), which indicates the quality of the retrieval. For the abovementioned period the cost function is large (Φ≫1), which means that the algorithm did not converge. This is partly due to the strong attenuation of the signal because of the optically thick cloud that was present during this period, but we will show in the following paragraphs that this is also related to an insufficient description of the microphysical properties of the ice crystals, in particular the phase function in the backscattering direction.

Our retrievals strongly depend on the phase function in the backscattering
direction, which defines the backscatter-to-extinction ratio. As explained in
Sect. 3.1, we calculate the backscatter-to-extinction
ratio for cirrus clouds from Eq. (3) and hence our
retrievals strongly depend on the single-scattering albedo, *ϖ*_{0}, and
the phase function in the exact backscattering direction, *P*_{11}(*π*),
which are obtained from the BV2015 microphysical model. The single-scattering
albedo is considered to be represented sufficiently accurately in this model.
Conversely, the phase function, especially in the exact backscattering
direction, is much more uncertain.

Figure 5 shows examples of phase functions of ice
crystals computed from the ensemble model and the BV2015 parametrization (see
Sect. 3.1), for a thin cirrus cloud with a small IWC and
a temperature of 250 K. The existence of a backscattering peak strongly
depends on the characteristics of the considered particles, especially their
heterogeneity represented respectively by their surface roughness and/or by
the presence of spherical inclusions
(Hess et al., 1998a; C.-Labonnote et al., 2001; Baran and Labonnote, 2006). The black line represents a phase
function obtained from the ensemble model, for a bulk ice composed of smooth
ice particles (e.g., smooth surface with no heterogeneity), and in this case
the phase function shows a strong increase in the backscattering direction.
Introducing particle heterogeneities, e.g., surface roughness and air
bubbles, leads to the disappearance of the backscattering peak (blue and
green lines in Fig. 5, computed from the same model
and bulk ice but by considering moderately and severely heterogeneous
particles, respectively). However, real ice clouds may consist of a mixture
of smooth and rough particles and different particle sizes, and their phase
functions in the backscattering direction have not yet been characterized
sufficiently accurately. The analytic phase function of Baran et al. (2001)
(represented by the red line in Fig. 5) that is
implemented in our algorithm does not include enhanced backscattering. For
scattering angles larger than 95^{∘} the parametrization assumes a
constant value. Recent publications of Zhou and Yang (2015) and Ding et al. (2016)
suggest that this assumption is not exact enough to realistically represent
the phase function of ice crystals, even for highly heterogeneous particles.
They found that a narrow backscattering peak also exists for ice particles
with rough surfaces and that the backscattering is generally underestimated.
Zhou and Yang (2015) showed that the phase function of real bulk ice crystals at
180^{∘} should be 1.5 to 2.0 times larger than the phase function at
175^{∘}, which is clearly not the case for the analytical phase
function used in this study.

Having that in mind, we tested the influence of the backscatter-to-extinction
ratio in our algorithm. Figure 6 shows the
retrieved IWC profiles for the lidar profiles measured on 30 November 2016
at 16:20 UTC (Fig. 6a) and 18:11 UTC (Fig. 6b) for different backscatter-to-extinction ratios
*k*^{′}, where ${k}^{\prime}=\mathit{\kappa}\cdot k$ with *k* the original
backscatter-to-extinction ratio computed from Eq. (3). The
blue line represents the retrieval result for a factor of *κ*=1.0, and the
red and green lines are for modified backscatter-to-extinction ratios with factors
of *κ*=1.5 and *κ*=2.0, respectively. When the
backscatter-to-extinction ratio is enhanced, the retrieved IWC decreases. However, the effect of modifying the backscatter-to-extinction
ratio is rather strong and integrating the IWC over the whole cloud results,
for the profile measured at 18:11 UTC, in an IWP of
4.22 ± 1.01 g m^{−2} for the backscatter-to-extinction ratio
modified by a factor of *κ*=2.0 and results in 5.98 ± 1.43 g m^{−2}
for a factor of *κ*=1.5, compared to
10.32 ± 2.47 g m^{−2} with the backscatter-to-extinction ratio
calculated directly from the BV2015 microphysical model. For the
geometrically thick cloud measured at 16:20 UTC, the use of the
backscatter-to-extinction ratio computed directly from Eq. (3) results in a strongly increasing IWC towards the cloud
top, which seems to be unrealistic. This peak of IWC at the cloud top is
reduced significantly for retrievals performed with the modified
backscatter-to-extinction ratios. Furthermore, the resulting IWP is reduced
by a factor of 4 comparing the retrievals assuming *κ*=1.0
(IWP = 32.21 ± 8.93 g m^{−2}) and *κ*=2.0
(IWP = 8.58 ± 2.25 g m^{−2}).

As discussed in Sect. 3.1, the challenge of inverting the lidar equation is to find ways to constrain the backscatter-to-extinction ratio, which is the major source of uncertainty in the lidar retrieval problem. Stephens et al. (2001) included a visible optical depth in the form of an additional measurement in the optimal estimation framework to retrieve both the backscatter-to-extinction ratio and the extinction profile together. Instead of relying on a retrieval product, such as optical depth, and because the integrated amount of ice depends strongly on the backscatter-to-extinction ratio, we use TIR radiometer measurements to constrain the backscatter-to-extinction ratio of cirrus clouds.

Since TIR radiances are sensitive to the integrated properties of the cloud, in particular the IWP, we can use them to constrain the amount of ice in the cloud and hence the backscatter-to-extinction ratio. CLIMAT radiometer measurements (see Sect. 2) in its three channels are available on 30 November 2016. To simulate these measurements, the linearized discrete ordinate radiative transfer (LIDORT) model (Spurr et al., 2001) has been used. This model requires as inputs profiles of atmospheric temperature, pressure and gases, especially water vapor, which are obtained from ECMWF reanalysis. Furthermore, the optical properties of aerosol and cloud particles deduced from the retrieved IWC and extinction profiles are used in the simulations. From these profiles, an optical thickness at the lidar wavelength for each model layer is calculated, which is linked via Mie theory for aerosols and via the BV2015 parametrization for ice clouds to the optical thickness at the wavelengths in the TIR. Other necessary inputs are the single-scattering albedo and the phase function coefficients for the representation of the Legendre polynomial. As mentioned in Sect. 3.1, the aerosol characteristics are obtained from the OPAC database (Hess et al., 1998b) for the urban aerosol type. Hence, the aerosol model between the lidar and the TIR wavelengths is coherent.

The normalized radiance (emitted by a source of brightness temperature *T*_{b})
received in channel C_{i} of the radiometer is characterized by the spectral
response *f*_{i}(*λ*) of the channel and can be expressed by

$$\begin{array}{}\text{(24)}& {L}_{i}={\displaystyle \frac{{\int}_{\mathrm{\Delta}{\mathit{\lambda}}_{i}}{B}_{\mathit{\lambda}}\left({T}_{\mathrm{b}}\right){f}_{i}\left(\mathit{\lambda}\right)\mathrm{d}\mathit{\lambda}}{{\int}_{\mathrm{\Delta}{\mathit{\lambda}}_{i}}{f}_{i}\left(\mathit{\lambda}\right)\mathrm{d}\mathit{\lambda}}},\end{array}$$

where *λ* is the wavelength, *B*_{λ}(*T*_{b}) the Planck function and
Δ*λ*_{i} the spectral band pass of channel C_{i}.

Figure 7 shows the normalized radiances
measured with CLIMAT in its three channels on 30 November 2016, between 15:00
and 20:00 UTC. The black lines represent the simulation with LIDORT for a
cloud-free atmosphere taking into account the aerosol extinction in the
layers below the cloud deduced from the lidar measurements and considering
water vapor and temperature profiles from an ECMWF reanalysis at 12:00 UTC. All
three channels show an increase in the signal due to the clouds present
between 16:00 and 18:30 as well as after 19:12. However, the signal for
the first cloudy period is much smaller than for the second period because of
the smaller COT in combination with a higher and hence colder cloud-base
altitude compared to the second period. Furthermore, Fig. 7 shows that the simulated radiances for the
cloud-free atmosphere are within the error range of the measured radiances
under cloud-free conditions (between 15:00 and 16:00 UTC as well as around 19:00 UTC)
for channels C11 and C12. Conversely, the measured radiance for channel
C09 is not reproduced by the radiative transfer simulations, which may be due
to an insufficient knowledge of the spectral response function
(*f*_{i}(*λ*) in Eq. 23) for this channel,
resulting in a convergence issue of our retrieval algorithm. We therefore
decided to not take this channel into account for the remainder of this
paper.

The aim of our method is to use these TIR radiances to constrain the
backscatter-to-extinction ratio correction factor (*κ*) and hence the
phase function in the backscattering direction. In the following, we aim to
show the potential of such an approach. Figure 8a shows the simulated normalized radiances (in
$\mathrm{W}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{\mathrm{sr}}^{-\mathrm{1}}\phantom{\rule{0.125em}{0ex}}\mathrm{\mu}{\mathrm{m}}^{-\mathrm{1}}$) for the lidar profile measured on
30 November 2016 at 18:11 UTC as a function of the correction factor
*κ* for the two CLIMAT channels C11 (represented in blue) and C12
(green). The dashed lines indicate the measurements and the shaded zones
around it the measurement error. Figure 8b
presents the corresponding COT computed from the retrieved extinction profile
(red crosses), while the black line represents the COT derived from the
transmission method with its corresponding error range (shaded grey zone). With
increasing correction factor for the backscatter-to-extinction ratio, the
retrieved COT decreases, which causes the simulated radiances to decrease as
well. Furthermore, one can see from Fig. 8a
that the simulated radiances for channel C12 are within the error range of
the measurements for a correction factor between *κ*=1.1 and
*κ*=1.5, whereas for channel C11 this is the case for correction factors
larger than *κ*=1.4. This leads to the conclusion that for this profile
the correction factor for the backscatter-to-extinction ratio should range
between 1.4 and 1.5 to find a retrieval of the IWC profile that would allow
both the lidar and TIR forward model to converge towards the corresponding
measurements. Additionally, the COT computed from the retrieved extinction
profile (see Sect. 3.3.1) agrees well
with the COT derived from the transmission method for this range of
correction factors (see Fig. 8b). These results
indicate that the TIR radiances can help to refine the phase function in the
backscattering direction and that the analytic phase function of
Baran et al. (2001) may not be exact enough to represent the phase function of
real ice crystals in the exact backscattering direction.

Figure 9 shows the same analysis for the lidar
profile measured at 16:20 UTC. The COT from the optimal estimation method
largely overestimates the COT from the transmission method for this profile
when the retrieval is performed with a correction factor of *κ*=1.0.
Moreover, Fig. 6a showed a rather unrealistic
increase in IWC at the cloud top in this case. Figure 9a suggests that the correction factor
constrained by the TIR radiances should range between *κ*=2.0 and
*κ*=2.3. A correction factor from this interval reduces the retrieved
COT considerably and slightly underestimates the COT obtained from the
transmission method.

However, the curves shown in Figs. 8 and 9 underline the importance of the quality of
the measurements in the TIR since small changes in the radiances may lead to
very different retrieved microphysics. Furthermore, Dubuisson et al. (2008)
showed that the atmosphere, especially the water vapor, has a very important
influence on ground-based TIR radiometer measurements and that the
sensitivity of these measurements to cloud properties is weaker for moist
atmospheres. However, the ECMWF reanalysis profile of water vapor indicates a
rather dry atmosphere for this day with a total amount of
0.62 g m^{−2} in the atmospheric column, so the water vapor as well
as the low aerosol optical depth have a rather small influence on the TIR
radiances measured during this case study.

The results presented in this section show that the ensemble of measurements
should be used to find a retrieval that corresponds best to all available
information. As mentioned above, the optimal estimation method is a
well-adapted tool to use different kinds of measurements in a common
retrieval framework. The results shown here confirm that the TIR radiances
provide an additional constraint for the amount of ice inside the cloud, which
strongly depends on the backscatter-to-extinction ratio. Therefore, the phase
function in the backscattering direction can be constrained by the TIR
radiometer measurements under the assumption that the single-scattering
albedo is accurately known (see Eq. 3). As a
consequence, we included the TIR radiometer measurements in the optimal
estimation framework of the lidar-only algorithm to retrieve, in addition to
the IWC and extinction profiles, the correction factor *κ* for the phase
function in the backscattering direction. This newly developed synergistic
algorithm is presented in the following section.

4 Synergy algorithm lidar – TIR

Back to toptop
The synergy algorithm is an expansion of the lidar-only algorithm, which
integrates the TIR radiometer measurements in the optimal estimation method.
The new state vector contains, in addition to the elements of the previous
state vector given by Eq. (9), the correction factor
*κ* for the phase function in the exact backscattering direction:

$$\begin{array}{ll}{\displaystyle}\mathit{x}& {\displaystyle}=\left[\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}\right({r}_{\mathrm{1}}),\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}({r}_{\mathrm{2}}),\phantom{\rule{0.125em}{0ex}}\mathrm{\cdots},\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}({r}_{j\mathrm{\_}\mathrm{bot}-\mathrm{1}}),\phantom{\rule{0.125em}{0ex}}\mathrm{IWC}({r}_{j\mathrm{\_}\mathrm{bot}}),\phantom{\rule{0.125em}{0ex}}\mathrm{\cdots},\\ \text{(25)}& {\displaystyle}& {\displaystyle}\mathrm{IWC}\left({r}_{j\mathrm{\_}\mathrm{top}}\right),\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}\left({r}_{j\mathrm{\_}\mathrm{top}+\mathrm{1}}\right),\phantom{\rule{0.125em}{0ex}}\mathrm{\cdots},\phantom{\rule{0.125em}{0ex}}{\mathit{\sigma}}_{\mathrm{p}}\left({r}_{N}\right),\phantom{\rule{0.125em}{0ex}}\mathit{\kappa}\phantom{\rule{0.125em}{0ex}}{]}^{T},\end{array}$$

where the new phase function ${P}_{\mathrm{11}}^{\prime}\left(\mathit{\pi}\right)$ in the backscattering direction is related to the previous one via ${P}_{\mathrm{11}}^{\prime}\left(\mathit{\pi}\right)=\mathit{\kappa}\cdot {P}_{\mathrm{11}}\left(\mathit{\pi}\right)$.

The measurement vector, initially containing the logarithm of the calibrated range-corrected lidar signal, is expanded by the measured radiances from the two channels of the CLIMAT instrument discussed above:

$$\begin{array}{ll}{\displaystyle}\mathit{y}& {\displaystyle}=\left[\phantom{\rule{0.125em}{0ex}}\mathrm{ln}\right(C\cdot P\left({r}_{\mathrm{1}}\right){r}_{\mathrm{1}}^{\mathrm{2}}),\phantom{\rule{0.125em}{0ex}}\mathrm{ln}(C\cdot P\left({r}_{\mathrm{2}}\right){r}_{\mathrm{2}}^{\mathrm{2}}),\phantom{\rule{0.125em}{0ex}}\mathrm{\cdots},\\ \text{(26)}& {\displaystyle}& {\displaystyle}\mathrm{ln}(C\cdot P({r}_{N}\left){r}_{N}^{\mathrm{2}}\right),\phantom{\rule{0.125em}{0ex}}{L}_{\mathrm{C}\mathrm{11}},\phantom{\rule{0.125em}{0ex}}{L}_{\mathrm{C}\mathrm{12}}\phantom{\rule{0.125em}{0ex}}{]}^{T}.\end{array}$$

The forward model for the lidar is the same as in the case of the lidar-only algorithm, given by the lidar equation in the form of Eq. (11), with the only modification that the backscatter-to-extinction ratio for ice cloud layers is now calculated by

$$\begin{array}{}\text{(27)}& {k}^{\prime}={\mathit{\varpi}}_{\mathrm{0}}\cdot {P}_{\mathrm{11}}^{\prime}\left(\mathit{\pi}\right)={\mathit{\varpi}}_{\mathrm{0}}\cdot \mathit{\kappa}\cdot {P}_{\mathrm{11}}\left(\mathit{\pi}\right).\end{array}$$

The forward model for the TIR radiances is the abovementioned radiative transfer model LIDORT (Spurr et al., 2001). The advantage of this model is that it provides not only radiances but also weighting functions for atmospheric and surface parameters. That means the Jacobians for surface parameters such as emissivity or temperature; profiles of Jacobians for the temperature, atmospheric gases, or IWC profiles; and column Jacobians for the integrated quantities, for example the (Jacobian about the) integrated water vapor in the whole atmospheric column, can be obtained together with the radiances from one single simulation. Therefore the use of this model considerably reduces the computation time of the algorithm in comparison to finite-difference calculations to obtain Jacobians. This numerical efficiency allows the use of a fine vertical resolution in the radiative transfer calculations without exceeding reasonable computation times. Hence, the radiative transfer calculations can be realized on a vertical grid corresponding to the lidar resolution inside the cirrus cloud (outside the cloud the vertical resolution is defined by the ECMWF reanalysis profiles on 137 levels). That means for the radiative transfer calculations in our algorithm, the extinction profile inside the cirrus cloud is calculated with the BV2015 microphysical model from the IWC given on the lidar resolution. As a consequence, thanks to the synergy with the lidar measurements and to the numerical efficiency of the radiative transfer model, our algorithm does not have to assume a homogeneous cloud like most inversion algorithms for cloud properties do.

The Jacobian of the synergistic algorithm contains, in addition to the
Jacobian of the lidar-only algorithm, two new rows for the sensitivity of the
TIR forward model to each state vector parameter and one new column for the
sensitivity of the forward model to the new state vector element *κ* (Eq. 28).

The sensitivities of the TIR radiances to the extinction profile outside the
cloud are set to zero since they are assumed to be small. The correction
factor *κ* does not have a direct influence on the TIR radiances; thus
the last two elements in the last column of the Jacobian matrix
(**K**) are also set to zero (see Eq. 28). The
sensitivities of the TIR radiances to the IWC profile inside the cloud are
calculated directly in LIDORT. Finally, the partial derivatives of the lidar
forward model with respect to *κ* are set to zero outside the cloud and
calculated analytically as a derivation of the forward model (Eq. 11) for the ice cloud layers,

$$\begin{array}{}\text{(29)}& {\displaystyle \frac{\partial {\mathrm{F}}_{j}}{\partial \mathit{\kappa}}}={\displaystyle \frac{{\mathit{\varpi}}_{\mathrm{0}}\cdot {P}_{\mathrm{11}}\left(\mathit{\pi}\right)\cdot {\mathit{\sigma}}_{\mathrm{p},j}}{{\mathit{\beta}}_{\mathrm{m},j}+{\mathit{\varpi}}_{\mathrm{0}}\cdot {P}_{\mathrm{11}}\left(\mathit{\pi}\right)\cdot \mathit{\kappa}\cdot {\mathit{\sigma}}_{\mathrm{p},j}}},\end{array}$$

where the layer extinction *σ*_{p,j} is calculated from the IWC with the
BV2015 parametrization.

As in the lidar-only algorithm, the variance–covariance matrices in the synergy algorithm are also considered to be diagonal. Concerning the lidar, they are defined in the same way as in the lidar-only algorithm (see Sect. 3.2) with the only difference that the error for the backscatter-to-extinction ratio in ice cloud layers is no longer considered since with the new algorithm we retrieve a correction factor for the phase function in the backscattering direction that is directly related to the backscatter-to-extinction ratio. Instead, an error of 1 % on the single-scattering albedo is integrated (compare with Eq. 26). For the variance–covariance matrix of the TIR forward model the considered non-retrieved parameters (and the errors attributed to them) are the following: surface emissivity (2 %), surface temperature (1 K), and the profiles of atmospheric temperature (1 K for each layer), water vapor (10 % for each layer) and ozone (2 % for each layer). The standard deviations are calculated via

$$\begin{array}{}\text{(30)}& {\mathit{\sigma}}_{{b}_{j}}={\displaystyle \frac{\partial \mathrm{F}}{\partial {b}_{j}}}\cdot {b}_{j}\cdot {\displaystyle \frac{{p}_{{b}_{j}}\left(\mathit{\%}\right)}{\mathrm{100}}},\end{array}$$

where *b*_{j} represents the considered non-retrieved parameter, ${p}_{{b}_{j}}\left(\mathit{\%}\right)$
its error in percent and $\frac{\partial \mathrm{F}}{\partial {b}_{j}}$ the
sensitivity of the forward model to this parameter. As mentioned above, the
latter can be calculated directly in LIDORT for all desired parameters (for a
detailed description of the calculation of Jacobians in LIDORT the reader is
referred to the LIDORT User's Guide; Spurr, 2012). The elements of the
diagonal variance–covariance matrix are then given by

$$\begin{array}{}\text{(31)}& {S}_{\mathit{\u03f5},i{i}_{\mathrm{TIR}}}={\mathit{\sigma}}_{y,{i}_{\mathrm{TIR}}}^{\mathrm{2}}+\sum _{j}{\mathit{\sigma}}_{{b}_{j},{i}_{\mathrm{TIR}}}^{\mathrm{2}},\end{array}$$

where ${\mathit{\sigma}}_{y,{i}_{\mathrm{TIR}}}^{\mathrm{2}}$ represents the measurement errors for each of the two channels of the TIR radiometer. This error depends on the calibration procedure and on the temperature of the instrument during the measurement because its sensitivity is a function of temperature. The calibration of the instrument is performed in the laboratory at room temperature. Unfortunately, our instrument does not have a thermal enclosure system and during field measurements it is exposed to atmospheric temperature influences. On 30 November 2016, the atmospheric temperature was low. Due to the poor knowledge of a coefficient to correct for the instrument's temperature, the assumed errors on the measured radiances are rather large. In particular, for channel C11 the error arising from this temperature correction ranges between 15 % and 30 % depending on the value of the radiance. The largest error percentages occur for small normalized radiances in combination with a cold instrument temperature. The radiances measured by channel C12 are larger and hence the error on the measurements of this channel is smaller and ranges between 5 % and 8 %.

This section presents some preliminary results of our new algorithm. Figure 10 shows the same example as given in Fig. 2 but obtained from the synergy algorithm. The a priori assumptions for the extinction and IWC profiles in the synergy algorithm are the same as in the lidar-only algorithm (see Sect. 3.3). Figure 10b shows that once the algorithm converged, the lidar forward model and the measured lidar signal overlay each other almost perfectly. As in the lidar-only algorithm, the relative difference between the forward model and the measurement is smaller than 1 % for all layers. Thus, the good convergence found in the lidar-only algorithm is confirmed in the synergy algorithm. Table 1 summarizes the radiometer measurements and the TIR forward model (LIDORT) corresponding to this profile before and after the algorithm's convergence (expressed in normalized radiances in $\mathrm{W}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{-\mathrm{2}}\phantom{\rule{0.125em}{0ex}}{\mathrm{sr}}^{-\mathrm{1}}\phantom{\rule{0.125em}{0ex}}\mathrm{\mu}{\mathrm{m}}^{-\mathrm{1}}$). Since the values of the TIR forward model after the iteration process are within the error range of the measurements, it can be concluded that the algorithm converged in the TIR as well.

The retrieved value for the correction factor *κ* for the phase function
in the backscattering direction is 1.48 ± 0.33, which is close to the
range of 1.5 to 2.0 found in the literature (Zhou and Yang, 2015) and confirms the
result shown in Fig. 8. The corresponding
retrieved IWC and extinction profiles are shown in Fig. 11. By comparing them to the result of the lidar-only algorithm (Fig. 3), it is obvious that
the IWC and the extinction are smaller for the synergy algorithm because the
backscatter-to-extinction ratio in the ice cloud is larger. The resulting IWP
is 6.13 ± 2.19 g m^{−2} compared to the initial IWP of
10.32 ± 2.47 g m^{−2} from the lidar-only algorithm. As a
consequence, the COT is also considerably reduced to a value of
0.239 ± 0.085 compared to 0.402 ± 0.096 from the lidar-only
algorithm and is in good agreement with the COT of 0.267 ± 0.126
derived from the transmission method considering the error ranges. This
indicates that the result of the synergy algorithm is more coherent than the
result of the lidar-only algorithm because the backscatter-to-extinction
ratio has been characterized more realistically.

The application of the synergy algorithm to the profile measured at 16:20 UTC
on 30 November 2016 results in a factor *κ* of 2.15 ± 0.33.
This value for the correction
factor corresponds to the region where the simulated TIR radiances converge
to the measurements in Fig. 9. Table 2 shows the radiometer measurements and the TIR
forward model after the iteration. As for the profile at 18:11 UTC, the
values of the TIR forward model after the iteration process are within the
error range of the measurements. Thus, it can be concluded that the algorithm
found a solution allowing both the lidar and the TIR forward model to
converge towards the corresponding measurements. For the retrieved correction
factor, the large IWC peak at the cloud top obtained from the lidar-only
algorithm for a correction factor of *κ*=1.0 is considerably reduced,
which results in a more realistic shape of the IWC profile. The IWP obtained
from the synergy algorithm is 7.79 ± 2.54 g m^{−2} compared to
the abovementioned value of 32.21 ± 8.93 g m^{−2} from the
lidar-only algorithm. Hence, the retrieval with the synergy algorithm results
in an important decrease in the IWP. However, the COT of 0.304 ± 0.099
obtained from the synergy algorithm underestimates the COT of
0.608 ± 0.186 obtained from the transmission method.

Finally, Fig. 12 presents the temporal evolution
of the retrieval results from the synergy algorithm for the time period from
15:00 to 20:00 UTC on 30 November 2016. Figure 12a
reiterates the measured lidar signal already shown in Fig. 4a, Fig. 12b shows
the retrieved IWC profiles, and Fig. 12c shows the
retrieved correction factor *κ* for the phase function in the
backscattering direction (blue) and the corresponding lidar ratio in steradians
(red),
which might be easier to interpret. Figure 12d
shows the comparison of the COT obtained from the synergy algorithm (blue)
and the transmission method (red), and Fig. 12e
presents the TIR radiometer measurements for channels C12 and C11 in green
and blue, respectively, and the forward model after convergence for channel
C12 in red and for channel C11 in violet, including their uncertainties. This
plot indicates that the majority of retrievals converge well in the TIR.
Furthermore, retrieval results are only shown in the other plots of this
panel if the normalized cost function is much smaller than unity. Hence, the
large number of results shown in this figure also indicates the overall good
convergence of our algorithm. The only retrievals that did not converge
correspond to either optically very thin clouds or to clouds that are thick
enough to attenuate the lidar signal completely (e.g., around 19:48 UTC). In
the second case, this could be related to the reduced size of the Jacobian
mentioned in Sect. 3.2. However, we believe that
this non-convergence is more likely due to physical reasons since the cloud
base at this time was located in low altitudes (around 6 km) and the
temperature at this altitude was rather warm for a cirrus cloud (between 245
and 250 K). In this temperature range, the presence of supercooled liquid
droplets is possible, which is not included in the BV2015 microphysical model.
Hence, this model probably does not represent the optical properties of this
cloud accurately enough. Nevertheless, compared to the lidar-only algorithm
the synergy algorithm converged for more profiles between 19:24 and 20:00 UTC and
the retrieved COT compares well to the COT from the transmission method
during this period. Hence, the TIR helped to constrain the
backscatter-to-extinction ratio through the IWP, allowing a better coherence
between the visible and the TIR forward model.

However, between 16:00 and 18:18 UTC the COT obtained from the synergy algorithm
underestimates the COT derived from the transmission method for most of the
retrievals (except around 18:12 UTC). It should be noted that the COT obtained
from the transmission method is an effective COT ^{*} and that the real
COT depends on the multiple-scattering factor for which
$\mathrm{COT}{}^{*}=\mathit{\eta}\phantom{\rule{0.125em}{0ex}}\cdot \phantom{\rule{0.125em}{0ex}}\mathrm{COT}$. Hence, for Fig. 12c, the effective optical thickness obtained from
the transmission method has been divided by the assumed multiple-scattering
factor for ice clouds (*η*=0.75) in order to be consistent with the
retrievals from the synergy algorithm. Thus, this corrected COT depends
strongly on the assumed multiple-scattering factor. Applying a larger value
of *η* (and hence reducing the effect of multiple scattering) would reduce
the optical thickness obtained from this method and would result in a value
that is closer to the result from the synergy algorithm. Conversely,
the COT retrieved with the synergy algorithm is constrained by the TIR
radiometer measurements and remains constant when applying another multiple-scattering factor. In the synergy algorithm, the retrieval of the correction
factor *κ* for the phase function in the backscattering direction would
change and hence the microphysics of the cirrus cloud would change. The influence of the
multiple-scattering factor on the retrievals of our synergy algorithm has to
be further investigated in future studies in order to draw more sophisticated
conclusions and the retrievals shown here should be understood as a first
test to show the potential of the algorithm.

However, the multiple-scattering factor alone cannot explain the inconsistency between the COT retrieved with the synergy algorithm and the COT derived from the transmission method. Another possible reason for this discrepancy may arise from the uncertainty in the transmission method itself because it depends on a good characterization of the molecular signal above the cloud and a good estimation of the cloud-top altitude. These parameters are related to rather large uncertainties due to the quite noisy micro-pulse lidar signal in the high altitudes of cirrus clouds. Furthermore, the discrepancy between the two COTs could also originate from a potential bias in the TIR radiometer measurements due to an inaccurate temperature correction as mentioned in Sect. 4.1 or from a potential bias in the TIR forward model due to an inaccurate description of the atmospheric water vapor profile since the TIR radiometer measurements are very sensitive to water vapor (Dubuisson et al., 2008). Finally, the difference in the COTs from the synergy algorithm and the transmission method could also originate from the microphysical model, which might not be perfect. The extinction at the lidar wavelength, which is calculated based on the IWP constrained by the TIR radiometer measurements, could be slightly underestimated. Figure 12c shows that if the IWP is larger, the difference between the COTs from the two methods becomes smaller. This can be explained by the fact that the contribution of the water vapor in the TIR radiometer measurements is more important for thin clouds than for thick clouds, leading to an underestimation of the IWP and consequently an underestimation of the extinction, especially in the case of thin cirrus clouds.

Despite these limitations, the retrievals of the lidar ratio shown in Fig. 12d are promising. For the geometrically and
optically thicker cloud between 16:00 and 17:12 UTC when a considerable increase
in the measured TIR radiances was observed (see Fig. 7), the average value of the retrieved lidar
ratio is 35.6 ± 4.5 sr, which is in agreement with the lidar ratios for
cirrus clouds reported in the literature ranging between 20 and 40 sr
(e.g., Chen et al., 2002; Giannakaki et al., 2007; Josset et al., 2012; Garnier et al., 2015). Between
17:12
and 18:12 UTC the retrieved lidar ratios are much higher (on average
52.2 ± 17.0 sr) but the cloud observed during this period is optically
very thin and the signal in the TIR radiances very small, so it is not
surprising that our algorithm reaches its limit here. For the optically
thicker cloud between 19:12 and 20:00 UTC the average of the retrieved lidar
ratio is 35.8 ± 9.1 sr, corresponding to the literature again.
However, as discussed above, the retrievals of the correction factor *κ*
and thus the lidar ratio depend on the assumed multiple-scattering factor.

Nevertheless, this new synergistic algorithm suggests that using information from both, active in the visible part of the electromagnetic spectrum and passive in the TIR part, allows us to obtain new information on bulk ice optical properties, especially on the amount of ice and its capability to backscatter the visible light. Moreover, it allows us to test existing microphysical models, particularly the BV2015 model and its original representation of bulk optical properties as a function of the in-cloud temperature and IWC. The results of this study point out the overall good coherence of the BV2015 model but also its limitations in representing all the different measured profiles, especially due to the poor representation of the exact backscattering characteristics of the bulk ice.

5 Conclusions

Back to toptop
In this paper a method to retrieve IWC profiles of cirrus clouds from the synergy of ground-based lidar and TIR radiometer measurements has been presented. The algorithm is based on optimal estimation theory and combines the visible lidar and TIR radiometer measurements in a common retrieval framework to retrieve profiles of IWC together with a correction factor for the backscatter intensity of bulk ice cloud particles.

As an initial step, an algorithm to retrieve IWC and extinction profiles (outside the cloud) from the lidar measurements alone was developed. Due to the backscatter-to-extinction ambiguity arising from the combination of scattering and absorption processes in the atmosphere, assumptions are required for the backscatter-to-extinction ratio, and the retrieval results strongly depend on these assumptions. As a consequence, the challenge is to find ways to reduce the uncertainties in the retrieval arising from insufficient knowledge of the backscatter-to-extinction ratio.

To overcome the backscatter-to-extinction ambiguity, we showed in a second step that it is possible to use TIR radiances to constrain the backscatter-to-extinction ratio defined as the product of the single-scattering albedo and the phase function in the backscattering direction. The latter has not yet been fully characterized and is associated with large uncertainties. Moreover, it strongly depends on the characteristics of the particles composing the cloud. However, the BV2015 microphysical model links the optical properties of cirrus clouds directly to the IWC without the need for assumptions about the particle shape and PSD. This model allows us to obtain the single-scattering albedo and the asymmetry parameter (from which the phase function is parametrized) as a function of IWC and in-cloud temperature alone. Our algorithm benefits from the fact that TIR radiances are sensitive to the integrated IWC over the whole cloud (IWP) and that the IWC of each layer governs the optical properties via the microphysical model. That means the backscatter intensity of the ice crystals is constrained by the TIR radiances under the assumption that the single-scattering albedo is represented sufficiently accurately in the microphysical model. Consequently, our synergy algorithm retrieves a profile of IWC together with a correction factor for the phase function of the ice crystals in the exact backscattering direction, which is assumed to be constant over the entire cloud profile. Hence, the integration of the TIR radiances into the optimal estimation framework allows us to retrieve the lidar ratio although we use backscattering profiles from a simple micro-pulse lidar.

It is important to note that the same microphysical model has been used to compute the bulk ice optical properties (i.e., the scattering and absorption coefficients as well as the asymmetry parameter and the phase function) for all wavelengths considered in this study. The consistency of this microphysical model over a large portion of the electromagnetic spectrum ranging from the visible to the infrared ranges has been tested in numerous studies. Nevertheless, the parametrization of these optical properties as a function of IWC and temperature may introduce some uncertainty. However, a personal communication from Anthony J. Baran (2018) suggests that the error introduced by such a parametrization is rather small (smaller than 5 %). Thus, we believe that the results presented in this paper are robust and mainly point out the misrepresentation of the phase function in the exact backscattering direction, which is a key result of this study.

Another achievement of our algorithm is the integration of information from the whole atmospheric profile, accessible thanks to the active lidar measurements, in the forward modeling of the TIR radiances. Most common retrieval algorithms for passive sensors assume a homogeneous cloud and include only information about the cloud altitude from active measurements in the radiative transfer calculations. The synergy between the lidar and the TIR radiometer measurements established in this paper allows us to account for the profile of IWC in the radiative transfer model. Furthermore, the extinction of aerosols that may be present in the atmosphere is included in the TIR forward model although further information on the aerosol type is required. In this study, the aerosol optical properties were fixed to a predefined aerosol model and an improvement of our method would be to better characterize the properties of the aerosols that are actually present during the measurement. It is worth noting that the high vertical resolution of the radiative transfer calculations in the TIR is possible thanks to the numerical efficiency of the radiative transfer model LIDORT discussed in Sect. 4.1, which allows us to obtain the radiances and Jacobians for different atmospheric parameters from a single simulation.

The results for the case study discussed in Sect. 4.2
show, for certain periods, a quite good agreement of the retrieved lidar ratios
from our synergy algorithm with the literature. When the cloud is optically
very thin, the signal in the TIR radiometer measurements is very small,
resulting in a large uncertainty in the retrieval, which seems to be rather
logical. However, it is important to keep in mind that the results depend on
several factors. All retrievals shown in this study were performed for a
multiple-scattering factor of *η*=0.75 for ice clouds, and the
backscatter-to-extinction ratio for aerosols was fixed to the value for an
urban water-soluble aerosol from the OPAC database. Changing those values may
change the retrievals, and further sensitivity studies with our algorithm are
necessary to evaluate the effects of (1) a varying multiple-scattering factor
and (2) using other aerosol models.

Furthermore, when regarding ground-based TIR radiometer measurements, a good
characterization of the surrounding atmosphere, especially the water vapor
profile, is crucial since the TIR radiances are very sensitive to water
vapor,
which is spatially and temporally highly variable. Hence, the ECMWF
reanalysis profiles used in this study, which are available for four time
steps at 00:00, 06:00, 12:00 and 18:00 UTC and have a spatial resolution of 1^{∘},
may not be accurate enough to characterize the local water vapor profile at
our measurement site during the measurement. It is certain that a better
characterization of the water vapor profile, e.g., from microwave radiometer
measurements, would help to reduce the uncertainties in our retrievals.

Finally, the quality of the measured TIR radiances plays an important role. For the case study presented here, the temperature correction of the sensitivity of the instrument results in a quite large uncertainty because of a large temperature difference between the temperatures during the measurement and the calibration. It was shown in Sect. 3.3.2 that it was not possible to simulate the clear-sky radiances measured with channel C09. Hence, this channel was not included in the analysis. This might be due to a bad characterization of the spectral response function of the instrument as mentioned above and/or to an insufficient temperature correction. Thus, another improvement of our method would be to isolate the instrument from atmospheric temperature influences.

Nevertheless, the first results obtained from this algorithm are promising and we showed that our method allows us to converge at the same time towards the measurements of two very different instruments. However, these results have to be confirmed in future studies for other measurement periods and measurement sites.

Data availability

Back to toptop
Data availability.

Data used in this paper are available upon request to the corresponding author.

Author contributions

Back to toptop
Author contributions.

FH, LCL, FP, and GB conceived the method, developed the retrieval algorithm and discussed the results. FH and LCL analyzed the data and prepared the figures. FH wrote the paper. BD conducted the calibration of the TIR radiometer. TP processed the lidar data. All co-authors reviewed the paper.

Competing interests

Back to toptop
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements

Back to toptop
Acknowledgements.

The authors thank the Région Hauts-de-France, the Ministère de
l'Enseignement Supérieur et de la Recherche (CPER Climibio) and the
European Fund for Regional Economic Development for their financial support.
The authors thank the CaPPA project (Chemical and Physical Properties of the
Atmosphere) funded by the French National Research Agency (ANR) through the
PIA (Programme d'Investissement d'Avenir) under contract
ANR-11-LABX-0005-01 and by the Regional Council “Hauts-de-France” and the
European Funds for Regional Economic Development (FEDER). The ACTRIS-FR
research infrastructure is acknowledged for financial support.

Edited by: Alexander Kokhanovsky

Reviewed by: two anonymous referees

References

Back to toptop
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., Wandinger, U., Riebesell, M., Weitkamp, C., and Michaelis, W.: Independent measurement of extinction and backscatter profiles in cirrus clouds by using a combined Raman elastic-backscatter lidar, Appl. Optics, 31, 7113–7131, https://doi.org/10.1364/AO.31.007113, 1992. a

Ansmann, A., Bösenberg, J., Brogniez, G., Elouragini, S., Flamant, P. H., Klapheck, K., Linn, H., Menenger, L., Michaelis, W., Riebesell, M., Senff, C., Thro, P.-Y., Wandinger, U., and Weitkamp, C.: Lidar network observations of cirrus morphological and scattering properties during the International Cirrus Experiment 1989: The 18 october 1989 case study and statistical analysis, J. Appl. Meteorol., 32, 1608–1622, https://doi.org/10.1175/1520-0450(1993)032<1608:LNOOCM>2.0.CO;2, 1993. a

Baran, A. J. and Francis, P. N.: On the radiative properties of cirrus cloud at solar and thermal wavelengths: A test of model consistency using high-resolution airborne radiance measurements, J. Quant. Spectrosc. Ra., 130, 763–778, https://doi.org/10.1256/qj.03.151, 2004. a

Baran, A. J. and Labonnote, L. C.: On the reflection and polarisation properties of ice cloud, J. Quant. Spectrosc. Ra., 100, 41–54, https://doi.org/10.1016/j.jqsrt.2005.11.062, 2006. a

Baran, A. J. and Labonnote, L. C.: A self-consistent scattering model for cirrus. I: The solar region, Q. J. Roy. Meteor. Soc., 133, 1899–1912, https://doi.org/10.1002/qj.164, 2007. a, b, c, d, e

Baran, A. J., Francis, P. N., Labonnote, L. C., and Doutriaux-Boucher, M.: A scattering phase function for ice cloud: Tests of applicability using aircraft and satellite multi-angle multi-wavelength radiance measurements of cirrus, Q. J. Roy. Meteor. Soc., 127, 2395–2416, https://doi.org/10.1002/qj.49712757711, 2001. a, b, c, d, e, f

Baran, A. J., Connolly, P. J., Heymsfield, A. J., and Bansemer, A.: Using in situ estimates of ice water content, volume extinction coefficient, and the total solar optical depth obtained during the tropical ACTIVE campaign to test an ensemble model of cirrus ice crystals, Q. J. Roy. Meteor. Soc., 137, 199–218, https://doi.org/10.1002/qj.731, 2011. a

Baran, A. J., Cotton, R., Furtado, K., Havemann, S., Labonnote, L. C., Marenco, F., Smith, A., and Thelen, J.-C.: A self-consistent scattering model for cirrus. II: The high and low frequencies, Q. J. Roy. Meteor. Soc., 140, 1039–1057, https://doi.org/10.1002/qj.2193, 2014a. a, b, c

Baran, A. J., Hill, P., Furtado, K., Field, P., and Manners, J.: A coupled cloud physics-radiation parameterization of the bulk optical properties of cirrus and its impact on the Met Office Unified Model Global Atmosphere 5.0 configuration, J. Climate, 27, 7725–7752, https://doi.org/10.1175/JCLI-D-13-00700.1, 2014b. a, b

Berthier, S., Chazette, P., Pelon, J., and Baum, B.: Comparison of cloud statistics from spaceborne lidar systems, Atmos. Chem. Phys., 8, 6965–6977, https://doi.org/10.5194/acp-8-6965-2008, 2008. a

Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V.-M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh, S. K., Sherwood, S., Stevens, B., and Zhang, X. Y.: Clouds and Aerosols, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 2013. a

Brogniez, G., Pietras, C., Legrand, M., Dubuisson, P., and Haeffelin, M.: A high-accuracy multiwavelength radiometer for in situ measurements in the thermal infrared. Part II: Behavior in field experiments, J. Atmos. Ocean. Tech., 20, 1023–1033, https://doi.org/10.1175/1520-0426(2003)20<1023:AHMRFI>2.0.CO;2, 2003. a

Campbell, J. R., Vaughan, M. A., Oo, M., Holz, R. E., Lewis, J. R., and Welton, E. J.: Distinguishing cirrus cloud presence in autonomous lidar measurements, Atmos. Meas. Tech., 8, 435–449, https://doi.org/10.5194/amt-8-435-2015, 2015. a, b

Campbell, J. R., Lolli, S., Lewis, J. R., Gu, Y., and Welton, E. J.: Daytime cirrus cloud top-of-the-atmosphere radiative forcing properties at a midlatitude site and their global consequences, J. Appl. Meteorol. Clim., 55, 1667–1679, https://doi.org/10.1175/JAMC-D-15-0217.1, 2016. a

Chen, W.-N., Chiang, C.-W., and Nee, J.-B.: Lidar ratio and depolarization ratio for cirrus clouds, Appl. Optics, 41, 6470–6476, https://doi.org/10.1364/AO.41.006470, 2002. a, b, c

Chiriaco, M., Chepfer, H., Noel, V., Delaval, A., Haeffelin, M., Dubuisson, P., and Yang, P.: Improving retrievals of cirrus cloud particle size coupling lidar and three-channel radiometric techniques, Mon. Weather Rev., 132, 1648–1700, https://doi.org/10.1175/1520-0493(2004)132<1684:IROCCP>2.0.CO;2, 2004. a

C.-Labonnote, L., Brogniez, G., Buriez, J.-C., Doutriaux-Boucher, M., Gayet, J.-F., and Macke, A.: Polarized light scattering by inhomogeneous hexagonal monocrystals: Validation with ADEOS-POLDER measurements, J. Geophys. Res., 106, 12139–12153, https://doi.org/10.1029/2000JD900642, 2001. a, b

Comstock, J. M. and Sassen, K.: Retrieval of cirrus cloud radiative and backscattering properties using combined lidar and infrared radiometer (LIRAD) measurements, J. Atmos. Ocean. Tech., 18, 1658–1673, https://doi.org/10.1175/1520-0426(2001)018<1658:ROCCRA>2.0.CO;2, 2001. a

Córdoba-Jabonero, C., Lopes, F. J. S., Landulfo, E., Cuevas, E., Ochoa, H., and Gil-Ojeda, M.: Diversity on subtropical and polar cirrus clouds properties as derived from both ground-based lidars and CALIPSO/CALIOP measurements, Atmos. Res., 183, 151–165, https://doi.org/10.1016/j.atmosres.2016.08.015, 2017. a

Delanoë, J. and Hogan, R. J.: A variational scheme for retrieving ice cloud properties from combined radar, lidar, and infrared radiometer, J. Geophys. Res., 113, D07204, https://doi.org/10.1029/2007JD009000, 2008. a

Delanoë, J. and Hogan, R. J.: Combined CloudSat-CALIPSO-MODIS retrievals of the properties of ice clouds, J. Geophys. Res., 115, D00H29, https://doi.org/10.1029/2009JD012346, 2010. a

Ding, J., Yang, P., Holz, R. E., Platnick, S., Meyer, K. G., Vaughan, M. A., Hu, Y., and King, M. D.: Ice cloud backscatter study and comparison with CALIPSO and MODIS satellite data, Opt. Express, 24, 620–636, https://doi.org/10.1364/OE.24.000620, 2016. a

Donovan, D. P. and van Lammeren, A. C. A. P.: Cloud effective particle size and water content profile retrievals using combined lidar and radar observations: 1. Theory and examples, J. Geophys. Res., 106, 27425–27448, https://doi.org/10.1029/2001JD900243, 2001. a

Dubuisson, P., Giraud, V., Pelon, J., Cadet, B., and Yang, P.: Sensitivity of thermal infrared radiation at the top of the atmosphere and the surface to ice cloud microphysics, J. Appl. Meteorol. Clim., 47, 2545–2560, https://doi.org/10.1175/2008JAMC1805.1, 2008. a, b, c

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, b

Field, P. R., Wood, R., Brown, P. R. A., Kaye, P. H., Hirst, E., Greenaway, R., and Smith, J. A.: Ice particle interarrival times measured with a fast FSSP, J. Atmos. Ocean. Tech., 20, 249–261, https://doi.org/10.1175/1520-0426(2003)020<0249:ipitmw>2.0.co;2, 2003. a

Field, P. R., Hogan, R. J., Brown, P. R. A., Illingworth, A. J., Choularton, T. W., and Cotton, R. J.: Parametrization of ice-particle size distributions for mid-latitude stratiform cloud, Q. J. Roy. Meteor. Soc., 131, 1997–2017, https://doi.org/10.1256/qj.04.134, 2005. a, b

Field, P. R., Heymsfield, A. J., and Bansemer, A.: Snow size distribution parameterization for midlatitude and tropical ice clouds, J. Atmos. Sci., 64, 4346–4365, https://doi.org/10.1175/2007JAS2344.1, 2007. a, b

Flamant, P. H., Cuesta, J., Denneulin, M.-L., Dabas, A., and Huber, D.: ADM-Aeolus retrieval algorithms for aerosol and cloud products, Tellus, 60A, 273–286, https://doi.org/10.1111/j.1600-0870.2007.00287.x, 2008. a

Garnier, A., Pelon, J., Dubuisson, P., Faivre, M., Chomette, O., Pascal, N., and Kratz, D. P.: Retrieval of cloud properties using CALIPSO Imaging Infrared Radiometer. Part I: Effective emissivity and optical depth, J. Appl. Meteorol. Clim., 51, 1407–1425, https://doi.org/10.1175/JAMC-D-11-0220.1, 2012. a

Garnier, A., Pelon, J., Dubuisson, P., Yang, P., Faivre, M., Chomette, O., Pascal, N., Lucker, P., and Murray, T.: Retrieval of cloud properties using CALIPSO Imaging Infrared Radiometer. Part II: Effective diameter and ice water path, J. Appl. Meteorol. Clim., 52, 2582–2599, https://doi.org/10.1175/JAMC-D-12-0328.1, 2013. a

Garnier, A., Pelon, J., Vaughan, M. A., Winker, D. M., Trepte, C. R., and Dubuisson, P.: Lidar multiple scattering factors inferred from CALIPSO lidar and IIR retrievals of semi-transparent cirrus cloud optical depths over oceans, Atmos. Meas. Tech., 8, 2759–2774, https://doi.org/10.5194/amt-8-2759-2015, 2015. a, b

Giannakaki, E., Balis, D. S., Amiridis, V., and Kazadzis, S.: Optical and geometrical characteristics of cirrus clouds over a Southern European lidar station, Atmos. Chem. Phys., 7, 5519–5530, https://doi.org/10.5194/acp-7-5519-2007, 2007. a, b

Hess, M., Koelemeijer, R. B. A., and Stammes, P.: Scattering matrices of imperfect hexagonal ice crystals, J. Quant. Spectrosc. Ra., 60, 301–308, https://doi.org/10.1016/S0022-4073(98)00007-7, 1998a. a

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/1520-0477(1998)079<0831:OPOAAC>2.0.CO;2, 1998b. a, b

Inoue, T.: On the temperature and effective emissivity determination of semi-transparent cirrus clouds by bi-spectral measurements in the 10 µm window region, J. Meteorol. Soc. Jpn., 63, 88–99, https://doi.org/10.2151/jmsj1965.63.1_88, 1985. a

Inoue, T.: A cloud type classification with NOAA 7 split-window measurements, J. Geophys. Res., 92, 3991–4000, https://doi.org/10.1029/JD092iD04p03991, 1987. a

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 CALIPSO-CloudSat observations using ocean surface echo, J. Geophys. Res., 117, D05207, https://doi.org/10.1029/2011JD016959, 2012. a

Keckhut, P., Borchi, F., Bekki, S., Hauchecorne, A., and Silaouina, M.: Cirrus classification at midlatitude from systematic lidar observations, J. Appl. Meteorol. Clim., 45, 249–258, https://doi.org/10.1175/JAM2348.1, 2006. a

King, M. D., Kaufman, Y. J., Menzel, W. P., and Tanré, D.: Remote sensing of cloud, aerosol, and water vapor properties from the Moderate Resolution Imaging Spectrometer (MODIS), IEEE T. Geosci. Remote, 30, 2–27, https://doi.org/10.1109/36.124212, 1992. a

King, M. D., Menzel, W. P., Kaufman, Y. J., Tanré, D., Gao, B.-C., Platnick, S., Ackerman, S. A., Remer, L. A., Pincus, R., and Hubanks, P. A.: Cloud and aerosol properties, precipitable water, and profiles of temperature and water vapor from MODIS, IEEE T. Geosci. Remote, 41, 442–458, https://doi.org/10.1109/TGRS.2002.808226, 2003. 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

Klett, J. D.: Lidar inversion with variable backscatter/extinction ratios, Appl. Optics, 24, 1638–1643, https://doi.org/10.1364/AO.24.001638, 1985. a

Legrand, M., Pietras, C., Brogniez, G., Haeffelin, M., Abuhassan, N. K., and Sicard, M.: A high-accuracy multiwavelength radiometer for in situ measurements in the thermal infrared. Part I: Characterization of the instrument, J. Atmos. Ocean. Tech., 17, 1203–1214, https://doi.org/10.1175/1520-0426(2000)017<1203:AHAMRF>2.0.CO;2, 2000. a, b

Levenberg, K.: A method for the solution of certain non-linear problems in least squares, Q. Appl. Math., 2, 164–168, https://doi.org/10.1090/qam/10666, 1944. a

Liou, K.-N.: Review: Influence of cirrus clouds on weather and climate processes: A global perspective, Mon. Weather Rev., 114, 1167–1199, https://doi.org/10.1175/1520-0493(1986)114<1167:IOCCOW>2.0.CO;2, 1986. a, b

Liu, J. J., Li, Z. Q., Zheng, Y. F., and Cribb, M.: Cloud-base distribution and cirrus properties based on micropulse lidar measurements at a site in southeastern China, Adv. Atmos. Sci., 32, 991–1004, https://doi.org/10.1007/s00376-014-4176-2, 2015. a

Mace, G. G., Zhang, Q., Vaughan, M., Marchand, R., Stephens, G., Trepte, C., and Winker, D.: A description of hydrometeor layer occurrence statistics derived from the first year of merged Cloudsat and CALIPSO data, J. Geophys. Res., 114, D00A26, https://doi.org/10.1029/2007JD009755, 2009. a

Macke, A., Mueller, J., and Raschke, E.: Single scattering properties of atmospheric ice crystals, J. Atmos. Sci., 53, 2813–2825, https://doi.org/10.1175/1520-0469(1996)053<2813:SSPOAI>2.0.CO;2, 1996. a

Marquardt, D. W.: An algorithm for least-squares estimation of nonlinear parameters, J. Soc. Ind. Appl. Math., 11, 431–441, https://doi.org/10.1137/0111030, 1963. a

McCormick, M. P., Winker, D. M., Browell, E. V., Coakley, J. A., Gardner, C. S., Hoff, R. M., Kent, G. S., Melfi, S. H., Menzies, R. T., Platt, C. M. R., Randall, D. A., and Reagan, J. A.: Scientific investigations planned for the Lidar In-Space Technology Experiment (LITE), B. Am. Meteorol. Soc., 74, 205–214, https://doi.org/10.1175/1520-0477(1993)074<0205:SIPFTL>2.0.CO;2, 1993. a

Mishchenko, M. I., Travis, L. D., Kahn, R. A., and West, R. A.: Modeling phase functions for dustlike tropospheric aerosols using a shape mixture of randomly oriented polydisperse spheroids, J. Geophys. Res., 102, 16831–16847, https://doi.org/10.1029/96JD02110, 1997. a

Nohra, R.: Étude des propriétés macrophysique et optiques de cirrus à l'aide d'un micro-lidar sur le site de Lille, PhD thesis, Université de Lille 1 Sciences et Technologies, Ecole Doctorale: Sciences de la Matière, du Rayonnement et de l'Environnement, 2016. a

Pandit, A. K., Gadhavi, H. S., Venkat Ratnam, M., Raghunath, K., Rao, S. V. B., and Jayaraman, A.: Long-term trend analysis and climatology of tropical cirrus clouds using 16 years of lidar data set over Southern India, Atmos. Chem. Phys., 15, 13833–13848, https://doi.org/10.5194/acp-15-13833-2015, 2015. a

Parol, F., Buriez, J. C., Brogniez, G., and Fouquart, Y.: Information content of AVHRR channels 4 and 5 with respect to the effective radius of cirrus cloud particles, J. Appl. Meteorol., 30, 973–984, https://doi.org/10.1175/1520-0450-30.7.973, 1991. a

Pelon, J., Mallet, M., Mariscal, A., Goloub, P., Tanré, D., Karam, D. B., Flamant, C., Haywood, J., Pospichal, B., and Victori, S.: Microlidar observations of biomass burning aerosol over Djougou (Benin) during African Monsoon Multidisciplinary Analysis Special Observation Period 0: Dust and Biomass-Burning Experiment, J. Geophys. Res., 113, D00C18, https://doi.org/10.1029/2008JD009976, 2008. a

Platt, C. M., Young, S. A., Carswell, A. I., Pal, S. R., McCormick, M. P., Winker, D. M., DelGuasta, M., Stefanutti, L., Eberhard, W. L., Hardesty, M., Flamant, P. H., Valentin, R., Forgan, B., Gimmestad, G. G., Jäger, H., Khmelevtsov, S. S., Kolev, I., Kaprieolev, B., ren Lu, D., Sassen, K., Shamanaev, V. S., Uchino, O., Mizuno, Y., Wandiger, U., Weitkamp, C., Ansmann, A., and Wooldridge, C.: The Experimental Cloud Lidar Pilot Study (ECLIPS) for cloud-radiation research, B. Am. Meteorol. Soc., 75, 1635–1654, https://doi.org/10.1175/1520-0477(1994)075<1635:TECLPS>2.0.CO;2, 1994. a

Platt, C. M. R.: Lidar and radiometric observations of cirrus clouds, J. Atmos. Sci., 30, 1191–1204, https://doi.org/10.1175/1520-0469(1973)030<1191:LAROOC>2.0.CO;2, 1973. a

Platt, C. M. R.: Remote sounding of high clouds: I. Calculation of visible and infrared optical properties from lidar and radiometer measurements, J. Appl. Meteorol., 18, 1130–1143, https://doi.org/10.1175/1520-0450(1979)018<1130:RSOHCI>2.0.CO;2, 1979. a

Platt, C. M. R., Scott, J. C., and Dilley, A. C.: Remote sounding of high clouds. Part VI: Optical properties of midlatitude and tropical cirrus, J. Atmos. Sci., 44, 729–747, https://doi.org/10.1175/1520-0469(1987)044<0729:RSOHCP>2.0.CO;2, 1987. a

Platt, C. M. R., Young, S. A., Austin, R. T., Patterson, G. R., Mitchell, D. L., and Miller, S. D.: LIRAD observations of tropical cirrus clouds in MCTEX. Part I: Optical properties and detection of small particles in cold cirrus, J. Atmos. Sci., 59, 3145–3162, https://doi.org/10.1175/1520-0469(2002)059<3145:LOOTCC>2.0.CO;2, 2002. a

Rodgers, C. D.: Retrieval of atmospheric temperature and composition from remote measurements of thermal radiation, Rev. Geophys. Space Ge., 14, 609–624, https://doi.org/10.1029/RG014i004p00609, 1976. a

Rodgers, C. D.: Characterization and error analysis of profiles retrieved from remote sounding measurements, J. Geophys. Res., 95, 5587–5595, https://doi.org/10.1029/JD095iD05p05587, 1990. a

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding: Theory and Practice, World Scientific Pub. Co. Inc., 2000. a, b, c, d

Saito, M., Iwabuchi, H., Yang, P., Tang, G., King, M. D., and Sekiguchi, M.: Ice particle morphology and microphysical properties of cirrus clouds inferred from combined CALIOP-IIR measurements, J. Geophys. Res.-Atmos., 122, 4440–4462, https://doi.org/10.1002/2016JD026080, 2017. a, b

Sassen, K.: The polarization lidar technique for cloud research: A review and current assessment, B. Am. Meteorol. Soc., 72, 1848–1866, https://doi.org/10.1175/1520-0477(1991)072<1848:TPLTFC>2.0.CO;2, 1991. a, b

Seifert, P., Ansmann, A., Müller, D., Wandinger, U., Althausen, D., Heymsfield, A. J., Massie, S. T., and Schmitt, C.: Cirrus optical properties observed with lidar, radiosonde, and satellite over the tropical Indian Ocean during the aerosol-polluted northeast and clean maritime southwest monsoon, J. Geophys. Res., 112, D17205, https://doi.org/10.1029/2006JD008352, 2007. a

Sicard, M., Spyak, P. R., Brogniez, G., Legrand, M., Abuhassan, N. K., Pietras, C., and Buis, J.-P.: Thermal-infrared field radiometer for vicarious cross-calibration: characterization and comparisons with other field instruments, Opt. Eng., 38, 345–356, https://doi.org/10.1117/1.602094, 1999. a

Sourdeval, O., Labonnote, L. C., Baran, A. J., and Brogniez, G.: A methodology for simultaneous retrieval of ice and liquid water cloud properties. Part I: Information content and case study, Q. J. Roy. Meteor. Soc., 141, 870–882, https://doi.org/10.1002/qj.2405, 2015. a

Sourdeval, O., Labonnote, L. C., Baran, A. J., Mülmenstädt, J., and Brogniez, G.: A methodology for simultaneous retrieval of ice and liquid water cloud properties. Part 2: Near-global retrievals and evaluation against A-Train products, Q. J. Roy. Meteor. Soc., 142, 3063–3081, https://doi.org/10.1002/qj.2889, 2016. a

Spurr, R. J. D.: User's Guide: LIDORT Version 3.6, RT Solutions, Inc., Cambridge, USA, 2012. a

Spurr, R. J. D., Kurosu, T. P., and Chance, K. V.: A linearized discrete ordinate radiative transfer model for atmospheric remote-sensing retrieval, J. Quant. Spectrosc. Ra., 68, 689–735, https://doi.org/10.1016/S0022-4073(00)00055-8, 2001. a, b

Stephens, G. L.: Cloud feedbacks in the climate system: A critical review, J. Climate, 18, 237–273, 2005. a

Stephens, G. L. and Webster, P. J.: Clouds and climate: Sensitivity of simple systems, J. Atmos. Sci., 38, 235–247, https://doi.org/10.1175/1520-0469(1981)038<0235:CACSOS>2.0.CO;2, 1981. a

Stephens, G. L., Engelen, R. J., Vaughan, M., and Anderson, T. L.: Toward retrieving properties of the tenuous atmosphere using space-based lidar measurements, J. Geophys. Res., 106, 28143–28157, https://doi.org/10.1029/2001JD000632, 2001. a, b, c, d, e, f, g, h

Stephens, G. L., Vane, D. G., Boain, R. J., Mace, G. G., Sassen, K., Wang, Z., Illingworth, A. J., O'Connor, E. J., Rossow, W. B., Durden, S. L., Miller, S. D., Austin, R. T., Benedetti, A., Mitrescu, C., and the CloudSat Science Team: The CloudSat mission and the A-Train: A new dimension of space-based observations of clouds and precipitation, B. Am. Meteorol. Soc., 83, 1771–1790, https://doi.org/10.1175/BAMS-83-12-1771, 2002. a

Stephens, G. L., Vane, D. G., Tanelli, S., Im, E., Durden, S., Rokey, M., Reinke, D., Partain, P., Mace, G. G., Austin, R., L'Ecuyer, T., Haynes, J., Lebsock, M., Suzuki, K., Waliser, D., Wu, D., Kay, J., Gettelman, A., Wang, Z., and Marchand, R.: CloudSat mission: Performance and early science after the first year of operation, J. Geophys. Res., 113, D00A18, https://doi.org/10.1029/2008JD009982, 2008. a

Strapp, J. W., Albers, F., Reuter, A., Korolev, A. V., Maixner, U., Rashke, E., and Vukovic, Z.: Laboratory measurements of the response of a PMS OAP-2DC, J. Atmos. Ocean. Tech., 18, 1150–1170, https://doi.org/10.1175/1520-0426(2001)018<1150:LMOTRO>2.0.CO;2, 2001. a

Stubenrauch, C. J., Rossow, W. B., Kinne, S., Ackerman, S., Cesana, G., Chepfer, H., Girolamo, L. D., Getzewich, B., Guignard, A., Heidinger, A., Maddux, B. C., Menzel, W. P., Minnis, P., Pearl, C., Platnick, S., Poulsen, C., Riedi, J., Sun-Mack, S., Walther, A., Winker, D., Zeng, S., and Zhao, G.: Assessment of global cloud datasets from satellites: Project and database initiated by the GEWEX radiation panel, B. Amer. Meteorol. Soc., 94, 1031–1049, https://doi.org/10.1175/BAMS-D-12-00117.1, 2013. a

Turner, D. D. and Eloranta, E. W.: Validating mixed-phase cloud optical depth retrieved from infrared observations with high spectral resolution lidar, IEEE Geosci. Remote S., 5, 285–288, https://doi.org/10.1109/LGRS.2008.915940, 2008. a

Vidot, J., Baran, A. J., and Brunel, P.: A new ice cloud parameterization for infrared radiative transfer simulation of cloudy radiances: Evaluation and optimization with IIR observations and ice cloud profile retrieval products, J. Geophys. Res.-Atmos., 120, 6937–6951, https://doi.org/10.1002/2015JD023462, 2015. a, b, c, d, e

Winker, D. M., Pelon, J. R., and McCormick, M. P.: The CALIPSO mission: spaceborne lidar for observation of aerosols and clouds, Proc. SPIE, 4893, Lidar Remote Sensing for Industry and Environment Monitoring III (21 March 2003), https://doi.org/10.1117/12.466539, 2003. a

Winker, D. M., Vaughan, M. A., Omar, A., Hu, Y., Powell, K. A., Liu, Z., Hunt, W. H., and Young, S. A.: Overview of the CALIPSO mission and CALIOP data processing algorithms, J. Atmos. Ocean. Tech., 26, 2310–2323, https://doi.org/10.1175/2009JTECHA1281.1, 2009. a

Winker, D. M., Pelon, J., Coakley Jr., J. A., Ackerman, S. A., Charlson, R. J., Colarco, P. R., Flamant, P., Fu, Q., Hoff, R. M., Kittaka, C., Kubar, T. L., Le Treut, H., McCormick, M. P., Mégie, G., Poole, L., Powell, K., Trepte, C., Vaughan, M. A., and Wielicki, B. A.: The CALIPSO mission: A global 3D view of aerosols and clouds, B. Am. Meteorol. Soc., 91, 1211–1230, https://doi.org/10.1175/2010BAMS3009.1, 2010. a

Wylie, D. P., Menzel, W. P., Woolf, H. M., and Strabala, K. I.: Four years of global cirrus cloud statistics using HIRS, J. Climate, 7, 1972–1986, https://doi.org/10.1175/1520-0442(1994)007<1972:FYOGCC>2.0.CO;2, 1994. a

Young, S. A.: Analysis of lidar backscatter profiles in optically thin clouds, Appl. Optics, 34, 7019–7031, https://doi.org/10.1364/AO.34.007019, 1995. a, b

Zhang, Y., Macke, A., and Albers, F.: Effect of crystal size spectrum and crystal shape on stratiform cirrus radiative forcing, Atmos. Res., 52, 59–75, https://doi.org/10.1016/S0169-8095(99)00026-5, 1999. a

Zhou, C. and Yang, P.: Backscattering peak of ice cloud particles, Opt. Express, 23, 11995–12003, https://doi.org/10.1364/OE.23.011995, 2015. a, b, c

Short summary

The paper presents a novel method to retrieve microphysical properties of cirrus clouds from the synergy of lidar and thermal infrared radiometer measurements. It highlights the advantages of combining two independent data sets resulting in a better characterization of the observed target. Our algorithm may help to improve the description of the backscattering features of the ice crystals composing the cloud and thereby improve our understanding of their interactions with atmospheric radiation.

The paper presents a novel method to retrieve microphysical properties of cirrus clouds from the...

Atmospheric Measurement Techniques

An interactive open-access journal of the European Geosciences Union