Depolarization Ratio Profiles Calibration and Observations of Aerosol and Cloud in the Tibetan Plateau Based on Polarization Raman Lidar

A brief description of the Water vapor, Cloud and Aerosol Lidar (WACAL) system is provided. To calibrate the volume linear depolarization ratio, the concept of “∆90◦-calibration” is applied in this study. This effective and accurate calibration method is adjusted according to the design of WACAL. Error calculations and analysis of the gain ratio, calibrated volume linear depolarization ratio and particle linear depolarization ratio are provided as well. In this method, the influences of the gain ratio, the rotation angle of the plane of polarization and the polarizing beam splitter are discussed in depth. Two groups of measurements with half wave plate (HWP) at angles of (0◦, 45◦) and (22.5◦, −22.5◦) are operated to calibrate the volume linear depolarization ratio. Then, the particle linear depolarization ratios measured by WACAL and CALIOP (the Cloud-Aerosol Lidar with Orthogonal Polarization) during the simultaneous observations were compared. Good agreements are found. The calibration method was applied in the third Tibetan Plateau Experiment of Atmospheric Sciences (TIPEX III) in 2013 and 2014 in China. Vertical profiles of the particle depolarization ratio of clouds and aerosol in the Tibetan Plateau were measured with WACAL in Litang (30.03◦N, 100.28◦E, 3949 m above sea level (a.s.l.)) in 2013 and Naqu (31.48◦N, 92.06◦E, 4508 m a.s.l.) in 2014. Then an analysis on the polarizing properties of the aerosol, clouds and cirrus over the Tibetan Plateau is provided. The particle depolarization ratio of cirrus clouds varies from 0.36 to 0.52, with a mean value of 0.44 ± 0.04. Cirrus clouds occurred between 5.2 and 12 km above ground level (a.g.l.). The cloud thickness ranges from 0.12 to 2.55 km with a mean thickness of 1.22 ± 0.70 km. It is found that the particle depolarization ratio of cirrus clouds become larger as the height increases. However, the increase rate of the particle depolarization ratio becomes smaller as the height increases.


Introduction
As a significant branch in the field of lidar, polarization lidar plays an important role in the detection of spatial and temporal distribution of aerosols and clouds.It is a well-established technique to measure the profile of the physical properties of clouds [1,2] and aerosol [3][4][5].Typically, polarization lidars use a linearly polarized laser and are equipped with a receiver that can separately measure the components of the return signal polarized perpendicular and parallel to the outgoing beam.The calibrated ratio of these two components is defined as depolarization ratio.Polarization-sensitive lidars have proven to be highly effective in discriminating spherical and nonspherical particles in the atmosphere.They are applied to distinguish ice clouds from water clouds and to identify layers with ice crystals in mixed-phase clouds [6][7][8].With polarization lidar measurements, the depolarization ratios of cirrus cloud are investigated.Sassen and Benson (2001) studied the midlatitude cirrus clouds by using 694 nm polarization lidar and they found the depolarization ratios of cirrus display a steady increase with height [9].Chen et al. (2002) performed the depolarization ratio measurements in Taiwan and they reported that the depolarization ratios of cirrus increase from approximately 0.2 to 0.5 in the height range of 11-15 km and are smaller than 0.3 between 16 and 17 km [10].Sunilkumar and Parameswaran (2005) measured tropical cirrus properties and they found the linear depolarization ratio within the clouds shows a small increase with decrease in temperature [11].Polarization lidar can also be used to do aerosol profiling [4,5,[12][13][14] and help to classify desert dust from other aerosols such as biomass burning aerosols, fresh smoke, and so forth.Gobbi et al. (2000) classified six particle types using depolarization ratio during the Saharan dust event observations in Crete [15].Sakai et al. (2002) reported a simple model of aerosol mixture based on the particle depolarization ratio and the Ångström exponent of aerosol backscatter coefficient [16].Murayama et al. (2004) observed Asian dust and Siberian smoke with multiwavelength polarization Raman lidar and large particle depolarization ratio difference between dust and smoke is found [17].For better understanding of the spectral dependence of the dust and smoke linear depolarization ratio on the size distribution, dual-wavelength aerosol polarization lidar [17][18][19] and triple-wavelength aerosol polarization lidar are established [20,21].With depolarization ratio, lidar ratio, Ångström exponent and other intensive optical properties of particles, the cloud phase, dust, smoke, urban, maritime and other types of aerosol are classified [22][23][24][25].In May 2006, the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) that performs global profiling of aerosols and clouds in the troposphere and lower stratosphere launched.CALIOP provides several optical properties of aerosol and clouds such as backscatter coefficient, particle linear depolarization ratio and color ratio.The global dataset allows a new view on earth's atmosphere.Algorithms based on integrated attenuated backscatter and depolarization ratio have been developed to classify aerosol and cloud [26].
In some papers, the linear depolarization ratio is defined as the ratio of the cross-polarized channel to the total channel [16,27].In the following text, the volume linear depolarization ratio is defined as the ratio of the cross-polarized to the parallel-polarized backscattered signal [28].The particle linear depolarization ratio is defined as the ratio of the perpendicular polarization component to the parallel component of aerosol (without molecular part) scattering [12,29,30].The planes of polarization of the two employed detectors are parallel and orthogonal to the plane of polarization of the transmitted linearly polarized laser.The accurate measurement of depolarization ratio depends on the determination of the system calibration constants including gain ratio, offset angle, cross-talk of polarizing beam splitter (PBS).Thus it needs careful calibration.In previous works, the Rayleigh calibration method was applied within the data analysis by assuming pure Rayleigh depolarization in an aerosol-free height range [8,31,32].Reichardt et al. (2003) presented a three-signal method for calibration of depolarization ratio.This method is applied to polar stratospheric cloud observations [33].In 2005, Alvarez et al. used a rotatable half-wave plate inserted into the optical path of the lidar receiver to obtain the system calibration constants of gain ratio and offset angle.The operational details of the technique were illustrated in the First International Satellite Cloud Climatology Project Regional Experiment of LaRC, NASA [34].In his work, by using a secondary PBS, the polarization purity of the perpendicular channel is assumed to be ensured.Therefore, the influence resulting from cross-talk of PBS are not discussed in detail.In 2009, Freudenthaler et al. described one accurate calibration method called "±45 • -calibration" (developed to ∆90 • -calibration method [35]), which solved the effect of PBS and gain ratio [28].In 2016, Bravo-Aranda et al. assessed lidar depolarization uncertainty by means of a polarimetric lidar simulator [36].
To conclude, it should be emphasized that, considering the actual lidar optical setup and measurement, the gain ratio, offset angle and the cross-talk of PBS are equally significant and need to be carefully calibrated.In the current paper, the existing "∆90 • -calibration method" is modified to fit the WAter vapor, Cloud and Aerosol Lidar (WACAL).
With a well calibrated polarization lidar, aerosol and clouds can be detected.In the current paper, Section 2 presents the WACAL system and the volume depolarization ratio calibration method.In Section 3, the results and an error analysis are provided.To cross check the calibrated depolarization ratio, the simultaneous measurements of particle depolarization ratio with WACAL and the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) were operated.The results are given as well in Section 3.
In Section 4, the application of the calibration method in the third Tibetan Plateau Experiment of Atmospheric Sciences (TIPEX III) in 2013 and 2014 in China is provided.The Tibetan Plateau is a vast elevated plateau in the middle of the Eurasian Continent with an average elevation above 4500 m above sea level, and has important roles in global and regional climate system [37].It lies at a critical and sensitive junction of four climatic systems: the Westerlies, the East Asian Monsoon, the Siberian cold polar airflow and the Indian monsoon.Consequently, the assessment of the cloud and aerosol over the Tibetan Plateau is crucial.By means of this calibration method, the volume and particle depolarization ratios measured in the Tibetan Plateau are calibrated.The polarizing properties of clouds, aerosol and cirrus clouds are presented in Section 4.

Optical Setup
A three-wavelength combined elastic-backscatter Raman lidar, Water vapor, cloud and aerosol lidar (WACAL) is established to perform continuous observations of aerosol, clouds and water vapor.The system is based on the second and third harmonic frequency of a compact, pulsed Nd:YAG laser, which emits pulses of 400, 120 and 710 mJ output energy at wavelengths of 355, 532 and 1064 nm, respectively, at a 30 Hz repetition rate.The optical receiver consists of four 308 mm diameter Newtonian telescopes.Five Hamamatsu 10721P-110 photomultipliers tubes (PMT) and one Hamamatsu G8931-20 APD are used to detect the lidar signals at wavelengths of 355, 387, 407, 532 (parallel-polarized), 532 (cross-polarized) and 1064 nm.The acquisition system employs a six-channel LICEL transient recorder including analog and photon counting modes.The vertical resolution of the signal is 3.75 m and the temporal resolution is 16 s.The details of WACAL are reported by Wu et al. in a separate paper [38].We provide the schematic diagram of WACAL in Figure 1 and the specifications of the polarization channel in Table 1, respectively.
In WACAL, depolarization ratio measurement is performed only at 532 nm.The polarization ratio of the laser at a wavelength of 532 nm is about 100:1.One telescope of the telescope assembly is utilized for the collection of elastic backscatter light in the polarization channel.One adjustable aperture is installed at the focal point of the telescope ahead of the polarization channel.The field of view (FOV) is determined by adjusting the pinhole aperture.In this study, the aperture is fixed as 2 mm and the FOV is determined as 1.3 mrad.The parallel and perpendicular backscattered light at a wavelength of 532 nm is collected to retrieve the depolarization ratio.The parallel-polarized and perpendicular-polarized light is separated by the polarizing beam splitter (PBS, CCM1-PBS-251, THORLABS) and detected by PMTs respectively.The transmissivity of the PBS for parallel-polarized signal is 95.5%, while the reflectivity for perpendicular-polarized signal is larger than 99.5%.The specifications of the optical elements of this channel are also listed in Table 1.To calibrate depolarization ratio measurements, one zero-order half wave plate of the type WPH10E-532 from THORLABS is installed in the receiver.To adjust its direction, the half wave plate is mounted in a manual continuous rotation mount of type CLR1 from THORLABS.Regardless of the calibration procedure, the half wave plate is always installed in the receiver.During the routine measurements, the direction of the half wave plate axis is set to be 0 • .

Methodology
In this paper, the volume linear depolarization ratio is defined as the ratio of the cross-polarized lidar return signal (P ⊥ ) to the parallel-polarized backscatter signal (P ) [28].
The backscatter signals can be determined according to the lidar equation Equation ( 1): where C is the system constant, β is the backscatter coefficient and α is the extinction coefficient.
The subscripts "m" and "a" indicate the backscatter and extinction coefficients caused by molecules and aerosols, respectively.P and P ⊥ can be described as Equations ( 2) and (3): and In the polarization lidar these two polarization components are separated using the polarizing beamsplitter cubes (PBS).However, the PBS is not accurate enough to discriminate P and P ⊥ .As a result, the cross-talk can decrease the accuracy of the measurement inevitably.
Here the volume linear depolarization ratio δ v [28] is calculated with P P and P S represent the power of parallel-and cross-polarized light after the half-wave plate respectively.φ is the transmitter-receiver misalignment angle (also called the rotation angle of the plane of polarization).

Calibration and Cross-Check of Depolarization Ratio
In Figure 2, the schematic diagrams for aligned and misaligned receivers are provided.Figure 2a presents a well aligned system which equips an ideal linearly polarized laser.Under this condition, the polarization plane of the laser is exactly parallel to the parallel axis of the PBS and the measured ratios of spherical particles should be zero.However, as shown in Figure 2b, misalignments between the optical axis of PBS and the transmitted beam are inevitable in general lidar systems.It leads to an erroneous calculation of the depolarization ratio [34].In this case, a half-wave plate can be used to realign the laser polarization vector with the parallel axis of the PBS.Following Figure 2c, by using a half-wave plate, the corresponding electric field vectors for the parallel and perpendicular components incident on the PBS are given by Equations ( 5) and ( 6) [34]: and where → E t represents the field transmitted through the half-wave plate, → E i represents the incident field, θ is the offset angle.The powers associated with the transmitted fields are and   (Reproduction from [34]).
After the half-wave plate, as shown in Figure 3, the polarization directions of the backscattered light have been rotated and the power of mutually perpendicular light can be rewritten as γ and φ denote the rotation angle of the HWP and the rotation of the plane of polarization respectively.Please note that angles on the right hand of PBS axis are positive and angles on the left size are negative.In the WACAL system, the PBS with type of CCM1-PBS-251 is from THORLABS and the cross talk cannot be ignored.The cross talk between the transmitted and reflected light needs to be taken into consideration.As a result, the power of transmitted and reflected light through the PBS detected by the PMTs and recorded by the data acquisition can be calculated by Equations (11)   vector with an offset angle of θ in a receiver without half-wave plate; and (c) the position of the laser polarization vector with an offset angle of θ in a receiver with a half-wave plate.(Reproduction from [34]).
After the half-wave plate, as shown in Figure 3, the polarization directions of the backscattered light have been rotated and the power of mutually perpendicular light can be rewritten as and γ and φ denote the rotation angle of the HWP and the rotation of the plane of polarization respectively.Please note that angles on the right hand of PBS axis are positive and angles on the left size are negative.
In the WACAL system, the PBS with type of CCM1-PBS-251 is from THORLABS and the cross talk cannot be ignored.The cross talk between the transmitted and reflected light needs to be taken into consideration.As a result, the power of transmitted and reflected light through the PBS detected by the PMTs and recorded by the data acquisition can be calculated by Equations ( 11) and ( 12): and The amplification factors T V and R V comprise the optical transmittance of the receiver and the electronic amplification in each channel.The exact values of the amplification are not necessary in this work.Furthermore, the measured ratio ( ) m r at height of r is defined as the ratio of R P and T P .It can be determined by Equation ( 13), which is derived by taking Equations ( 9) and ( 10) into Equations ( 11) and ( 12): where is the electro-optical gain ratio between the two channels.According to the Equations ( 11) and ( 12), the rotation angle of the plane of polarization φ can be derived corresponding to the minimum power of the perpendicular channel R P and the maximum power of the parallel channel T P .The first derivative of equations of perpendicular and parallel channels can be calculated by the equations below: and 2 2 0 2 sin( 42 The amplification factors V T and V R comprise the optical transmittance of the receiver and the electronic amplification in each channel.The exact values of the amplification are not necessary in this work.Furthermore, the measured ratio m(r) at height of r is defined as the ratio of P R and P T .It can be determined by Equation ( 13), which is derived by taking Equations ( 9) and ( 10) into Equations ( 11) and ( 12): where G = V R /V T is the electro-optical gain ratio between the two channels.According to the Equations ( 11) and ( 12), the rotation angle of the plane of polarization φ can be derived corresponding to the minimum power of the perpendicular channel P R and the maximum power of the parallel channel P T .
The first derivative of equations of perpendicular and parallel channels can be calculated by the equations below: and dP R /dγ is bigger than 0 and dP T /dγ is smaller than 0.
Consequently, as long as 2γ − φ ⊆ [0, 90], the power of the parallel channel P T is a monotonically decreasing function while the power of the perpendicular channel P R is a monotonically increasing function.
As a conclusion, the P T reaches to the maximum and P R to the minimum conversely in the case of 2γ − φ = 0 • .When 2γ − φ = 90 • , the P R achieves the maximum and P T is the minimum.In WACAL system, by rotating the half wave plate, we can derive the angle 2γ − φ corresponding to the minimum power of the perpendicular channel P R and the maximum power of the parallel channel P T .And the angle φ = 2γ is the rotation of the plane of polarization.
The retrieval of the gain ratio G is introduced in the following steps.During the process of calibration, two calibration angles of 0 • and 45 • are chosen to retrieve the gain ratio.The measured ratios can be calculated by Equations ( 16) and ( 17) respectively: and The gain ratio G is calculated by the method mentioned above.According to the Equation ( 18): the gain ratio depends on the measured ratios, the specification of the PBS (i.e., the reflectivities and the transmittances) and the rotation angle of the plane of polarization.
Similarly, the gain ratio can also be determined by means of the measurement group with HWP angles of (22.5 • , −22.5 • ): . (19) When operating the calibration process, we have to select an atmospheric region where the atmospheric condition is clean.This clean air region must be stable and free of aerosol and cloud layers.After the calibration measurements, the gain ratio and the rotation angle of the plane of polarization can be obtained.

Results of Gain Ratio, Volume Linear Depolarization Ratio and Particle Linear Depolarization Ratio
The calibration experiments were operated in Qingdao on 21 September 2015, 8, 18 and 19 January 2016.By rotating the half-wave plate, the rotation of the plane of polarization φ can be determined as 5 • .Since the calibration regions were free of aerosol, the δ v (r) can be replaced by the volume linear depolarization ratio of air molecules δ v m (r) at the wavelength of 532 nm (δ v (r) = δ v m (r)).Once φ and δ v m (r) were known, based on Equations ( 18) and (19), the gain ratio G can be obtained.
The gain ratios obtained during the calibration experiments are provided in Figures 4 and 5.The black lines in these figures indicate the gain ratio determined with measurement group of (0 • , 45 • ) while the blue lines denote the gain ratio determined with measurement group of (22.5 • , −22.5 • ).The mean value and standard deviation of gain ratio are provided as well in the panels.The red lines are the gain ratio differences between these two measurement groups.
From these calibration measurements, the standard deviation of the gain ratio from (0 • , 45 • ) group is about 0.05 while the standard deviation of gain ratio from (22.5 • , −22.5 • ) group is about 0.13.The slight deviation of gain ratios derived from these two measurement groups may result from the uncertainties of the rotation angles of the HWP.From these calibration measurements, the standard deviation of the gain ratio from (0 °, 45 °) group is about 0.05 while the standard deviation of gain ratio from (22.5 °, −22.5 °) group is about 0.13.The slight deviation of gain ratios derived from these two measurement groups may result from the uncertainties of the rotation angles of the HWP.
Based on Equation (4), the calibrated volume depolarization ratio can be represented as a function of component variables including the gain ratio G, the rotation of the plane of polarization φ , the reflectivities and the transmittances of PBS via:   From these calibration measurements, the standard deviation of the gain ratio from (0 °, 45 °) group is about 0.05 while the standard deviation of gain ratio from (22.5 °, −22.5 °) group is about 0.13.The slight deviation of gain ratios derived from these two measurement groups may result from the uncertainties of the rotation angles of the HWP.
Based on Equation (4), the calibrated volume depolarization ratio can be represented as a function of component variables including the gain ratio G, the rotation of the plane of polarization φ , the reflectivities and the transmittances of PBS via: Based on Equation (4), the calibrated volume depolarization ratio can be represented as a function of component variables including the gain ratio G, the rotation of the plane of polarization φ, the reflectivities and the transmittances of PBS via: Once the volume depolarization ratio was calibrated, the particle linear depolarization ratio δ p (r) can be retrieved via [28] where R is the backscatter ratio which can be retrieved by Fernald/Klett inversion [39,40].

WACAL and CALIOP Measurements Comparison
To cross check the performance of the calibration method, the comparison of particle depolarization ratios at wavelength of 532 nm from simultaneous measurements with WACAL and CALIOP were made.The CALIOP measurement data are downloaded from the Atmospheric Science Data Center, NASA (https://eosweb.larc.nasa.gov/project/calipso/calipso_table).During the routine atmospheric observation in Qingdao, the simultaneous observation experiments were performed for 5 days including 30 September 2013, 19 December 2013 and 4 January 2014, 12 and 28 May 2014. Figure 6 shows the measurement tracks of CALIOP and the WACAL site.In these measurements, the range between the footprints of CALIOP and WACAL site should be less than 50 km (inside the black dashed circle in Figure 6).
Once the volume depolarization ratio was calibrated, the particle linear depolarization ratio ( ) p δ r can be retrieved via [28] (1 where R is the backscatter ratio which can be retrieved by Fernald/Klett inversion [39,40].

WACAL and CALIOP Measurements Comparison
To cross check the performance of the calibration method, the comparison of particle depolarization ratios at wavelength of 532 nm from simultaneous measurements with WACAL and CALIOP were made.The CALIOP measurement data are downloaded from the Atmospheric Science Data Center, NASA (https://eosweb.larc.nasa.gov/project/calipso/calipso_table).During the routine atmospheric observation in Qingdao, the simultaneous observation experiments were performed for 5 days including 30 September 2013, 19 December 2013 and 4 January 2014, 12 and 28 May 2014. Figure 6 shows the measurement tracks of CALIOP and the WACAL site.In these measurements, the range between the footprints of CALIOP and WACAL site should be less than 50 km (inside the black dashed circle in Figure 6).From the simultaneous observations, the particle linear depolarization ratio profiles measured with WACAL and CALIOP were obtained.For the purpose of comparison, in every measurement case, the particle depolarization ratio profile from CALIOP is averaged from the profiles (blue solid dots) inside the black dashed circle in Figure 6.Comparisons of these five cases are provided in Figure 7. From these profiles, good agreements are found.Discrepancies of the particle depolarization ratio From the simultaneous observations, the particle linear depolarization ratio profiles measured with WACAL and CALIOP were obtained.For the purpose of comparison, in every measurement case, the particle depolarization ratio profile from CALIOP is averaged from the profiles (blue solid dots) inside the black dashed circle in Figure 6.Comparisons of these five cases are provided in Figure 7. From these profiles, good agreements are found.Discrepancies of the particle depolarization ratio between these two instruments mainly result from the difference between the footprint of CALIOP and the WACAL site.Considering the meteorological condition in Qingdao, the combination and interaction of atmosphere from ocean and continent make the load in the air unstable spatially and temporally.Consequently, the difference of measurement sites between WACAL and CALIOP contributes to the main deviation of particle depolarization ratio.Furthermore, because of the inhomogeneous atmosphere in the boundary layer, WACAL and CALIOP may measure different air mass during the simultaneous measurements, which result in the particle depolarization ratio differences in the lower heights.
Remote Sens. 2018, 10, x FOR PEER REVIEW 12 of 22 between these two instruments mainly result from the difference between the footprint of CALIOP and the WACAL site.Considering the meteorological condition in Qingdao, the combination and interaction of atmosphere from ocean and continent make the load in the air unstable spatially and temporally.Consequently, the difference of measurement sites between WACAL and CALIOP contributes to the main deviation of particle depolarization ratio.Furthermore, because of the inhomogeneous atmosphere in the boundary layer, WACAL and CALIOP may measure different air mass during the simultaneous measurements, which result in the particle depolarization ratio differences in the lower heights.

Error of Gain Ratio
With the calibration method proposed in this paper, the gain ratio can be determined.In order to evaluate its accuracy, the error should also be calculated.Take the measurement group of (0 °, 45 °) as an example for error analysis.The total error of the gain ratio includes systematic and random errors and can be expressed by where Δ is the uncertainty associated with the various quantities and the propagation factor ( ) f x x F is defined as

Error of Gain Ratio
With the calibration method proposed in this paper, the gain ratio can be determined.In order to evaluate its accuracy, the error should also be calculated.Take the measurement group of (0 • , 45 • ) as an example for error analysis.The total error of the gain ratio includes systematic and random errors and can be expressed by where ∆ is the uncertainty associated with the various quantities and the propagation factor The propagation factor is a factor by which the relative error in f (x) is magnified with respect to the relative error in the component variables x.
In this study, the random error is introduced by measured ratio and it is calculated via [41] The systematic error is caused by the bias of φ and it is determined by Since the F G φ is far smaller than 1, the error introduced by offset angle φ can be ignored.Therefore, the Equation ( 22) is simplified as Following Equations ( 22) and ( 23), the error of the gain ratio determined with measurement group of (22.5 • , −22.5 • ) can also be obtained via The error of measured ratio is calculated by where P background is the background of the signal.

Error of Calibrated Volume Depolarization Ratio
In this section, the error analysis of the calibrated volume depolarization ratio is provided.The total error of the volume depolarization ratio is determined with Assuming δ v = A/B, A x = ∂A/∂x and B x = ∂B/∂x, the propagation factor can be represented as From (30), this factor varies with the gain ratio G, the rotation of the plane of polarization φ and the measured ratio m.To illustrate the behavior of the volume depolarization ratio typical error, Table 2 gives some examples of the volume depolarization ratio and its propagated error from mean values of the measured ratio, the gain ratio and the rotation angle of the plane of polarization.Please note that the measured ratio error depends on measurements.The value of 5% is only used for instance in this table.It is variable in routine measurements.
From Table 2 and Figure 8, it is evident that the propagation factor of the gain ratio F δ v G is always small and is typically between 0 and 1.The propagation factors of the polarization plane F δ v φ and of the measured ratio F δ v m vary depending on the measured ratio.It can be found that in the case of the large values of the measured ratio, the volume linear depolarization ratio is mainly effected by errors of the measured ratio and the gain ratio.And in the case of the small values of the measured ratio, the volume linear depolarization ratio is mainly effected by the errors of the measured ratio and the offset angle.Thus, the error of the offset angle φ has little influence in detection of aerosols with large depolarization ratio, e.g., dust.However, it cannot be ignored since its pronounced impact on the detection of aerosols with small depolarization ratio, e.g., marine aerosol and biomass burning.From this perspective, the calibration for the offset angle is critical.Different from the influence of the error of φ, the error of the gain ratio has distinct impact on the detection of aerosols with large depolarization ratio but causes little impact on detection of aerosols with small depolarization ratio.
Table 2. Error of volume depolarization ratio propagated from the errors in the gain ratio, the measured ratio and the rotation angle of the plane of polarization.Note that the error of the measured ratio m, depends on measurements.The value 5% is only used for instance in this table.It is variable in routine measurements.and the measured ratio m .To illustrate the behavior of the volume depolarization ratio typical error, Table 2 gives some examples of the volume depolarization ratio and its propagated error from mean values of the measured ratio, the gain ratio and the rotation angle of the plane of polarization.Please note that the measured ratio error depends on measurements.The value of 5% is only used for instance in this table.It is variable in routine measurements.
From Table 2 and Figure 8, it is evident that the propagation factor of the gain ratio vary depending on the measured ratio.It can be found that in the case of the large values of the measured ratio, the volume linear depolarization ratio is mainly effected by errors of the measured ratio and the gain ratio.And in the case of the small values of the measured ratio, the volume linear depolarization ratio is mainly effected by the errors of the measured ratio and the offset angle.Thus, the error of the offset angle φ has little influence in detection of aerosols with large depolarization ratio, e.g., dust.However, it cannot be ignored since its pronounced impact on the detection of aerosols with small depolarization ratio, e.g., marine aerosol and biomass burning.From this perspective, the calibration for the offset angle is critical.Different from the influence of the error of φ , the error of the gain ratio has distinct impact on the detection of aerosols with large depolarization ratio but causes little impact on detection of aerosols with small depolarization ratio.

Error of Particle Depolarization Ratio
Once the volume depolarization ratio was calibrated, the particle depolarization ratio can be determined via Equation (21).The error of the particle depolarization ratio can be calculated by The propagation factors F δ p R and F δ p δ v are determined by Equation (23).

Observation Campaigns-TIPEX III 2013 and TIPEX III 2014
The calibration method was applied in TIPEX III in 2013 and 2014 in China.The two depolarization observation field experiments were performed in the Tibetan Plateau in Litang (30.03 • N, 100.28 • E, 3949 m above sea level (a.s.l.)) in 2013 and in Naqu (31.48 • N, 92.06 • E, 4508 m a.s.l.) in 2014.Figure 9 shows the locations and Table 3 presents the measurement time of these two campaigns.
The propagation factors

Observation Campaigns-TIPEX III 2013 and TIPEX III 2014
The calibration method was applied in TIPEX III in 2013 and 2014 in China.The two depolarization observation field experiments were performed in the Tibetan Plateau in Litang (30.03 °N, 100.28 °E, 3949 m above sea level (a.s.l.)) in 2013 and in Naqu (31.48 °N, 92.06 °E, 4508 m a.s.l.) in 2014.Figure 9 shows the locations and Table 3 presents the measurement time of these two campaigns.

Aerosol and Mid-Level Clouds Observations and Discussions
Two case studies during the campaigns listed in Table 3 were selected for analysis related to the measurements of the polarization state of the mid-level cloud systems and boundary layer aerosol.
During the TIPEX III 2013 in Litang, the depolarization ratio calibration process was performed at the night of 30 July 2013.It was rainy in that afternoon and the atmosphere was very clean after the rainfall.In Figure 10a, the gain ratio is presented.Then it is applied in the depolarization ratio measurements.The observations on 28 July during the TIPEX III 2013 in Litang are chosen as a measurement case.From the profile and temporal development of particle depolarization ratio shown in Figure 10, the particle depolarization ratio of the clouds located at a height of 3.7-4.5 km above ground level (a.g.l.) is 0.19 ± 0.03.The clouds can be classified as mixed phase clouds.From all the measurements in this campaign, most mid-level clouds over Litang have similar characteristics in summertime and the clouds in Figure 10c are proposed to be typical cases.Another layer was suspended between 2.5 km and 3.6 km a.g.l.below the clouds.The particle depolarization ratio of this layer is 0.08 ± 0.001.The backscatter coefficient of this layer is 0.031 (km sr) −1 ± 0.005 (km sr) −1 , which is much smaller than the backscatter coefficient of the clouds.Consequently, this layer is an

Aerosol and Mid-Level Clouds Observations and Discussions
Two case studies during the campaigns listed in Table 3 were selected for analysis related to the measurements of the polarization state of the mid-level cloud systems and boundary layer aerosol.
During the TIPEX III 2013 in Litang, the depolarization ratio calibration process was performed at the night of 30 July 2013.It was rainy in that afternoon and the atmosphere was very clean after the rainfall.In Figure 10a, the gain ratio is presented.Then it is applied in the depolarization ratio measurements.The observations on 28 July during the TIPEX III 2013 in Litang are chosen as a measurement case.From the profile and temporal development of particle depolarization ratio shown in Figure 10, the particle depolarization ratio of the clouds located at a height of 3.7-4.5 km above ground level (a.g.l.) is 0.19 ± 0.03.The clouds can be classified as mixed phase clouds.From all the measurements in this campaign, most mid-level clouds over Litang have similar characteristics in summertime and the clouds in Figure 10c are proposed to be typical cases.Another layer was suspended between 2.5 km and 3.6 km a.g.l.below the clouds.The particle depolarization ratio of this layer is 0.08 ± 0.001.The backscatter coefficient of this layer is 0.031 (km sr) −1 ± 0.005 (km sr) −1 , which is much smaller than the backscatter coefficient of the clouds.Consequently, this layer is an aerosol layer.From Figure 10c, a thin clean atmosphere region between the clouds and the aerosol layer can be found.Furthermore, the particle depolarization ratio near the ground is 0.04 ± 0.01.Consequently, the aerosol load in Litang near the ground is rare and free of non-spherical particles such as biomass burning aerosols (0.1 < δ p < 0.2) and pure dust (0.2 < δ p < 0.4) [23].aerosol layer.From Figure 10c, a thin clean atmosphere region between the clouds and the aerosol layer can be found.Furthermore, the particle depolarization ratio near the ground is 0.04 ± 0.01.Consequently, the aerosol load in Litang near the ground is rare and free of non-spherical particles such as biomass burning aerosols (0.1 < p δ < 0.2) and pure dust (0.2 < p δ < 0.4) [23].During the TIPEX III 2014 in Naqu, the depolarization ratio calibration process was operated at the night of 26 July 2014.On that day, the atmospheric condition was clean and cloudless.In Figure 11a, the gain ratio is provided and then it is applied to the volume linear depolarization ratio and particle linear depolarization ratio measurements.The observation on 12 August is chosen as a case study.In Figure 11b,c, temporal changes of the particle depolarization ratio of the clouds and ground aerosols are shown.Classified by the volume depolarization ratio (0.084 ± 0.01), the clouds locating at the range of 2.5-3.3 km a.g.l. are water clouds.The particle depolarization ratio of the atmosphere near ground is 0.16 ± 0.02, which indicates the load of nonspherical aerosols.Since the solar background light was strong at this moment in Naqu, the ground heating was still ongoing.We propose that the nonspherical aerosols may be the mixture of anthropogenic aerosol (such as biomass burning aerosols) and local dust.In Figure 11b, the particle depolarization ratios measured at 19:30-19:40 and 19:40-19:50 are presented.We can find that the particle depolarization ratios inside clouds are stable while the particle depolarization ratios near the ground are slight larger over time.Please note that the white areas in Figure 11c indicate too low signal-to-noise ratios (SNR) caused by the strong sun light.So that a good retrieval of depolarization ratio was not possible.During the TIPEX III 2014 in Naqu, the depolarization ratio calibration process was operated at the night of 26 July 2014.On that day, the atmospheric condition was clean and cloudless.In Figure 11a, the gain ratio is provided and then it is applied to the volume linear depolarization ratio and particle linear depolarization ratio measurements.The observation on 12 August is chosen as a case study.In Figure 11b,c, temporal changes of the particle depolarization ratio of the clouds and ground aerosols are shown.Classified by the volume depolarization ratio (0.084 ± 0.01), the clouds locating at the range of 2.5-3.3 km a.g.l. are water clouds.The particle depolarization ratio of the atmosphere near ground is 0.16 ± 0.02, which indicates the load of nonspherical aerosols.Since the solar background light was strong at this moment in Naqu, the ground heating was still ongoing.We propose that the nonspherical aerosols may be the mixture of anthropogenic aerosol (such as biomass burning aerosols) and local dust.In Figure 11b, the particle depolarization ratios measured at 19:30-19:40 and 19:40-19:50 are presented.We can find that the particle depolarization ratios inside clouds are stable while the particle depolarization ratios near the ground are slight larger over time.Please note that the white areas in Figure 11c indicate too low signal-to-noise ratios (SNR) caused by the strong sun light.So that a good retrieval of depolarization ratio was not possible.
Remote Sens. 2018, 10, x FOR PEER REVIEW 16 of 22 aerosol layer.From Figure 10c, a thin clean atmosphere region between the clouds and the aerosol layer can be found.Furthermore, the particle depolarization ratio near the ground is 0.04 ± 0.01.Consequently, the aerosol load in Litang near the ground is rare and free of non-spherical particles such as biomass burning aerosols (0.1 < p δ < 0.2) and pure dust (0.2 < p δ < 0.4) [23].During the TIPEX III 2014 in Naqu, the depolarization ratio calibration process was operated at the night of 26 July 2014.On that day, the atmospheric condition was clean and cloudless.In Figure 11a, the gain ratio is provided and then it is applied to the volume linear depolarization ratio and particle linear depolarization ratio measurements.The observation on 12 August is chosen as a case study.In Figure 11b,c, temporal changes of the particle depolarization ratio of the clouds and ground aerosols are shown.Classified by the volume depolarization ratio (0.084 ± 0.01), the clouds locating at the range of 2.5-3.3 km a.g.l. are water clouds.The particle depolarization ratio of the atmosphere near ground is 0.16 ± 0.02, which indicates the load of nonspherical aerosols.Since the solar background light was strong at this moment in Naqu, the ground heating was still ongoing.We propose that the nonspherical aerosols may be the mixture of anthropogenic aerosol (such as biomass burning aerosols) and local dust.In Figure 11b, the particle depolarization ratios measured at 19:30-19:40 and 19:40-19:50 are presented.We can find that the particle depolarization ratios inside clouds are stable while the particle depolarization ratios near the ground are slight larger over time.Please note that the white areas in Figure 11c indicate too low signal-to-noise ratios (SNR) caused by the strong sun light.So that a good retrieval of depolarization ratio was not possible.

High-Level Clouds Observations and Discussions
In July and August, the high-altitude cirrus occurrence probability is in Naqu high.High altitude cirrus clouds have been identified as one important regulator of the radiance balance of the earth-atmosphere system [42].During the TIPEX III, several cirrus conditions were observed and all of the measurements are presented in Figure 12.Based on their heights, they are divided into 21 separate cirrus systems.For further analysis, the cirrus clouds particle depolarization ratio and base/top heights are derived.In this study, the cirrus cloud base heights and top heights are determined with an algorithm that combined "differential zero-cross algorithm" and "threshold algorithm" [38].The particle depolarization ratios and heights of cirrus are summarized in Table 4.In Table 4, the mean values and standard deviations of particle depolarization ratios and heights are provided.
In Figure 13a-c, the cirrus structures, the occurrence frequency of cirrus cloud base/top heights and cirrus middle height are presented.Note that in Figure 13a, the x-axis label corresponds to the cirrus serial number in Table 4.The results show that cirrus clouds are occurring between 5.2 km and 12 km a.g.l.The cloud thickness ranges from 0.12 to 2.55 km with a mean thickness of 1.22 ± 0.70 km.The particle depolarization ratio differs from 0.36 to 0.52, with a mean value of 0.44 ± 0.04.Then the height dependence of particle depolarization ratio is investigated.In Figure 13d, the correlation between cirrus middle height and particle depolarization ratio is provided.The vertical bars present the statistical uncertainties of cirrus particle depolarization ratio and the black line is the fitted curve.Similar with the results of mid-latitude cirrus clouds from Sassen and Benson, 2001 [9] and Chen et al., 2002 [10], It is found the particle depolarization ratios of cirrus clouds become larger with the height increase.Furthermore, a two-order polynomial function is found between heights and particle depolarization ratios of cirrus clouds over the Tibetan Plateau, which means that the increase rate of particle depolarization ratio decrease as the height increase.With the increase of cirrus middle height, the temperature gradually decreases.At these lower temperatures, the irregular types of ices with larger size and greater distortion are formed, which result in the increase of particle depolarization ratio of the cirrus clouds.Additionally, the height dependence of cirrus particle depolarization ratio at lower height is stronger than that at higher heights.This should be explained by the evolution of the cirrus particles at different heights and temperatures.Under the condition of extremely high height and low temperature, the cirrus particles become large and stable, which weaken their height dependence and temperature dependence.

Conclusions
In this study, an application of the " 90 Δ °-calibration" method modified to fit the particular system of WACAL is presented, as well as the error analysis.We also presented atmospheric observations during TIPEX III 2013 and TIPEX III 2014.Using WACAL, we observed the atmosphere over the Tibetan Plateau and obtained information of the aerosol and cloud polarizing properties.
The key findings of this study are summarized as below: (1) The calibration method for volume linear depolarization ratio is presented and the error analysis is provided.In this method, the influences of the gain ratio, the rotation of the plane of polarization and the cross talk of polarizing beam splitter are determined by using a half wave plate.Two groups of calibration measurements are operated and the results are compared.The calibrated particle depolarization ratio is also compared with the particle depolarization ratio from simultaneous measurements with CALIOP and good agreements are found.(2) Based on the error analysis, the gain ratio determined in this paper is mainly effected by the random error of measured ratios at (0 °, 45 °) and (22.5 °, −22.5 °).The offset angle φ has a small effect on the determination of gain ratio.However, the polarized backscatter light must be strictly rotated by 90 ° (equivalent to half wave plate strictly rotated by 45 °) in this method.The errors of the volume depolarization ratio are introduced by the errors of the measured ratio, the offset angle and the gain ratio.The influences of them vary depending on the polarization properties of the measured air mass.
(3) The calibration method is applied in TIPEX 2013 and TIPEX 2014.By mean of this method, two particle depolarization ratio measurement cases in Litang and Naqu are provided.The investigations on the polarizing properties of the aerosol and clouds over the Tibetan Plateau are provided.The cloud phases, spherical particles and nonspherical particles are identified by using depolarization ratios.

Conclusions
In this study, an application of the "∆90 • -calibration" method modified to fit the particular system of WACAL is presented, as well as the error analysis.We also presented atmospheric observations during TIPEX III 2013 and TIPEX III 2014.Using WACAL, we observed the atmosphere over the Tibetan Plateau and obtained information of the aerosol and cloud polarizing properties.The key findings of this study are summarized as below: (1) The calibration method for volume linear depolarization ratio is presented and the error analysis is provided.In this method, the influences of the gain ratio, the rotation of the plane of polarization and the cross talk of polarizing beam splitter are determined by using a half wave plate.Two groups of calibration measurements are operated and the results are compared.The calibrated particle depolarization ratio is also compared with the particle depolarization ratio from simultaneous measurements with CALIOP and good agreements are found.(2) Based on the error analysis, the gain ratio determined in this paper is mainly effected by the random error of measured ratios at (0 • , 45 • ) and (22.5 • , −22.5 • ).The offset angle φ has a small effect on the determination of gain ratio.However, the polarized backscatter light must be strictly rotated by 90 • (equivalent to half wave plate strictly rotated by 45 • ) in this method.The errors of the volume depolarization ratio are introduced by the errors of the measured ratio, the offset angle and the gain ratio.The influences of them vary depending on the polarization properties of the measured air mass.(3) The calibration method is applied in TIPEX 2013 and TIPEX 2014.By mean of this method, two particle depolarization ratio measurement cases in Litang and Naqu are provided.The investigations on the polarizing properties of the aerosol and clouds over the Tibetan Plateau are provided.The cloud phases, spherical particles and nonspherical particles are identified by using depolarization ratios.(4) During TIPEX III in Naqu, all the measurement cases of cirrus clouds are presented.In total, 21 separate cirrus clouds are found and the particle depolarization ratio, cirrus cloud base/top heights and cirrus middle height are retrieved.According to the results, the particle depolarization ratio of cirrus clouds differs from 0.36 to 0.52, with a mean value of 0.44 ± 0.04.Cirrus clouds occurred between 5.2 km and 12 km a.g.l.The cloud thickness ranges from 0.12 to 2.55 km with a mean thickness of 1.22 ± 0.70 km.The occurrence frequencies of cirrus cloud base/top/middle heights are provided as well.Finally, the correlation between cirrus middle height and particle depolarization ratio is studied.It indicates that the particle depolarization ratio of cirrus clouds increases with the increasing height.While the increase rate of particle depolarization ratio decreases as the height increases.

Figure 2 .
Figure 2. Receiver geometries: (a) a well aligned configuration; (b) the position of the laser polarization vector with an offset angle of θ in a receiver without half-wave plate; and (c) the position of the laser polarization vector with an offset angle of θ in a receiver with a half-wave plate.

Figure 2 .
Figure 2. Receiver geometries: (a) a well aligned configuration; (b) the position of the laser polarizationvector with an offset angle of θ in a receiver without half-wave plate; and (c) the position of the laser polarization vector with an offset angle of θ in a receiver with a half-wave plate.(Reproduction from[34]).

Figure 3 .
Figure 3. Signal power components in a receiver of WACAL with a PBS with the reflectivities of P R and S R and the transmittances of P T and S T for linear polarized light parallel and perpendicular to the plane of the PBS. P , ⊥ P , P P and S P represent the power of backscatter light at the corresponding polarization direction.(Developed from [28]).

dP dγ is bigger thanFigure 3 .
Figure3.Signal power components in a receiver of WACAL with a PBS with the reflectivities of R P and R S and the transmittances of T P and T S for linear polarized light parallel and perpendicular to the plane of the PBS.P , P ⊥ , P P and P S represent the power of backscatter light at the corresponding polarization direction.(Developed from[28]).

Figure 6 .
Figure 6.The track of the CALIOP and the site of the WACAL during the simultaneous measurements.The colored lines indicate the footprints of CALIOP and the purple star denotes the WACAL site.The blue solid dots indicate the footprint of CALIOP and CALIOP provides one product profile every dot.The black dashed circle represents a 50-km radius area which is centered on WACAL.Note that the tracks of CALIOP on 19 December 2013 (purple) and 28 May 2014 (green) overlapped each other.The particle depolarization ratios are at 532 nm and the CALIOP measurement data are downloaded from the Atmospheric Science Data Center, NASA.

Figure 6 .
Figure 6.The track of the CALIOP and the site of the WACAL during the simultaneous measurements.The colored lines indicate the footprints of CALIOP and the purple star denotes the WACAL site.The blue solid dots indicate the footprint of CALIOP and CALIOP provides one product profile every dot.The black dashed circle represents a 50-km radius area which is centered on WACAL.Note that the tracks of CALIOP on 19 December 2013 (purple) and 28 May 2014 (green) overlapped each other.The particle depolarization ratios are at 532 nm and the CALIOP measurement data are downloaded from the Atmospheric Science Data Center, NASA.
typically between 0 and 1.The propagation factors of the polarization plane v

Figure 8 .
Figure 8.(a) Propagation factor of the gain ratio

Figure 8 .
Figure 8.(a) Propagation factor of the gain ratio F δ v G ; (b) propagation factor of the offset angle F δ v φ ; (c) propagation factor of the measured ratio F δ v m and (d) correlation between volume depolarization ratio and the measured ratio.

F
are determined by Equation (23).

Figure 10 .
Figure 10.(a) Gain ratio determined during TIPEX III 2013 in Litang; (b) Profile of particle depolarization ratio (dark green line) with error bar (horizontal line) and (c) temporal development of particle depolarization ratio measured on 28 July 2013 in Litang.

Figure 11 .
Figure 11.(a) Gain ratio determined during TIPEX III 2014 in Naqu; (b) Profile of particle depolarization ratio (dark green line) with error bar (horizontal line) and (c) temporal development of particle depolarization ratio measured on 12 August 2014 in Naqu.

Figure 10 .
Figure 10.(a) Gain ratio determined during TIPEX III 2013 in Litang; (b) Profile of particle depolarization ratio (dark green line) with error bar (horizontal line) and (c) temporal development of particle depolarization ratio measured on 28 July 2013 in Litang.

Figure 10 .
Figure 10.(a) Gain ratio determined during TIPEX III 2013 in Litang; (b) Profile of particle depolarization ratio (dark green line) with error bar (horizontal line) and (c) temporal development of particle depolarization ratio measured on 28 July 2013 in Litang.

Figure 11 .
Figure 11.(a) Gain ratio determined during TIPEX III 2014 in Naqu; (b) Profile of particle depolarization ratio (dark green line) with error bar (horizontal line) and (c) temporal development of particle depolarization ratio measured on 12 August 2014 in Naqu.

Figure 11 .
Figure 11.(a) Gain ratio determined during TIPEX III 2014 in Naqu; (b) Profile of particle depolarization ratio (dark green line) with error bar (horizontal line) and (c) temporal development of particle depolarization ratio measured on 12 August 2014 in Naqu.

Figure 12 .
Figure 12.Particle depolarization ratio of cirrus measured during TIPEX III in Naqu.Figure 12. Particle depolarization ratio of cirrus measured during TIPEX III in Naqu.

Figure 12 .
Figure 12.Particle depolarization ratio of cirrus measured during TIPEX III in Naqu.Figure 12. Particle depolarization ratio of cirrus measured during TIPEX III in Naqu.

Figure 13 .
Figure 13.Measurements of cirrus clouds during TIPEX III 2014 in Naqu.(a) Cirrus structure and (b) the frequency occurrence of cloud base/top heights; (c) Frequency occurrence of cirrus cloud middle height and (d) correlation between particle depolarization ratio and cirrus middle height.

Figure 13 .
Figure 13.Measurements of cirrus clouds during TIPEX III 2014 in Naqu.(a) Cirrus structure and (b) the frequency occurrence of cloud base/top heights; (c) Frequency occurrence of cirrus cloud middle height and (d) correlation between particle depolarization ratio and cirrus middle height.

Table 1 .
System specification of the depolarization channel of WACAL.

Table 1 .
System specification of the depolarization channel of WACAL. ) )tan

Table 2 .
Error of volume depolarization ratio propagated from the errors in the gain ratio, the measured ratio and the rotation angle of the plane of polarization.Note that the error of the measured ratio m , depends on measurements.The value 5% is only used for instance in this table.It is variable in routine measurements.

Table 3 .
Summary of the measurement of the water vapor, cloud and aerosol lidar.

Table 4 .
Heights and depolarization ratio of cirrus measured in Naqu.