Characteristics of Cloud and Aerosol Derived from Lidar Observations during Winter in Lhasa, Tibetan Plateau

: In order to investigate the variations of cloud and aerosol vertical profiles over the Tibetan Plateau (TP) in winter, we performed ground-based lidar observations in Lhasa, a city on the TP, from November 2021 to January 2022. The profiles of extinction coefficient, depolarization ratio, and signal-to-noise ratio (SNR) were retrieved using the atmospheric echo signals collected by the lidar. Clouds were identified by the range-correction echo signals and classified into water clouds, mixed clouds, horizontally oriented ice crystal clouds (HOICC), and ice clouds by the depolarization ratio and the hourly temperature from the European Centre for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5). The clouds mainly appeared at a height of 3~5 km from 14:00–22:00 Beijing Time throughout the field campaign. The height and frequency (~30%) for cloud appearance were significantly lower than that reported in previous studies in summer. The cloud categories were dominated by mixed clouds and ice clouds during the observation period. The proportions of ice clouds gradually increased with increasing heights. After eliminating profiles influenced by clouds, the aerosol extinction coefficient and depolarization ratio were obtained, and the atmospheric boundary layer height (ABLH) was calculated. The aerosol extinction coefficient decreased with increasing height in the ABLH, and there were no obvious changes for the aerosol extinction coefficient above the ABL. The aerosol extinction coefficients near the Earth’s surface presented two peaks, appearing in the morning and evening, respectively. The high aerosols at the surface in the morning continually spread upward for 4–5 h and finally reached an altitude of 1 km with the development of ABLH. In addition, the depolarization ratio of aerosols decreased slowly with increasing altitudes. There was no obvious diurnal variation for depolarization ratios, indicating partly that the source of aerosols did not change significantly. These results are beneficial in understanding the evolution of cloud and aerosol vertical profiles over the TP.


Introduction
Clouds and aerosols are important for many physical and chemical processes in the atmosphere.Clouds determine the surface precipitation and atmospheric radiation balance [1,2].Aerosol particles can be used as cloud condensation nuclei to change the microphysical properties of clouds and then affect the radiative forcing of clouds.They Remote Sens. 2024, 16, 2074 2 of 16 also affect the energy balance of the Earth-atmosphere system by scattering and absorbing short-wave solar radiation and emitting and capturing long-wave radiation [3][4][5].Tibetan Plateau (TP) is located in the central and eastern parts of the Eurasian continent and serves as a role of "the world water tower" [6].Due to its unique landform, fragile ecosystem, and particular monsoon circulation, TP is considered a sensitive area for human activities [7].The dust aerosol over the TP can heat the middle troposphere by absorbing solar radiation.The heating produces atmospheric dynamical feedback, i.e., the so-called elevated-heatpump (EHP) effect, which enhances the rate of snow melt in the Himalayas and TP [8].The increase of cloud cover on the lee slope of the TP increases the surface temperature, which leads to the destabilization and desaturation in the boundary layer [9].In addition, the increase of low-level cloud cover at night is greater than the decrease of low-level cloud cover and total cloud cover in daytime, which influences surface temperature and its range of diurnal variation over the TP [10].Investigating the spatiotemporal variations of cloud and aerosol on the TP can help us to understand the effects of their radiative feedback on the climate and environment in this region.
Lhasa is one of the major cities on the TP, with an average elevation of more than 3600 m [11].Due to urbanization and tourism development, fossil fuels have become the primary energy consumption in Lhasa [12].Primitive biomass burning and religious burning also occur frequently [13].There have been some studies focusing on the variation characteristics of clouds and aerosols over the TP.It was shown that the cloud cover in July over the TP accounted for a high proportion throughout the day due to the influence of monsoons, especially in the southern and northwestward areas of the TP, and the semidiurnal (2:00-10:00 UTC) cloud-cover cycle was weaker in July than in January, April, and October [14].The dependences of cloud cover on altitude over the TP were found to show an obvious seasonal variation, with a single peak (7~11 km) from January to April, double peaks (5~8 km and 11~17 km, respectively) after mid-June, and a single peak (5~10 km) after mid-August [15].According to the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) aerosol observation, the dust plume in spring over the TP is mainly from the desert and the Gobi region in the northeastern TP [16].In addition, the local religious burning in Lhasa in summer is a significant potential source of aerosols and accounts for 80% of total source emissions on the religious activity days [13].However, the vertical distribution of clouds and aerosols over the TP in winter are still not fully understood at present.
In recent years, lidar has been widely used for remote sensing detection of clouds and aerosols [17].Lidar has the advantages of high spatial and temporal resolution in obtaining profile information of atmospheric features.Based on polarization detection technology, the parallel and perpendicular polarization components of the received backscattered signal are separated and then the corresponding depolarization ratios are calculated [18].The depolarization ratio is used to distinguish the scattering of spherical particles from non-spherical particles and identify aerosols, cloud droplets, and ice crystals [18][19][20].The joint observations of lidar and weather radar in summer in Naqu, Tibet indicated that the frequency of low clouds in summer is about 50%, the frequency of liquid clouds and mixed phase clouds is higher in the morning and afternoon, and ice clouds mainly occur in the afternoon and midnight [2].The polarization lidar observations in Litang and Naqu over the TP showed that the average depolarization ratio is 0.44 for cirrus clouds, and the height of cirrus clouds is between 5.2-12 km above ground level.The depolarization ratio of cirrus particles increases with increasing height, but the situation for the particle depolarization ratio is opposite [21].The observations by satellite-borne lidar CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) from 2007 to 2014 found that the dust on the TP mainly came from the Taklimakan Desert, in addition to the Gurbantunggut Desert and the Great Indian Thar Desert in some cases [22].These studies have mainly concentrated on the characteristics of clouds and aerosol in summer over the TP through lidar observations.The lidar observations of clouds and aerosols over the TP are scarce in winter, and their spatiotemporal evolutions are not yet fully understood.
In this study, we performed ground-based lidar observations in Lhasa from 2 November 2021 to 25 January 2022.Clouds were identified by the range-corrected lidar signal and classified by the depolarization ratios and hourly temperature profiles from ERA5.The vertical distributions and diurnal variations of the extinction coefficient and depolarization ratio were also investigated.Section 2 introduces the observation site and instrument, as well as the methods for the extinction coefficient and depolarization ratio retrieval, cloud identification, and classification.Section 3 presents the vertical distribution and its temporal evolution of cloud, extinction coefficient, and depolarization ratio during the observation period.Finally, a discussion and conclusions are provided in Sections 4 and 5, respectively.

Site and Instrument
The lidar observations were performed at the Lhasa Meteorological Station (91.14 • E, 29.66 • N, altitude 3552.5 m, a.s.l) in the downtown of Lhasa, the capital of the Tibet Autonomous Region, from 2 November 2021 to 25 January 2022.Within a radius of 3 km around the site, there are mainly low buildings, including several temples (such as Jokhang Temple and Norbulingka Temple), where spices and biomass fuels are burned for religious activities.In addition, the world-famous Potala Palace is located about 1.5 km southwest of the site.The site belongs to the temperate monsoon semi-arid plateau climate, characterized by strong solar radiation, long sunny days, low rainfall, and cold and dry winters.The Barkhor Street ambient air quality monitoring station is about 500 m away, southwest of the site, where hourly concentrations of surface particulate matter 2.5 (PM 2.5 ) are operationally released in real time.More details of the natural and human environment around the observation site can be found in previous studies [12,23,24].
The lidar system used in the study was developed by the Anhui Institute of Optics and Fine Mechanics, Chinese Academy of Sciences (AIOFM, CAS).The lidar consists of a laser emitting system, a receiving system, a data acquisition system, and a control system [25].The laser emitting system consists of a laser and an optical unit, which emits laser pulses at a wavelength of 532 nm into the atmosphere.The receiving system consists of a telescope and a subsequent optical system.The telescope belongs to a Cassegrain telescope with a field of view angle of less than 1 mrad and a diameter of 150 mm.The atmospheric backscattered light received by the telescope is split by the polarizing prism, filtered on each channel, and then detected by photomultiplier tubes (PMT).The echo signals can be obtained via the data acquisition system.The main technical specifications of the lidar system are given in Table 1.During the 85-day observation period, the lidar operated stably except for maintenance and power outages, and the data acquisition rate was approximately 95%.The original vertical resolution of the lidar was set to 7.5 m, with a temporal resolution of 5 min.During data pre-processing, the original signals were smoothed using a 7-point window, and the background noises were removed.In addition, data quality control was carried out based on the principle of 3 times standard deviation.

Retrieval of Extinction Coefficient and Depolarization Ratio
With the assumption of only single scattering, the lidar equation can be expressed as the following [26]: where P(z, λ) is the power of the echo signal at distance z; C is the lidar system constant; P 0 (λ) is the initial laser emission power.The atmospheric backscatter coefficient β(z, λ) is the sum of the atmospheric molecule backscatter coefficient β m (z, λ) and the aerosol backscatter coefficient β a (z, λ).The atmospheric extinction coefficient α(z, λ) is the sum of the atmospheric molecule extinction coefficient α m (z, λ) and the aerosol extinction coefficient α a (z, λ).
The aerosol extinction coefficients can be retrieved by the Collis method [27], the Klett method [28], and the Fernald method [29]; wherein the Fernald method is widely used with high inversion accuracy.In this study, we used the Fernald method to retrieve the extinction coefficient of aerosol.The Fernald method assumes that the backscattering coefficient of particles is proportional to the extinction coefficient, as shown in Equations ( 2) and (3): Therefore, the aerosol backscatter coefficient can be written as follows [29]: where z c is the height of reference point, and S a is the ratio of extinction coefficient and backscatter coefficient for atmospheric aerosol (usually referred to as the lidar ratio).The lidar ratio is one of the important parameters for the retrieval of the aerosol extinction coefficient, which is related to the physical and chemical properties of the aerosol, such as composition, particle size distribution, and shape [30,31].S m = 8π/3 is the ratio of extinction coefficient and backscatter coefficient for atmospheric molecules.β m is calculated from the vertical profile of air molecule density provided by the American standard atmospheric model according to the Rayleigh scattering theory [32], where z c is the height of the reference point.β a (z c ) at wavelength 532 nm is calculated by Equation ( 5) [33,34]: The depolarization ratio is defined as the ratio of the cross-polarized signal to the parallelly polarized signal at the same wavelength, which is usually used as the indicator of particle shape.Commonly, a larger depolarization ratio corresponds to a more irregular shape of the particles [35].In this study, the depolarization ratio was calculated by the 532 nm echo signal, as shown in Equation ( 6): where δ(z) is the depolarization ratio; β 532⊥ (z) and β 532∥ (z) are the backscatter coefficients at the wavelength of 532 nm for cross and parallel polarization channels, respectively.
Figure A1 shows the range-corrected echo signal, extinction coefficient and depolarization ratio profile on 10 November 2021, respectively.

Cloud Identification and Classification
Based on the range-corrected lidar signal, we use the method of value distribution equalization (VDE) to identify clouds [36].The VDE method has several advantages, such as preventing the expansion of noises at far distances, making aerosol/cloud signals stand out clearly from noises, and making SNR at a far distance comparable to that at a near distance [36].
In the first part, the semi-discretization method (SDP) is used to reduce the influence of noise on the power of the echo signal P [37].The background noise (signals at 14-15 km height above the surface in this study) is subtracted from P and the standard deviation (σ n ) of the background noise is calculated.Then smoothing P with a window of 0.06 km obtains the processed echo signal P s (z).The processing of P s (z i )(i = 1, 2, 3, . . ., N) is shown in Equations ( 7)-( 9): where P d (z) denotes the average of P d1 (z) and P d2 (z), which retains the characteristics of sharp signals and minimizes the impact of atmospheric turbulence.
The second part is to treat P d (z) using the VDE method to identify clouds and aerosol layers.First, P d (z) is sorted in ascending order to obtain the intermediate vari- able R(i), (i = 1, 2, 3, . . ., N).The index of an element of R(i) in P d (z) is denoted as I(i), (i = 1, 2, 3, . . ., N).The minimum and maximum values of P d (z) are expressed as M I and M A , respectively.Then, the scale mapping coefficients P e (i) are shown as follows: with P e (i) and the signal ranges, new ascending data values can be calculated: Finally, a new dataset P N (z) can be constructed via the aforementioned steps: The baseline B(z) is constructed from points (z 1 , M A ) and (z N , M I ).The position where P N (z) is greater than B(z) is regarded as the appearance of the cloud and aerosol layer.We removed layers with a thickness of less than 50 m, probably caused by random fluctuations [36].
The third part distinguishes the cloud layer from the aerosol layer by calculating the slope of the rising part and the falling part of the aforementioned identification layer.The following threshold function F(z) should be constructed: where P b (z) is the echo signal after subtracting the background noise.The maximum and minimum values of F(z) are denoted as K and D, respectively.The aerosol layer and cloud layer can be distinguished by Equation ( 15) [36,38,39]: In addition to the clouds identified by the above method, the cloud phase can be further classified by using depolarization ratios and the ERA5 hourly temperature profile [40].The temperature profile is gridded into 37 layers from 1000 hPa to 1 hPa and interpolated to match the temporal and vertical resolution of the lidar profile.According to the temperature and depolarization ratio data, the cloud phases are divided into four categories, i.e., water cloud, mixed phase cloud, horizontally oriented ice crystal cloud (HOICC), and ice cloud [2,41].The classification basis is shown in Table 2: Table 2. Classification basis of cloud phase.

Cloud Phase
where β a (z) is the aerosol backscattering coefficient, T is the temperature by ERA5, and δ(z) is the depolarization ratio.

Characteristics of Cloud Vertical Distribution
According to the classification algorithm described in Section 2.3, the vertical distribution of clouds in Lhasa during the effective observation period was obtained.Because the strong attenuation of clouds can cause a sharp drop in the SNR, we eliminated data with an SNR below 10.In addition, the receiving field of view of the lidar is less than 1 mrad, which can effectively reduce the receiving of multiple scattered signals.Figure 1 shows the average diurnal variations of cloud appearance frequency in different months.We marked the locations identified as clouds in the observation data matrix as 1 and the other locations as 0 to obtain the marker matrix.The sum of the daily marker matrix is divided by the number of days to obtain the average diurnal variations of cloud appearance frequency.In November 2021, the height of the cloud was concentrated at 3-6 km, and the time of appearance was concentrated at 2:00~10:00 Beijing Time (BJT, Universal Time Coordinated + 8 h) and 16:00~22:00 BJT, with the highest appearance frequency of ~23%.In December 2021, the height of cloud appearance was concentrated at the height of 1.5~3 km during the period of 6:00~16:00, but it was at 3 km~6 km during time segments of 4:00~8:00 BJT and 16:00~22:00 BJT, with the highest appearance frequency of ~29%.In January 2022, the height of the cloud appearance was concentrated between 2~7 km during the time period of 16:00~22:00 BJT.During the observation period, the maximum of cloud appearance frequency was 41% in January 2022.The clouds occurred almost all day in December 2021, presenting a different diurnal pattern of cloud appearance frequency from that in the other two months.Multi-layer clouds appeared in the early morning and afternoon in December 2021.
January 2022, the height of the cloud appearance was concentrated between 2~7 km during the time period of 16:00~22:00 BJT.During the observation period, the maximum of cloud appearance frequency was 41% in January 2022.The clouds occurred almost all day in December 2021, presenting a different diurnal pattern of cloud appearance frequency from that in the other two months.Multi-layer clouds appeared in the early morning and afternoon in December 2021.On average, almost no clouds appeared around noon during the field campaign (Figure 2a).The heights of the cloud layer were primarily in the range of 2~7 km, although On average, almost no clouds appeared around noon during the field campaign (Figure 2a).The heights of the cloud layer were primarily in the range of 2~7 km, although there were clouds below 2 km.The frequency of cloud appearance between 14:00~22:00 was significantly higher than that at the other time of the day, which is similar to the results reported in previous research studies [2,34].The peak frequency (~29%) of cloud appearance was significantly lower than that at the same region in summer [2,32].This is probably related to the cold and dry climate conditions in winter.The water clouds were not identified, probably due to the fact that the temperature was below zero at the height of the cloud appearance in winter of Lhasa.HOICC appeared below 5 km, with an occurrence probability below 1.5% during the observation period (Figure 2b).The proportion of HOICC in ice clouds was close to the global research result reported by Noel et al. [41].The occurrence probability of mixed cloud was less than 10% in the height range of 1~6 km, mainly around 20:00 BJT (Figure 2c).It should be mentioned that the mixed clouds probably included the supercooled water clouds due to limitations in cloud classification methods.The ice clouds mainly occurred from 16:00 BJT to 22:00 BJT, with maximum frequency at a height of 3-5 km (Figure 2d).The occurrence time of ice clouds is similar to that observed in summer in the Naqu region [2], which is about 200 km away from Lhasa.However, the mean height for the maximum frequency of ice cloud appearance is lower by 3 km in winter in Lhasa, probably due to the seasonal difference in surface radiation.
methods.The ice clouds mainly occurred from 16:00 BJT to 22:00 BJT, with maximum frequency at a height of 3-5 km (Figure 2d).The occurrence time of ice clouds is similar to that observed in summer in the Naqu region [2], which is about 200 km away from Lhasa.However, the mean height for the maximum frequency of ice cloud appearance is lower by 3 km in winter in Lhasa, probably due to the seasonal difference in surface radiation.Figure 3a shows the frequencies of cloud phases at different heights in three months.The frequencies for each cloud phase first increased and then decreased with the increasing heights, and the peak frequency appeared at 5 km.The frequencies of clouds also increased significantly from November 2021 to January 2022.Figure 3b shows the proportions of different cloud phases in the total cloud at different heights in three months.The proportion of mixed clouds showed a decreasing trend with increasing height from November 2021 to January 2022, while the situation was the opposite for the proportion of Figure 3a shows the frequencies of cloud phases at different heights in three months.The frequencies for each cloud phase first increased and then decreased with the increasing heights, and the peak frequency appeared at 5 km.The frequencies of clouds also increased significantly from November 2021 to January 2022.Figure 3b shows the proportions of different cloud phases in the total cloud at different heights in three months.The proportion of mixed clouds showed a decreasing trend with increasing height from November 2021 to January 2022, while the situation was the opposite for the proportion of ice clouds.On average, the maximum of total cloud frequency was ~11% at the 5 km height, dominated by mixed clouds and ice clouds (Figure 3c).At the 1 km height, mixed clouds accounted for 90%, with the remaining 10% for ice clouds and HOICC (Figure 3d).In addition, the maximum proportion (~70%) of ice clouds appeared at 10 km during the field campaign.Only mixed clouds and ice clouds were found above a height of 6 km.
ice clouds.On average, the maximum of total cloud frequency was ~11% at the 5 km height, dominated by mixed clouds and ice clouds (Figure 3c).At the 1 km height, mixed clouds accounted for 90%, with the remaining 10% for ice clouds and HOICC (Figure 3d).In addition, the maximum proportion (~70%) of ice clouds appeared at 10 km during the field campaign.Only mixed clouds and ice clouds were found above a height of 6 km.

Monthly and Diurnal Variations of Vertical Distribution of Aerosol Optical Properties
The extinction coefficient and depolarization ratio are affected by different physicochemical properties and aerosol sources [5].The relationship between aerosol extinction coefficient and PM2.5 concentration can be simply expressed as a pm 2.5 KC  =  [42].The specific coefficient K is related to aerosol size distribution, refractive index, and atmospheric relative humidity.When the relative humidity is below 70%, the aerosol hydrophilic growth could be negligible.So, the specific coefficient K could be considered as constant under the condition of less than 70% relative humidity in the planetary boundary layer [42].In order to verify the reliability of lidar observation in Lhasa, we performed a linear regression analysis between the extinction coefficient in the near-surface layer (i.e., 200-250 m in this study) and the PM2.5 mass concentration measured at the adjacent station during the field campaign.As shown in Figure 4, the correlation coefficient is 0.77 between them, partly implying that the lidar observations are reliable.

Monthly and Diurnal Variations of Vertical Distribution of Aerosol Optical Properties
The extinction coefficient and depolarization ratio are affected by different physicochemical properties and aerosol sources [5].The relationship between aerosol extinction coefficient and PM 2.5 concentration can be simply expressed as α a = K × C pm2.5 [42].The specific coefficient K is related to aerosol size distribution, refractive index, and atmospheric relative humidity.When the relative humidity is below 70%, the aerosol hydrophilic growth could be negligible.So, the specific coefficient K could be considered as constant under the condition of less than 70% relative humidity in the planetary boundary layer [42].In order to verify the reliability of lidar observation in Lhasa, we performed a linear regression analysis between the extinction coefficient in the near-surface layer (i.e., 200-250 m in this study) and the PM 2.5 mass concentration measured at the adjacent station during the field campaign.As shown in Figure 4, the correlation coefficient is 0.77 between them, partly implying that the lidar observations are reliable.Given the low SNR of the echo signal in the vertical channel above 3 km, we mainly analyzed the vertical evolutions of aerosol optical parameters below 3 km.The data of extinction coefficient and depolarization ratio were rejected once they were affected by the clouds or outside three times the standard deviation in the following statistical analysis.Figure 5 shows the vertical profiles of the extinction coefficient and depolarization ratio in each month as well as during the whole observation period.Both the extinction coefficient and depolarization ratio decreased with increasing altitude.For example, the extinction coefficients were 0.059 ± 0.03 in November 2021, 0.067 ± 0.03 in December 2021, and 0.065 ± 0.08 in January 2022 at 200 m, respectively.The corresponding extinction coefficients at 1 km were 0.043 ± 0.004, 0.047 ± 0.01 and 0.048 ± 0.03, respectively.This indicates that the rather low aerosol concentrations in Lhasa decreased with increasing altitude.In addition, there were no significant monthly trends for the extinction coefficient and depolarization ratio, indicating no significant differences in particle shapes and aerosol sources.Given the low SNR of the echo signal in the vertical channel above 3 km, we mainly analyzed the vertical evolutions of aerosol optical parameters below 3 km.The data of extinction coefficient and depolarization ratio were rejected once they were affected by the clouds or outside three times the standard deviation in the following statistical analysis.Figure 5 shows the vertical profiles of the extinction coefficient and depolarization ratio in each month as well as during the whole observation period.Both the extinction coefficient and depolarization ratio decreased with increasing altitude.For example, the extinction coefficients were 0.059 ± 0.03 in November 2021, 0.067 ± 0.03 in December 2021, and 0.065 ± 0.08 in January 2022 at 200 m, respectively.The corresponding extinction coefficients at 1 km were 0.043 ± 0.004, 0.047 ± 0.01 and 0.048 ± 0.03, respectively.This indicates that the rather low aerosol concentrations in Lhasa decreased with increasing altitude.In addition, there were no significant monthly trends for the extinction coefficient and depolarization ratio, indicating no significant differences in particle shapes and aerosol sources.Given the low SNR of the echo signal in the vertical channel above 3 km, we mainly analyzed the vertical evolutions of aerosol optical parameters below 3 km.The data of extinction coefficient and depolarization ratio were rejected once they were affected by the clouds or outside three times the standard deviation in the following statistical analysis.Figure 5 shows the vertical profiles of the extinction coefficient and depolarization ratio in each month as well as during the whole observation period.Both the extinction coefficient and depolarization ratio decreased with increasing altitude.For example, the extinction coefficients were 0.059 ± 0.03 in November 2021, 0.067 ± 0.03 in December 2021, and 0.065 ± 0.08 in January 2022 at 200 m, respectively.The corresponding extinction coefficients at 1 km were 0.043 ± 0.004, 0.047 ± 0.01 and 0.048 ± 0.03, respectively.This indicates that the rather low aerosol concentrations in Lhasa decreased with increasing altitude.In addition, there were no significant monthly trends for the extinction coefficient and depolarization ratio, indicating no significant differences in particle shapes and aerosol sources.

November
December January Total  Figure 6 shows the average diurnal variations of the aerosol extinction coefficient and the depolarization ratio in the lower layer (i.e., between 200-400 m) in each month and during the whole observation period.The diurnal variations of the extinction coefficient presented a "double-peak" pattern, around 12:00 BJT and 20:00 BJT, respectively.The averaged diurnal variation of ground PM 2.5 concentration measured at the adjacent monitoring station also presented double peaks, which appeared around 10:00 BJT and 20:00 BJT, respectively (Figure A2).The diurnal cycles were partly influenced by the local source emissions, such as production activities and religious activities [43].The time differences for the first peak between aerosol extinction and PM 2.5 concentrations were partly related to the particle transport from the lower layer to the higher layer with the development of the boundary layer.The change in the depolarization ratio and extinction coefficient in November 2021 (Figure 6a) and December 2021 (Figure 6b) is consistent, but this consistency is not obvious in January 2022 (Figure 6c).Overall, the diurnal patterns of depolarization ratio were similar to that of aerosol extinction coefficient.Considering the depolarization ratio variation in the range of 0.06-0.1, the rather small fluctuation indicates that the type of aerosol particles did not change significantly during the observation period.
Remote Sens. 2024, 16, x FOR PEER REVIEW 11 of 16 Figure 6 shows the average diurnal variations of the aerosol extinction coefficient and the depolarization ratio in the lower layer (i.e., between 200-400 m) in each month and during the whole observation period.The diurnal variations of the extinction coefficient presented a "double-peak" pattern, around 12:00 BJT and 20:00 BJT, respectively.The averaged diurnal variation of ground PM2.5 concentration measured at the adjacent monitoring station also presented double peaks, which appeared around 10:00 BJT and 20:00 BJT, respectively (Figure A2).The diurnal cycles were partly influenced by the local source emissions, such as production activities and religious activities [43].The time differences for the first peak between aerosol extinction and PM2.5 concentrations were partly related to the particle transport from the lower layer to the higher layer with the development of the boundary layer.The change in the depolarization ratio and extinction coefficient in November 2021 (Figure 6a) and December 2021 (Figure 6b) is consistent, but this consistency is not obvious in January 2022 (Figure 6c).Overall, the diurnal patterns of depolarization ratio were similar to that of aerosol extinction coefficient.Considering the depolarization ratio variation in the range of 0.06-0.1, the rather small fluctuation indicates that the type of aerosol particles did not change significantly during the observation period.The atmospheric boundary layer height (ABLH) was also retrieved from lidar echo signals by the gradient-based method [44][45][46].The aerosol lidar-derived ABLH is actually the height of the surface aerosol layer [47].Assuming that the distribution of aerosol concentration during the daytime is primarily influenced by turbulent mixing, the boundary layer height measured by lidar is considered as the mixed layer height (MLH), which is equivalent to the convective boundary layer height obtained from radiosonde observations [47,48].For the stable nocturnal boundary layer (NBL), the height determined by aerosol lidar is either the top of the residual layer or the top of the surface mechanical driven ML (MLH, or the height of ground inversion layer) [47].Figure 7 shows the diurnal variations of ABLH overlaid onto that of the extinction coefficient profile.It is clearly shown that the aerosols were basically distributed within the atmospheric boundary layer (ABL).The ABLH gradually increased from hundreds of meters at ~10:00 BJT to 1 km at ~15:00 BJT.Such ABL development driven by the increase of surface thermal radiation The atmospheric boundary layer height (ABLH) was also retrieved from lidar echo signals by the gradient-based method [44][45][46].The aerosol lidar-derived ABLH is actually the height of the surface aerosol layer [47].Assuming that the distribution of aerosol concentration during the daytime is primarily influenced by turbulent mixing, the boundary layer height measured by lidar is considered as the mixed layer height (MLH), which is equivalent to the convective boundary layer height obtained from radiosonde observations [47,48].For the stable nocturnal boundary layer (NBL), the height determined by aerosol lidar is either the top of the residual layer or the top of the surface mechanical driven ML (MLH, or the height of ground inversion layer) [47].Figure 7 shows the diurnal variations of ABLH overlaid onto that of the extinction coefficient profile.It is clearly shown that the aerosols were basically distributed within the atmospheric boundary layer (ABL).The ABLH gradually increased from hundreds of meters at ~10:00 BJT to 1 km at ~15:00 BJT.Such ABL development driven by the increase of surface thermal radiation during the day led to the vertical diffusion of aerosols and played an important role in the dilution effect of air pollution.For another air pollution period (around 20:00 BJT), the ABLHs were about 0.6 km and remained nearly at the same level for several hours, which were not conducive to the diffusion of air pollution.In addition, the ABLH diurnal variation pattern did not change significantly over the months during the field campaign.during the day led to the vertical diffusion of aerosols and played an important role in the dilution effect of air pollution.For another air pollution period (around 20:00 BJT), the ABLHs were about 0.6 km and remained nearly at the same level for several hours, which were not conducive to the diffusion of air pollution.In addition, the ABLH diurnal variation pattern did not change significantly over the months during the field campaign.

Discussion
Here we report for the first time long continuous lidar observations in Lhasa during winter.The lidar observation shows that the cloud height in winter of Lhasa is significantly lower than that in summer, which is partly related to the weaker thermal force caused by the weakening surface radiation in winter [2].The clouds always appear except around noon, which is probably related to the local boundary layer processes [49].The frequency of cloud appearance is also significantly lower in winter than that in summer, which is probably connected with the colder and drier climate conditions in winter.In addition, the clouds in winter of Lhasa are dominated by ice clouds and mixed clouds.Due to the limitations of the cloud classification method and uncertainties of temperature data, almost no water clouds were identified.Probably, there were super-cooled water clouds that were not distinguishable.In addition, the uncertainty of the inversion results caused by multiple scattering and the strong extinction effect of clouds needs further study.
The consistency between aerosol extinction coefficient and PM2.5 concentrations partially indicates that the lidar observation was reliable during the field campaign.Overall, the extinction coefficients decrease with increasing altitude and the depolarization ratios vary over a small range, indicating that aerosol particles are mainly concentrated in the lower layer and the aerosol type does not change significantly.Thus, the primary sources of aerosol in Lhasa are dominated by local emissions [50].The clear bimodal structure of aerosol diurnal variation reveals that we should pay more attention to the reduction of local source emissions in the two key periods.In addition, the diurnal variation of the aerosol vertical profile presents a clear diffusion process of aerosols with increasing ABLH, which implies that the ABL evolutions play an important role in air pollution in Lhasa.

Discussion
Here we report for the first time long continuous lidar observations in Lhasa during winter.The lidar observation shows that the cloud height in winter of Lhasa is significantly lower than that in summer, which is partly related to the weaker thermal force caused by the weakening surface radiation in winter [2].The clouds always appear except around noon, which is probably related to the local boundary layer processes [49].The frequency of cloud appearance is also significantly lower in winter than that in summer, which is probably connected with the colder and drier climate conditions in winter.In addition, the clouds in winter of Lhasa are dominated by ice clouds and mixed clouds.Due to the limitations of the cloud classification method and uncertainties of temperature data, almost no water clouds were identified.Probably, there were super-cooled water clouds that were not distinguishable.In addition, the uncertainty of the inversion results caused by multiple scattering and the strong extinction effect of clouds needs further study.
The consistency between aerosol extinction coefficient and PM 2.5 concentrations partially indicates that the lidar observation was reliable during the field campaign.Overall, the extinction coefficients decrease with increasing altitude and the depolarization ratios vary over a small range, indicating that aerosol particles are mainly concentrated in the lower layer and the aerosol type does not change significantly.Thus, the primary sources of aerosol in Lhasa are dominated by local emissions [50].The clear bimodal structure of aerosol diurnal variation reveals that we should pay more attention to the reduction of local source emissions in the two key periods.In addition, the diurnal variation of the aerosol vertical profile presents a clear diffusion process of aerosols with increasing ABLH, which implies that the ABL evolutions play an important role in air pollution in Lhasa.

Conclusions
Lidar observations were performed from November 2021 to January 2022 in Lhasa, one of the major cities on the TP.We identified clouds based on the distance-corrected echo signal of lidar and classified clouds into water clouds, HOICC, mixed clouds, and ice clouds.
Then, we investigated the temporal variations of the vertical distribution for each cloud phase.We found that the distribution height and frequency of clouds were significantly different from those reported in summer.We also presented temporal variations of the vertical distribution of aerosol optical properties as well as ABLH retrieved from lidar observation.The main conclusions are as follows: 1.
Clouds always appeared except around noon during the field campaign, concentrated at a height of about 5 km.The height and frequency of cloud appearance in this study were significantly lower than that reported in summer.The cloud types are dominated by mixed clouds and ice clouds in the winter of Lhasa.The relative proportions of mixed clouds gradually decrease with increasing altitude, but the situation is opposite for ice clouds.

2.
The aerosol extinction coefficients decrease with increasing altitude.The depolarization ratios vary in the range of 0.06-0.1,implying that the major sources of aerosol particles did not change significantly during the observation period.There were two peaks for aerosol diurnal variation, appearing in the morning and late afternoon, respectively.The aerosols were basically distributed within the ABL.With the ABL development in the morning, the aerosol particles continued to diffuse upward for about 4 h until reaching a height of 1 km at ~15:00 BJT.For another air pollution period (around 20:00 BJT), the ABLHs were about 0.6 km and remained basically at the same level, which is not conducive to the diffusion of air pollution.In addition, the ABLH diurnal variation pattern changed slightly over the months during the field campaign.
As a whole, the results above enhance the understanding the characteristics of clouds and aerosols in winter over the TP.In the future, the comparisons made by different observation methods will be significant.Also, it is necessary to undertake long-term continuous observations by lidar in different regions of the TP in order to increase the spatial-temporal representativeness.

Figure 2 .
Diurnal variations of the vertical distribution of cloud frequencies during the field campaign for (a) all clouds, (b) HOICCs, (c) mixed clouds, and (d) ice clouds (BJT).

Figure 2 .
Figure 2. Diurnal variations of the vertical distribution of cloud frequencies during the field campaign for (a) all clouds, (b) HOICCs, (c) mixed clouds, and (d) ice clouds (BJT).

Figure 3 .
Figure 3. (a) The frequency of cloud phases at different heights in three months; (b) the proportion of different cloud phases in the total cloud at different heights in three months; (c)the mean frequency of cloud phases at different heights during the observation period; (d) the averaged proportion of different cloud phases in the total cloud at different heights during the observation period (HOICC: cyan; mixed cloud: light blue; ice cloud: blue).

Figure 3 .
Figure 3. (a) The frequency of cloud phases at different heights in three months; (b) the proportion of different cloud phases in the total cloud at different heights in three months; (c) the mean frequency of cloud phases at different heights during the observation period; (d) the averaged proportion of different cloud phases in the total cloud at different heights during the observation period (HOICC: cyan; mixed cloud: light blue; ice cloud: blue).

Figure 4 .
Figure 4. Linear regression analysis between the hourly extinction coefficient at the height range of 200~250 m and corresponding PM2.5 mass concentration measured at the adjacent station during the field campaign.

Figure 5 .
Figure 5. Mean vertical profiles of extinction coefficient (upper panel) and depolarization ratio (lower panel) (a,e) in November 2021, (b,f) in December 2021, (c,g) in January 2022, and (d,h) during the whole observation period (BJT).The error bars indicate the standard deviations.

Figure 4 .
Figure 4. Linear regression analysis between the hourly extinction coefficient at the height range of 200~250 m and corresponding PM 2.5 mass concentration measured at the adjacent station during the field campaign.

Figure 4 .
Figure 4. Linear regression analysis between the hourly extinction coefficient at the height range of 200~250 m and corresponding PM2.5 mass concentration measured at the adjacent station during the field campaign.

Figure 5 .
Figure 5. Mean vertical profiles of extinction coefficient (upper panel) and depolarization ratio (lower panel) (a,e) in November 2021, (b,f) in December 2021, (c,g) in January 2022, and (d,h) during the whole observation period (BJT).The error bars indicate the standard deviations.

Figure 5 .
Figure 5. Mean vertical profiles of extinction coefficient (upper panel) and depolarization ratio (lower panel) (a,e) in November 2021, (b,f) in December 2021, (c,g) in January 2022, and (d,h) during the whole observation period (BJT).The error bars indicate the standard deviations.

Table 1 .
Main technical specifications of lidar system.