An algorithm for correcting instrumental effects in polarization lidar studies is discussed. Cross-talk between the perpendicular and parallel polarization channels and imperfect polarization of the transmitted laser beam are taken into account. On the basis of the Mueller formalism it is shown that - with certain assumptions - the combined effects of imperfect polarization of the transmitted laser pulse, non-ideal properties of transmitter and receiver optics and cross-talk between parallel and perpendicular polarization channels can be described by a single parameter, which is essentially the overall system depolarization.
© Optical Society of America
For more than 20 years polarization sensitive backscatter lidars are widely used for atmospheric remote sensing (e.g. [1, 2, 3, 4]). A polarization lidar detects the change of the polarization due to scattering in the atmosphere with respect to the transmitted linearly polarized light. The Mueller formalism represents a suitable tool to describe these changes. The 4-component Stokes vector of the transmitted laser light, =(I, Q, U, V), is defined as
E ‖,⊥ denote the two perpendicular components of the electric field with respect to the plane of polarization of the transmitted light, and 〈…〉 the temporal average. For simplicity normalized Stokes vectors and Mueller matrices are used except for the scattering matrix. For a scattering process the linear depolarization factor δ is defined 
where i ‖,⊥ are the components of the scattered light with linear polarization parallel or perpendicular to the scattering plane. For backscattering the reference plane is not defined by the scattering geometry and is chosen to be parallel to the plane of polarization of the emitted light. It is now easy to express the depolarization factor by the Stokes parameters:
In polarization lidar studies the molecular scattering is used as a reference for data evaluation and therefore the value of the molecular depolarization factor is an important parameter when analyzing the depolarization of the backscattered light. In addition to the minor wavelength dependence, there are two limiting cases for the value of the molecular depolarization factor, which are given by the scattering process [6, 7]. For scattering on the central Cabannes line the depolarization factor is δ=0.00365 at λ=532 nm and for Rayleigh scattering , which, in addition to the Cabannes scattering, also includes the rotational Raman bands, the depolarization factor is δ=0.0144 at λ=532 nm [6, 7].
The lidar equation describes the relationship between the backscatter coefficient β and the backscatter signal intensity i
where c is a system constant, z the scattering altitude, and T the atmospheric transmission.
In the following, the correction algorithm will be described in terms of backscatter ratios. The backscatter ratio S is defined as the ratio of total (molecular plus aerosol) backscatter coefficient and molecular backscatter coefficient
Similarly the parallel, perpendicular and total backscatter ratios
are expressed in terms of the parallel, perpendicular and total backscatter coefficients, β ‖, β ⊥ and β T≡β ‖+β ⊥, respectively.
The volume depolarization δ V(≡δ) and aerosol depolarization δ A are given by the ratio of perpendicular and parallel backscatter signals
In a lidar application a scattering process transforms the Stokes vector of the transmitted light in a new Stokes vector describing the backscattered light. The relation between and is given by the 4×4 Mueller matrix F,
In general, the Mueller matrix F can be calculated according to
Fs : Mueller matrix of the atmospheric scattering process
Fp : Mueller matrix of the analyzer optics.
N.B.: We use this rather complex and general mathematical treatment here to facilitate the reader’s own adoptions and modifications.
Now here, in order to reduce the number of free parameters we make the following assumptions:
1. The transmitted laser light is linearly polarized with a small unpolarized component α, i.e. =(1,1 - α, 0,0) where (|α|≪1). This includes to first order the effects of emitter optics, e.g. reflections on mirrors with polarization dependent reflectivity.
2. Only scattering in backward direction is considered. It is assumed that there is complete overlap between the transmitter and receiver beams in the altitude range of interest.
3. The scattering takes place at molecules and aerosols, where scattering particles are randomly orientated in space and for each particle with an asymmetric shape there is a corresponding mirror-symmetric particle. Then, the resulting scattering matrix Fs in backward direction is diagonal. It may be written as:
4. There is no multiple scattering.
5. The polarization analyzers are aligned exactly parallel or perpendicular to the plane of reference but are subject to cross-talk, i.e.,
where B ‖,⊥ describes the cross-talk factor.
6. The quantum efficiency of the photodetectors is independent of polarization. Thus the detected signal is proportional to total intensity after the analyzers, which is given by the first component of the final Stokes vector.
For an ideal lidar instrument (α=B ‖,⊥=0) the intensities in the two polarization channels are given by the first component of the resulting Stokes vector
Note that i ⊥+i ‖=1. For a non-ideal instrument the measured intensities are
For simplicity of the further treatment we combine α with B ⊥ resp. B ‖ into two independent parameters, C and δ C, defined as:
Similarly, the backscatter coefficients are given by,
The receiver and transmitter optics imperfections are assumed to be only a few percent (δ C, C, α, B ‖,⊥≪1). For nonchiral scatterers depolarization in the backward direction does not exceed 100% and thus β ⊥ β ‖ (see e.g. ). It follows that
(1- C)β ‖≫ C β ⊥,
and we obtain
The combined effects of imperfect polarization of the transmitted laser pulse, non-ideal properties of the transmitter and receiver optics and cross-talk between parallel and perpendicular polarization channels are now expressed in terms of one parameter δ C.
2.1 Correction formula for cross-polarized backscatter ratio
In the following a relationship between measured, uncorrected backscatter ratios , and cross-talk corrected backscatter ratios S ⊥, S ‖ is derived . In terms of the measured backscatter coefficients the measured perpendicular backscatter ratio is given by
Neglecting terms involving the product δ C δ R and solving for S ⊥ and S ‖ we arrive at the final result
In order to determine the correction factor δ C we assume S ⊥=1, i.e. a cloud consisting of liquid (assumed to be spherical, i.e. non-depolarizing) particles is observed. From eqn. 18 it follows
again ignoring terms involving the product δ C δ R. The correction factor δ C is obtained by linear regression of the measured, uncorrected backscatter ratios and . We note that δ C can be calculated from both, slope δ C/(δ C+δ R) and intercept δ R/(δ C+δ R). Alternatively, eqn. 21 can be solved for δ C,
The relation for the correction of aerosol depolarization follows directly from eqn. 8. While cross-talk corrected volume depolarization δ V could be obtained from S ‖,⊥ and eqn. 7 as well, in practice δ V should be calculated from the definition
The measured volume depolarization is given by
using eqn. 17. The normalization constant k is determined by imposing the condition =δ R at aerosol-free altitudes. Thus, k=δ R/(δ C+(1-δ C)δ R). Inserting this result
in eqn. 24 and solving for δ V yields
Note that there is no limitation in the application of eqns. 20 ff. to subsequent aerosol measurements provided δ C is suficiently small, i.e. of the same order of magnitude as δ R.
Here, we present a lidar measurement of a liquid polar stratospheric cloud (PSC) with the Alfred Wegener Institute’s multiwavelength and polarization aerosol lidar at the Primary Arctic Station of the Network for the Detection of Stratospheric Change (NDSC) in Ny-Ålesund, Spitsbergen (79°N, 12°E). For a comprehensive description see [12, 13] and citations therein. The polarization measurements utilize the second harmonic of a Nd:YAG laser at a wavelength of λ=532 nm. Typically, two polarizing beamsplitters in series are used to separate the two polarization directions parallel (‖) and perpendicular (⊥) with respect to the plane of polarization of the transmitted beam. For the winter 1996/1997 only one beam splitter was used with consequences for instrumental bias as we will show below. Interference filters with a full width at half maximum (FWHM) of 10 nm are used for background suppression. Since the detector bandwidth is wide enough to detect the full rotational Raman band the molecular depolarization factor is δ R=1.44% as discussed above.
We determine the correction parameter δ C using atmospheric observations as follows: data with a strictly linear correlation between minimum and indicating the presence of liquid, non-depolarizing particles were selected. The correction factor δ C is then calculated from the instrumental contribution to the perpendicular channel. As an example, the (,) dependence is shown in Fig. 1 for one day in winter 1996/97, an observation with an almost purely liquid PSC. A linear least-squares fit of vs. (weighted with the standard deviation) yields the correction parameter δ C. Then, knowledge of δ C is used to calculate S ⊥ and δ A for all data of a measurement period according to eqn. 20. Now all those data points which are consistent with S ⊥=1, δ A=0 within the uncertainty limits are selected and the linear regression is applied to the restricted data set, yielding an improved value of δ C; the procedure is iterated once more to check convergence. From this example we calculate a overall system depolarization of δ C=2.17%. We conclude that manufacturers’ specifications for laser depolarization and imperfections of the polarizers should in general not be relied upon. Possible sources of uncertainty are (a) degradation of dielectric coatings leading to changes in optical properties, (b) unknown polarization-dependent properties of transmitter and detector optics (c) misalignment of the detector polarization plane with respect to the transmitter polarization plane. Note that other possible instrumental error sources (e.g. lidar signal contamination from other lidar beams) can be present in lidar polarization measurements that cannot be treated with the same correction factor such as done here, but can be avoided experimentally.
4 Summary and conclusions
We have shown how the backscatter ratios, volume and aerosol depolarization can be affected by instrumental effects. A correction algorithm is described and discussed in light of polarization lidar observations of polar stratospheric clouds. In practice, only one parameter is needed to describe the instrumental effect; it can be determined experimentally. Application of this correction technique can be found in [11, 12, 13].
Helpful discussions with our colleagues at AWI and the University of Bonn are gratefully acknowledged.
References and links
2. C. M. R. Platt, “Lidar observations of a mixed-phase altostratus cloud,” J. Appl. Meteorol. , 16, 339–345 (1977) [CrossRef]
3. K. Sassen, “Depolarization of laser light backscattered by artificial clouds,” J. Appl. Meteorol. , 13, 923–933 (1973) [CrossRef]
4. R. M. Schotland, R. M. Sassen, K. Stone, and R. J., “Observations by lidar of linear depolarization ratios by hydrometeors,” J. Appl. Meteorol. , 10, 1011–1017 (1971) [CrossRef]
5. H. C. van de Hulst, “Light Scattering by Small Particles,” Dover Publications, New York(1957)
6. D.R. Bates, “Rayleigh Scattering by Air,” Planet. Space Sci. 32, 785–790 (1984) [CrossRef]
9. G. Beyerle, “Untersuchungen polarer stratosphärischer Aerosole vulkanischen Ursprungs und polarer stratosphärischer Wolken mit einem Mehrwellenlängen-Lidar auf Spitzbergen (79° N, 12° O),” Berichte zur Polarforschung 138/’94, Alfred-Wegener-Institut für Polar- und Meeresforschung, Bremerhaven (1994)
11. G. Baumgarten, “Erste Messungen des Bonner Rayleigh/Mie/Raman-Lidar auf Esrange, Schweden, zur Untersuchung von dynamisch induzierten polaren Stratosphärenwolken im Januar 1997,” Diploma thesis, IB-97-26 University of Bonn, Germany (1997)
12. J. Biele, “Polare stratosphärische Wolken: Lidar-Beobachtungen, Charakterisierung von Entstehung und Entwicklung,” Berichte zur Polarforschung 03/’99, Alfred-Wegener-Institut für Polar- und Meeresforschung, Bremerhaven (1999)
13. J. Biele, A. Tsias, B. P. Luo, K. S. Carslaw, R. Neuber, G. Beyerle, and Th. Peter, “Non-equilibrium co-existence of Solid and Liquid Particles in Arctic Stratospheric Clouds,” J. Geophy. Res. submitted (2000)