# Determination and climatology of the diurnal cycle of the atmospheric mixing layer height over Beijing 2013-2018: lidar measurements and implications for air pollution

Haofei Wang, Zhengqiang Li, Yang Lv, Ying Zhang, Hua Xu, Jianping Guo, Philippe Goloub

► To cite this version:

Haofei Wang, Zhengqiang Li, Yang Lv, Ying Zhang, Hua Xu, et al. Determination and climatology of the diurnal cycle of the atmospheric mixing layer height over Beijing 2013-2018: lidar measurements and implications for air pollution. Atmospheric Chemistry and Physics, 2020, 20, pp.8839-8854. 10.5194/acp-20-8839-2020. insu-03686298

HAL Id: insu-03686298

https://insu.hal.science/insu-03686298

Submitted on 3 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Determination and climatology of the diurnal cycle of the atmospheric mixing layer height over Beijing 2013–2018: lidar measurements and implications for air pollution

Haofei Wang<fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Zhengqiang Li<fnref n="1" />, Yang Lv<fnref n="1" />,<fnref n="2" />, Ying Zhang<fnref n="1" />, Hua Xu<fnref n="1" />, Jianping Guo<fnref n="4" />, and Philippe Goloub<fnref n="5" />

**Correspondence:** Zhengqiang Li (lizq@radi.ac.cn)

Received: 23 February 2020 – Discussion started: 16 March 2020

Revised: 8 June 2020 – Accepted: 20 June 2020 – Published: 27 July 2020

**Abstract.** The atmospheric mixing layer height (MLH) determines the space in which pollutants diffuse and is thus conducive to the estimation of the pollutant concentration near the surface. The study evaluates the capability of lidar to describe the evolution of the atmospheric mixing layer and then presents a long-term observed climatology of the MLH diurnal cycle. Detection of the mixing layer heights (MLH<sub>L</sub> and MLH'<sub>L</sub>) using the wavelet method based on lidar observations was conducted from January 2013 to December 2018 in the Beijing urban area. The two dataset results are compared with radiosonde as case studies and statistical forms. MLH<sub>L</sub> shows good performance in calculating the convective layer height in the daytime and the residual layer height at night. While MLH'<sub>L</sub> has the potential to describe the stable layer height at night, the performance is limited due to the high range gate of lidar. A nearly 6-year climatology for the diurnal cycle of the MLH is calculated for convective and stable conditions using the dataset of MLH<sub>L</sub> from lidar. The daily maximum MLH<sub>L</sub> characteristics of seasonal change in Beijing indicate that it is low in winter (1.404±0.751 km) and autumn (1.445±0.837 km) and high in spring (1.647±0.754 km) and summer (1.526±0.581 km). A significant phenomenon is found from 2014 to 2018: the magnitude of the diurnal cycle of MLH<sub>L</sub> increases year by year, with peak values of 1.291±0.646 km, 1.435±0.755 km,

1.577 ± 0.739 km, 1.597 ± 0.701 km and 1.629 ± 0.751 km, respectively. It may partly benefit from the improvement of air quality. As to converting the column optical depth to surface pollution, the calculated PM<sub>2.5</sub> using MLH<sub>L</sub> data from lidar shows better accuracy than that from radiosonde compared with observational PM<sub>2.5</sub>. Additionally, the accuracy of calculated PM<sub>2.5</sub> using MLH<sub>L</sub> shows a diurnal cycle in the daytime, with the peak at 14:00 LST. The study provides a significant dataset of MLH<sub>L</sub> based on measurements and could be an effective reference for atmospheric models of surface air pollution calculation and analysis.

## 1 Introduction

The height of the mixing layer (MLH) is a crucial parameter for near-surface air quality forecasts, pollutant dispersion and the quantification of pollutant emissions (Haeffelin et al., 2012; Seibert et al., 2000; Baars et al., 2008; Liu and Liang, 2010; de Bruine et al., 2017). The pollutants discharged in the boundary layer diffuse vertically under the drive of turbulence (Gan et al., 2011; Monks et al., 2009; Guo et al., 2016) and finally become completely mixed over this layer if sufficient time is given (Emeis et al., 2008). The MLH determines the space in which pollutants diffuse and is thus con-

<footnote n="1">¹State Environmental Protection Key Laboratory of Satellite Remote Sensing, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing, 100101, China</footnote>

<footnote n="2">²University of Chinese Academy of Sciences, Beijing, 100101, China</footnote>

<footnote n="3">³Key Laboratory of Radiometric Calibration and Validation for Environmental Satellites, National Satellite Meteorology Center, Beijing, 100081, China</footnote>

<footnote n="4">⁴State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing, 100081, China</footnote>

<footnote n="5">⁵Laboratoire d'Optique Atmospherique, UMR8518, CNRS – Université de Lille 1, Villeneuve d'Ascq, Lille, 59000, France</footnote>
ducive to the estimation of the pollutant concentration near the surface, which might be detrimental to the health of humans and ecosystems (Emeis et al., 2007; Collaud Coen et al., 2014; Singh et al., 2016; Mues et al., 2017).

Within the planetary boundary layer (PBL) the height of the mixing layer (ML) is defined as the height up to which vertical dispersion by turbulent mixing of air pollutants takes place due to the thermal structure of the PBL (Seibert et al., 2000; Schafer et al., 2006; Emeis et al., 2007). The MLH depends largely on the synoptic weather situation (Emeis et al., 2008). The MLH can be estimated by the detection of variance in the mechanical turbulence, the temperature enabling convection or the substance content in the low troposphere. Singh et al. (2016) investigated the evolution of the local boundary layer in the central Himalayan region using a radar wind profiler detecting wind components based on the signal-to-noise ratio profile. Collaud Coen et al. (2014) compared MLH measurements of a microwave radiometer from an atmospheric temperature profile with other measurements in the Swiss Plateau. Mues et al. (2017) used a ceilometer to retrieve the MLH based on the aerosol backscatter signal in the Kathmandu Valley. These measurements are based on different atmospheric parameters, different measuring instruments and various analysis algorithms, leading to MLH results obtained by different methods being inconsistent (Collaud Coen et al., 2014).

In order to move from the general definition to practical measurements, it is necessary to separately consider the structure of the convective boundary layer (CBL) and the stable boundary layer (SBL). In the case of fair weather days, the PBL height has a well-defined structure and diurnal cycle (Collaud Coen et al., 2014; see Fig. S1 in the Supplement). The PBL development under strong convection driven mainly by solar heating is called the CBL (Collaud Coen et al., 2014). The nocturnal SBL shows a more complex internal structure, including a stable layer caused by radiative cooling from the ground, and gradually merges into a neutral layer called the residual layer (RL) (Stull, 1988; Mahrt et al., 1998; Salmond and McKendry, 2005; Collaud Coen et al., 2014). The RL height is the top of the neutral layer and the beginning of the stable free troposphere. The pollutants discharged from the surface at night are restricted to the SBL, while the pollutant emissions from past days tend to stay in the RL. In addition to the dominance of the CBL in the afternoon, the SBL and neutral boundary layer may be formed under certain weather conditions (Stull, 1988; Poulos et al., 2002; Medeiros et al., 2005; Zhang et al., 2018).

Since most atmospheric column aerosol particles are usually present in atmosphere below the MLH, the MLH can be used to convert the aerosol optical thickness of the column observed by sun photometers and satellites to the concentration of near-surface pollutants (Sifakis et al., 1998; Dandou et al., 2002; Emeis et al., 2007). Particulate can be used as an important indicator of atmospheric layering because its vertical distribution is strongly affected by the thermal struc-

ture of the atmosphere (Neff and Coulter, 1986). Provided the vertical aerosol distribution adapts rapidly to the variational thermal dynamics of the boundary layer, the MLH can be retrieved from the analysis of this aerosol distribution (Emeis et al., 2008).

By the measurement of the profile of aerosol, lidar offers a direct and continuous way to monitor the diurnal cycle of the different layers constituting the PBL (Seibert et al., 2000; De Haij et al., 2006; Emeis et al., 2008; Liu and Liang, 2010; Tang et al., 2016; Su et al., 2017, 2019). With recent upgrades of hardware, a ceilometer, also known as an automated low-power lidar or an automated lidar-ceilometer (ALC), has been demonstrated to capably determine the MLH (Wiegner et al., 2006; Martucci et al., 2007, 2010; Geiß et al., 2017; Kotthaus and Grimmond, 2018; Mues et al., 2017). Recent studies have compared remote sensing measurements (lidar, radar wind profiler, microwave radiometer) with radiosonde (RS) (Wiegner et al., 2006; Milroy et al., 2012; Sawyer and Li, 2013; Cimini et al., 2013; Tang et al., 2016; Singh et al., 2016; Mues et al., 2017; Su et al., 2019), with which convection weather cases have good correlation with differences of 100–300 m, while non-convective weather conditions lead to much larger differences in the MLH estimations if the approaches are intended to measured different structures of the ML such as the CBL, SBL or RL (Collaud Coen et al., 2014). Meteorological radiosondes usually acquire the MLH in the morning (08:00 LT) and at night (20:00 LT), when the diurnal cycle of the ML combined with the stable and convective PBL cannot be well characterized.

In existing studies, numerical simulations, ground-based remote sensing or meteorological radiosondes are used to obtain the characterization of the MLH during short time periods in Beijing, mainly focusing on heavy pollution events (Yang et al., 2005; Quan et al., 2013; Hu et al., 2014; Zhang et al., 2015), which underscores the scarcity of continuous high-resolution measurements for a long time period. Depending on the measured atmospheric parameters and observational uncertainties, different measurement approaches may reveal different aspects of PBL structure (Seibert et al., 2000; Seidel et al., 2010; Beyrich and Leps, 2012). Thus, it is of great significance to apply consistent algorithms to consistent types of atmospheric structure parameters when comparing the MLH from different times.

The main aim of this study, therefore, is to present a long-term observed climatology of the MLH diurnal cycle based on lidar observations. For that, the capability of lidar to describe the diurnal evolution of the mixing layer height is evaluated first. The data and methods used are described in Sect. 2. Section 3 contains the results and discussion, which consist of the comparison of lidar-derived MLH with radiosonde measurements, the climatology of the MLH in Beijing, and implications for surface pollution retrieval. Then, it is concluded in Sect. 4.
## 2 Data and methods

### 2.1 Site and lidar measurements

Beijing is the capital of China with about 20 000 000 citizens. Beijing is located on flat terrain in the North China Plain (altitude of 20–60 m), with the Taihang Mountains in the west and the Yan Mountains in the north (altitude of 1000–1500 m). Similar to many other metropolitan areas, Beijing suffers from episodes of poor air quality, in particular fine particulate matter (PM2.5). In the study, the observatory (116.379° E, 40.005° N) in the metropolitan area of Beijing is located on the building roof (59 m a.s.l.) of the Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences. A micro-pulse lidar (CE370, CIMEL, France) was used to detect the atmospheric aerosol structure at this site. The laser used is a frequency-doubled Nd:YAG with pulse repetition frequency 4.7 kHz and energy 8–20 µJ. The CE370 operates at a wavelength of 532 nm with the vertical resolution of 15 m and can detect a long-range profile up to 30 km every second. For the enhancement of the signal-to-noise ratio, 60 profiles are averaged to restore as one, with the time resolution of 1 min. The signal received from the lidar is processed by the subtraction of atmospheric background using the averaged value of the signal received from the height between 22 and 30 km. There are remnants of the previous signal in the system in the absence of optical signal reception, and these signals are called after-pulses. They are monitored every week and removed from the receive signal. Due to the design of the lidar, the received view close to the ground does not completely coincide with the transmitted view. There is a detection blind area in the lidar, and a geometric overlap factor is used to correct the mismatch of the field of view. Then, the correction of the range (range-corrected signal, RCS) is used to retrieve the MLH. The profile of the RCS is expressed as $f(z)$, with $z$ the measurement height. Logarithm calculation of the RCS (expressed in $\ln(\beta'_{532})$) is presented in the lidar image (Campbell et al., 2003; Yan et al., 2014; Su et al., 2019). MLH estimation from lidar systems is based on the measurement of the sudden drop in aerosol backscatter at the top of the mixing layer (Seibert et al., 2000). The period of lidar measurements is from 2013 to 2018, nearly 6 years. Except for the lidar data for 2013 mainly in winter and spring (months 1–4 and 11–12), the measurements for 2014–2018 are all annual continued observations.

### 2.2 MLH derived from lidar

Wavelet transforms are commonly used in many studies for MLH determination from lidar observations (Cohn and Angevine, 2000; Davis et al., 2000; Brooks, 2003; De Haij et al., 2006; Baars et al., 2008; Su et al., 2019). When the maximum value of the attenuated backscattering profile convolved with the Haar function is reached, the corresponding height is the MLH. The equation of the wavelet is defined as

follows:

$$h\left(\frac{z-b}{a}\right) = \begin{cases} 1, & b - \frac{a}{2} \le z \le b; \\ -1, & b < z \le b + \frac{a}{2}, \\ 0, & \text{else} \end{cases} \quad (1)$$

where $b$ is the transformation of the equation, the equation is cantered and $a$ is the expansion of the equation. The equation of the wavelet covariance transformation $W_f(a, b)$, namely the convolution of $f(z)$ with the wavelet function, is defined as follows:

$$W_f(a, b) = \frac{1}{a} \int_{z_b}^{z_t} f(z)h\left(\frac{z-b}{a}\right) dz, \quad (2)$$

where $f(z)$ represents the RCS at different heights, $z_b$ denotes the lower limit of the height of the profile and $z_t$ represents the upper limit of the height. A valid MLH$_L$ is detected, corresponding to the value $b$ when $W_f(a, b)$ reaches the biggest local maximum with a coherent scale of $a$ (Brooks, 2003; De Haij et al., 2006; Emeis et al., 2008). In this study, the expansion of $a$ is selected as 420, 435 and 450 m, and the final $W_f(a, b)$ is calculated from the averaged corresponding values. Another layer, called MLH'$_L$, is detected simultaneously by the first local maximum $W_f(a, b)$ from $z_b$, which is assumed to be smaller than or equal to MLH$_L$ (De Haij et al., 2006; Mues et al., 2017). Actually, every local maximum corresponds to an aerosol layer and several internal layers appear, making the allocation of a local maximum to an atmospheric feature very difficult (Morille et al., 2007; Geiß et al., 2017; Poltera et al., 2017; Kotthaus and Grimmond, 2018). Since the absolute maximum in the vertical gradient of the lidar profiles is characterized by a rapid decrease in pollutant concentrations, the MLH$_L$ can be associated with the CBL height during daytime (Haeffelin et al., 2012; Poltera et al., 2017) and the RL height during nighttime (Collaud Coen et al., 2014). However, the interpretation of the first local maximum (MLH'$_L$) is critical.

To form a diurnal cycle of the MLH from these several layers, a geodesic approach was applied to pathfinderTURB (Poltera et al., 2017), while COBOLT (Geiß et al., 2017) uses a time-height tracking approach with moving windows. Nevertheless, these methods are all based on the selection of the lowest detected aerosol layer. The height of the lowest detected aerosol layer was regarded as the daytime MLH and the nocturnal stable boundary layer, as reported by Mues et al. (2017) and Kotthaus and Grimmond (2018), respectively. Su et al. (2019) developed a depend time and depend stability (DTDS) algorithm, started with the lowest point and tracked depending on the time and stability, but the nocturnal MLH with SBL height is not evaluated. The detection of nocturnal boundary layer heights, in contrast to the residual layer, is a major challenge (Haeffelin et al., 2012; Lotteraner and Piringer, 2016; de Bruine et al., 2017). Thus, one of the objectives of this study is to investigate the usefulness of MLH'$_L$ from CE-370 to capture the SBL height over Beijing.
MLH retrievals are eliminated if a cloud flag is marked when the cloud base is found within 6 km of the surface, and a threshold is selected to distinguish between clouds and aerosol layers. To improve the retrieval, a Gaussian filter is applied to retrievals to smooth the temporal variability, and unrealistic outliers are deleted. Due to the limitation of the algorithm and insufficient lidar overlap, the minimum range of the MLH calculation from CE-370 is on the order of 250 m, which is higher than the order of 50 m with a ceilometer. This is due to the optical design of the ceilometer using the same lens for the emitter and the receiver optical paths, which suffers from a low signal-to-noise ratio and provides a lower overlap with the limited power transmitted from the optical design. Detecting significant vertical gradients of attenuated backscatter can be challenging on a clear day when aerosol content is low (Eresmaa et al., 2012; Haeffelin et al., 2012). Compared to CE-370, a ceilometer usually needs a large-scale temporal and vertical average at the cost of a reduction of retrieval in relatively clean atmospheric conditions (de Bruine et al., 2017; Kotthaus and Grimmond, 2018).

## 2.3 MLH from radiosonde

Radiosonde (RS) measurements are one of most widely used methods, especially in China, to derive the SBL height and CBL height due to the ability to characterize the thermodynamic and dynamic states of the boundary layer (Piringer et al., 2007; Seidel et al., 2010; Guo et al., 2019). Meteorological radiosondes perform measurements at the international standard weather station (39.484°N, 116.282°E), which was located nearly 11 km from the lidar station. It includes two categories: conventional observations around the year, which are performed at 00:00 UTC (08:00 LST) in the morning and at 12:00 UTC (20:00 LST) in the evening each day, and intensified observations only in summer, which are taken at 06:00 UTC (14:00 LT) in the afternoon. The observed meteorological parameters include atmospheric pressure (P), temperature (T), relative humidity (RH), wind speed (WS) and wind direction (WD).

The bulk Richardson number ($Ri_b$) is a dimensionless parameter combining the thermal energy and the vertical wind shear; it is widely used in MLH climatology (Seidel et al., 2012; Collaud Coen et al., 2014). $Ri_b$ is defined as the ratio of turbulence associated with buoyancy to that induced by mechanical shear, which is expressed as

$$Ri_b = \frac{gz(\theta(z) - \theta(z_s))}{\theta_{z_s}(U^2(z) + V^2(z))}, \quad (3)$$

where z is the height ($z > z_s$), the subscript “s” denotes the surface), $\theta$ characterizes the virtual potential temperature, U and V indicate the two horizontal wind velocity components, and g represents the Earth’s gravitational constant. The MLH$_{RS}$ corresponds to the first elevation z, with $Ri_b$ greater than a critical threshold taken as 0.25 (Stull, 1988; Seidel et al., 2012; Guo et al., 2016, 2019). In most cases, the exact

threshold value has only a small impact on the PBL height due to the large slope of $Ri_b$ in this interval (Collaud Coen et al., 2014).

## 2.4 Air pollution model

Data on the MLH are usually combined within an atmospheric model to obtain the surface air pollutant concentration. For example, PM$_{2.5}$ remote sensing (PMRS) models, such as that derived by Zhang and Li (2015), have the ability to calculate the mass concentration of PM$_{2.5}$ above the ground. The PMRS method is designed to employ currently available remote sensing parameters, including aerosol optical depth (AOD), fine-mode fraction (FMF), planetary boundary layer height (PBL height) and atmospheric relative humidity (RH), to derive PM$_{2.5}$ from instantaneous remote sensing measurements under different pollution levels (Zhang and Li, 2015; Li et al., 2016; Yan et al., 2017). PM$_{2.5}$ is calculated following the PMRS model as

$$PM_{2.5} = \frac{AOD}{PBLH} \cdot \frac{FMF \cdot VE_f(FMF) \cdot \rho_{2.5,dry}}{f_0(RH)}, \quad (4)$$

where AOD indicates aerosol optical depth and FMF represents the fine-mode fraction; $VE_f$ is the ratio of the volume and extinction of fine-mode aerosol, which can be calculated from FMF as $VE_f(FMF) = 0.2887FMF^2 - 0.4663FMF + 0.356$. The parameter $\rho_{2.5,dry}$ indicates the density of dry PM$_{2.5}$, while $f_0(RH)$ represents the particle hydroscopic growth function, which is $f_0(RH) = (1 - RH/100)^{-1}$. The PBL height can be derived from remote sensing and radiosonde measurements.

All the parameters are observed by the instruments employed in the same observatory as the lidar. The optical parameters of the column aerosols (AOD and FMF) are obtained by a sky-sun photometer (CE318-DP, CIMEL, France), which is affiliated with the Aerosol RObotic Network (AERONET) (Holben et al., 1998; Dubovik et al., 2000). Measurements are automatically scheduled with direct sun irradiance measurements of about 15 min each and angular sky radiance scanning of about 1 h each (Li et al., 2015; Che et al., 2014; Wang et al., 2019). Atmospheric meteorological data (relative humidity, wind speed, wind direction, etc.) are obtained by an automatic meteorological monitoring station (BLJW-4). The PM$_{2.5}$ mass concentration is obtained by a PM$_{2.5}$ monitor (BAM-1020; MetOne, USA), which shows good agreement with the measurements of the national monitoring network near the observatory. All the data are quality-controlled and calculated as 1 h averages, and the measurement period is from 2014 to 2018. The MLH obtained from both lidar and radiosonde within the period is used in the model to calculate the surface PM$_{2.5}$. For convenient comparison with air quality and meteorological parameters, all MLH results are 1 h averaged.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (a) Daily backscatter profiles from lidar for 2–4 March 2017. The lines connected by black dots in the figure represent the retrieved MLH<sub>L</sub>, while white dots line indicate the MLH'<sub>L</sub>, and the top of the purple triangle indicates the MLH<sub>RS</sub> identified from radiosonde. The horizontal axis represents the local standard time (LST), and the vertical axis represents the height. The color bar denotes the logarithm of the attenuated backscattering coefficient (ln(β'<sub>532</sub>)). (b, c, d) The vertical profile of RCS (orange curve) from lidar and the wavelet coefficient (red curve) of RCS, as well as the vertical profile of temperature (T) (blue curve) and relative humidity (RH) (cyan curve) for (b) 20:00 LST on 3 March 2017, (c) 08:00 LST on 4 March 2017 and (d) 20:00 LST on 4 March 2017, as indicated by the white-outlined triangle in panel (a).</figcaption></figure>

# 3 Results and discussion

## 3.1 MLH operational measurement

A selection of typical atmospheric conditions included in the dataset of lidar measurements is plotted in Figs. 1–4. The heights of the mixing layer, MLH<sub>L</sub> and MLH'<sub>L</sub>, are obtained from different criteria using the wavelet covariance transform method. As shown in Fig. 1, the development of a convective mixing layer could clearly be observed, with a sharp decrease in aerosol backscatter between the mixing layer and the free atmosphere. MLH<sub>RS</sub> is also presented, accompanied by the evolution of MLH<sub>L</sub> and MLH'<sub>L</sub>.

As shown in Fig. 1a, MLH<sub>L</sub> and MLH'<sub>L</sub> increase at sunrise. On 2 March 2017, MLH<sub>L</sub> and MLH'<sub>L</sub> show an obvious diurnal cycle, with a maximum up to 1.0 km. During the evening of 3 March 2017 and the early morning of 4 March 2017, the aerosol layers present two visually obvious layers; MLH'<sub>L</sub> characterizes the first layer height and

MLH<sub>L</sub> represents the upper-layer top. On the next day, due to the existence of cloud, the MLH results are discrete. In the evening and early morning, MLH'<sub>L</sub> deviated from MLH<sub>L</sub> and approached MLH<sub>RS</sub>. As shown as the vertical profile from the lidar (RCS, wavelet) and radiosonde (RH, T) at 20:00 (LST) on 3 March 2017 (Fig. 1b), both MLH'<sub>L</sub> and MLH<sub>RS</sub> demonstrate 0.53 km, while MLH<sub>L</sub> shows 1.22 km. Figure 1c indicates that MLH'<sub>L</sub> (0.86 km) approaches MLH<sub>RS</sub> (0.61 km), although it is a little higher (0.25 km) and much lower than MLH<sub>L</sub> (1.62 km). In these cases, the result of MLH'<sub>L</sub>, indicated by the first local maximal aerosol gradient, agrees well with MLH<sub>RS</sub>. However, this is not always true, as shown in Fig. 1d. MLH'<sub>L</sub> (0.752 km) is much higher than MLH<sub>RS</sub> (0.243 km) but equal to MLH<sub>L</sub>. It is related to the fact that the stable layer height obtained from radiosonde in this case is out of the range of lidar detection (0.255 km), so MLH'<sub>L</sub> from lidar is disabled to determine the stable layer height.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Similar to Fig. 1, but for (a) 25–27 August 2014, and with vertical profiles for (b) 14:00 LST on 25 August 2014, (c) 14:00 LST on 27 August 2014 and (d) 20:00 LST on 27 August 2014.</figcaption></figure>

In the summertime (JJA) when the radiosonde is additionally launched at 14:00 LST to detect the convective boundary layer, it can provide the comparison between lidar and RS measurements in the afternoon. As shown in Fig. 2a, $MLH_L$ undergoes a rapid increase in the morning and reaches the peak in the afternoon, while $MLH'_L$ grows with a smaller magnitude. On 25 and 26 August 2014, the aerosol load is relatively low, and $MLH_L$ reaches the peak at around 3.0 km in the afternoon, while a lower $MLH_L$ peak on 27 August 2014 is observed with a high aerosol content. In the afternoon on these three days, $MLH_L$ is consistent with $MLH_{RS}$, while $MLH'_L$ is frequently under $MLH_{RS}$. The measurement of RS in the evening and early morning presents very low values at 0.2–0.3 km. The detailed information represented in Fig. 2b shows that $MLH_L$ is equal to $MLH_{RS}$, which reaches up to 2.95 km, while $MLH'_L$ is only at 1.24 km. Under clear convective conditions on 25 and 26 August 2014, when there are weak vertical gradients in the aerosol load indicated by a weak RCS, lidar can still obtain good MLH results compared to radiosonde, as shown in Fig. 2c and d with $MLH_L$ (2.25 and 2.18 km) approaching $MLH_{RS}$ (2.96 and 2.50 km). The slightly lower value of $MLH_{RS}$ is associated

with the fact that aerosol within the mixing layer needs some time to adjust to the thermal structure, and there is a delay in reaching the thermodynamic PBL height (Stull, 1988; Collaud Coen et al., 2014).

As shown in Fig. 3a, the peaks of the three days gradually increase, with values of 1.0, 2.0 and 2.5 km. Due to the high temporal variability of the distribution of aerosol, $MLH_L$ presents discontinuity on 17 June 2017. $MLH_L$ presents good evolution of the mixing layer height on 15 and 17 June 2017 compared to $MLH_{RS}$. $MLH'_L$ corresponds to $MLH_{RS}$ for most of the time in the morning and evening from 15 to 17 June 2017. But when the stable layer height from radiosonde is around or below 0.25 km, for example at 08:00 LST on 16 June 2017 with the RS measurement of 0.27 km (Fig. 3b), $MLH'_L$ misses the height of the SBL but points to the height of the residual layer (0.61 km). When the stable layer height is higher than 0.25 km, $MLH'_L$ (0.62 km) tends to approach to $MLH_{RS}$ (0.62 km) (Fig. 3d). However, in the afternoon $MLH_L$ was closer to $MLH_{RS}$ compared to $MLH'_L$ (Fig. 3c). As shown in Fig. 4a, $MLH_L$ ($MLH'_L$) presents the diurnal cycle with the maximum of 1.2 km,
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Similar to Fig. 1, but for (a) 15–17 June 2017, and with vertical profiles for (b) 08:00 LST on 16 June 2017, (c) 14:00 LST on 16 June 2017 and (d) 20:00 LST on 16 June 2017.</figcaption></figure>

while the next two days stay stable the whole day, and the height of the SBL is missed by $MLH'_L$ (Fig. 4b–d).

## 3.2 Intercomparison of different MLH approaches

A comparison of MLH estimated by lidar and radiosonde at 08:00 and 20:00 LST, is shown in Fig. 5. The same observation period of nearly 6 years (2013–2018) is considered, during which the data are continuous except for 2013. As shown in the histograms in Fig. 5, the total column is the annual relative frequency, and the different colors indicate the contribution of each season to the total. There is a wide discrepancy between $MLH_L$ and $MLH_{RS}$ at 08:00 and 20:00 LST. The frequency of MLH from a radiosonde lower than 0.25 km is nearly 35%, whereas there are no data for $MLH_L$ from lidar due to the limited detection range. This lower value mainly occurs in winter and autumn, when MLH tends to be lower (Tang et al., 2016). Specifically, the rate of $MLH_L$ from lidar lower than 0.5 km is nearly 18% and 12% at 08:00 and 20:00 LST, respectively, while the corresponding frequency for radiosonde is beyond 75% and 66%. The frequency of the larger $MLH_L$ value at 20:00 LST is bigger

than that of 08:00 LST from both lidar and radiosonde. It is reasonable that the residual layer has not yet collapsed entirely at 20:00 LST, while the CBL has not developed well in the early morning. As for the $MLH'_L$, its distribution trend is more similar to $MLH_{RS}$ than $MLH_L$ (see Fig. S2), and the correlation between $MLH'_L$ and $MLH_{RS}$ is a little higher than that between $MLH_L$ and $MLH_{RS}$, in spite of the fact that it is still not good (see Figs. S3 and S4). This indicates that $MLH'_L$ has the potential to determine the SBL height as a radiosonde does.

As for the seasonal variation of both lidar and RS measurements at 08:00 LST, the frequency of the larger $MLH_L$ value in summer is minimal, indicating that the summer MLH is lower than in other seasons. As for radiosonde, $MLH_L$ lower than 0.25 km is mostly distributed in winter, with a rate of around 15% for both 08:00 and 20:00 LST, and the frequency decreases rapidly when $MLH_L$ becomes higher than 0.25 km.

The poor agreement between MLH from lidar and $MLH_{RS}$ is also reported in the study of Su et al. (2019), in which it is shown that the correlation of the PBL height measurement between the lidar and radiosonde is 0.14 at 06:30 LST. The
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Similar to Fig. 1, but for (a) 16–18 October 2014, and with vertical profiles for (b) 08:00 LST on 17 October 2014, (c) 20:00 LST on 17 October 2014 and (d) 08:00 LST on 18 October 2014.</figcaption></figure>

significant scatter in the morning and evening is associated with the complicated structure of the boundary layer, as indicated by the existence of a stable boundary layer and residual layer (Su et al., 2019; Tang et al., 2016). In this study, regardless of MLH<sub>L</sub> and MLH'<sub>L</sub>, more than 35 % of the measurements of the SBL height are not within the scope of the lidar detection. Additionally, in the evening and early morning in some cases, a sufficiently clear variety cannot be found in the backscatter profile at the top of the SBL within the previously well-mixed layer (Russell et al., 1974; Seibert et al., 2000).

The comparison of MLH<sub>L</sub> and MLH<sub>RS</sub> at 14:00 LST in summer is presented in Fig. 6, with both mainly indicating the CBL height. MLH<sub>L</sub> shows very good agreement with MLH<sub>RS</sub>, with a correlation coefficient of 0.692 and a root mean square error (RMSE) of 0.573 km. It is noted that the slope of the linear fitting line is smaller than the 1 : 1 line, indicating that MLH<sub>RS</sub> tends to be larger than MLH<sub>L</sub> in the afternoon, which is consistent with the case study. Although the comparison is only for summer, it can be generally concluded that MLH<sub>L</sub> from lidar in the afternoon characterizes the CBL height with good accuracy. As shown in Fig. S4, the correlation of MLH'<sub>L</sub> and MLH<sub>RS</sub> at 14:00 LST is 0.330, and the value of MLH'<sub>L</sub> is generally lower than MLH<sub>RS</sub>, in-

dicating that it is overall unsuitable for MLH'<sub>L</sub> to describe the CBL height in the afternoon.

In fact, complete agreement is not expected between MLH<sub>L</sub> (MLH'<sub>L</sub>) derived from lidar and MLH<sub>RS</sub> from radiosonde for several reasons. First, the two systems measure different atmospheric parameters (aerosol for lidar and temperature, humidity, wind for radiosonde) with varying height resolution and accuracy, and these parameters are influenced in different ways by the processes occurring within the PBL (Seibert et al., 2000). Additionally, it is difficult to identify a clear upper boundary of the mixing layer because the measured parameter is actually not a fixed point but rather a transition layer between two atmospheric states (Stull, 1988; Garratt, 1992; Collaud Coen et al., 2014).

A dataset containing nearly 6 years of measurements in Beijing is used for assessment of the overall performance of the wavelet MLH algorithm (MLH<sub>L</sub> and MLH'<sub>L</sub>) with respect to the diurnal availability, as shown in Fig. 7. It can be seen that the expected shape representing the growth of a convective mixing layer is observed. Owing to solar heating of the surface, when the convective layer begins to rise due to upward convection in the early morning and the nocturnal residual layer tends to collapse, MLH<sub>L</sub> from lidar presents the minimal value. After that, MLH<sub>L</sub> grows contin-
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Comparison of the frequency distribution of all MLH<sub>L</sub> (2013–2018) retrieved from lidar and MLH<sub>RS</sub> from radiosonde with supplementary information on seasonal variation. The MLHs from (a) lidar and (b) radiosonde at 08:00 (LST) as well as (c) lidar and (d) radiosonde at 20:00 (LST) are presented. Noted that for presenting the detailed distribution, MLH<sub>L</sub> adds up to 20%, while MLH<sub>RS</sub> adds up to 45%.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Comparisons between MLH<sub>L</sub> in summer derived from lidar and MLH<sub>RS</sub> from radiosonde at 14:00 (LST). The red line indicates the linear fitting of 321 samples, while the black dashed line represents the 1 : 1 line.</figcaption></figure>

uously and reaches its maximum height around 15:00 LST, with the value of 1.449, similar to results found for Vienna (Lotteraner and Piringer, 2016) and Berlin (Geiß et al., 2017). The shaded areas indicate the temporal variability as calculated from the standard deviation of MLH<sub>L</sub> (MLH'<sub>L</sub>). It is on the order of 600 m for MLH<sub>L</sub> (300 m for MLH'<sub>L</sub>) during 09:00–15:00 LST, while the other period is on the order of 700 km (400 m for MLH'<sub>L</sub>). The larger standard deviation is attributed to the variability of the residual layer.

The diurnal cycles derived from MLH<sub>L</sub> match the RS results well in the afternoon, but they are larger than MLH<sub>RS</sub> in the early morning and evening. Contrarily, MLH'<sub>L</sub> tends to approach MLH<sub>RS</sub> in the early morning and evening but

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Diurnal cycles of the mixing layer height. The red line indicates the MLH<sub>L</sub> retrieved from lidar, and the blue line represents the MLH'<sub>L</sub> from lidar. The shaded areas show the standard deviation of MLH<sub>L</sub> and MLH'<sub>L</sub>. Purple triangles indicate the MLH<sub>RS</sub> averaged from routine RS data at 08:00 and 20:00 (LST) as well as from summer radiosonde at 14:00 (LST). The purple line indicates the standard deviation of MLH<sub>RS</sub>.</figcaption></figure>

stay far from MLH<sub>RS</sub> in the afternoon. The difference between MLH<sub>L</sub> (1.405 ± 0.675 km, mean ± standard deviation of the mean) and MLH<sub>RS</sub> (1.524 ± 0.582 km) at 14:00 LST is around 0.120 km. This is reasonable considering that RS data are acquired at 14:00 LST only in summertime when MLH<sub>L</sub> is usually larger throughout the year. However, the discrepancy between MLH'<sub>L</sub> (0.912 ± 0.315 km) and MLH<sub>RS</sub> (1.524 ± 0.582 km) at 14:00 is around 0.612 km. Actually, MLH<sub>L</sub> is nearly 0.46 km larger than MLH'<sub>L</sub> throughout the day, with a bigger standard deviation. As for the measurement at 08:00 LST, the difference between MLH<sub>L</sub> (1.196 ± 0.710 km) and MLH<sub>RS</sub> (0.434 ± 0.364 km) is around 0.762 km, which is larger than the difference between MLH'<sub>L</sub> (0.755 ± 0.334 km) and MLH<sub>RS</sub> (0.434 ± 0.364 km). The discrepancy between lidar and RS measurements at 20:00 LST is similar.

Overall, MLH'<sub>L</sub> can capture the high SBL height in the nocturnal time, when it is larger than 300 m. The stable layer height detected by MLH<sub>L</sub> in the nighttime is the layer in which ground-emitted atmospheric pollutants are trapped; it contributes to the assessment of the surface pollutant concentration when there are emissions in the nocturnal time using numerical models (Collaud Coen et al., 2014). Due to incomplete optical overlap, in some cases the point derived from MLH'<sub>L</sub> is the residual layer height rather than the low nocturnal SBL height. And in the daytime, MLH'<sub>L</sub> tends to be lower than the CBL height. In the studies of Mues et al. (2017) and Kotthaus and Grimmond (2018), the MLH in the daytime is usually assigned as the lowest layer detected by a ceilometer. Using higher-power lasers (CE-370) with an increasing signal-to-noise ratio (SNR) and a small gradient detected,
<table><caption>**Table 1.** Statistics of boundary layer height seasonal change.</caption><thead><tr><th>MLH<sub>L</sub> km<sup>-1</sup></th><th>Spring</th><th>Summer</th><th>Autumn</th><th>Winter</th></tr></thead><tbody><tr><th>Mean</th><td>1.409</td><td>1.261</td><td>1.297</td><td>1.228</td></tr><tr><th>Maximum</th><td>1.647</td><td>1.526</td><td>1.445</td><td>1.404</td></tr><tr><th>Minimum</th><td>1.126</td><td>0.932</td><td>1.117</td><td>1.098</td></tr></tbody></table>

the attribution of the lowest layer in the daytime may remain open, since the first local maximum gradient (MLH<sub>L</sub>) does not always correspond to the biggest local maximum (MLH<sub>L</sub>). Our study indicates that MLH<sub>L</sub> retrieves the consistent RL height during the night following the CBL diurnal maximal. The RL height corresponds to trapped atmospheric constituents discharged some hours before, which can be employed to convert column-mean optical depths into near-surface air quality information from remote sensing. And the SBL height provided by radiosonde at 08:00 and 20:00 LST can be considered complementary to the lidar approaches.

## 3.3 Climatology of MLH in Beijing

### 3.3.1 Seasonal variation

The seasonal mean diurnal cycle of the MLH<sub>L</sub> from lidar is shown in Fig. 8. An evident seasonal variation in the magnitude of the diurnal cycle is observed. As shown in Table 1 and Fig. 8, the smallest MLH<sub>L</sub> magnitude is found in winter, with the peak value of 1.404 ± 0.751 km at 15:00 LST, whereas spring demonstrates the maximum magnitude with 1.647 ± 0.754 km at 15:00 LST. The maximum in summer is 1.526 ± 0.581 km, and the maximum in autumn is 1.445 ± 0.837 km. From the all-day average of the four seasons, the averages in spring, summer, autumn and winter are 1.409, 1.261, 1.297 and 1.228 km, respectively. In summer, MLH<sub>L</sub> acquired by lidar at 14:00 and 15:00 LST is 1.430 and 1.507 km, respectively, while MLH<sub>RS</sub> at 14:00 LST is 1.524 km. The measurement of MLH<sub>L</sub> at 15:00 LST is closer to MLH<sub>RS</sub> at 14:00 LST. This is consistent with the case study in that it takes some time for aerosol to diffuse upward with the drive of thermal turbulence. As for the statistical variation, the values of autumn MLH<sub>L</sub> vary most at each hour with the bigger standard deviation, indicating great fluctuations in the long measurement period, while the variation of summer MLH<sub>L</sub> values for most hours is relatively stable.

It should be noted that summer exhibits the biggest amplitude of the diurnal variation of MLH<sub>L</sub>, with the deepest drop (0.93 km) increasing to the peak value of 1.51 km. Tang et al. (2016) indicate that the lower MLH<sub>L</sub> value for summer nights and early mornings is attributed to effect of the mountain plain wind. When the local mountain breeze from the northeast in the summer night superimposes the surface cooling, leading to an increase in the thickness of the inversion layer, the height of the mixed layer gradually decreases. After sunrise, with the drive of thermal turbulence,

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Seasonal variation of the diurnal cycles of MLH<sub>L</sub> retrieved from lidar (dot lines) and the standard deviation of MLH<sub>L</sub> (histograms). The red triangle indicates the MLH<sub>RS</sub> measured at 14:00 LST by radiosonde.</figcaption></figure>

the residual layer height observed by lidar is gradually replaced by a convective boundary layer height, with MLH<sub>L</sub> increasing rapidly; after 12:00 LT, the plain wind from the southwesterly direction gradually dominates. Previous studies have suggested that the seasonal variation in the MLH<sub>L</sub> may be associated with radiation flux (Stull, 1988; Kamp and McKendry, 2010; Munoz and Undurraga, 2010), which is consistent with our results. The observational data from Tang et al. (2016) indicate that radiation flux in spring is more than that in summer. The relatively low values in autumn and winter are likely related to the low radiation flux.

### 3.3.2 Interannual variation

Interannual variations of the MLH<sub>L</sub> diurnal cycle are investigated in Beijing from 2013 to 2018, as shown in Fig. 9. Diurnal variations of MLH<sub>L</sub> in different years all have the same patterns but with a different magnitude. Clearly, from 2013 to 2018, the values of the diurnal circle of MLH<sub>L</sub> increase year by year, including both the RL height at night and the CBL height in the daytime. Since the data for 2013 are mainly from winter and spring, the MLH<sub>L</sub> seems stable, unlike the amplitude in other years. As shown in Fig. 9 and Table 2, from 2014 to 2018, the MLH<sub>L</sub> all-day maximum values around 15:00 LST grow year by year, with values of 1.291 ± 0.646 km, 1.435 ± 0.755 km, 1.577 ± 0.739 km, 1.597 ± 0.701 km and 1.629 ± 0.751 km, respectively. The all-day averages of MLH<sub>L</sub> are 1.110, 1.216, 1.352, 1.391 and 1.502 km, respectively, also showing an increasing trend. This indicates that the volume available for the dispersion of pollutants is extending, which is beneficial to the mitigation of surface pollution. As shown in Fig. S5, from 2014 to 2018, the cumulative increase in the mean MLH<sub>L</sub> for the whole
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Interannual variation of the diurnal cycles of averaged MLH<sub>L</sub> retrieved from lidar (dot lines) and the standard deviation of MLH<sub>L</sub> (histograms). Due to the incomplete data for 2013, the MLH<sub>L</sub> data for 2013 are presented as a dotted line.</figcaption></figure>

<table><caption>**Table 2.** Statistics of boundary layer height interannual change.</caption><thead><tr><th>MLH<sub>L</sub> km<sup>-1</sup></th><th>2013</th><th>2014</th><th>2015</th><th>2016</th><th>2017</th><th>2018</th></tr></thead><tbody><tr><th>Mean</th><td>1.118</td><td>1.110</td><td>1.216</td><td>1.352</td><td>1.391</td><td>1.502</td></tr><tr><th>Maximum</th><td>1.207</td><td>1.291</td><td>1.435</td><td>1.577</td><td>1.597</td><td>1.629</td></tr><tr><th>Minimum</th><td>1.025</td><td>0.975</td><td>0.989</td><td>1.121</td><td>1.158</td><td>1.287</td></tr></tbody></table>

day was 0.392 km, and the total increase in the maximum is 0.338 km. As for the annual increase in MLH<sub>L</sub>, the average of the all-day increments in 2016 is the largest (0.136 km), while the average of the all-day increments in 2017 is the smallest (0.039 km).

In particular, the interannual variation of MLH<sub>L</sub> in the period from 10:00 to 15:00 LST is calculated, when the PBL is characterized by an obvious convective boundary layer. From 2014 to 2018, the average CBL height shows a significant increasing trend of 1.075, 1.212, 1.324, 1.351 and 1.533 km, respectively. The total increase in the average CBL height is 0.458 km. As for the annual increase in CBL height, in 2018 it is the largest (0.182 km), while the average increment for the whole day in 2017 is the smallest (0.027 km).

It is found that, based on measurements from 2014 to 2017, MLH<sub>L</sub> has a strong negative correlation with AOD ($R = -0.41$) and with relative humidity ($R = -0.21$), while MLH<sub>L</sub> presents a positive correlation with wind speed ($R = 0.43$) and shows no correlation with temperature (Fig. 10). As shown by the study of Wang et al. (2019), from 2014 to 2017 in Beijing, AOD and surface PM<sub>2.5</sub> have a tendency to decrease year by year, while MLH<sub>L</sub> increases gradually. The reduction of AOD and surface PM<sub>2.5</sub> is revealed to relate to pollution emission control in recent years (Zhang et al., 2019). Compared with 2016, relative humidity increased in 2017 and wind speed weakened, which is not beneficial for the development of the MLH but is consistent with the small

<figure><img src="image_12.png" /><figcaption>**Figure 10.** The correlation between MLH<sub>L</sub> and AOD, relative humidity, wind speed and temperature for measurements from 2014 to 2017.</figcaption></figure>

increase in MLH<sub>L</sub> (0.027 km) in 2017. In addition to the effects of meteorological conditions, the increase in MLH<sub>L</sub> benefits from the improvement of air quality in Beijing in recent years (Wang et al., 2019). Due to the scattering and absorbing of aerosol, solar radiation received from the ground decreases. It is thermal buoyancy generated from surface radiation that drives the PBL to develop. Thus, the development of the MLH is suppressed under a high aerosol load. Hence, with the relief of the radiation effect by aerosol during these years, the turbulence increases, thus leading to a larger PBL height (Ding et al., 2016; Li et al., 2017; Wang et al., 2019). The MLH affects the concentration of pollutants near the surface, while the total radiation of aerosol within the column atmosphere in turn influences the MLH.

## 3.4 Implications for surface pollution retrieval

The vertical structure of the ML is important for pollution concentrations at the surface due to its impact on the volume into which pollutants are mixed (Miao et al., 2018). Mues et al. (2017) reported that black carbon concentrations show a clear anticorrelation with MLH measurements. Hu et al. (2014) found a negative correlation between near-surface O<sub>3</sub> and MLH for seven cities in the North China Plain. In the study, as shown in Fig. 11, the correlation between MLH<sub>L</sub> and observed PM<sub>2.5</sub> data from the same observatory shows a high negative correlation ($R = -0.569$) with 4 years of measurement (2014–2017). Actually, the pollutant concentration near the surface is affected by the overall effect of local emissions and meteorological conditions, with variation in different spatiotemporal distributions. The MLH is just one of these influencing factors. Geiß et al. (2017) indicated that
<figure><img src="image_13.png" /><figcaption>**Figure 11.** The correlation between MLH<sub>L</sub> and PM<sub>2.5</sub>_OBS for measurements from 2014 to 2017.</figcaption></figure>

when the MLH and near-surface concentrations are linked, it
is necessary to take the locations and the details, i.e., meteo-
rological conditions and local sources, of the MLH retrieval
into account. In fact, all the data used in our study are from
the same observatory. And the PMRS model used to calcu-
late the surface PM<sub>2.5</sub> concentration includes the parameters
for emissions (AOD) and meteorological conditions (RH).

Due to the difference in the sources of the MLH, lidar
and radiosonde, the comparison of derived PM<sub>2.5</sub>_lidar and
PM<sub>2.5</sub>_RS with in situ observational PM<sub>2.5</sub> data at 08:00 LST
is presented in Fig. 12a and b. MLH from lidar shows rea-
sonably good performance for the retrieval of PM<sub>2.5</sub> in the
morning, with a correlation coefficient of 0.741 and RMSE
46.69 µg m<sup>-3</sup>. However, the calculated PM<sub>2.5</sub> from MLH<sub>RS</sub>
obviously overestimates the surface pollution, with lower
correlation coefficients and a larger standard deviation. The
large overestimation should be attributed to the underestima-
tion of the aerosol layer height. In the morning when the PBL
is not well developed, above MLH<sub>RS</sub> there is still a large
amount of aerosol, as seen in the lidar images in Figs. 1–4.
The discrepancy makes sense given the method using the ob-
served total amount of pollutant in the column atmosphere,
including emissions from the surface and the residual aerosol
from the day before. Therefore, MLH<sub>L</sub> from lidar, as a good
indicator of the aerosol layer height, is more suitable for
estimating surface air pollution from column-mean optical
depths.

As presented in Fig. 12c, the calculated PM<sub>2.5</sub>_lidar data
for the daytime period (08:00–17:00 LST) show higher cor-
relation (0.846) than those for only the early morning, and
the slightly larger RMSE (55.58 µg m<sup>-3</sup>) is associated with
the larger number of samples in the statistics. Considering the
uncertainty of the series of parameters used in the model, the
agreement between calculated PM<sub>2.5</sub>_lidar and in situ mea-
surements is reasonably good. Actually, the accuracy also
shows a diurnal cycle, with the peak of correlation coeffi-

<figure><img src="image_14.png" /><figcaption>**Figure 12.** Comparisons between observed PM<sub>2.5</sub> and PM<sub>2.5</sub> calculated from the PMRS model using (a) MLH<sub>L</sub> and (b) MLH<sub>RS</sub> at 08:00 (LST) as well as (c) MLH<sub>L</sub> for the period of 08:00–17:00 (LST). (d) The correlation coefficient and RMSE between observed PM<sub>2.5</sub> and PM<sub>2.5</sub> calculated from the PMRS model using MLH<sub>L</sub> for each hour from 08:00 to 17:00 LST.</figcaption></figure>

cients (0.927) at 14:00 LST (Fig. 12d). The correlations at
12:00, 13:00, 14:00 and 15:00 LST were 0.894, 0.922, 0.927
and 0.900, respectively. The higher accuracy may be due to
the completed mixing of the aerosol at noon and the vertical
distribution of the aerosol tending to be uniform. The corre-
lation between 08:00, 09:00 and 17:00 LST is less than 0.8,
which is related to the complex boundary layer structure in
the morning and at nightfall. It is difficult to achieve full mix-
ing of the aerosol in the stable boundary layer or the residual
layer. The smaller RMSE is related to the limited number
of samples. Therefore, the daily variation of the accuracy in
the calculated surface pollutant using MLH retrieval by lidar
varies with the daily variation of aerosol mixing uniformity
at different times during the daytime. Based on the observa-
tional data, PM<sub>2.5</sub> tends to peak in the morning and evening.
In contrast, afternoon usually witnessed lower mass concen-
trations due to the rapid vertical diffusion of aerosols (Guo et
al., 2016, 2017). Thus, MLH<sub>L</sub> from lidar can offer a signifi-
cant contribution to retrieving the diurnal circle of surface air
pollution.

4 Summary and conclusions

To acquire high-resolution observations of MLH diurnal
variation, a study using lidar was performed from Jan-
uary 2013 to December 2018 in the Beijing urban area. The
detection of the MLH based on two wavelet methods (MLH<sub>L</sub>
and MLH′~L~ applied to lidar observations was conducted. The two data results are compared with radiosonde as case studies and statistical forms. The temporal resolution (two or three measurements per day) of PBL detection by RS is not able to provide the mixing layer height diurnal cycle, regardless of its good precision. The MLH shows good performance for the convective layer height in the daytime and the residual layer height at night. MLH~′~L has the potential to describe the stable layer height at night sometimes, even though the capability is limited due to the highly incomplete overlap with the lidar used in the study. The stable layer height detected by MLH~′~L in the nighttime is the layer in which ground-emitted atmospheric pollutants are trapped; it contributes to the assessment of the surface pollutant concentration when there are emissions in the nocturnal time using numerical models. While the residual layer height corresponds to trapped atmospheric constituents discharged some hours before, it can be employed to convert column-mean optical depths into near-surface air quality information from remote sensing. And MLH~′~L does not always capture the convective layer height as MLH~L~ in the afternoon. Nevertheless, MLH~′~L could be useful as complementary information for the stable layer height in datasets of MLH~L~.

A nearly 6-year climatology for the MLH~L~ diurnal cycle is calculated for convective and stable conditions. It is true that the height of the mixing layer obtained by different approaches may be different. We focus on the temporal change in the aerosol layer height with a consistent method using a dataset of MLH~L~. The maximum MLH~L~ characteristics of seasonal change in Beijing indicate that it is low in winter (1.404±0.751 km) and autumn (1.445±0.837 km) and high in spring (1.647±0.754 km) and summer (1.526±0.581 km). A significant phenomenon is found from 2014 to 2018: the magnitude of the diurnal cycle of MLH~L~ increases year by year. The cumulative increase in the mean MLH~L~ for the whole day is 0.392 km, and the total increase in the maximum is 0.338 km. It may partly benefit from the improvement of air quality. As for converting the column optical depth to surface pollution, the calculated PM~2.5~ using MLH~L~ data from lidar shows better accuracy than that from radiosonde compared with observational PM~2.5~. Additionally, the accuracy of calculated PM~2.5~ using MLH~L~ shows a diurnal cycle in the daytime, with the peak at 14:00 LST. For the operational measurement of PBL height, the MLH from lidar has the capability to mark the diurnal circle of the mixing layer height and can be used as an effective parameter for the vertical distribution of aerosols, providing an important reference to obtain near-ground pollutant concentrations with remote sensing.

Actually, interpreting data from aerosol lidar is often not straightforward because the detected aerosol layers are not always the result of ongoing vertical mixing; they may originate from advective transport or past accumulation processes (Russell et al., 1974; Coulter, 1979; Baxter, 1991; Batchvarova et al., 1999; Tang et al., 2016). Each detection

method has good performances only for defined ML struc-
tures and under specific meteorological conditions. There-
fore, the combination of several methods and instruments
may contribute to characterizing the complete diurnal cy-
cle of the complex ML structure (Wiegner et al., 2006; de
Bruine et al., 2017; Morille et al., 2007; Kotthaus and Grim-
mond, 2018).

*Data availability.* The lidar data used in this study can be requested from the corresponding author (lizq@radi.ac.cn).

**Supplement.** The supplement related to this article is available on-line at: https://doi.org/10.5194/acp-20-8839-2020-supplement.

**Author contributions.** ZL conceived and designed the study. JG collected and processed the radiosonde data. YL, HX and PG collected the lidar data. HW improved the MLH retrieval algorithm from lidar. HW processed and analyzed the data and prepared the paper.

***Competing interests.*** The authors declare that they have no conflict of interest.

***Special issue statement.*** This article is part of the special issue "Satellite and ground-based remote sensing of aerosol optical, physical, and chemical properties over China". It is not associated with a conference.

***Acknowledgements.*** The authors acknowledge assistance with the acquisition of radiosonde data by Guo Jianping. The authors would like to acknowledge the China Meteorological Administration for obtaining the radiosonde data in Beijing. The authors thank Tang Guiqian for helpful discussions on algorithm improvement.

***Financial support.*** This research has been supported by the National Natural Science Foundation of China (grant nos. 41671367 and 41925019; program no. 305090306).

***Review statement.*** This paper was edited by Sachin S. Gunthe and reviewed by two anonymous referees.

**References**

Baars, H., Ansmann, A., Engelmann, R., and Althausen, D.: Continuous monitoring of the boundary-layer top with lidar, Atmos. Chem. Phys., 8, 7281–7296, https://doi.org/10.5194/acp-8-7281-2008, 2008.
Batchvarova, E., Cai, X., Gryning, S. E., and Steyn, D.: Modelling internal boundary layer development in a region with complex coastline, Bound.-Lay. Meteorol., 90, 1–20, https://doi.org/10.1023/a:1001751219627, 1999.

Baxter, R.: Determination of mixing heights from data collected during the 1985 SCCCAMP field program, J. Appl. Meteor., 30, 598–606, 1991.

Brooks, I. M.: Finding boundary layer top: application of a wavelet covariance transform to lidar backscatter profiles, J. Atmos. Ocean. Technol., 20, 1092–1105, https://doi.org/10.1175/1520-0426(2003)020<1092:FBLTAO>2.0.CO;2, 2003.

Beyrich, F. and Leps, J. P.: An operational mixing height data set from routine radiosoundings at Lindenberg: Methodology, Meteorol. Z., 21, 337–348, https://doi.org/10.1127/0941-2948/2012/0333, 2012.

Campbell, J. R., Welton, E. J., Spinhirne, J. D., Ji, Q., Tsay, S. C., Piketh, S. J., Barenbrug, M., and Holben, B. N.: Micropulse lidar observations of tropospheric aerosols over northeastern South Africa during the ARREX and SAFARI 2000 dry season experiments, J. Geophys. Res.-Atmos., 108, 8497, https://doi.org/10.1029/2002JD002563, 2003.

Che, H., Xia, X., Zhu, J., Li, Z., Dubovik, O., Holben, B., Goloub, P., Chen, H., Estelles, V., Cuevas-Agulló, E., Blarel, L., Wang, H., Zhao, H., Zhang, X., Wang, Y., Sun, J., Tao, R., Zhang, X., and Shi, G.: Column aerosol optical properties and aerosol radiative forcing during a serious haze-fog month over North China Plain in 2013 based on ground-based sunphotometer measurements, Atmos. Chem. Phys., 14, 2125–2138, https://doi.org/10.5194/acp-14-2125-2014, 2014.

Cimini, D., De Angelis, F., Dupont, J.-C., Pal, S., and Haeffelin, M.: Mixing layer height retrievals by multichannel microwave radiometer observations, Atmos. Meas. Tech., 6, 2941–2951, https://doi.org/10.5194/amt-6-2941-2013, 2013.

Cohn, S. A. and Angevine, W. M.: Boundary Layer Height and Entrainment Zone Thickness Measured by Lidars and Wind-Profiling Radars, J. Appl. Meteor., 39, 1233–1247, 2000.

Collaud Coen, M., Praz, C., Haefele, A., Ruffieux, D., Kaufmann, P., and Calpini, B.: Determination and climatology of the planetary boundary layer height above the Swiss plateau by in situ and remote sensing measurements as well as by the COSMO-2 model, Atmos. Chem. Phys., 14, 13205–13221, https://doi.org/10.5194/acp-14-13205-2014, 2014.

Coulter, R. L.: A comparison of three methods for measuring mixing layer height, J. Appl. Meteor., 18, 1495–1499, https://doi.org/10.1175/1520-0450(1979)018<1495:ACOTMF>2.0.CO;2, 1979.

Dandou, A., Bossioli, E., Tombrou, M., Sifakis, N., Paronis, D., Soulakellis, N., and Sarigiannis, D.: The importance of mixing height in characterizing pollution levels from aerosol optical thickness derived by satellite, Water Air Soil Poll., 2, 17–28, 2002.

Davis, K. J., Gamage, N., Hagelberg, C. R., Kiemle, C., Lenschow, D. H., and Sullivan, P. P.: An objective method for deriving atmospheric structure from airborne lidar observations, J. Atmos. Ocean. Technol., 17, 1455–1468, https://doi.org/10.1175/1520-0426(2000)017<1455:AOMFDA>2.0.CO;2, 2000.

de Bruine, M., Apituley, A., Donovan, D. P., Klein Baltink, H., and de Haij, M. J.: Pathfinder: applying graph theory to consistent tracking of daytime mixed layer height with backscatter lidar, At-

mos. Meas. Tech., 10, 1893–1909, https://doi.org/10.5194/amt-10-1893-2017, 2017.

De Haij, M., Wauben, W., and Baltink, H. K.: Determination of mixing layer height from ceilometer backscatter profiles, in: Proc. SPIE 6362, Remote Sensing of Clouds and the Atmosphere XI, Stockholm, Sweden, 11–14 September 2006, 63620R, https://doi.org/10.1117/12.691050, 2006.

Ding, A. J., Huang, X., Nie, W., Sun, J. N., Kerminen, V. M., Petäjä, T., Su, H., Cheng, Y. F., Yang, X. Q., Wang, M. H., Chi, X. G., Wang, J. P., Virkkula, A., Guo, W. D., Yuan, J., Wang, S. Y., Zhang, R. J., Wu Y. F., Song, Y., Zhu, T., Zilitinkevich, S., Kulmala, M., and Fu, C. B.: Enhanced haze pollution by black carbon in megacities in China, Geophys. Res. Lett., 43, 2873–2879, https://doi.org/10.1002/2016GL067745, 2016.

Dubovik, O., Smirnov, A., Holben, B. N., King, M. D., Kaufman, Y. J., Eck, T. F., and Slutsker, I.: Accuracy assessments of aerosol optical properties retrieved from Aerosol Robotic Network (AERONET) Sun and sky radiance measurements, J. Geophys. Res.-Atmos., 105, 9791–9806, https://doi.org/10.1029/2000JD900040, 2000.

Emeis, S., Jahn, C., Munkel, C., Munsterer, C., and Schäfer, K.: Multiple atmospheric layering and mixing-layer height in the Inn valley observed by remote sensing, Meteorol. Z., 16, 415–424, 10.1127/0941-2948/2007/0203, 2007.

Emeis, S., Schäfer, K., and Münkel, C.: Surface-based remote sensing of the mixing-layer height – a review, Meteorol. Z., 17, 621–630, https://doi.org/10.1127/0941-2948/2008/0312, 2008.

Eresmaa, N., Härkönen, J., Joffre, S. M., Schultz, D. M., Karpinen, A., and Kukkonen, J.: A three-step method for estimating the mixing height using ceilometer data from the Helsinki test bed, J. Appl. Meteorol. Clim., 51, 2172–2187, https://doi.org/10.1175/JAMC-D-12-058.1, 2012.

Gan, C. M., Wu, Y. H., Madhavan, B. L., Gross, B., and Moshary, F.: Application of active optical sensors to probe the vertical structure of the urban boundary layer and assess anomalies in air quality model PM₂.₅ forecasts, Atmos. Environ. 45, 6613–6621, 2011.

Garratt, J. R.: The Atmospheric Boundary Layer, Cambridge University Press, England, 316 pp., 1992.

Geiß, A., Wiegner, M., Bonn, B., Schäfer, K., Forkel, R., von Schneidemesser, E., Münkel, C., Chan, K. L., and Nothard, R.: Mixing layer height as an indicator for urban air quality?, Atmos. Meas. Tech., 10, 2969–2988, https://doi.org/10.5194/amt-10-2969-2017, 2017.

Guo, J., Miao, Y., Zhang, Y., Liu, H., Li, Z., Zhang, W., He, J., Lou, M., Yan, Y., Bian, L., and Zhai, P.: The climatology of planetary boundary layer height in China derived from radiosonde and reanalysis data, Atmos. Chem. Phys., 16, 13309–13319, https://doi.org/10.5194/acp-16-13309-2016, 2016.

Guo, J., Xia F., Zhang Y., Liu, H., Li, J., Lou, M., He, J., Yan, Y., Wang, F., Min, M., and Zhai, P.: Impact of diurnal variability and meteorological factors on the PM₂.₅ – AOD relationship: Implications for PM₂.₅ remote sensing, Environ. Pollut., 221, 94–104, https://doi.org/10.1016/j.envpol.2016.11.043, 2017.

Guo, J., Li, Y., Cohen, J. B., Li, J., Chen, D., Xu, H., Liu, L., Yin, J., Hu, K., and Zhai, P.: Shift in the temporal trend of boundary layer height in China using long-term (1979–2016) radiosonde data, Geophys. Res. Lett., 46, 6080–6089, https://doi.org/10.1029/2019GL082666, 2019.
Haeffelin, M., Angelini, F., Morille, Y., Martucci, G., Frey, S., Gobbi, G. P., Lolli, S., O'Dowd, C. D., Sauvage, L., Xueref-Remy, I., Wastine, B., and Feist, D. G.: Evaluation of mixing-height retrievals from automatic profiling lidars and ceilometers in view of future integrated networks in Europe, Bound.-Lay. Meteorol. 143, 49–75, https://doi.org/10.1007/s10546-011-9643-z, 2012.

Holben, B. N., Eck, T. F., Slutsker, I., Tanré, D., Buis, J. P., Setzer, A., Vermote, E., Reagan, J. A., Kaufman, Y. J., and Nakajima, T.: AERONET – A federated instrument network and data archive for aerosol characterization, Remote Sens. Environ., 66, 1–16, https://doi.org/10.1016/s0034-4257(98)00031-5,1998.

Hu, X., Ma, Z., Lin, W., Zhang, H., Hu, J., Wang, Y., Xu, X., Fuentes, J. D., and Xue, M.: Impact of the Loess Plateau on the atmospheric boundary layer structure and air quality in the North China Plain: a case study, Sci. Total Environ., 499, 228–237, https://doi.org/10.1016/j.scitotenv.2014.08.053, 2014.

Kamp, D. and McKendry, I.: Diurnal and seasonal trends in convec-tive mixed-layer heights estimated from two years of continuous ceilometer observations in Vancouver, BC, Bound.-Lay. Meteorol., 137, 459–475, https://doi.org/10.1007/s10546-010-9535-7, 2010.

Kotthaus, S. and Grimmond, C. S. B.: Atmospheric Boundary Layer Characteristics from Ceilometer Measurements Part 1: A new method to track mixed layer height and classify clouds, Q. J. Roy. Meteor. Soc., 144, 1525–1538, https://doi.org/10.1002/qj.3299, 2018.

Li, Z. Q., Li, D. H., Li, K. T., Xu, H., Chen, X. F., Chen, C., Xie, Y. S., Li, L., Li, L., and Li, W.: Sun-sky radiometer observation network with the extension of multi-wavelength polarization measurements, J. Remot. Sens. 19, 495–519, 2015.

Li, Z., Zhang, Y., Shao, J., Li, B., Hong, J., Liu, D., Li, D., Wei, P., Li, W., Li, L., Zhang, F., Guo, J., Deng, Q., Wang, B., Cui, C., Zhang, W., Wang., Z., Lv, Y., and Qie, L.: Remote sensing of atmospheric particulate mass of dry PM<sub>2.5</sub> near the ground: Method validation using ground-based measurements, Remote Sens. Environ., 173, 59–68, https://doi.org/10.1016/j.rse.2015.11.019, 2016.

Li, Z. Q., Guo, J. P., Ding, A. J., Liao, H., Liu, J. J., Sun, Y. L., Wang, T., Xue, H., Zhang, H., and Zhu, B.: Aerosol and boundary-layer interactions and impact on air quality, Natl. Sci. Rev., 4, 810–833, https://doi.org/10.1093/nsr/nwx117, 2017.

Liu, S. and Liang, X. Z.: Observed diurnal cycle climatology of planetary boundary layer height, J. Climate, 23, 5790–5809, https://doi.org/10.1175/2010JCLI3552.1, 2010.

Lotteraner, C. and Piringer, M.: Mixing-Height Time Series from Operational Ceilometer Aerosol-Layer Heights, Bound.-Lay. Meteorol., 161, 265–287, https://doi.org/10.1007/s10546-016-0169-2, 2016.

Mahrt, L., Sun, J., Blumen, W., Delany, T., and Oncley, S.: Noc-turnal boundary-layer regimes, Bound.-Lay. Meteorol., 88, 255–278, 1998.

Martucci, G., Matthey, R., Mitev, V., and Richner, H.: Com-parison between backscatter lidar and radiosonde measure-ments of the diurnal and nocturnal stratification in the lower troposphere, J. Atmos. Ocean. Technol., 24, 1231–1244, https://doi.org/10.1175/JTECH2036.1, 2007.

Martucci, G., Milroy, C. and O'Dowd, C. D.: Detection of cloud-base height using Jenoptik CHM15K and Vaisala

CL31 ceilometers, J. Atmos. Ocean. Technol., 27, 305–318, https://doi.org/10.1175/2009JTECHA1326.1, 2010.

Medeiros, B., Hall, A. and Stevens, B.: What controls the mean depth of the PBL?, J. Climate, 18, 3157–3172, https://doi.org/10.1175/JCLI3417.1, 2005.

Miao Y., Liu S., Guo J., Huang, S., Yan, Y., and Lou, M.: Unraveling the relationships between boundary layer height and PM<sub>2.5</sub> pollution in China based on four-year ra-diosonde measurements, Environ. Pollut., 243, 1186–1195, https://doi.org/10.1016/j.envpol.2018.09.070, 2018.

Milroy, C., Martucci, G., Lolli, S., Loaec, S., Sauvage, L., Xueref-Remy, I., Lavric, J. V., Ciais, P., Feist, D. G., Biavati, G., and O'Dowd, C. D.: An Assessment of Pseudo-Operational Ground-Based Light Detection and Ranging Sensors to Determine the Boundary-Layer Structure in the Coastal Atmosphere, Adv. Meteorol., 2012, 929080, https://doi.org/10.1155/2012/929080, 2012.

Monks, P. S., Granier, C., Fuzzi, S., Stohl, A., Williams, M. L., Akimoto, H., Amann, M., Baklanov, A., Baltensperger, U., Bey, I., Blake, N., Blake, R. S., Carslaw, K., Cooper, O. R., Dentener, F., Fowler, D., Fragkou, E., Frost, G. J., and von Glasow, R.: Atmospheric composition change – global and regional air quality, Atmos. Environ. 43, 5268–5350, https://doi.org/10.1016/j.atmosenv.2009.08.021, 2009.

Morille, Y., Haeffelin, M., Drobinski, P., and Pelon, J.: STRAT: An automated algorithm to retrieve the vertical structure of the atmo-sphere from single-channel lidar data, J. Technol., 24, 761–775, https://doi.org/10.1175/JTECH2008.1, 2007.

Mues, A., Rupakheti, M., Münkel, C., Lauer, A., Bozem, H., Hoor, P., Butler, T., and Lawrence, M. G.: Investigation of the mix-ing layer height derived from ceilometer measurements in the Kathmandu Valley and implications for local air quality, At-mos. Chem. Phys., 17, 8157–8176, https://doi.org/10.5194/acp-17-8157-2017, 2017.

Munoz, R. C. and Undurraga, A. A.: Daytime mixed layer over the Santiago Basin: Description of two years of observations with a lidar ceilometer, J. Appl. Meteorol. Clim., 49, 1728–1741, https://doi.org/10.1175/2010JAMC2347.1, 2010.

Neff, W. D. and Coulter, R. L.: Acoustic remote sensing, in: Probing the Atmospheric Boundary Layer, edited by: Lenschow, D. H., Amer. Meteor. Soc. Boston, MA, 201–239, 1986.

Piringer, M., Joffre, S., Baklanov, A., Christen, A., Deserti, M., Ridder, K., Emeis, S., Mestayer, P., Tombrou, M., Mid-dleton, D., Baumann-Stanzer, K., Dandou, A., Karppinen, A. and Burzynski, J.: The surface energy balance and the mixing height in urban areas – activities and recommenda-tions of COST-Action 715, Bound.-Lay. Meteorol., 124, 3–24, https://doi.org/10.1007/s10546-007-9170-0, 2007.

Poltera, Y., Martucci, G., Collaud Coen, M., Hervo, M., Emmeneger, L., Henne, S., Brunner, D., and Haefele, A.: Pathfinder- TURB: an automatic boundary layer algorithm. Development, validation and application to study the impact on in situ mea-surements at the Jungfraujoch, Atmos. Chem. Phys., 17, 10051–10070, https://doi.org/10.5194/acp-17-10051-2017, 2017.

Poulos, G. S., Blumen, W., Fritts, D. C., Lundquist, J. K., Sun, J., Burns, S. P., Nappo, C., Banta, R., Newsom, R., Cuxart, J., Ter-radellas, E., Balsley, B., and Jensen, M.: CASES-99: A compre-hensive investigation of the stable nocturnal boundary layer, B.
Am. Meteorol. Soc., 83, 555-581, https://doi.org/10.1175/1520-0477(2002)083<0555:CACIOT>2.3.CO;2, 2002.

Quan, J., Gao, Y., Zhang, Q., Tie, X., Cao, J., Han, S., Meng, J., Chen, P., and Zhao, D.: Evolution of planetary boundary layer under different weather conditions, and its impact on aerosol concentrations, Particuology, 11, 34-40, https://doi.org/10.1016/j.partic.2012.04.005, 2013.

Russell, P. B., Uthe, E. E., Ludwig, F. L., and Shaw, N. A.: A comparison of atmospheric structure as observed with monostatic acoustic sounder and lidar techniques, J. Geophys. Res.-Atmos., 79, 5555-5566, 1974.

Salmond, J. A. and McKendry, J. A.: A review of turbulence in the very stable nocturnal boundary layer and its implications for air quality, Prog. Phys. Geog., 29, 171-188, 2005.

Sawyer, V. and Li, Z.: Detection, variations and intercomparison of the planetary boundary layer depth from radiosonde, lidar and infrared spectrometer, Atmos. Environ., 79, 518-528, https://doi.org/10.1016/j.atmosenv.2013.07.019, 2013.

Schafer, K., Emeis, S., Hoffmann, H., and Jahn, C.: Influence of mixing layer height upon air pollution in urban and sub-urban area, Meteorol. Z., 15, 647-658, https://doi.org/10.1127/0941-2948/2006/0116, 2006.

Seibert, P., Beyrich, F., Gryning, S. E., Joffre, S., Rasmussen, A., and Tercier, P.: Review and intercomparison of operational methods for the determination of the mixing height, Atmos. Environ., 34, 1001-1027, https://doi.org/10.1016/S1474-8177(02)80023-7, 2000.

Seidel, D. J., Ao, C. O., and Li, K.: Estimating climatological planetary boundary layer heights from radiosonde observations: Comparison of methods and uncertainty analysis, J. Geophys. Res.-Atmos., 115, D16113, https://doi.org/10.1029/2009JD013680, 2010.

Seidel, D. J., Zhang, Y., Beljaars, A., Golaz, J. C., Jacobson, A. R., and Medeiros, B.: Climatology of the planetary boundary layer over the continental United States and Europe, J. Geophys. Res.-Atmos., 117, D17106, https://doi.org/10.1029/2012JD018143, 2012.

Sifakis, N., Soulakellis, N., and Paronis, D.: Quantitative mapping of air pollution density using Earth observations: a new processing method and application to an urban area, Int. J. Rem. Sens., 19, 3289-3300, 1998.

Singh, N., Solanki, R., Ojha, N., Janssen, R. H. H., Pozzer, A., and Dhaka, S. K.: Boundary layer evolution over the central Himalayas from radio wind profiler and model simulations, Atmos. Chem. Phys., 16, 10559-10572, https://doi.org/10.5194/acp-16-10559-2016, 2016.

Stull, R. B.: An Introduction to Boundary Layer Meteorology, Kluwer Academic, 666 pp., 1988.

Su, T., Li, J., Li, C., Xiang, P., A. K.-H. Lau, Guo, J., Yang, D., and Miao, Y.: An intercomparison of long-term planetary boundary layer heights retrieved from CALIPSO, ground-based lidar, and radiosonde measurements over Hong Kong, J. Geophys. Res.-Atmos., 122, 3929-3943, https://doi.org/10.1002/2016JD025937, 2017.

Su, T., Li, Z., and Kahn, R.: A new method to retrieve the diurnal variability of planetary boundary layer height from lidar under different thermodynamic stability conditions, Remote Sens. Environ., 237, 111519, https://doi.org/10.1016/j.rse.2019.111519, 2019.

Tang, G., Zhang, J., Zhu, X., Song, T., Münkel, C., Hu, B., Schäfer, K., Liu, Z., Zhang, J., Wang, L., Xin, J., Suppan, P., and Wang, Y.: Mixing layer height and its implications for air pollution over Beijing, China, Atmos. Chem. Phys., 16, 2459-2475, https://doi.org/10.5194/acp-16-2459-2016, 2016.

Wang, H. F., Li, Z. Q., Lv, Y., Xu, H., Li, K., Li, D., Hou, W., Zheng, F., Wei, Y., and Ge, B.: Observational study of aerosol-induced impact on planetary boundary layer based on lidar and sunphotometer in Beijing, Environ. Pollut., 252, 897-906, https://doi.org/10.1016/j.envpol.2019.05.070, 2019.

Wiegner, M., Emeis, S., Freudenthaler, V., Heese, B., Junkermann, W., Münkel, C., Schäfer, K., Seefeldner, M., and Vogt, S.: Mixing Layer Height over Munich, Germany: Variability and comparisons of different methodologies, J. Geophys. Res.-Atmos., 111, D13201, https://doi.org/10.1029/2005JD006593, 2006.

Yan, H., Li, Z., Huang, J., Cribb, M., and Liu, J.: Long-term aerosol-mediated changes in cloud radiative forcing of deep clouds at the top and bottom of the atmosphere over the Southern Great Plains, Atmos. Chem. Phys., 14, 7113-7124, https://doi.org/10.5194/acp-14-7113-2014, 2014.

Yan, X., Shi, W., Li, Z., Li, Z., Luo, N., Zhao, W., Wang, H., and Yu, X.: Satellite-based PM₂.₅ estimation using fine-mode aerosol optical thickness over China, Atmos. Environ., 170, 290-302 https://doi.org/10.1016/j.atmosenv.2017.09.023, 2017.

Yang, H., Liu, W., Lu, Y., Xie, P., Xu, L., Zhao, X., Yu, T., and Yu, J.: PBL observations by lidar at Peking, Opt. Technol., 31, 221-226, 2005.

Zhang, Q., Quan, J., Tie, X., Li, X., Liu, Q., Gao, Y., and Zhao, D.: Effects of meteorology and secondary particle formation on visibility during heavy haze events in Beijing, China, Sci. Total Environ., 502, 578-584, https://doi.org/10.1016/j.scitotenv.2014.09.079, 2015.

Zhang, Q., Zheng, Y., Tong, D., Tong, D., Shao, M., Wang, S., Zhang, Y., Xu, X., Wang, J., He, H., Liu, W., Ding, Y., Lei, Li., J., Wang, Z., Zhang, X., Wang, Y., Cheng, J., Liu, Y., Shi, Q., Yan, L., Geng, G., Hong, C., Li, M., Liu, F., Zheng, B., Cao, J., Ding, A., Gao, J., Fu, Q., Huo, J., Liu, B., Liu, Z., Yang, F., He, K., and Hao, J.: Drivers of improved PM₂.₅ air quality in China from 2013 to 2017, P. Natl. Acad. Sci. USA, 116, 24463-24469, https://doi.org/10.1073/pnas.1907956116, 2019.

Zhang, W. C., Guo, J. P., Miao, Y. C., Liu, H., Song, Y., Fang, Z., He, J., Lou, M. Y., Yan, Y., Li, Y., and Zhai, P. M.: On the summertime planetary boundary layer with different thermodynamic stability in China: a radiosonde perspective, J. Climate, 31, 1451-1465, https://doi.org/10.1175/JCLI-D-17-0231.1, 2018.

Zhang, Y. and Li, Z.: Remote sensing of atmospheric fine particulate matter (PM₂.₅) mass concentration near the ground from satellite observation, Remote Sens. Environ., 160, 252-262, https://doi.org/10.1016/j.rse.2015.02.005, 2015.