Lidar at 1064 nm and Ka-band millimetre-wave cloud radar (MMCR) are powerful tools for detecting the height distribution of cloud boundaries and can monitor the entire life cycle of cloud layers. In this study, lidar and MMCR are employed to jointly detect cloud boundaries under different conditions. By enhancing the echo signal of lidar at 1064 nm and combining its signal-to-noise ratio (SNR), the cloud signal can be accurately extracted from the aerosol signals and background noise. The interference signal is eliminated from Doppler spectra of the MMCR by using the noise ratio of the smallest measurable cloud signal (SNR_min⁡) and the spectral point continuous threshold (N_ts). Moreover, the quality control of the reflectivity factor of MMCR obtained by the inversion is conducted, which improves the detection accuracy of the cloud signal. We analysed three typical cases studies; case one presents two interesting phenomena: (a) at 19:00–20:00 CST (China standard time), the ice crystal particles at the cloud top boundary are too small to be detected by MMCR, but they are well detected by lidar. (b) At 19:00–00:00 CST, the cirrus cloud changes to altostratus where the cloud particles eventually grow into large sizes, producing precipitation. Further, MMCR has more advantages than lidar in detection of the cloud top boundary within this period. Considering the advantages of the two devices, the change characteristics of the cloud boundary in Xi'an from December 2020 to November 2021 were analysed, with MMCR detection data as the main data and lidar data as the assistant data. The seasonal variation characteristics of clouds show that, in most cases, high clouds often occur in summer and autumn, and the low clouds are usually in winter. The normalized cloud cover shows that the maximum and minimum cloud cover occur in summer and winter, respectively. Furthermore, the cloud boundary frequency distribution results for the whole of the observation period show that the cloud bottom boundary below 1.5 km is more than 1 %, the frequency within the height range of 3.06–3.6 km is approximately 0.38 %, and the frequency above 8 km is less than 0.2 %. The cloud top boundary frequency distribution exhibits the characteristics of a bimodal distribution. The first narrow peak lies at approximately 1.0–3.1 km, and the second peak appears at 6.4–9.8 km.

The use of Raman and high-spectral lidars enables measurements of a stratospheric aerosol extinction profile independent of backscatter, and a multi-wavelength (MW) lidar can obtain additional information that can aid in retrieving the microphysical characteristics of the sampled aerosol. The inversion method for retrieving aerosol particle size distributions and microphysical particle parameters from MW lidar data was studied. An inversion algorithm for retrieving aerosol particle size distributions based on the regularization method was established. Based on the inversion of regularization, the inversion method was optimized by choosing the base function closest to the aerosol distribution. The logarithmic normal distribution function was selected over the triangle function as the base function for the inversion. The averaging procedure was carried out for three main types of aerosol. The 1% averaging result near the minimum of the discrepancy gave the best estimate of the particle parameters. The accuracy and stabilization of the optimized algorithm for microphysical parameters were tested by scores of aerosol size distributions. The systematic effects and random errors impacting the inversion were also considered, and the algorithm was tested by the data, showing 10% systematic error and 15% random error. At the same time, the reliability of the proposed algorithm was also verified by using the aerosol particle size distribution data of the aircraft. The inversion results showed that the algorithm was reliable in retrieving the aerosol particle size distributions at vertical heights using lidar data.

Telescopes in polarization lidar often modify the input polarization of the return signal, such that the telescope may significantly impact the depolarization estimates of aerosol and introduce error to the polarization lidar measurements. The error cannot be corrected by a traditional calibration constant. We present a method to correct the polarization effect of the telescope. We analyze the polarization effect of a telescope on the basis of the Mueller formalism, and introduce an algorithm for correcting the depolarization parameter of aerosol. A Newton telescope and a Cassegrain telescope are often chosen as the receiver in lidar. Their polarization models are established, and the Mueller matrices are calculated. The components of these matrices are dependent on wavelength, incident angle of the incoming light, and surface properties. The polarization impact of the telescope in lidar can be calibrated by a parameter, and the effects of different telescopes are discussed. The polarization crosstalk induced by the Newton telescope is obvious. The depolarization parameters change greatly with coating and wavelength, and they are calculated and presented. Whereas the crosstalk of a Cassegrain telescope is much smaller, the error can reach the level of ${10}^{-3}$ and can be negligible. The method presented in this paper could also be upgraded by taking into account all of the optical devices instead of only the telescope.

A macro-vertical structure is closely related to weather evolution and the energy budget balance of the atmospheric system of the Earth. In this study, radiosonde data were used to identify a cloud vertical structure (CVS) using the adjusted relative humidity threshold method. To evaluate the reliability and stability of this method, the results obtained based on the spatiotemporal matching criteria established in this study were compared with Ka-band millimetre-wave cloud radar (MMCR) observation data. This comparison showed that both devices exhibit high consistency in low-level cloud detection. With the increase in the cloud height, the frequency of the cloud appearance detection by the radiosonde became higher than that by the MMCR. In spring, the results of the CVS detection by the two devices were in good agreement. Specifically, the determination coefficients of the modified degrees of freedom (adjusted R-square) of the cloud base height (CBH) and cloud top height (CTH) detected by the two devices were 0.934 and 0.879, respectively. The horizontal drift of the radiosonde was the smallest in summer, and the adj. R-square values of the CBH and CTH were 0.814 and 0.852, respectively. The CVS observation results by the radiosonde and the MMCR were significantly different in autumn (the adj. R-Square values of the CBH and CTH were 0.715 and 0.629, respectively). In winter, the adj. R-Square values of the CBH and CTH observed by the radiosonde and the MMCR were 0.958 and 0.710, respectively. The statistics and analysis of the results of the distribution characteristics of the CVSs using radiosonde data from 2019 to 2021 from Xi’an showed that the average CTH and CBH were at 7–10 km and 3–5 km, respectively. The frequencies of the cloud absence, rainfall, and two- and three-layer clouds were the highest in the winter (34.36%), autumn (12.99%), and summer, respectively.

A complex optical system used in polarization lidars often modifies the input polarization of the return signal so that it may significantly impact depolarization estimates and introduce errors to polarization lidar measurements. In most cases, retardation, depolarization, and misalignment of the system exist at the same time and interact with each other. Polarization effects of the system cannot be represented by a simple correction coefficient, so they cannot be removed using a traditional calibration method. Detailed analysis and correction technologies were provided to remove systematic biases in estimating depolarization values from a polarization lidar owing to multiple optical components. The Mueller matrices from an emitter to a receiver were calculated, and the expression for an aerosol depolarization parameter including system polarization effects was derived and obtained. In addition, the correction algorithm based on the Mueller matrix was introduced and provided. A polarization lidar was established, and the polarization characteristics of its optical components were measured with a laboratory ellipsometer; then, the Mueller matrix of the receiver was calculated and obtained. Lidar observations were performed, and our correction algorithm was applied to lidar field data. The results show that the correction method can significantly remove systematic polarization effects.