Long-term variation in aerosol lidar ratio in Shanghai based on Raman lidar measurements

Accurate lidar ratio (LR) and better understanding of its variation characteristics can not only improve the retrieval accuracy of parameters from elastic lidar, but also play an important role in assessing the impacts of aerosols on climate. Using the observational data of a Raman lidar in Shanghai from 2017 to 2019, LRs at 355 nm were retrieved and their variations and influence factors were analyzed. Within the height range of 0.5–5 km, about 90 % of the LRs were distributed in 10–80 sr with an average value of 41.0± 22.5 sr, and the LR decreased with the increase in height. The volume depolarization ratio (δ) was positively correlated with LR, and it also decreased with the increase in height, indicating that the vertical distribution of particle shape was one of the influence factors of the variations in LR with height. LR had a strong dependence on the original source of air masses. Affected by the aerosols transported from the northwest, the average LR was the largest, 44.2± 24.7 sr, accompanied by the most irregular particle shape. The vertical distribution of LR was affected by atmospheric turbidity, with the greater gradient of LR under clean conditions. The LR above 1 km could be more than 80 sr, when Shanghai was affected by biomass burning aerosols.


Introduction
Aerosols in the atmosphere can affect the earth's climate by absorbing and scattering solar radiation (direct effect of aerosols) (Huang et al., 2014;Wang et al., 2013) or acting as cloud condensation nuclei, which can affect cloud physical properties and precipitation (indirect effect of aerosols) (Huang et al., 2006;Liu et al., 2019a, b;Yan and Wang, 2020). In general, the vertical distribution information of aerosols is required to improve our understanding of aerosol climate effects (Ferrare et al., 2001;Sicard et al., 2011). For example, Wang et al. (2020b) found that dust-forced radiative heating decreased significantly as the transport height of dust aerosols decreased. A study by Lu et al. (2020) showed that anomalous elevated aerosol layers above 2 km led to warming in the upper atmosphere (+0.32 K/d) and cooling on the surface (−0.04 K/d). In addition, the results of Liu et al. (2019b) showed that the aerosols mixing with the cloud layer in the vertical direction can significantly reduce the effective radius of ice particles and prolong the life of clouds.
As an active remote-sensing instrument, the elastic scattering lidar can obtain vertical distribution information of aerosols; however, it is necessary to assume an aerosol extinction-to-backscattering ratio (i.e., lidar ratio, LR) in the retrieval process (Fernald, 1984;Welton et al., 2001), which can result in significant errors for the extinction coefficient, followed by aerosol optical depth (AOD). To our knowledge, the LRs at 355, 532, and 1064 nm are usually assumed to be 50 sr in China (Fan et al., 2018; Published by Copernicus Publications on behalf of the European Geosciences Union.

5378
T. Liu et al.: Long-term variation in aerosol lidar ratio in Shanghai 2015;Lv et al., 2020;Ma et al., 2019). In addition, CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) onboard CALIPSO (Cloud-Aerosol Lidar Infrared Pathfinder Satellite Observations) can observe the vertical distribution of global aerosol optical properties, and its data products have been widely used around the world (Kim et al., 2018). The CALIOP algorithm first determines the type of aerosol according to the aerosol classification algorithm and then uses lookup table of multiple types of aerosol to determine the LR (Kim et al., 2018;Omar et al., 2009). Therefore, the quality of the CALIOP aerosol products depends on the accuracy of aerosol type identification and the consistency between actual LR and that in the lookup table (Painemal et al., 2019). LR is a complicated function of time and space, which depends on aerosol size distribution and particle composition (Reagan et al., 1988). LR is also affected by meteorological elements (Salemink et al., 1984), such as relative humidity (RH), which can change the aerosol particle size distribution and refractive index (Young et al., 1993). Therefore, good knowledge of accurate LR and its variation characteristics can not only improve the retrieval accuracy of parameters from elastic lidar, but also obtain information on aerosol types to trace the source of pollutants (Franke et al., 2001).
According to the definition of LR in Müller's (2003) study, . (1) LR was negatively correlated with phase function at 180 • (P (π )) and single-scattering albedo (ω). P (π ) is related to the sphericity of particles, and the sphericity information of particles can be obtained from polarization lidar. ω is indicative of aerosol absorption properties. The aerosol absorption properties and their vertical distributions play a crucial role in evaluating direct radiative forcing of aerosols (McComiskey et al., 2008;Zarzycki and Bond, 2010). Moreover, absorbing aerosols increase atmospheric stability by reducing the solar radiation reaching the surface during the day. By contrast, absorbing aerosols on the near-surface layer heat the surface and increase the atmospheric instability during night (Jacobson, 1998;Jacobson and Kaufman, 2006). The vertical profile of LR can reflect vertical variation in aerosol absorption properties, which can provide a basis for studying aerosol radiative forcing and causes of pollution (Mishchenko et al., 2004). LR can be obtained by a variety of methods, such as highspectrum-resolution lidar (HSRL), Raman lidar, and joint retrieval using sun photometer and elastic lidar (Zhao et al., 2018). Raman lidar can independently retrieve the extinction coefficient and backscatter coefficient of aerosols and obtain LR by combining elastic backscatter and Raman backscatter signals (Ansmann et al., 1992), which is the most widely used independent measurement method at present. Moreover, the LR measured by Raman lidar is a useful index to study the variations in aerosol physical properties (Ferrare et al., 2001).
A large number of observations and analysis of LR have been carried out based on Raman lidar all over the world. Since the establishment of the European Aerosol Research Lidar Network (EARLINET) in 2000, long time series observation data of vertical distribution and LR for various types of aerosol have been obtained on the European continent (Müller et al., 2007;Wandinger et al., 2016). In South Korea and Japan, the LR of Asian dust and biomass burning aerosols has also been studied based on Raman lidar (Murayama et al., 2004;Noh et al., 2007Noh et al., , 2008. The LR observed around the world usually shows different values due to different types of aerosols. However, long-term observations and research of LR in China are limited  due to the limitation of observation instruments. In particular, the observations and studies of LR are still rare in east China; however, range-resolved LR profiles based on independent measurement on a regional scale are very important. On one hand, the range-resolved LR obtained from groundbased Raman lidar can not only be used for comparison with 355 nm LR obtained from ATLID (Atmospheric LIDar) on EarthCARE (Earth Clouds and Radiation Explorer) planned to be launched by ESA (European Space Agency) (Liu et al., 2020a;Nicolae et al., 2018), but can also provide a reliable basis for the inversion hypothesis of elastic lidar in Shanghai and surrounding areas and improve product reliability for elastic lidar networks such as the Asian dust and aerosol lidar observation network. On the other hand, vertical distribution of aerosol absorption properties reflected by LR can be used as an input parameter for regional climate models (Mehta et al., 2018), which can further improve the calculation accuracy of radiative forcing.
In addition, studying the influence factors of LR in Shanghai can be conducive to understanding the LR variation characteristics and determining the source of pollutants. With these motivations, the vertical and temporal variations in LR and its influence factors in Shanghai were analyzed using the results retrieved from Raman lidar, which laid a solid foundation for the quantitative study of pollution and its causes in the future. The Raman depolarization lidar (Raymetrics S.A., Athens, Greece, model LR231-D300) used in this study is deployed on the roof of a building (31.1916 • N, 121.4333 • E; 67 m above sea level, a.s.l.) in downtown Shanghai. The site is 35 km from the East China Sea coastline, surrounded by populated residential and commercial areas. The laser of the lidar system is an Nd: YAG pulse laser equipped with a water cooling device. The system can emit 355, 532, and 1064 nm laser pulses with a pulse width of 5.4 ns and a repetition rate of 20 Hz. The lidar detection system can receive 355, 532, and 1064 nm elastic scattering signals. The 532 nm elastic scattering channel has two polarization channels: 532P (parallel) and 532S (perpendicular). In addition, the lidar can also receive the vibrational Raman scattering signals of nitrogen (387 nm) and water vapor molecules (408 nm) with an incident wavelength of 355 nm. Therefore, 355 nm LRs were obtained and discussed in this study. The transient recorder can detect signals in two modes: analog and photon counting. The analog mode is suitable for detecting strong signals at low altitude, and the photon-counting mode is suitable for detecting weak signals at high altitude. In order to better combine the advantages of the two modes to retrieve aerosol optical properties, it is necessary to glue signals of the two modes together (Newsom et al., 2009;Walker et al., 2014). The range resolution of the raw signal is 7.5 m, and the temporal resolutions are 1 min (before 25 January 2019) and 2 min (after 25 January 2019). Unless otherwise stated, the time in the paper was UTC (universal time coordinated).

HYSPLIT-4
HYSPLIT-4 (Hybrid Single Particle Lagrangian Integrated Trajectory Model, Version 4) is a professional model jointly developed by the National Oceanic and Atmospheric Administration (NOAA) and the Australian Bureau of Meteorology for calculating and analyzing transport and diffusion trajectories of atmospheric pollutants and has been widely used in many studies around the world (Huang et al., 2012;Noh et al., 2007). It supports the input of a variety of meteorological data, and meteorological reanalysis data provided by NOAA were used in this study.

MERRA-2
MERRA-2 (The Modern-Era Retrospective Analysis for Research and Applications, Version 2) is an atmospheric reanalysis dataset provided by the National Aeronautics and Space Administration (NASA) and the Global Modeling and Assimilation Office (GMAO) (Gelaro et al., 2017). The aerosol optical property data used in this study were derived from the 1 h average product of the MERRA-2 tavg1_2d_aer_Nx dataset and CO column concentration from the MERRA-2 tavg1_2d_chm_Nx. The spatial resolution of the two datasets was 0.625 • × 0.5 • . During data processing, the 24 h data of each day were averaged to obtain daily average data.

ERA5
ERA5 is a global atmospheric reanalysis dataset provided by the European Centre for Medium-range Weather Forecasting (ECMWF) . In recent years, some studies have evaluated the accuracy of reanalysis data provided by the ECMWF based on radiosonde data. For example, Luo et al. (2020) found that the average RH discrep- ancy between the ERA-Interim and radiosonde was within 10 % below 500 hPa. Song et al. (2020) found that the root mean square error (RMSE) of ERA5 RH was 3.85 % compared with the RH profile of the radiosonde. The above results show that RH from reanalysis data has good accuracy, and it has been widely used in various research fields (Sajadi et al., 2020;Tzanis et al., 2019;Xiao et al., 2020). The RH data of ERA5 used in this paper were divided into 37 layers vertically (1-1000 hPa). The temporal resolution was 1 h and the spatial resolution was 0.125 • × 0.125 • . In order to ensure the spatial consistency, the RH data (31.25 • N, 121.375 • E) closest to the location of the Raman lidar were used.

Lidar data processing
Original signals need to be pre-processed before retrieval, including background subtraction, photon-counting signal dead-time correction, signal gluing, and overlap correction . The calculation of the glue coefficients in this study used the methods proposed by Newsom et al. (2009). In order to reduce the influences of lidar incompletely overlapping detection areas on retrieved results, only signals in the complete overlap area were used for retrieval. In addition, affected by the location altitude of the Raman lidar and the least square method used in the retrieval process, the lowest height of LR obtained by the Raman method was 569.5 m (a.s.l.). Since the Raman lidar used in this study can detect the Raman scattering signal of 387 nm nitrogen and signal-to-noise ratios of Raman signals in daytime are much lower than that in the nighttime, the 355 nm LR at night can be obtained through retrieval. The retrieval results of raw signals were counted by hour, and the hours with more than 15 minutes of retrieval results were regarded as effective observation hours. The retrieval results within the effective observation hours were averaged to obtain hourly average data. During the observation period, data of 667 effective observation hours were obtained through retrieval and statistics. The monthly distribution is shown in Fig. 1.
Previous studies indicated that there are some sources of errors in the retrieval. The relative errors of particle extinction coefficients caused by assumed air density profile are 1.5 % (Masonis, 2002), and the relative errors of particle 5380 T. Liu et al.: Long-term variation in aerosol lidar ratio in Shanghai backscatter coefficients caused by reference height can be 10 % (Ansmann et al., 1992). The mean deviations of particle extinction coefficients caused by signal detection are within 15 % in the 350-2000 m height range and within 20 % in the 3000-4000 m height range (Pappalardo et al., 2004). The difference is caused by different signal-to-noise ratios at low altitude and high altitude. Due to the low signal-to-noise ratio, there were usually more missing values at high altitudes.
3 Results and discussions 3.1 LR temporal and vertical variations 3.1.1 Variation in the averaged LR Figure 2a shows the averaged LR profile for 667 h. Because of the variability of aerosol particle size and microphysical properties with height (Singh et al., 2005), the averaged LR was characterized by large variability, ranging from 17 to 82 sr. LR reached a maximum at the height of 600 m and decreased with the increase in height. The averaged LR in the height range of 0.5-5 km was 41.0 ± 22.5 sr, and it was relatively small (24.8 ± 13.7 sr) above 2 km. These discrepancies were in good agreement with the results, which were mostly less than 25 sr in the altitude range of 2 km to 3 km, that Hee et al. (2016) found. LR depends on the aerosol size distribution and refractive index (Takamura et al., 1994;Young et al., 1993). The small LR may be caused by the shape effect of aerosol particles (not obviously nonspherical) and relatively low absorption efficiency .
In order to investigate the variations in LR at different altitude ranges, Fig. 2b presents the averaged LR for different altitude ranges. The averaged LR from 0.5 to 1 km was 68.2 ± 19.5 sr, which was in good agreement with the 355 nm LR in Oklahoma, USA, observed by Ferrare et al. (2001). The mean values of LR were between 40 sr and 50 sr in the altitude range of 1-2 km, and mean values of LR were usually less than 40 sr above 2 km, which was related to low aerosol concentration and the low absorption efficiency of aerosols (Hänel et al., 2012;Hee et al., 2016). Furthermore, the slope of LR for the different height ranges as shown in Fig. 2b also gradually decreased with the increased height. Below 3 km, LR decreased rapidly with the increase in altitude, and the largest slope was −17.83 below 1.5 km. The reason for that was that temperature inversion in the low layer of the planetary boundary layer (PBL) at night weakened the vertical movement of the atmosphere and inhibited the diffusion of pollutants emitted by human activities such as vehicles and fossil fuel combustion. The accumulation of pollutants in the low layer of the PBL resulted in significant differences in aerosol vertical distribution and rapidly decreased aerosol extinction coefficients in the PBL (Liu et al., 2017;Wang et al., 2020a). However, above 3 km, the low aerosol concentration and homogeneous vertical distribution of aerosol led to small differences in LR at different height ranges. Figure 3 shows the frequency distribution of LR for different altitude ranges. Overall, LRs were widely distributed in the altitude range of 0.5-5 km. In most cases (about 90 %) LR ranged from 10 to 80 sr with the highest frequency of 17.3 % between 40 and 50 sr. It should be noted that the number of observations trailed off at larger LR and the frequency of abnormally large LR (>90 sr) was about 4 %. LR also had a wide distribution range within 0.5-2 km, and the frequency of 40-50 sr was the highest (24.6 %), which was similar to the range of 0.5-5 km. Large LRs (>60 sr) were mainly distributed in the range of 0.5-2.0 km, suggesting that aerosols in this height range had a strong absorbing ability. Although there were a few large LRs (>60 sr) above 2 km, LRs were mainly distributed between 0-40 sr with the highest frequency of 34 % between 10 and 20 sr. Figure 4a presents seasonal variation in LR over Shanghai during the observation period. The seasonal average LR was the largest in autumn with 47.6 ± 25.1 sr and the lowest in spring with 39.1 ± 19.6 sr. Generally, the LR of aerosol particles with stronger absorption ability is larger (Müller et al., 2007). Black carbon (BC) has a strong ability to absorb visible light (Chow et al., 2009), and aerosols rich in BC tend to show a larger LR, such as biomass burning aerosols (Giannakaki et al., 2016). Smoke produced by the burning of crop residues during the harvest season in autumn often affects Shanghai (Xu et al., 2018). According to the results found by Wang et al. (2014), the BC concentration in Shanghai was the lowest in spring and higher in autumn and winter. It indicated that the aerosol absorption ability in spring was weaker compared with that in autumn and winter.

Temporal variations in LR
The statistics of the averaged LR for a different height range in each month were shown in Fig. 4b. LR for all months decreased with the increase in altitude. The averaged LR below 2 km was the largest in October, which was attributed to smoke aerosols produced by biomass burning in the surrounding cities and rural areas during the harvest season (Nie et al., 2015). In view of LR vertical variations in different months, aerosols with LR > 40 sr were confined to within 1.5 km from March to August, while in the other months it was 2 km, especially in October and November (within 2.5 km). This revealed the point that the diffusion heights of aerosols with absorbing properties were characterized by strong seasonality. In spring and summer, the diffusion height of absorbing aerosols was lower; by contrast, it was higher in autumn and winter. The discrepancies in the seasonal height distribution of aerosols can be attributed to the monsoon climate (He et al., 2006;Liu et al., 2020b;Wang et al., 2016). Shanghai, in the southeast of China, was affected by the subtropical monsoon climate, and the prevailing winds were southeast and northwest in spring and southeast in summer   (Cai et al., 2010). Clean air from the sea could reduce air pollution to a large extent by diluting pollutant concentration . In particular, in March and April the LR above 2.5 km was higher than that in the other months, which may be due to the influence of dust aerosols brought by the prevailing northwest wind in spring. Liu et al. (2020b) pointed out that dust aerosols from Inner Mongolia and the Gobi desert had a high frequency in spring at the altitude of 3-5 km in Shanghai. The prevailing winds in autumn and winter were northeast and northwest, respectively. Absorbing aerosols originating from north China resulted in relatively large LRs at high altitude in Shanghai. Liu et al. (2012) reported that vast majority of aerosol particles in the Yangtze Delta region (including Shanghai) were below 2 km. In order to precisely analyze the variation characteristics of the LR in Shanghai, Fig. 5 shows LRs of 667 effective observation hours below 2 km. The abnormally large LRs (>80 sr) were usually distributed in the PBL, meaning that local emissions were the main source of strongly absorbing aerosols. From the figure, one can also conclude that the number of larger LRs within the PBL followed a decreasing trend with the passage of observation time, alluding to gradually reduced emissions of absorbing aerosols in Shanghai. The reduction was in good agreement with the reductions of BC and particulate matter (PM) concentrations caused by a series of energy-saving and emission-reduction measures such as the Shanghai Clean Air Action Plan (2018-2022) implemented by the Shanghai government in recent years (Wei et al., 2020).

Reasons for LR variation with height
From Eq. (1), we found that LR was negatively correlated with P (π ), and the nonspherical geometry of particles can cause a reduction in the P (π ) (Müller, 2003). δ can reflect the regularity of particle shape (Novitsky and Philbrick, 2005) and is defined as the ratio of vertical backscatter coefficient and parallel backscatter coefficient (Behrendt and Nakamura, 2002). The smaller value of δ is, the closer the particle shape is to spherical (Gobbi, 1998). Since the Raman depolarization lidar used in this study can detect the 532 nm polarization signals, 532 nm δ was used to analyze the sphericity of particles. It can be seen from Fig. 6a that LR was proportional to δ with a correlation coefficient of 0.86, which was consistent with Reagan's (1988) study that  LR increased with the increase in particle asphericity. LR is the ratio of the extinction coefficient and the backscatter coefficient. Therefore, there are two reasons for the increase in LR with the increasing particle asphericity. On one hand, the backscatter coefficient decreases significantly with the increase in particle asphericity; on the other hand, the extinction coefficient is sensitive to the particle cross section and is less affected by particle shape.
It is worth noting that LR responding to large δ was in the range of 100-120 sr. Generally, the LR of dust aerosols with a large δ was between 40 and 60 sr (Murayama et al., 2004;Noh et al., 2007). Hee et al. (2016) found that the 355 nm LR of aged forest fire aerosols was relatively large, ranging from 80 to 120 sr. And previous studies have found that some aged forest fire aerosols also showed a large depolarization ratio (Hu et al., 2019;Murayama et al., 2004). There might be two reasons for this phenomenon. One is that dust aerosols on the surface are lifted into the biomass burning plume (Müller et al., 2007) and the other is the nonsphericity of particles due to the coagulation of smoke particles during the aging process (Reid et al., 1998). Therefore, the LR in the range of 100-120 sr may correspond to large δ because of aged forest fire aerosols. In order to explore whether the decrease in LR with increasing height is affected by the particle shape effect, Fig. 6b shows the average profile of δ at lidar observation time. The average and median of δ gradually decreased with the increase in height, indicating that the particle shape became more regular with the increase in height. The result evidenced the inference by Tesche et al. (2007) that the regular particle shape was one of the reasons for the smaller LR at high altitudes.

Influences of aerosol sources on LR
In order to further understand the influences of wind directions on LR and its vertical distribution, a cluster analysis of back trajectories was used to study the transport of atmospheric aerosols. Based on the HYSPLIT-4 model (Franke et al., 2001;Noh et al., 2007), the 72 h backward trajectories at the height of 1000 m were shown in Fig. 7a. The cluster analysis resulted in four main air mass directions (Hänel et al., 2012;Pietruczuk and Podgorski, 2009). Backward trajectory cluster analysis based on the HYSPLIT model is widely used in atmospheric aerosol research Xu et al., 2018;Zhang et al., 2020). We performed a significance test on the cluster analysis results, and the one-way ANOVA showed that P <0.05, indicating that the LR of the four clusters referred to as 1, 2, 3, and 4 was significantly different. Similarly, there were significant differences in δ among the four clusters.
The averaged LR and δ between 0.5 and 5 km for the four clusters are shown in Fig. 7b, and Fig. 7c presents the distribution of LR and δ at different heights for the four clusters. The averaged LR (38.7 ± 24.2 sr) and δ (0.030 ± 0.021) of air mass 1 were the lowest in all clusters. As air mass 1 came from the western Pacific Ocean, it could bring abundant marine aerosols. Sea salt particles are characterized by a coarse mode and are spherical in wet conditions. In addition, marine aerosols have lower LRs than dust aerosols because of their weak absorption ability (Papagiannopoulos et al., 2018). Interestingly, the 355 nm LR of clean marine aerosols was usually between 20 and 35 sr (Alexander and Protat, 2019;Hee et al., 2016), which was smaller than that of aerosols brought by air mass 1 in the study. It was observed that the averaged LR in the range of 0.5-1 km was more than 60 sr, which suggested that clean marine aerosols from the sea might be mixed with local absorbing aerosols in Shanghai (Franke et al., 2001;Müller et al., 2007).   The averaged LR affected by the aerosols brought by air mass 2 was approximately equivalent to the LR affected by air mass 1, with an averaged value of 39.4 ± 19 sr at 0.5-5 km altitude. The source region of air mass 2 was located in Inner Mongolia, which could bring dust aerosols. During the transport over the ocean, clean marine aerosols mixed with dust aerosols led to LRs larger than that of the clean marine aerosols.
The averaged LR of air mass 3 was 44.2 ± 24.7 sr, the largest in the four clusters. As air mass 3 passed through the northern Chinese region where pollution levels and pollutant amounts were relatively high, aerosol particles brought by the air mass had a strong absorbing ability. In addition, it was evident that δ corresponding to air mass 3 below 2.5 km was larger than that of the other three air masses in Fig. 7c. The larger δ hinted at a high contribution of irregularly shaped aerosol particles. Air mass 3 passing through dust source areas of Mongolia and Inner Mongolia could bring abundant dust aerosols to Shanghai (Huang et al., 2012). Dust aerosols usually show larger particle depolarization ratios than other types of aerosols due to their irregular particle shape (Fu et al., 2009;Kai et al., 2008;Murayama et al., 1999). For example, Huang et al. (2012) found that the aerosol 532 nm depolarization ratio in the case of dust pollution was significantly greater than that in the cases of secondary inorganic pollution and biomass burning pollution. As a consequence, the larger δ corresponding to air mass 3 was attributed to dust aerosols, and long-range-transported dust aerosols could reach a height of about 2.5 km.
The averaged LR affected by aerosols from air mass 4 was 42.6 ± 21.8 sr. It should be noted that the averaged δ affected by aerosols brought by air mass 4 was small, which was comparable to that of air mass 1, indicating a high contribution of spherical aerosol particles, but LR was larger. By observing the trajectory of air mass 4, we found that it passed through Hubei Province with a high industrial level  and Anhui Province with heavy pollution from biomass burning (Wu et al., 2020). The industrial and smoke aerosols accompanying air mass 4 were approximately spherical (Giannakaki et al., 2016;Müller et al., 2007) and could be responsible for the smaller δ.
In summary, the variations and vertical distributions of LR and δ in Shanghai were caused by the synthetic impacts of long-range transport from different source areas and local emissions. Large LR and δ were the results of smoke and dust aerosols from the northwest. And the mixing of aerosols from the sea and locally emitted absorbing aerosols resulted in a smaller δ and slightly larger LR than that of clean marine aerosols.

Influence of atmospheric turbidity on LR
As mentioned previously, the vertical variations in absorbing aerosols and their influence factors played an important role in evaluating the aerosol radiation effect and studying the cause of pollution (Mishchenko et al., 2004). The LR vertical variation under different atmospheric turbidity has rarely been discussed; however, previous studies have analyzed the vertical profiles of LR in different pollution degree cases and their main concern was the averaged LR of the aerosol layer Wang et al., 2016). AOD is an important parameter to characterize aerosol optical properties, which can reflect aerosol content in the atmosphere and is also an important index to evaluate atmospheric quality and visibility (Cheng et al., 2015;Hess et al., 1998;Qi et al., 2013). Previous studies have shown a positive correlation between AOD and LR by analyzing averaged LR for different AOD ranges (Ferrare et al., 2001;He et al., 2006) due to the increase in aerosol absorption and extinction caused by the increase in small particles (Takamura et al., 1994).
AOD was obtained by integrating 355 nm extinction coefficients in the range of 0.5-2 km. The averaged profile of LR below 2 km in different AOD ranges was drawn as shown in Fig. 8. Under clean conditions, LR decreased more dramatically with the increase in height. By contrast, the lack of significant vertical variability of LR under high atmosphericturbidity conditions reflected the homogeneous vertical distribution of absorbing aerosols. The result that the vertical slope of LR presented a decreasing trend with increasing atmospheric turbidity can be explained by aerosol radiative effects on thermal structure and atmospheric stability. Under high atmospheric-turbidity conditions, aerosol particles that absorb a large amount of solar radiation during the day radiatively warm the surface at night but radiatively cool the air above the surface (Jacobson and Kaufman, 2006;Ramanathan et al., 2005). The decrease in the atmosphere stability due to the temperature difference increases vertical turbulence and results in the homogeneous vertical distribution of aerosols. By contrast, in the clear and pollution-free nights, the surface radiation cooling results in temperature inversion near the ground. The stable atmosphere is not conducive to the lifting of absorbing aerosols, resulting in a significant vertical variation in LR.

The main aerosol types causing the abnormal variation in LR
As shown in Fig. 5, abnormally large LR occurred occasionally in relatively high locations approximately above the top of PBL in spite of a usual decay trend in LR with height. To investigate the reasons, we selected 5 d with LR > 80 sr over 1 km. According to previous reports, biomass burning aerosols are relatively small and spherical, and their strong absorption ability makes them have a large LR (Papagiannopoulos et al., 2018). For example, Amiridis et al. (2009) observed smoke plumes from biomass burning over Greece and found that the 355 nm LR ranged from 40 to 100 sr. Giannakaki et al. (2016) used a Raman lidar in South Africa and found that the 355 nm LR of biomass burning aerosols was 92 ± 10 sr. Generally, the determina- tion of aerosol type with a large LR observed by Raman lidar was based on fire data and backward trajectory models such as HYSPLIT (Hee et al., 2016;Noh et al., 2008). In our study, the spatial distribution of biomass burning tracers was used to determine whether the abnormally large LRs were related to biomass burning aerosols (Huang et al., 2012), which could lay the foundation for future research on aerosol threedimensional spatial distribution and pollution causes. Biomass burning is one of the important sources of PM, organic carbon (OC) and BC in the atmosphere (Wu et al., 2020). It also emits pollutant gases such as CO, SO 2 , NO x , and HCN (Andreae and Merlet, 2001;Kalluri et al., 2020;Randel et al., 2010). CO can be used as a tracer for biomass burning. For example, Huang et al. (2012) found that CO column concentration in the biomass burning zone was significantly different from that in the non-biomass burning zone. In addition, biomass burning would also result in high AOD and AAOD (absorbing aerosol optical depth) in the region due to strong absorption of biomass burning aerosols (Shaik et al., 2019). For example, He et al. (2015) found that 500 nm AOD increased from 0.73 to 1.00 when analyzing the smoke plume of biomass burning in Shanghai. Similarly, Vadrevu et al. (2011) found that the AOD of the wheat stalks burning season and the rice stalk burning season were both high by using satellite data in India: 0.598 and 0.58, respectively. CAOD is the optical depth of black carbon aerosols and organic carbon aerosols. Hence, 550 nm AOD, AAOD, CAOD, and CO column concentrations were used as tracers to determine whether these five cases were affected by biomass burning. Figure 9 depicts the spatial distribution of four tracers for the five cases. AOD in Shanghai during these 5 d ranged from 0.45 to 1.05, pointing to heavy aerosol pollution. The average value of AAOD was between 0.0375 and 0.1, and CAOD was between 0.12 and 0.24, which reasonably demonstrated the presence of absorbing aerosols and carbon aerosols in Shanghai (Shaik et al., 2019). Additionally, the CO column concentrations in Shanghai for these 5 d were relatively high, all more than 9 × 10 −4 kg/m 2 , which reflected significant possibility of smoke advection. Therefore, it could be inferred that Shanghai was affected by biomass burning aerosols from local rural areas or surrounding provinces in the five cases.
Although the abnormally large LR above 1 km was mainly relevant for the advection of biomass burning aerosols, it should be noted that the increasing aerosol extinction caused by the increase in RH could also result in large LR (Salemink et al., 1984). For example, Ackerman (1998) found that the LR of continental aerosols increased from 40 to 80 sr with RH. Figure 10 presents LR and RH profiles for three cases and shows that LR was a function of RH. The abnormally large LRs above 1 km had a good corresponding relationship with high RH, which demonstrated that the abnormally larger LRs were also related to high RH.

Conclusions
For the first time, a long-term (2017-2019) observation based on Raman lidar was carried out in Shanghai. The aerosol 355 nm LR was retrieved, and the variations in LR and their influence factors were analyzed. In the height range of 0.5-5 km, about 90 % of LRs were distributed in 10-80 sr, with an average of 41.0 ± 22.5 sr, and LR decreased with the increase in height. The averaged LR in autumn was the largest, which was 47.6 ± 25.1 sr. The averaged LRs in summer and winter were close (41.0 ± 21.6 and 42.0 ± 27.3 sr, respectively), and the averaged LR in spring was the smallest. The seasonal variation in LR was closely related to the seasonal variations in BC concentration. In the height range of 0.5-2 km, the monthly average LR was the largest in October, which was related to biomass burning aerosols produced by burning straw in the surrounding areas during the harvest season. In addition, affected by the prevailing winds in spring and summer, the aerosols with LRs > 40 sr were confined within 1.5 km from March to August.
LR and δ were positively correlated, meaning that the more regular the particle shape is, the smaller the LR is. δ decreased with the increase in height, which proved that particle shape was one of the factors affecting the vertical distribution of LR. LR had a strong dependence on the source directions of air masses. Large LR coincided with the air masses from the northwest, while the air masses from the east led to small LR. In addition, the shape of aerosol particles was the most irregular due to aerosols brought by air masses from the northwest. The vertical distribution of LR was affected by atmospheric turbidity, and the slope of LR vertical variation decreased with the increase in AOD.  We analyzed the spatial distribution of 500 nm AOD, AAOD, CAOD, and CO column concentrations of five cases with LRs > 80 sr and found that large LRs above 1 km in Shanghai were related to biomass burning aerosols. In addition, the large LRs above 1 km could be also related to high RH.
Data availability. The data presented in this paper are available from the corresponding authors upon request.
Author contributions. TL retrieved the data and wrote the paper. QH and YC formulated the project goals and edited and reviewed the paper. JL, QL, and WG downloaded and analyzed the reanalysis data. GH, WS, and XY revised the paper.
We also gratefully acknowledge the ECMWF for the provision of the ERA5 dataset. Review statement. This paper was edited by Jianping Huang and reviewed by two anonymous referees.