# Dust Aerosols Detected Using a Ground-Based Polarization Lidar and CALIPSO over Wuhan (30.5°N, 114.4°E), China

Authors: Yun He; Fan Yi
Journal: Advances in Meteorology
 (2015)

---

## Abstract

The vertical distribution, horizontal range, and optical properties of Asian dust were obtained using a ground-based depolarization lidar and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) over a two-year measurement period (2010–2012) in Wuhan (30.5°N, 114.4°E), China. The depolarization lidar registered 13 dust events, most of which occurred in the spring (5 events) and winter (6 events). The dust layers occurred at heights of approximately 1.4–3.5 km. The horizontal ranges of the dust plumes were approximately 750–2400 km, based on the CALIPSO data. The average volume depolarization ratio (δ), particle depolarization ratio (δ
p), extinction and optical depth (AOD) of the dust layers were 0.12, 0.22, 0.19 km−1, and 0.32, respectively. The dust layers observed in the winter occurred at a lower height and had larger mean extinction and AOD, and smaller mean δ and δ
p than the spring dust layers. These wintertime features may result from a lower troposphere temperature inversion, the mixing of local aerosols, and hygroscopic growth under suitable relative humidity conditions. A dust event in April 2011 spanned 9 days. Compared with the observations at other sites, the dust layers over Wuhan exhibited more turbid along with suppressed nonspherical particle shape.


---

## Body

## 1. Introduction

Water loss and soil erosion induce dust emissions, which affect the large-scale atmospheric environment. Dust aerosols primarily originate in arid and semiarid areas. Dust particles entrained from desert surfaces by strong air convection are transported with the wind in the atmosphere [1–5].There are three vast desert areas on Earth: the Sahara Desert, the East Asian Desert, and the Arabian Peninsula Desert. These deserts contribute 40% of the total aerosol loading in the atmosphere [6]. Specifically, Asian dust contributes 800 Tg of dust aerosols to the atmosphere which occupies 40–80% of the total dust emission (~1-2 Pg) [7] every year, and 30% settles back onto the desert areas of Asia, 20% is transported on an intercontinental scale, and the remaining 50% is transported in the atmosphere over the Pacific or farther [8]. Asian dust affects areas far from the source regions via long-range transport [9, 10]. During this transport, dust particles undergo complex physical and chemical processes and play an important role in the atmospheric radiation budget [4, 11, 12].Lidar, due to its high temporal and spatial resolution, is an important tool for measuring the vertical distribution and optical properties of dust, and ground-based lidar has been widely used in recent studies to observe Asian dust [13–17]. Sakai et al. [16] observed a dust event on 23 April 2001 in Tsukuba. The aerosol depolarization ratios at a height of approximately 5 km exceeded 0.2, and the optical depth between heights of 4 km and 7 km was 0.18. Kim et al. [18] observed an enhanced aerosol extinction of ~0.015 km−1 at a height of 3–5 km in Seoul, Korea, in the spring seasons during the period 1997–2004, which indicates the effect of the dust. Since its launch in 2006, CALIPSO lidar has been widely used to obtain both local and global perspectives of dust optical properties and horizontal ranges [19–24]. Using CALIPSO data, Kim et al. [22] observed a dust plume that developed over South Korea on 29–31 May 2008, finding that a majority of the plume traveled 250–300 km across South Korea. Also using CALIPSO data, Huang et al. [20] observed that dust layers over the Tibetan Plateau during the summer of 2006 appeared frequently at heights of approximately 4–7 km and had a mean volume depolarization ratio and color ratio of 0.21 and 0.83, respectively. Most lidar observations have focused on northwest and north China, Japan, and Korea; however, few such studies have focused on southeast China [23, 25, 26]. Southeast China, which is not located in the path of the prevailing westerlies, may exhibit certain unique dust features after transport. Long-term joint observations using ground-based lidar and CALIPSO lidar in Wuhan provide an opportunity for studying the vertical distribution, horizontal range, and optical properties of dust in this region of China. In this paper, polarization lidar was used to identify nonspherical dust particles and measure the dust’s optical properties, and the CALIPSO data were used to assess the horizontal range of the dust plume.The paper is organized as follows. The current situation of dust research is first introduced. In Section2, the instruments and data processing are described. In Section 3, the results of the dust observations over Wuhan are presented with two case studies and statistical analysis. In Section 4, the differences between the dust characteristics in the spring and winter over Wuhan are discussed. Special dust optical properties are also discussed by making comparisons with other sites in this section. Finally, the conclusions are provided in Section 5.


## 2. Instruments and Data

### 2.1. Polarization Lidar System

Wuhan (30.5°N, 114.4°E) is located in southeast China. The polarization lidar system that was used is installed at Wuhan University. The lidar system parameters are presented in Table1. The polarization lidar transmitter used an Nd:YAG laser to produce an emission of ~25–40 mJ per pulse at 532 nm with a repetition rate of 20 Hz and a pulse width of 6 ns. The polarization purity of the transmitting laser was approximately 1/10000. On the receiver side, the backscattered photons were collected using a 300 mm Cassegrain telescope that was followed by a field-stop iris and a collimating lens. Double cubic polarizing beam splitters was employed to split the beam into two independent polarization channels in parallel and perpendicular directions relative to the outgoing laser beam. The filters for the two channels have a band width of ~10.6 cm−1 (0.3 nm at 532 nm). The two closest rotational Raman lines (N2) are ±11.9 cm−1 away from the elastic Cabannes line. Furthermore the pure rotational Raman lines are ~3 orders of magnitude weaker in intensity than the elastic Cabannes line. Thus the cross-talk of pure rotation Raman line was negligible [27]. The cross-talk between the two polarization channels is considerably suppressed by using double cubic polarizing beam splitters at each channel. The beams of the two channels were then detected using two photomultiplier tubes (PMT). The output pulses from the PMT were gathered using a signal acquisition instrument (Licel).Table 1
Specification of the depolarization lidar system at Wuhan University.
Transmitter
Receiver
Laser
Continuum Inlite II-20
Telescope
300 mm Cassegrain
Wavelength
532 nm
Diameter
300 mm
Energy/pulse
25–40 mJ
Field of view
1 mrad
Repetition rate
20 Hz
PMT
Hamamatsu 5783P
Beam divergence
6 ns
Acquisition instrument
Licel TR40-160The spatial and temporal resolutions of the lidar system were 30 m and 1 minute, respectively. The lower lidar observation limit was 0.3 km, which was determined based on the overlap of the laser and the field of view of the telescope. Our calibration for the polarization lidar was made by observing a fully depolarized light from thick-cloud-covered sky [28]. The resulting gain ratio (k) between the two channels was 0.065. The aerosol extinction was calculated using the combined parallel and perpendicular channels [29]. The volume depolarization ratio δ was calculated using the equations
(1)
δ
=
k
I
⊥
I
∥
,
where I
∥ and I
⊥ are the signals from the parallel and perpendicular channels, respectively. δ is the volume depolarization ratio caused by molecules and aerosols.The particle depolarization ratioδ
p was derived from the total-to-Rayleigh backscattering ratio BR and the volume depolarization ratio δ from the following equation [30]:(2)
δ
p
z
=
δ
z
BR
z
+
BR
z
δ
m
-
δ
m
-
δ
m
BR
z
-
1
+
BR
z
δ
m
-
δ
z
,
where δ
m is the molecular depolarization ratio for the lidar system. The value of molecular depolarization ratio measured by lidar depends on the wavelength and spectrum width of the detected light. The band widths of two filters used in our lidar system are both 0.3 nm. On the basis of calculation by Behrendt and Nakamura [27], the molecular depolarization ratio computed for Cabannes line is slightly larger than the theoretical value of 0.00363 for linear polarized incident light and 180-degree backscattering (see Figure 4 of [27]). Noting that the cross-talk of pure rotational Raman is negligible and the molecular depolarization ratio does not change with the variation of atmospheric temperature, the value of 0.4% was used here. The value of BR is defined as BR
z
=
1
+
β
a
(
z
)
/
β
m
(
z
) [31], where β
m and β
a are the molecular and aerosol backscattering coefficients, respectively.The aerosol extinction was calculated from the Mie backscattering signal using the method of Fernald [32]. Simultaneous measurement with lidar and sun photometer indicates that the average lidar ratio in the troposphere is ~50 sr (via comparing the aerosol optical depths derived from the lidar (at 532 nm) and the sun photometer (at 500 nm) [33, 34]). In practice, Murayama et al. [35] derived the extinction coefficient of aerosols by using the algorithm of Fernald [32], the lidar ratio of 50 sr was chosen to reproduce the optical depth derived from the sun photometer. The lidar ratio of 50 sr was also used in Asian dust observation in Taiwan [36] and Hefei, China [26]. Noh et al. [37] measured the lidar ratio of 51
±
6 sr for Asian dust with a Raman lidar at Gwangju, Korea, in spring. Therefore, the assumed lidar ratio of 50 sr is reasonable in this study. The molecular extinction was calculated using the U.S. Standard Atmosphere Model temperature and pressure data (U.S. GPO, Washington, D.C., 1976). The aerosol optical depth (AOD) and the centroid height (z
c) of a dust layer are given by
(3)
AOD
=
∫
Z
b
Z
t
α
a
(
z
)
d
z
,
z
c
=
∫
Z
b
Z
t
z
·
α
a
z
·
d
z
∫
Z
b
Z
t
α
a
z
·
d
z
,
where z
t and z
b are the dust layer top and base heights, respectively, and α
a is the aerosol extinction coefficient.


### 2.2. Satellite Data

The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument is carried on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite and used to acquire vertical profiles of elastic backscatter at 532 nm and 1064 nm near nadir during both daytime and nighttime. The nearest CALIOP orbits pass over Wuhan at approximately 02:00 LST and 14:00 LST every day. Profiles of linear depolarization at 532 nm are provided. The ratio of 1064 nm backscattering to 532 nm backscattering, which is defined as the “color ratio,” is also provided. This ratio represents the aerosol particle size. The depolarization ratio of dust is large because of its nonspherical shape, which is a good indicator for separating dust from other aerosols [38]. Liu et al. [23] derived the probability density functions of volume depolarization ratio (VDR) for four types of aerosol (desert dust, biomass burning, continental, and maritime). The maritime and continental aerosols had a VDR distribution peaking at zero. The biomass burning aerosols had a nonzero peak at ~0.02. The VDR distribution of the dust aerosols was centered at ~0.17. Based on these probability density functions, when a VDR threshold is set to 0.06, only ~1.4% of dust aerosol particles were left out and only ~2.3% of non-dust aerosol particles were mistakenly classified as dust aerosol. Hence, in practical lidar data analysis (CALIPSO), the volume depolarization ratio value of 0.04 or 0.06 was employed as a threshold value to distinguish between dust and nondust aerosols by Liu et al. [39] and Huang et al. [21], respectively.In this study, the CALIPSO level 1B data product was used to validate the performance of the polarization lidar and study the horizontal range of the dust plumes. The level 2 5 km aerosol profile data was used to exhibit the vertical feature mask and aerosol subtypes mask. In this product, bits 1–3 of theAtmospheric_Volume_Description represent the feature mask (1: clean air, 2: cloud, and 3: aerosol) and bits of 10–12 of theAtmospheric_Volume_Description represent the aerosol subtypes (0: not determined, 1: clean marine, 2: dust, 3: polluted continental, 4: clean continental, 5: polluted dust, 6: smoke, and 7: other). The Level 2 333 m and 5 km cloud layer data products were used to screen the clouds, as mentioned in Section 4.1.Aura OMI (Ozone monitoring Instrument) AI (Aerosol Index) data were used to validate the presence of dust outbreak in source region. The AI indicates presence of ultraviolet- (UV-) absorbing aerosols in the atmosphere and is derived from a residual of the measured UV reflectance. The AI data is useful to detect dust aerosols [21, 25]. The daily averaged data with a horizontal resolution of 1° × 1° was used in this paper.


### 2.3. Sun Photometer

A sun photometer (CE-318) was installed at our lidar site in 2008 and has been in continuous operation since then. The instrument measures the direct solar radiance and sky radiance in 8 wavelength bands (340, 380, 440, 500, 670, 870, 940, and 1020 nm). The level 1.5 version data were used to provide spectral distribution of aerosol optical depth (AOD) after calibration and cloud screening. The uncertainties of optical depth were of the order of 0.02 (at 440–1020 nm) and 0.04 (at 340–380 nm) at around noon [40]. The level 1.5 AOD data were used to derive the particle size distribution (PSD) using the algorithm of King et al. [41]. The equation used to retrieve PSD can be written as
(4)
τ
λ
=
∫
r
min
⁡
r
max
⁡
Q
ext
r
,
λ
,
m
π
r
2
n
r
d
r
,
where τ
λ is the aerosol optical depth at wavelength λ; r is the radius of aerosol particles; r
max
⁡ and r
min
⁡ are, respectively, the maximum and minimum radius confining the integration range; m is the complex refractive exponent; Q
ext
(
r
,
λ
,
m
) is the extinction efficiency factor from Mie theory, which is the function of particle radius, wavelength, and complex refractive exponent; and n
(
r
) is the unknown particle size distribution. In this study, the integration range (radius) is 0.05–10 μm. The AODs at 6 wavelengths (380, 440, 500, 670, 870, and 1020 nm) are used in the calculation.


### 2.4. HYSPLIT Backward Trajectory Analysis and Meteorological Data

The Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model was used to back-calculate the trajectories of the air mass based on atmospheric wind field data from the National Center for Environmental Prediction (NCEP) analyses [42]. It was thus possible to estimate an approximate source of each Asian dust plume.A radiosonde was launched at 08:00 LST and 20:00 LST every day from a location approximately 24 km from the lidar site. The relatively humidity (RH) and temperature data were used in this study.


## 2.1. Polarization Lidar System

Wuhan (30.5°N, 114.4°E) is located in southeast China. The polarization lidar system that was used is installed at Wuhan University. The lidar system parameters are presented in Table1. The polarization lidar transmitter used an Nd:YAG laser to produce an emission of ~25–40 mJ per pulse at 532 nm with a repetition rate of 20 Hz and a pulse width of 6 ns. The polarization purity of the transmitting laser was approximately 1/10000. On the receiver side, the backscattered photons were collected using a 300 mm Cassegrain telescope that was followed by a field-stop iris and a collimating lens. Double cubic polarizing beam splitters was employed to split the beam into two independent polarization channels in parallel and perpendicular directions relative to the outgoing laser beam. The filters for the two channels have a band width of ~10.6 cm−1 (0.3 nm at 532 nm). The two closest rotational Raman lines (N2) are ±11.9 cm−1 away from the elastic Cabannes line. Furthermore the pure rotational Raman lines are ~3 orders of magnitude weaker in intensity than the elastic Cabannes line. Thus the cross-talk of pure rotation Raman line was negligible [27]. The cross-talk between the two polarization channels is considerably suppressed by using double cubic polarizing beam splitters at each channel. The beams of the two channels were then detected using two photomultiplier tubes (PMT). The output pulses from the PMT were gathered using a signal acquisition instrument (Licel).Table 1
Specification of the depolarization lidar system at Wuhan University.
Transmitter
Receiver
Laser
Continuum Inlite II-20
Telescope
300 mm Cassegrain
Wavelength
532 nm
Diameter
300 mm
Energy/pulse
25–40 mJ
Field of view
1 mrad
Repetition rate
20 Hz
PMT
Hamamatsu 5783P
Beam divergence
6 ns
Acquisition instrument
Licel TR40-160The spatial and temporal resolutions of the lidar system were 30 m and 1 minute, respectively. The lower lidar observation limit was 0.3 km, which was determined based on the overlap of the laser and the field of view of the telescope. Our calibration for the polarization lidar was made by observing a fully depolarized light from thick-cloud-covered sky [28]. The resulting gain ratio (k) between the two channels was 0.065. The aerosol extinction was calculated using the combined parallel and perpendicular channels [29]. The volume depolarization ratio δ was calculated using the equations
(1)
δ
=
k
I
⊥
I
∥
,
where I
∥ and I
⊥ are the signals from the parallel and perpendicular channels, respectively. δ is the volume depolarization ratio caused by molecules and aerosols.The particle depolarization ratioδ
p was derived from the total-to-Rayleigh backscattering ratio BR and the volume depolarization ratio δ from the following equation [30]:(2)
δ
p
z
=
δ
z
BR
z
+
BR
z
δ
m
-
δ
m
-
δ
m
BR
z
-
1
+
BR
z
δ
m
-
δ
z
,
where δ
m is the molecular depolarization ratio for the lidar system. The value of molecular depolarization ratio measured by lidar depends on the wavelength and spectrum width of the detected light. The band widths of two filters used in our lidar system are both 0.3 nm. On the basis of calculation by Behrendt and Nakamura [27], the molecular depolarization ratio computed for Cabannes line is slightly larger than the theoretical value of 0.00363 for linear polarized incident light and 180-degree backscattering (see Figure 4 of [27]). Noting that the cross-talk of pure rotational Raman is negligible and the molecular depolarization ratio does not change with the variation of atmospheric temperature, the value of 0.4% was used here. The value of BR is defined as BR
z
=
1
+
β
a
(
z
)
/
β
m
(
z
) [31], where β
m and β
a are the molecular and aerosol backscattering coefficients, respectively.The aerosol extinction was calculated from the Mie backscattering signal using the method of Fernald [32]. Simultaneous measurement with lidar and sun photometer indicates that the average lidar ratio in the troposphere is ~50 sr (via comparing the aerosol optical depths derived from the lidar (at 532 nm) and the sun photometer (at 500 nm) [33, 34]). In practice, Murayama et al. [35] derived the extinction coefficient of aerosols by using the algorithm of Fernald [32], the lidar ratio of 50 sr was chosen to reproduce the optical depth derived from the sun photometer. The lidar ratio of 50 sr was also used in Asian dust observation in Taiwan [36] and Hefei, China [26]. Noh et al. [37] measured the lidar ratio of 51
±
6 sr for Asian dust with a Raman lidar at Gwangju, Korea, in spring. Therefore, the assumed lidar ratio of 50 sr is reasonable in this study. The molecular extinction was calculated using the U.S. Standard Atmosphere Model temperature and pressure data (U.S. GPO, Washington, D.C., 1976). The aerosol optical depth (AOD) and the centroid height (z
c) of a dust layer are given by
(3)
AOD
=
∫
Z
b
Z
t
α
a
(
z
)
d
z
,
z
c
=
∫
Z
b
Z
t
z
·
α
a
z
·
d
z
∫
Z
b
Z
t
α
a
z
·
d
z
,
where z
t and z
b are the dust layer top and base heights, respectively, and α
a is the aerosol extinction coefficient.


## 2.2. Satellite Data

The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument is carried on the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite and used to acquire vertical profiles of elastic backscatter at 532 nm and 1064 nm near nadir during both daytime and nighttime. The nearest CALIOP orbits pass over Wuhan at approximately 02:00 LST and 14:00 LST every day. Profiles of linear depolarization at 532 nm are provided. The ratio of 1064 nm backscattering to 532 nm backscattering, which is defined as the “color ratio,” is also provided. This ratio represents the aerosol particle size. The depolarization ratio of dust is large because of its nonspherical shape, which is a good indicator for separating dust from other aerosols [38]. Liu et al. [23] derived the probability density functions of volume depolarization ratio (VDR) for four types of aerosol (desert dust, biomass burning, continental, and maritime). The maritime and continental aerosols had a VDR distribution peaking at zero. The biomass burning aerosols had a nonzero peak at ~0.02. The VDR distribution of the dust aerosols was centered at ~0.17. Based on these probability density functions, when a VDR threshold is set to 0.06, only ~1.4% of dust aerosol particles were left out and only ~2.3% of non-dust aerosol particles were mistakenly classified as dust aerosol. Hence, in practical lidar data analysis (CALIPSO), the volume depolarization ratio value of 0.04 or 0.06 was employed as a threshold value to distinguish between dust and nondust aerosols by Liu et al. [39] and Huang et al. [21], respectively.In this study, the CALIPSO level 1B data product was used to validate the performance of the polarization lidar and study the horizontal range of the dust plumes. The level 2 5 km aerosol profile data was used to exhibit the vertical feature mask and aerosol subtypes mask. In this product, bits 1–3 of theAtmospheric_Volume_Description represent the feature mask (1: clean air, 2: cloud, and 3: aerosol) and bits of 10–12 of theAtmospheric_Volume_Description represent the aerosol subtypes (0: not determined, 1: clean marine, 2: dust, 3: polluted continental, 4: clean continental, 5: polluted dust, 6: smoke, and 7: other). The Level 2 333 m and 5 km cloud layer data products were used to screen the clouds, as mentioned in Section 4.1.Aura OMI (Ozone monitoring Instrument) AI (Aerosol Index) data were used to validate the presence of dust outbreak in source region. The AI indicates presence of ultraviolet- (UV-) absorbing aerosols in the atmosphere and is derived from a residual of the measured UV reflectance. The AI data is useful to detect dust aerosols [21, 25]. The daily averaged data with a horizontal resolution of 1° × 1° was used in this paper.


## 2.3. Sun Photometer

A sun photometer (CE-318) was installed at our lidar site in 2008 and has been in continuous operation since then. The instrument measures the direct solar radiance and sky radiance in 8 wavelength bands (340, 380, 440, 500, 670, 870, 940, and 1020 nm). The level 1.5 version data were used to provide spectral distribution of aerosol optical depth (AOD) after calibration and cloud screening. The uncertainties of optical depth were of the order of 0.02 (at 440–1020 nm) and 0.04 (at 340–380 nm) at around noon [40]. The level 1.5 AOD data were used to derive the particle size distribution (PSD) using the algorithm of King et al. [41]. The equation used to retrieve PSD can be written as
(4)
τ
λ
=
∫
r
min
⁡
r
max
⁡
Q
ext
r
,
λ
,
m
π
r
2
n
r
d
r
,
where τ
λ is the aerosol optical depth at wavelength λ; r is the radius of aerosol particles; r
max
⁡ and r
min
⁡ are, respectively, the maximum and minimum radius confining the integration range; m is the complex refractive exponent; Q
ext
(
r
,
λ
,
m
) is the extinction efficiency factor from Mie theory, which is the function of particle radius, wavelength, and complex refractive exponent; and n
(
r
) is the unknown particle size distribution. In this study, the integration range (radius) is 0.05–10 μm. The AODs at 6 wavelengths (380, 440, 500, 670, 870, and 1020 nm) are used in the calculation.


## 2.4. HYSPLIT Backward Trajectory Analysis and Meteorological Data

The Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model was used to back-calculate the trajectories of the air mass based on atmospheric wind field data from the National Center for Environmental Prediction (NCEP) analyses [42]. It was thus possible to estimate an approximate source of each Asian dust plume.A radiosonde was launched at 08:00 LST and 20:00 LST every day from a location approximately 24 km from the lidar site. The relatively humidity (RH) and temperature data were used in this study.


## 3. Observation Results

To validate the performance of our lidar system, a direct profile comparison of the lidar-derived optical parameters was first performed using the ground-based polarization lidar and almost simultaneous and nearly colocated CALIPSO lidar (Figure1). The 20 profiles closest to our lidar site obtained from the CALIPSO data were averaged to compare with the 1 min accumulated polarization lidar data at 02:24 LST on 23 April 2011. The horizontal distances between the satellite footprint and our lidar site are less than 50 km. Individual profiles from CALIPSO have an integration time of ~1 s (equivalent to a horizontal distance of ~6.7 km). They are then smoothed with a window length of 420 m (14 range bins at altitudes below 8.2 km and 7 range bins at altitudes above 8.2 km). The 1 min averaged profile from our ground-based lidar represents an accumulation of photon counts from 1000 laser shots. It is then smoothed with a window length of 210 m (7 range bins). The vertical volume depolarization ratio distributions measured using the two lidar systems quantitatively agree. The attenuated backscatter profiles exhibit a similar double-layer vertical structure of the dust layers. This comparison confirms the reliability of the vertical optical properties measured using our polarization lidar.Comparison of the 532 nm (a) volume depolarization ratio profiles and (b) attenuated backscatter coefficient obtained by the depolarization lidar (blue curve) and simultaneous and nearly colocated CALIPSO measurements (red curve). The polarization lidar profile represents a 1 min integration from 02:24 to 02:25 LST on 23 April 2011, while the CALIPSO profile stands for a ~1 s integration where the horizontal distances between the satellite footprint and our lidar site were less than 50 km.
(a)
(b)In this study, 13 dust events (41 observation days) observed with the ground-based polarization lidar at Wuhan were studied to reveal the vertical distribution and optical properties of dust transported over long distances from October 2010 to September 2012. Most of the dust events occurred in the spring (5 cases) and winter (6 cases). The remaining two events occurred in the summer and fall. Huang et al. [21] defined a dust plume as having a physical aerosol layer thickness exceeding 0.9 km and a volume depolarization ratio exceeding 0.06 throughout the layer. This definition was used in this study to identify dust layers. The CALIPSO level 1B data during the same two years were used to estimate the horizontal range of the dust plumes. The aerosol feature mask from CALIPSO level 2 5 km aerosol profile data was shown in case studies to validate the presence of dust aerosol. Out of 13 dust events registered by ground-based lidar, 10 events were simultaneously observed with CALIPSO lidar. The remaining 3 events were missed by CALIPSO observation because of the presence of cloud. First, two case studies are discussed, representing the typical dust features in the spring and winter. Thereafter, the overall statistical characteristics of the dust plumes are discussed.


### 3.1. Spring Dust Case: 22–29 April 2011

The longest dust event over Wuhan was observed on 22–29 April 2011 (Figure2). This figure presents the temporal variation of the extinction coefficient (a) and volume depolarization ratio (b) at 532 nm. A persistent, thick dust plume of this nature has never been observed over Wuhan before. Dust layers with large δ were clearly observed at heights of approximately 0–4 km and 5.5–8.0 km beginning on the afternoon of 22 April. The upper layer steadily rose until 23 April before descending and dissipating on 25 April. The lower layer filled the lower troposphere (ranging from the surface to approximately 6 km) during the entire event.Temporal variation in the extinction coefficient (a) and volume depolarization ratio (b) at 532 nm. These data were obtained using the depolarization lidar on 22–29 April 2011 over Wuhan.
(a)
(b)To examine the source of the dust plume, a three-day back-calculated air mass trajectory simulation from Wuhan that started at 1000UTC on 22, 25, and 28 April was performed using the HYSPLIT model, respectively. The air masses of five heights were simulated for each of three days. The Aura OMI AI distribution provided likely source region for this dust event (shown in Figures3(a)–3(c)). As seen in Figure 3, some calculated backward trajectories appeared to link to the dust-outbreak source regions in a reasonable time sequence, while on some calculated backward trajectories the dust source (shown by Figures 3(a)–3(c)) is invisible. The air masses on 22 April at the different heights originated in three regions. The air mass near the surface was tracked backward to the surface of north China. The lower dust layer below 4 km was from the Gobi desert 1-2 days before its appearance over Wuhan. This finding is similar to the result of Kwon et al. [5], who reported that the dust from the Gobi desert always affects the lower troposphere (<4 km) in downstream regions. The upper layer (approximately 6 km) originated from the Tibet Plateau (with an altitude of 5-6 km). However, the dust outbreak in corresponding source region is invisible because of the miss of OMI AI data. The air masses on 25 April showed that the dust layer at about 3-4 km was tracked backward to the Gobi near the border of China and Mongolia. These results suggest that the dust particles for the event are likely from Gobi of Mongolia. Additionally, another dust layer at about 3–5 km on 28 April was tracked backward to the Taklimakan region.Aura OMI Aerosol Index (AI) distribution over East China on (a) 20 April 2011, (b) 23 April 2011, and (c) 26 April 2011. Three-day back-calculated trajectories starting at Wuhan (30.5°N, 114.4°E) from 1000UTC on (d) 22 April 2011 (at 0.5 km, 2 km, 3 km, 6 km, and 7 km), (e) 25 April 2011 (at 1 km, 2 km, 3 km, 4 km, and 6 km), and (f) 28 April 2011 (at 1 km, 2 km, 3 km, 4 km, and 5 km), using the HYSPLIT model. The time-height display for the calculation of (d), (e), and (f) is shown in (g), (h), and (i), respectively.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)As shown in Table2, the mean top height, base height, thickness, and centroid height of the dust layers were 5.02 km, 1.52 km, 3.50 km, and 3.04 km, respectively. The standard deviations, as well as the maximum and the minimum value of the optical properties, were also shown in the table. The mean δ was 0.12. In addition, extremely large δ of 0.2-0.3 were observed below 5 km from midday on 26 April to the morning of 27 April. The δ
p (0.24) value was much larger than the δ value, which indicated the dominance of dust particles within the layer. The mean dust layer extinction was 0.11 km−1. As the mean extinction exceeded 0.1 km−1 in the observations, the atmosphere was very turbid. The mean dust layer AOD was 0.38.Table 2
Statistical characteristics of the dust events observed from October 2010 to September 2012 at Wuhan measured using the depolarization lidar and CALIPSO. The standard deviation of the mean values is shown. The maximum values and the minimum values of the optical properties are also shown in the parenthesis. Here, (H) and (L) represent the upper and lower layers, respectively.
Date
Top (km)
Base (km)
Thickness (km)
Centroid (km)
δ
δ
p
Extinction (km−1)
AOD
Horizontal range
Vertical width
8–10 November 2010
2.62 ± 0.99
0.73 ± 1.03
1.89 ± 0.63
1.52 ± 0.98
0.09 ± 0.02(0.06–0.32)
0.18 ± 0.09(0.09–0.23)
0.09 ± 0.04(0.01–0.21)
0.18 ± 0.09(0.01–0.43)
2000 km
Surface-3 km
2–5 December 2010
2.43 ± 1.29
0.61 ± 1.33
1.82 ± 0.48
1.04 ± 0.77
0.11 ± 0.04(0.06–0.41)
0.17 ± 0.05(0.07–0.42)
0.19 ± 0.21(0.01–3.03)
0.34 ± 0.27(0.01–3.09)
2400 km
Surface-2 km
7–9 December 2010
2.22 ± 0.24
⩽0.3
1.92 ± 0.23
1.02 ± 0.11
0.15 ± 0.02(0.11–0.20)
0.27 ± 0.04(0.18–0.33)
0.13 ± 0.03(0.07–0.28)
0.26 ± 0.06(0.13–0.61)
1300 km
Surface-2 km
26–28 December 2010
2.08 ± 0.82
0.64 ± 0.72
1.44 ± 0.38
1.16 ± 0.73
0.08 ± 0.03(0.06–0.37)
0.19 ± 0.07(0.08–0.41)
0.09 ± 0.10(0.02–0.84)
0.13 ± 0.17(0.02–0.37)
750 km
Surface-5 km
20-21 February 2011
1.74 ± 0.73
0.61 ± 0.75
1.12 ± 0.14
1.04 ± 0.84
0.10 ± 0.01(0.08–0.13)
0.14 ± 0.05(0.08–0.53)
0.40 ± 0.17(0.05–0.78)
0.46 ± 0.19(0.05–0.83)
1900 km (L) 900 km (H)
Surface-2 km (L) 2–8 km (H)
12-13 March 2011
3.44 ± 0.60
2.03 ± 0.20
1.41 ± 0.52
2.63 ± 0.23
0.13 ± 0.01(0.07–0.17)
0.28 ± 0.02(0.25–0.31)
0.06 ± 0.01(0.04–0.09)
0.09 ± 0.02(0.06–0.14)
—
—
22–29 April 2011
5.02 ± 1.99
1.52 ± 2.37
3.50 ± 1.72
3.04 ± 2.07
0.12 ± 0.03(0.07–0.21)
0.24 ± 0.09(0.09–0.40)
0.11 ± 0.08(0.01–0.58)
0.38 ± 0.30(0.01–1.84)
950 km (L) 750 km (H)
1–4 km (L) 3–9 km (H)
17–20 May 2011
4.74 ± 1.59
0.76 ± 0.93
3.98 ± 1.83
2.44 ± 0.93
0.14 ± 0.04(0.06–0.24)
0.27 ± 0.11(0.06–0.60)
0.14 ± 0.16(0.02–0.89)
0.48 ± 0.30(0.02–0.99)
1100 km
1–7 km
08-09 July 2011
2.73 ± 0.87
0.73 ± 1.06
2.00 ± 0.41
1.38 ± 1.02
0.12 ± 0.03(0.07–0.56)
0.23 ± 0.05(0.15–0.26)
0.15 ± 0.07(0.02–0.50)
0.30 ± 0.13(0.01–0.67)
—
—
23–25 December 2011
2.12 ± 1.38
0.81 ± 1.29
1.31 ± 0.31
1.41 ± 1.52
0.10 ± 0.01(0.08–0.21)
0.17 ± 0.11(0.09–0.40)
0.26 ± 0.16(0.02–0.80)
0.35 ± 0.23(0.02–1.07)
1600 km (L) 800 km (H)
Surface-2 km (L) 3–7 km (H)
6 January 2012
1.46 ± 0.09
⩽0.3
1.16 ± 0.09
0.75 ± 0.12
0.09 ± 0.01(0.08–0.11)
0.13 ± 0.10(0.09–0.16)
0.38 ± 0.14(0.13–0.70)
0.47 ± 0.19(0.14–0.76)
—
—
31 March-1 April 2012
3.65 ± 1.98
1.51 ± 1.28
2.14 ± 0.39
2.36 ± 1.61
0.09 ± 0.01(0.07–0.12)
0.22 ± 0.11(0.12–0.25)
0.08 ± 0.06(0.01–0.30)
0.15 ± 0.09(0.01–0.43)
>1360 km (L) 2300 km (H)
Surface-3 km (L) 3–10 km (H)
25–28 April 2012
4.00 ± 2.34
2.13 ± 2.24
1.87 ± 0.62
2.95 ± 2.40
0.14 ± 0.05(0.07–0.33)
0.17 ± 0.13(0.17–0.44)
0.11 ± 0.09(0.01–0.29)
0.20 ± 0.17(0.01–0.64)
1800 km
Surface-6 km
Average
3.51
1.45
2.06
2.34
0.12
0.22
0.19
0.32
—
—The vertical profiles of extinction, volume depolarization ratio, particle depolarization ratio, backscattering ratio, relative humidity, and temperature on 25 April are shown in Figure4. The error bars represent the standard deviations. A double-layer structure with variable extinctions and δ was observed above 2 km. The extinction and δ exhibited similar height distributions, indicating the dominance of nonspherical particles in the dust layer. The lower layer (2.0–4.5 km) δ ranged from 0.06 to 0.2, and the extinction ranged from 0.01 km−1 to 0.08 km−1. The respective maximum δ  (0.2) and extinction (0.08 km−1) both occurred at a height of approximately 3 km. The upper layer (5–8 km) exhibited a smaller extinction and δ. The extinction was extremely large below 2 km, which was likely due to the combined contributions of dust and local aerosols. The values of δ
p were very large, ranging from 0.2 to 0.4 within the lower dust layer. The BRs within the dust layer ranged from 1.1 to 2.2. The lower layer exhibited a slightly larger BR with a peak of 2.2 around 3.2 km. The dust layer was extremely dry; that is, the RH was less than 10% at a height of 2–8 km.Profiles of (a) aerosol extinction, (b) volume depolarization ratio, (c) particle depolarization ratio, (d) backscattering ratio measured using the depolarization lidar, (e) relative humidity, and (f) temperature measured using radiosonde during the half hour of 20:00–20:30 LST on 25 April 2011. The error bars represent the standard deviations.
(a)
(b)
(c)
(d)
(e)
(f)Figure5 shows plots of the attenuated backscatter and volume depolarization ratio at 532 nm, the color ratio, the vertical feature mask, the aerosol subtype mask, and the trajectory of CALIPSO observed over southeast China on 22 April. The dust plumes are denoted by the black rectangle. The dust plumes are represented by large δ (0.2-0.3) and color ratios (0.7–1.0) due to the nonspherical shape and large particle size of the dust. The double-layer structure observed using the depolarization lidar was also seen in the CALIPSO data. The aerosol subtype mask indicates that the lower layer is the mixture of dust and polluted dust and the upper layer contains only pure dust. Moreover, the CALIPSO data suggest a horizontal range of ~950 km and ~750 km of the layers at heights of approximately 1–4 km and 4–9 km, respectively. Specially, the height of the upper layer decreases from north to south.CALIPSO altitude-orbit cross-section measurements of (a) 532 nm total attenuated backscatter coefficient, (b) volume depolarization ratio, and (c) 1064 nm/532 nm backscatter color ratio; the figure shows the (d) vertical feature mask, (e) aerosol subtype mask, and (f) CALIPSO trajectory in the map (the corresponding orbit is 2011-04-22T18-16-45) over southeast China between 20° and 50°N on 22 April 2011. The red line (see (a)) marks the location of Wuhan. The locations of Taklimakan, Gobi, and Wuhan are marked in (f).
(a)
Attenuated backscatter coefficient (km−1 sr−1)
(b)
Volume depolarization ratio
(c)
Color ratio
(d)
Vertical feature mask
(e)
Aerosol subtypes
(f)
Trajectory of CALIPSO


### 3.2. Winter Dust Case: 7–9 December 2010

Figure6 shows the temporal variation in the extinction (a) and volume depolarization ratio (b) at 532 nm on 7–9 December 2010, which represents a typical wintertime dust plume. Dust was observed for at least two days; however, the beginning and end of the dust event were not observed. Large δ (ranging from 0.1 to 0.25) and extinctions exceeding 0.1 km−1 were observed below 3 km during the event, which are typical optical property values in the winter.Same as Figure2, except for the date, 7–9 December 2010.
(a)
(b)A three-day back-calculated air mass trajectory simulation from Wuhan starting at 1000UTC on 7 December and 9 December was performed using the HYSPLIT model. As shown in Figure7, three starting heights were simulated with the model. The results indicate that the dust plume observed below 3 km over Wuhan was due to the long-range transport of Gobi dust 2 days before the dust event began. The large values of OMI AI (up to 2.0) indicated the occurrence of dust outbreak in Gobi region.Aura OMI Aerosol Index (AI) distribution over East Asia on (a) 6 December 2010 and (b) 7 December 2010. Three-day back-calculated trajectories starting at Wuhan (30.5°N, 114.4°E) from 1000UTC on (c) 07 December 2010 (at 0.5 km, 1 km, 2 km, 3 km, and 5 km) and (d) 09 December 2010 (at 0.5 km, 1 km, 2 km, 3 km, and 4.5 km), using the HYSPLIT model. The time-height display for the calculation of (c) and (d) is shown in (e) and (f), respectively.
(a)
(b)
(c)
(d)
(e)
(f)As shown in Table2, the mean top height, base height, thickness, and centroid height of the dust layers are 2.22 km, <0.30 km, 1.92 km, and 1.02 km, respectively. The mean δ is 0.15. Specifically, extremely large δ that exceed 0.2 were observed inside the dust layer from the morning of 8 December to midday on 9 December. The mean δ
p is 0.27. The mean extinction (0.13 km−1) and the mean AOD (0.26) indicate that dust layers in winter are always very turbid.Figure8 shows extinction, volume depolarization ratio, particle depolarization ratio, backscattering ratio, relative humidity, and temperature vertical profiles on 7 December from the polarization lidar and radiosonde. A pronounced single dust layer was observed below 2 km with extinctions ranging from 0.03 km−1 to 0.20 km−1 and peaking at 0.6 km. The δ ranged from 0.07 to 0.16 within the dust layer and rapidly decreased above 1.5 km to nearly 0 above 2.5 km. The δ
p within the dust layer ranged from 0.20 to 0.33. The peak value of the δ
p appeared at 1.7 km, which was different from the height of 0.6 km for aerosol extinction; thus, anthropogenic aerosols within the PBL contributed greatly below 1 km. The RH decreased with height from 40% at the surface to only 2% above 2.5 km. RHs of 20–50% were typical in wintertime dust layers.Same as Figure4, except for the dates, 7 December 2010.
(a)
(b)
(c)
(d)
(e)
(f)Figure9 shows plots of the attenuated backscatter and volume depolarization ratio at 532 nm, the color ratio, the vertical feature mask, the aerosol subtype mask, and the trajectory of CALIPSO observed over southeast China on 7 December. Dust plumes with large δ (0.2-0.3) and color ratios (0.7–1.0) are denoted with a black rectangle. The single dust layer observed from the surface to 2 km displays the same characteristics as those seen in the ground-based lidar observations. Dust plumes below a height of 3 km frequently occur between 30°N and 40°N over mainland China in the winter. The dust layers typically have a horizontal range of ~1300 km. The CALIPSO data confirm that the high δ observed at this low height using ground-based lidar was not caused by local floating dust but rather by a long-range transport dust plume.Same as Figure5, except for the date, 8 December 2010. The corresponding CALIPSO orbit is 2010-12-08T18-00-13.
(a)
Attenuated backscatter coefficient (km−1 sr−1)
(b)
Volume depolarization ratio
(c)
Color ratio
(d)
Vertical feature sask
(e)
Aerosol subtypes
(f)
Trajectory of CALIPSO


### 3.3. Statistical Characteristics of Dust over Wuhan

Table2 lists the characteristics of the 13 dust events extracted from the depolarization lidar measurements collected between October 2010 and September 2012. All of the selected events spanned at least 10 hours and exhibited no complex mixing between dust and clouds. Table 2 lists the mean top height, base height, thickness, centroid height, volume depolarization ratio, particle depolarization ratio, extinction, and AOD of these dust events. The standard deviations, maximum and minimum values, are also shown in the table. The average dust layer height was 1.45–3.51 km; the maximum and minimum heights were approximately 5.02 km and less than 0.3 km, respectively. The average centroid height was 2.34 km, which was the altitude of highest dust loading. The δ ranged from 0.09 to 0.15; the mean was 0.12. The δ
p ranged from 0.13 to 0.33; the mean was 0.22. The extinctions ranged from 0.06 km−1 to 0.40 km−1; the mean was 0.19 km−1. Moreover, the AODs ranged from 0.09 to 0.48; the mean was 0.32. The vertical and horizontal ranges of the dust plumes obtained from the CALIPSO data are also shown in Table 2. The retrieval failures in 4 cases were due to missing CALIPSO data or the failure to identify a dust plume due to the presence of clouds. Among the remaining 9 cases, 4 cases exhibit a double-layer structure. The horizontal ranges of the dust plumes were approximately 750 km to 2400 km. The dust plume thicknesses ranged from approximately 2 km to 7 km.Dust emissions over Wuhan are frequent during the spring and winter. Comparisons of the vertical distribution and optical properties of dust between the two seasons are listed in Table3. The dust layer was higher in the spring (approximately 1.59–4.17 km) than winter (approximately 0.54–2.01 km). The mean centroid height of 2.68 km in the spring was higher than in winter (1.07 km). Moreover, the dust layer was thicker in the spring (2.58 km) than in winter (1.46 km). The mean δ was slightly larger in the spring (0.12) than in winter (0.10). The mean δ
p was much larger in spring (0.27) than in winter (0.18), indicating the greater nonsphericity of dust particles in the spring. In contrast, the mean extinction in the winter (0.24 km−1) was 1.4 times larger than in spring (0.10 km−1). Furthermore, the mean AOD in the winter (0.34) was 31% larger than in spring (0.26). Based on the CALIPSO data, the dust plumes were thicker and higher in the spring, and the plumes were generally thinner and closer to the surface in the winter.Table 3
Comparison of dust properties over Wuhan during spring and winter. The standard deviations of the mean values are shown. The maximum values and the minimum values of the optical properties are also shown in the parenthesis.
Time
Top (km)
Base (km)
Thickness (km)
Centroid (km)
δ
δ
p
α
a (km−1)
AOD
Spring
4.17 ± 0.68
1.59 ± 0.54
2.58 ± 0.97
2.68 ± 0.66
0.12 ± 0.02(0.09–0.14)
0.27 ± 0.04(0.22–0.33)
0.10 ± 0.03(0.06–0.14)
0.26 ± 0.13(0.09–0.48)
Winter
2.01 ± 0.35
0.54 ± 0.20
1.46 ± 0.31
1.07 ± 0.62
0.10 ± 0.02(0.08–0.15)
0.18 ± 0.05(0.13–0.27)
0.24 ± 0.14(0.09–0.40)
0.34 ± 0.15(0.13–0.47)To obtain more information on the microphysical character of dust aerosols, the particle size distribution (PSD) of the aerosols was computed using sun photometer (CE-318) measurements [41]. Figure 10 shows the retrieved aerosol size distributions for the dust cases in spring (blue) and winter (red). The PSD was bimodal with the peak of the coarse model at 2.0 μm in both seasons, indicating the typical particle radius of dust aerosols. The peak of the fine model was approximately 0.13 μm in winter and 0.07 μm in spring.Figure 10
Particle size distribution of aerosols over Wuhan during the dust cases in spring and winter, derived from sun photometer data (CE-318).Similar bimodal distributions were also observed in spring and winter in Wuhan by Gong et al. [43]. The fine and coarse modes in both spring and winter peaked at 0.2 μm and 4 μm, respectively. The particle size distribution (PSD) value was larger in coarse mode in winter than in spring. Liu et al. [44] observed two spring dust events at Taihu, China, in 2009. The fine and coarse modes in the PSD observed by sun photometer peaked at 0.1 μm and 2 μm, respectively. Liu et al. [45] presented the seasonal mean PSD based on the sun photometer observations from 2006 to 2009 at Taihu. They noticed a similar bimodal distribution during winter and spring with nearly the same peaks in both fine and coarse modes. The fine and coarse modes were observed at 0.1-0.2 μm and 3.0-4.0 μm, respectively. These observational results (bimodal distribution with a fine mode at 0.1-0.2 μm and a coarse mode at 2.0–4.0 μm) are similar to that given in the current work. A main reason for the similarity is that Wuhan has weather condition and geographical latitude similar to the sites where the earlier observations were made. However, the observations near the dust source region of Asian dust showed rather different results. Zheng et al. [46] observed the PSD at Dunhuang, China, from January 1999 to March 2001. The mean PSD in spring and winter had a fine mode with a central radius of ~0.25 μm, and a coarse mode with a central radius of ~7.7 μm. The central radius of coarse mode was remarkably larger than the downstream regions of observations. It is seemed that extremely large dust particles were removed during the dust plume transport.


## 3.1. Spring Dust Case: 22–29 April 2011

The longest dust event over Wuhan was observed on 22–29 April 2011 (Figure2). This figure presents the temporal variation of the extinction coefficient (a) and volume depolarization ratio (b) at 532 nm. A persistent, thick dust plume of this nature has never been observed over Wuhan before. Dust layers with large δ were clearly observed at heights of approximately 0–4 km and 5.5–8.0 km beginning on the afternoon of 22 April. The upper layer steadily rose until 23 April before descending and dissipating on 25 April. The lower layer filled the lower troposphere (ranging from the surface to approximately 6 km) during the entire event.Temporal variation in the extinction coefficient (a) and volume depolarization ratio (b) at 532 nm. These data were obtained using the depolarization lidar on 22–29 April 2011 over Wuhan.
(a)
(b)To examine the source of the dust plume, a three-day back-calculated air mass trajectory simulation from Wuhan that started at 1000UTC on 22, 25, and 28 April was performed using the HYSPLIT model, respectively. The air masses of five heights were simulated for each of three days. The Aura OMI AI distribution provided likely source region for this dust event (shown in Figures3(a)–3(c)). As seen in Figure 3, some calculated backward trajectories appeared to link to the dust-outbreak source regions in a reasonable time sequence, while on some calculated backward trajectories the dust source (shown by Figures 3(a)–3(c)) is invisible. The air masses on 22 April at the different heights originated in three regions. The air mass near the surface was tracked backward to the surface of north China. The lower dust layer below 4 km was from the Gobi desert 1-2 days before its appearance over Wuhan. This finding is similar to the result of Kwon et al. [5], who reported that the dust from the Gobi desert always affects the lower troposphere (<4 km) in downstream regions. The upper layer (approximately 6 km) originated from the Tibet Plateau (with an altitude of 5-6 km). However, the dust outbreak in corresponding source region is invisible because of the miss of OMI AI data. The air masses on 25 April showed that the dust layer at about 3-4 km was tracked backward to the Gobi near the border of China and Mongolia. These results suggest that the dust particles for the event are likely from Gobi of Mongolia. Additionally, another dust layer at about 3–5 km on 28 April was tracked backward to the Taklimakan region.Aura OMI Aerosol Index (AI) distribution over East China on (a) 20 April 2011, (b) 23 April 2011, and (c) 26 April 2011. Three-day back-calculated trajectories starting at Wuhan (30.5°N, 114.4°E) from 1000UTC on (d) 22 April 2011 (at 0.5 km, 2 km, 3 km, 6 km, and 7 km), (e) 25 April 2011 (at 1 km, 2 km, 3 km, 4 km, and 6 km), and (f) 28 April 2011 (at 1 km, 2 km, 3 km, 4 km, and 5 km), using the HYSPLIT model. The time-height display for the calculation of (d), (e), and (f) is shown in (g), (h), and (i), respectively.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)As shown in Table2, the mean top height, base height, thickness, and centroid height of the dust layers were 5.02 km, 1.52 km, 3.50 km, and 3.04 km, respectively. The standard deviations, as well as the maximum and the minimum value of the optical properties, were also shown in the table. The mean δ was 0.12. In addition, extremely large δ of 0.2-0.3 were observed below 5 km from midday on 26 April to the morning of 27 April. The δ
p (0.24) value was much larger than the δ value, which indicated the dominance of dust particles within the layer. The mean dust layer extinction was 0.11 km−1. As the mean extinction exceeded 0.1 km−1 in the observations, the atmosphere was very turbid. The mean dust layer AOD was 0.38.Table 2
Statistical characteristics of the dust events observed from October 2010 to September 2012 at Wuhan measured using the depolarization lidar and CALIPSO. The standard deviation of the mean values is shown. The maximum values and the minimum values of the optical properties are also shown in the parenthesis. Here, (H) and (L) represent the upper and lower layers, respectively.
Date
Top (km)
Base (km)
Thickness (km)
Centroid (km)
δ
δ
p
Extinction (km−1)
AOD
Horizontal range
Vertical width
8–10 November 2010
2.62 ± 0.99
0.73 ± 1.03
1.89 ± 0.63
1.52 ± 0.98
0.09 ± 0.02(0.06–0.32)
0.18 ± 0.09(0.09–0.23)
0.09 ± 0.04(0.01–0.21)
0.18 ± 0.09(0.01–0.43)
2000 km
Surface-3 km
2–5 December 2010
2.43 ± 1.29
0.61 ± 1.33
1.82 ± 0.48
1.04 ± 0.77
0.11 ± 0.04(0.06–0.41)
0.17 ± 0.05(0.07–0.42)
0.19 ± 0.21(0.01–3.03)
0.34 ± 0.27(0.01–3.09)
2400 km
Surface-2 km
7–9 December 2010
2.22 ± 0.24
⩽0.3
1.92 ± 0.23
1.02 ± 0.11
0.15 ± 0.02(0.11–0.20)
0.27 ± 0.04(0.18–0.33)
0.13 ± 0.03(0.07–0.28)
0.26 ± 0.06(0.13–0.61)
1300 km
Surface-2 km
26–28 December 2010
2.08 ± 0.82
0.64 ± 0.72
1.44 ± 0.38
1.16 ± 0.73
0.08 ± 0.03(0.06–0.37)
0.19 ± 0.07(0.08–0.41)
0.09 ± 0.10(0.02–0.84)
0.13 ± 0.17(0.02–0.37)
750 km
Surface-5 km
20-21 February 2011
1.74 ± 0.73
0.61 ± 0.75
1.12 ± 0.14
1.04 ± 0.84
0.10 ± 0.01(0.08–0.13)
0.14 ± 0.05(0.08–0.53)
0.40 ± 0.17(0.05–0.78)
0.46 ± 0.19(0.05–0.83)
1900 km (L) 900 km (H)
Surface-2 km (L) 2–8 km (H)
12-13 March 2011
3.44 ± 0.60
2.03 ± 0.20
1.41 ± 0.52
2.63 ± 0.23
0.13 ± 0.01(0.07–0.17)
0.28 ± 0.02(0.25–0.31)
0.06 ± 0.01(0.04–0.09)
0.09 ± 0.02(0.06–0.14)
—
—
22–29 April 2011
5.02 ± 1.99
1.52 ± 2.37
3.50 ± 1.72
3.04 ± 2.07
0.12 ± 0.03(0.07–0.21)
0.24 ± 0.09(0.09–0.40)
0.11 ± 0.08(0.01–0.58)
0.38 ± 0.30(0.01–1.84)
950 km (L) 750 km (H)
1–4 km (L) 3–9 km (H)
17–20 May 2011
4.74 ± 1.59
0.76 ± 0.93
3.98 ± 1.83
2.44 ± 0.93
0.14 ± 0.04(0.06–0.24)
0.27 ± 0.11(0.06–0.60)
0.14 ± 0.16(0.02–0.89)
0.48 ± 0.30(0.02–0.99)
1100 km
1–7 km
08-09 July 2011
2.73 ± 0.87
0.73 ± 1.06
2.00 ± 0.41
1.38 ± 1.02
0.12 ± 0.03(0.07–0.56)
0.23 ± 0.05(0.15–0.26)
0.15 ± 0.07(0.02–0.50)
0.30 ± 0.13(0.01–0.67)
—
—
23–25 December 2011
2.12 ± 1.38
0.81 ± 1.29
1.31 ± 0.31
1.41 ± 1.52
0.10 ± 0.01(0.08–0.21)
0.17 ± 0.11(0.09–0.40)
0.26 ± 0.16(0.02–0.80)
0.35 ± 0.23(0.02–1.07)
1600 km (L) 800 km (H)
Surface-2 km (L) 3–7 km (H)
6 January 2012
1.46 ± 0.09
⩽0.3
1.16 ± 0.09
0.75 ± 0.12
0.09 ± 0.01(0.08–0.11)
0.13 ± 0.10(0.09–0.16)
0.38 ± 0.14(0.13–0.70)
0.47 ± 0.19(0.14–0.76)
—
—
31 March-1 April 2012
3.65 ± 1.98
1.51 ± 1.28
2.14 ± 0.39
2.36 ± 1.61
0.09 ± 0.01(0.07–0.12)
0.22 ± 0.11(0.12–0.25)
0.08 ± 0.06(0.01–0.30)
0.15 ± 0.09(0.01–0.43)
>1360 km (L) 2300 km (H)
Surface-3 km (L) 3–10 km (H)
25–28 April 2012
4.00 ± 2.34
2.13 ± 2.24
1.87 ± 0.62
2.95 ± 2.40
0.14 ± 0.05(0.07–0.33)
0.17 ± 0.13(0.17–0.44)
0.11 ± 0.09(0.01–0.29)
0.20 ± 0.17(0.01–0.64)
1800 km
Surface-6 km
Average
3.51
1.45
2.06
2.34
0.12
0.22
0.19
0.32
—
—The vertical profiles of extinction, volume depolarization ratio, particle depolarization ratio, backscattering ratio, relative humidity, and temperature on 25 April are shown in Figure4. The error bars represent the standard deviations. A double-layer structure with variable extinctions and δ was observed above 2 km. The extinction and δ exhibited similar height distributions, indicating the dominance of nonspherical particles in the dust layer. The lower layer (2.0–4.5 km) δ ranged from 0.06 to 0.2, and the extinction ranged from 0.01 km−1 to 0.08 km−1. The respective maximum δ  (0.2) and extinction (0.08 km−1) both occurred at a height of approximately 3 km. The upper layer (5–8 km) exhibited a smaller extinction and δ. The extinction was extremely large below 2 km, which was likely due to the combined contributions of dust and local aerosols. The values of δ
p were very large, ranging from 0.2 to 0.4 within the lower dust layer. The BRs within the dust layer ranged from 1.1 to 2.2. The lower layer exhibited a slightly larger BR with a peak of 2.2 around 3.2 km. The dust layer was extremely dry; that is, the RH was less than 10% at a height of 2–8 km.Profiles of (a) aerosol extinction, (b) volume depolarization ratio, (c) particle depolarization ratio, (d) backscattering ratio measured using the depolarization lidar, (e) relative humidity, and (f) temperature measured using radiosonde during the half hour of 20:00–20:30 LST on 25 April 2011. The error bars represent the standard deviations.
(a)
(b)
(c)
(d)
(e)
(f)Figure5 shows plots of the attenuated backscatter and volume depolarization ratio at 532 nm, the color ratio, the vertical feature mask, the aerosol subtype mask, and the trajectory of CALIPSO observed over southeast China on 22 April. The dust plumes are denoted by the black rectangle. The dust plumes are represented by large δ (0.2-0.3) and color ratios (0.7–1.0) due to the nonspherical shape and large particle size of the dust. The double-layer structure observed using the depolarization lidar was also seen in the CALIPSO data. The aerosol subtype mask indicates that the lower layer is the mixture of dust and polluted dust and the upper layer contains only pure dust. Moreover, the CALIPSO data suggest a horizontal range of ~950 km and ~750 km of the layers at heights of approximately 1–4 km and 4–9 km, respectively. Specially, the height of the upper layer decreases from north to south.CALIPSO altitude-orbit cross-section measurements of (a) 532 nm total attenuated backscatter coefficient, (b) volume depolarization ratio, and (c) 1064 nm/532 nm backscatter color ratio; the figure shows the (d) vertical feature mask, (e) aerosol subtype mask, and (f) CALIPSO trajectory in the map (the corresponding orbit is 2011-04-22T18-16-45) over southeast China between 20° and 50°N on 22 April 2011. The red line (see (a)) marks the location of Wuhan. The locations of Taklimakan, Gobi, and Wuhan are marked in (f).
(a)
Attenuated backscatter coefficient (km−1 sr−1)
(b)
Volume depolarization ratio
(c)
Color ratio
(d)
Vertical feature mask
(e)
Aerosol subtypes
(f)
Trajectory of CALIPSO


## 3.2. Winter Dust Case: 7–9 December 2010

Figure6 shows the temporal variation in the extinction (a) and volume depolarization ratio (b) at 532 nm on 7–9 December 2010, which represents a typical wintertime dust plume. Dust was observed for at least two days; however, the beginning and end of the dust event were not observed. Large δ (ranging from 0.1 to 0.25) and extinctions exceeding 0.1 km−1 were observed below 3 km during the event, which are typical optical property values in the winter.Same as Figure2, except for the date, 7–9 December 2010.
(a)
(b)A three-day back-calculated air mass trajectory simulation from Wuhan starting at 1000UTC on 7 December and 9 December was performed using the HYSPLIT model. As shown in Figure7, three starting heights were simulated with the model. The results indicate that the dust plume observed below 3 km over Wuhan was due to the long-range transport of Gobi dust 2 days before the dust event began. The large values of OMI AI (up to 2.0) indicated the occurrence of dust outbreak in Gobi region.Aura OMI Aerosol Index (AI) distribution over East Asia on (a) 6 December 2010 and (b) 7 December 2010. Three-day back-calculated trajectories starting at Wuhan (30.5°N, 114.4°E) from 1000UTC on (c) 07 December 2010 (at 0.5 km, 1 km, 2 km, 3 km, and 5 km) and (d) 09 December 2010 (at 0.5 km, 1 km, 2 km, 3 km, and 4.5 km), using the HYSPLIT model. The time-height display for the calculation of (c) and (d) is shown in (e) and (f), respectively.
(a)
(b)
(c)
(d)
(e)
(f)As shown in Table2, the mean top height, base height, thickness, and centroid height of the dust layers are 2.22 km, <0.30 km, 1.92 km, and 1.02 km, respectively. The mean δ is 0.15. Specifically, extremely large δ that exceed 0.2 were observed inside the dust layer from the morning of 8 December to midday on 9 December. The mean δ
p is 0.27. The mean extinction (0.13 km−1) and the mean AOD (0.26) indicate that dust layers in winter are always very turbid.Figure8 shows extinction, volume depolarization ratio, particle depolarization ratio, backscattering ratio, relative humidity, and temperature vertical profiles on 7 December from the polarization lidar and radiosonde. A pronounced single dust layer was observed below 2 km with extinctions ranging from 0.03 km−1 to 0.20 km−1 and peaking at 0.6 km. The δ ranged from 0.07 to 0.16 within the dust layer and rapidly decreased above 1.5 km to nearly 0 above 2.5 km. The δ
p within the dust layer ranged from 0.20 to 0.33. The peak value of the δ
p appeared at 1.7 km, which was different from the height of 0.6 km for aerosol extinction; thus, anthropogenic aerosols within the PBL contributed greatly below 1 km. The RH decreased with height from 40% at the surface to only 2% above 2.5 km. RHs of 20–50% were typical in wintertime dust layers.Same as Figure4, except for the dates, 7 December 2010.
(a)
(b)
(c)
(d)
(e)
(f)Figure9 shows plots of the attenuated backscatter and volume depolarization ratio at 532 nm, the color ratio, the vertical feature mask, the aerosol subtype mask, and the trajectory of CALIPSO observed over southeast China on 7 December. Dust plumes with large δ (0.2-0.3) and color ratios (0.7–1.0) are denoted with a black rectangle. The single dust layer observed from the surface to 2 km displays the same characteristics as those seen in the ground-based lidar observations. Dust plumes below a height of 3 km frequently occur between 30°N and 40°N over mainland China in the winter. The dust layers typically have a horizontal range of ~1300 km. The CALIPSO data confirm that the high δ observed at this low height using ground-based lidar was not caused by local floating dust but rather by a long-range transport dust plume.Same as Figure5, except for the date, 8 December 2010. The corresponding CALIPSO orbit is 2010-12-08T18-00-13.
(a)
Attenuated backscatter coefficient (km−1 sr−1)
(b)
Volume depolarization ratio
(c)
Color ratio
(d)
Vertical feature sask
(e)
Aerosol subtypes
(f)
Trajectory of CALIPSO


## 3.3. Statistical Characteristics of Dust over Wuhan

Table2 lists the characteristics of the 13 dust events extracted from the depolarization lidar measurements collected between October 2010 and September 2012. All of the selected events spanned at least 10 hours and exhibited no complex mixing between dust and clouds. Table 2 lists the mean top height, base height, thickness, centroid height, volume depolarization ratio, particle depolarization ratio, extinction, and AOD of these dust events. The standard deviations, maximum and minimum values, are also shown in the table. The average dust layer height was 1.45–3.51 km; the maximum and minimum heights were approximately 5.02 km and less than 0.3 km, respectively. The average centroid height was 2.34 km, which was the altitude of highest dust loading. The δ ranged from 0.09 to 0.15; the mean was 0.12. The δ
p ranged from 0.13 to 0.33; the mean was 0.22. The extinctions ranged from 0.06 km−1 to 0.40 km−1; the mean was 0.19 km−1. Moreover, the AODs ranged from 0.09 to 0.48; the mean was 0.32. The vertical and horizontal ranges of the dust plumes obtained from the CALIPSO data are also shown in Table 2. The retrieval failures in 4 cases were due to missing CALIPSO data or the failure to identify a dust plume due to the presence of clouds. Among the remaining 9 cases, 4 cases exhibit a double-layer structure. The horizontal ranges of the dust plumes were approximately 750 km to 2400 km. The dust plume thicknesses ranged from approximately 2 km to 7 km.Dust emissions over Wuhan are frequent during the spring and winter. Comparisons of the vertical distribution and optical properties of dust between the two seasons are listed in Table3. The dust layer was higher in the spring (approximately 1.59–4.17 km) than winter (approximately 0.54–2.01 km). The mean centroid height of 2.68 km in the spring was higher than in winter (1.07 km). Moreover, the dust layer was thicker in the spring (2.58 km) than in winter (1.46 km). The mean δ was slightly larger in the spring (0.12) than in winter (0.10). The mean δ
p was much larger in spring (0.27) than in winter (0.18), indicating the greater nonsphericity of dust particles in the spring. In contrast, the mean extinction in the winter (0.24 km−1) was 1.4 times larger than in spring (0.10 km−1). Furthermore, the mean AOD in the winter (0.34) was 31% larger than in spring (0.26). Based on the CALIPSO data, the dust plumes were thicker and higher in the spring, and the plumes were generally thinner and closer to the surface in the winter.Table 3
Comparison of dust properties over Wuhan during spring and winter. The standard deviations of the mean values are shown. The maximum values and the minimum values of the optical properties are also shown in the parenthesis.
Time
Top (km)
Base (km)
Thickness (km)
Centroid (km)
δ
δ
p
α
a (km−1)
AOD
Spring
4.17 ± 0.68
1.59 ± 0.54
2.58 ± 0.97
2.68 ± 0.66
0.12 ± 0.02(0.09–0.14)
0.27 ± 0.04(0.22–0.33)
0.10 ± 0.03(0.06–0.14)
0.26 ± 0.13(0.09–0.48)
Winter
2.01 ± 0.35
0.54 ± 0.20
1.46 ± 0.31
1.07 ± 0.62
0.10 ± 0.02(0.08–0.15)
0.18 ± 0.05(0.13–0.27)
0.24 ± 0.14(0.09–0.40)
0.34 ± 0.15(0.13–0.47)To obtain more information on the microphysical character of dust aerosols, the particle size distribution (PSD) of the aerosols was computed using sun photometer (CE-318) measurements [41]. Figure 10 shows the retrieved aerosol size distributions for the dust cases in spring (blue) and winter (red). The PSD was bimodal with the peak of the coarse model at 2.0 μm in both seasons, indicating the typical particle radius of dust aerosols. The peak of the fine model was approximately 0.13 μm in winter and 0.07 μm in spring.Figure 10
Particle size distribution of aerosols over Wuhan during the dust cases in spring and winter, derived from sun photometer data (CE-318).Similar bimodal distributions were also observed in spring and winter in Wuhan by Gong et al. [43]. The fine and coarse modes in both spring and winter peaked at 0.2 μm and 4 μm, respectively. The particle size distribution (PSD) value was larger in coarse mode in winter than in spring. Liu et al. [44] observed two spring dust events at Taihu, China, in 2009. The fine and coarse modes in the PSD observed by sun photometer peaked at 0.1 μm and 2 μm, respectively. Liu et al. [45] presented the seasonal mean PSD based on the sun photometer observations from 2006 to 2009 at Taihu. They noticed a similar bimodal distribution during winter and spring with nearly the same peaks in both fine and coarse modes. The fine and coarse modes were observed at 0.1-0.2 μm and 3.0-4.0 μm, respectively. These observational results (bimodal distribution with a fine mode at 0.1-0.2 μm and a coarse mode at 2.0–4.0 μm) are similar to that given in the current work. A main reason for the similarity is that Wuhan has weather condition and geographical latitude similar to the sites where the earlier observations were made. However, the observations near the dust source region of Asian dust showed rather different results. Zheng et al. [46] observed the PSD at Dunhuang, China, from January 1999 to March 2001. The mean PSD in spring and winter had a fine mode with a central radius of ~0.25 μm, and a coarse mode with a central radius of ~7.7 μm. The central radius of coarse mode was remarkably larger than the downstream regions of observations. It is seemed that extremely large dust particles were removed during the dust plume transport.


## 4. Discussion

### 4.1. Different Characteristics in Spring and Winter

As mentioned above, the dust layers are higher in spring than in winter. Two likely explanations are presented below. First, Asian dust is generated by cold air outbreaks and driven by westerly jets. Asian dust outbreaks are more intense in the spring than in winter due to the higher frequency of strong winds. Kurosaki and Mikami [47] found good agreement in the yearly variations and spatial distributions between the surface wind and dust outbreaks over East Asia during the period 2000–2002. The stronger surface winds in the spring effectively carry the desert dust particles to higher altitudes in the free troposphere. The elevated dust aerosols are then transported over long distances. As a result, the spring dust plume over downstream regions (such as Wuhan) is higher than the winter plume. Second, the temperature inversion that occurs frequently in the lower troposphere over Wuhan during the winter restricts the vertical convective transportation of dust and other aerosols in the lower troposphere. Zhang et al. [48] observed two lower tropospheric inversion layers over Yichang (111°18′E, 30°42′N) in January 2007. Their observations reveal the inversion layers are closely related to the intensive inertial gravity waves propagation in winter tropospheric jets. The temperature profiles from the radiosondes were statistically analyzed over Wuhan. Of the 27 radiosondes launched during the 6 winter cases, a temperature inversion (a layer thicker than 0.12 km with a lapse rate greater than −2 K/km) at heights of 1.5–3 km was observed in 60% of the radiosonde measurements. Of the 32 radiosondes launched during the 5 spring cases, a temperature inversion in the same region was observed in only 28% of the radiosonde measurements. The probability of a temperature inversion in the winter was approximately twice that of spring. Therefore, it is likely that the dust particles are inhibited from leaving the lower troposphere via convective transport due to the restriction caused by the temperature inversion in the winter.As shown in Table3, the mean dust layer extinction and AOD in the spring are 0.10 km−1 and 0.26, respectively. These values increase by 140% and 31% in the winter, respectively. What caused the larger extinction and AOD in the winter? Our conjecture is that the lower dust plume height in the winter provides more opportunities for effective mixing of dust particles and local pollution aerosols near the surface. Consequently, the large extinction and AOD observed in the winter are caused by the combined contributions of dust particles and other aerosol types near the surface. In contrast, the mean base height of 1.59 km in the spring reduces the probability of the mixing of dust and near-ground aerosols.The meanδ of 0.10 in winter is slightly smaller than that in spring (0.12). Particularly, the value of δ
p in winter (0.18) is much smaller than that in spring (0.27). This finding is potentially attributed to the mixing of dust particles and water vapor in the lower troposphere. Figure 11 shows the correlation between δ
p and RH in the spring and winter. The δ
p and RH are weakly linearly correlated in the spring (correlation coefficient of −0.19); however, the δ
p and RH exhibit a strong, negative linear correlation in the winter (correlation coefficient of −0.71), which indicates the strong effect of RH on the δ
p. In the winter, plenty of dust aerosols transported at a relatively low horizontal velocity and low height are determined by regional weather systems and topography. Therefore, dust particles have abundant time to mix with the local water vapor and pollution aerosols. Dust particles may interact with acidic pollutants (such as SO2 and NO
x) under high RH conditions. The neutralization process (such as the reaction between CaCO3 and H2SO4) promotes the hygroscopic growth of dust particles [16, 25]. Under high relative humidity conditions, activated dust particle would experience hygroscopic growth to produce enhanced scattering and extinction [25, 26, 49], as well as reducing the depolarization ratio [16, 50, 51]. Murayama et al. [50] observed the δ
p decreased with the increasing of RH. Ikegami et al. [51] found that dust particles could be covered with sulfuric acid solution in the troposphere. The internal mixture would lead to the change of depolarization ratio and its dependence to RH.Scatter diagrams of the correlation between particle depolarization ratio and relative humidity measured during winter (a) and spring (b) dust events over Wuhan from October 2010 to September 2012. The relative humidity and particle depolarization ratio were simultaneously measured using the radiosondes and ground-based depolarization lidar. Here,R represents the linear correlation coefficient.
(a)
(b)With the prevailing westerlies, most Asian dust events occur in the spring [14–17]; few events occur in the winter [52, 53]. Kim et al. [53] and Kim [54] observed that more than 80% of dust emissions occur in the spring, based on long-term observations made in Korea, which indicates that dust affects Japan and Korea mostly in the spring. In contrast, approximately half of the dust events over Wuhan are observed in the winter, indicating that winter is another season when plenty of intrusive dust plumes occur. Two years of CALIPAO level 1B data (the same period of the polarization lidar observations) were employed to validate the frequency of seasonal dust events (Table 4). Two regions were chosen for the comparison: southeast China (28°–33°N, 109°–119°E) and the area of Korea and Japan (33°–38°N, 127°–137°E). Wuhan is located in the middle of the selected southeast China region. Cloud screening was first performed using CALIPSO level 2 333 m and 5 km cloud layer products following the method of Liu et al. [23]. After screening for clouds, there remained 6230 and 8537 5 km profiles for southeast China and the Japan-Korea area, respectively. In southeast China, the dust occurrence is highest in the spring and slightly lower in winter; the frequency is lowest in the summer. The Korea-Japan region exhibits similar seasonal variations; however, for every season, the occurrence of dust is relatively less there than in southeast China. There is another difference between the two regions. The dust occurrence probability over southeast China in the winter (0.68) is nearly the same as that as in the spring (0.73). However, the dust occurrence probability in the Korea-Japan region in the winter (0.26) is slightly less than half of that in the spring (0.56). These findings based on two years of CALIPSO data agree with other reports from Korea and Japan and the results from our ground-based depolarization lidar. In the winter, dust from Mongolia and northern China is always entrained by more southeastward cold frontal systems.Table 4
Comparison of seasonal dust occurrences in southeast China (25–35°N, 110–118°E) and Korea-Japan (33°–38°N, 127°–137°E), based on CALIPSO data from October 2010 to September 2012.
Dust occurrence
Spring
Summer
Fall
Winter
Southeast China
0.73
0.37
0.56
0.68
Korea-Japan
0.56
0.13
0.25
0.26


### 4.2. Comparisons of Optical Properties with Other Observations

The volume depolarization ratios listed in Table2 are smaller than those of other Asian dust events reported in previous studies. The dust aerosol volume depolarization ratios are 0.21 on the Tibet Plateau [20], 0.17–0.25 in Japan [15], 0.18 in Korea [55], 0.35 in Beijing [56], and 0.19 in Taiwan [25]. The dust plumes are always transported at a low horizontal velocity as they leave the westerly jet and move southward to Wuhan. Enough time is provided for the dust particles to externally and internally mix with other urban aerosols. Cottle et al. [57] observed the maximum volume depolarization ratio within the dust layers reduced substantially from 0.27 over Vancouver to 0.15 over Egbert by the time the dust layers travelled across North America during April 2010. They explained this finding that was the result of continued loss of dust, external mixing with the background aerosol layer or more likely a combination of both. Large quantity of pollutant (such as soot, sulfate, and nitrate) is released due to the heating from November to next March over northern China. The atmospheric diffusion condition in winter is usually unfavorable, which causes the persistence of aerosols aggregation [58]. As a result, a plenty of anthropogenic aerosols intrude into the transported dust plume and weaken its overall nonspherical shape property. Local emissions from vehicles and coal-fired power plants, industrial waste gas, and city construction plant dust yield abundant local pollution aerosols over Wuhan and nearby areas. Thus the lidar-observed dust plumes over Wuhan were actually mixed with those aerosols that intruded into during transport and were produced locally. Chemical processes may occur between the dust particles and the acidic pollution under suitable RH conditions. Then, the hydrophobic dust particles may become hydrophilic [25]. These mixing and chemical processes may make the dust particles more spherical (reducing the particle depolarization ratio).The extinctions measured over Wuhan are larger than those measured at other sites. The mean extinction measured on 29–31 May 2008 in Seoul, Korea, was 0.08–0.12 km−1. Further statistical dust characteristics were examined based on 14 dust events occurring between April 2006 and March 2009. The extinctions ranged from 0.10 km−1 to 0.20 km−1; the mean was 0.14 km−1 [22]. A similar extinction range (0.02–0.23 km−1) was measured using Raman lidar in Thessaloniki, Greece, based on 12 cases observed during the period 2001-2002 [59]. The mean AOD (0.32) in the dust layer over Wuhan is relatively large compared to other downstream regions. Liu et al. [39] observed an extensive dust storm originating on 17 August 2006 in North Africa using CALIPSO lidar. Dust layer AODs (0.08-0.09) were measured on 28 August 2006 in the Gulf of Mexico using the NASA high spectral resolution lidar. Papayannis et al. [60] measured dust layer AODs (0.10–0.25) using 8 lidar stations in Europe from May 2000 to December 2002. The extinctions and AODs over Wuhan were larger than those reported in other studies, which was caused by high local aerosol loadings. Pure dust particles are removed over time with continual sedimentation during transport. Hence, dust particles contribute increasingly less to the total extinction and AOD, while other aerosols inside the dust layer contribute increasingly more. Hänel et al. [61] generally observed an elevated aerosol layer at 2.5–5 km height caused by the complex mixture of aged desert dust, biomass burning smoke, and industrial pollution over eastern Asia at the Global Atmospheric Watch station of Shangdianzi near Beijing, China. In conclusion, the mixture of dust aerosol and locally emitted aerosols cause the large extinction and large AOD in Wuhan.


## 4.1. Different Characteristics in Spring and Winter

As mentioned above, the dust layers are higher in spring than in winter. Two likely explanations are presented below. First, Asian dust is generated by cold air outbreaks and driven by westerly jets. Asian dust outbreaks are more intense in the spring than in winter due to the higher frequency of strong winds. Kurosaki and Mikami [47] found good agreement in the yearly variations and spatial distributions between the surface wind and dust outbreaks over East Asia during the period 2000–2002. The stronger surface winds in the spring effectively carry the desert dust particles to higher altitudes in the free troposphere. The elevated dust aerosols are then transported over long distances. As a result, the spring dust plume over downstream regions (such as Wuhan) is higher than the winter plume. Second, the temperature inversion that occurs frequently in the lower troposphere over Wuhan during the winter restricts the vertical convective transportation of dust and other aerosols in the lower troposphere. Zhang et al. [48] observed two lower tropospheric inversion layers over Yichang (111°18′E, 30°42′N) in January 2007. Their observations reveal the inversion layers are closely related to the intensive inertial gravity waves propagation in winter tropospheric jets. The temperature profiles from the radiosondes were statistically analyzed over Wuhan. Of the 27 radiosondes launched during the 6 winter cases, a temperature inversion (a layer thicker than 0.12 km with a lapse rate greater than −2 K/km) at heights of 1.5–3 km was observed in 60% of the radiosonde measurements. Of the 32 radiosondes launched during the 5 spring cases, a temperature inversion in the same region was observed in only 28% of the radiosonde measurements. The probability of a temperature inversion in the winter was approximately twice that of spring. Therefore, it is likely that the dust particles are inhibited from leaving the lower troposphere via convective transport due to the restriction caused by the temperature inversion in the winter.As shown in Table3, the mean dust layer extinction and AOD in the spring are 0.10 km−1 and 0.26, respectively. These values increase by 140% and 31% in the winter, respectively. What caused the larger extinction and AOD in the winter? Our conjecture is that the lower dust plume height in the winter provides more opportunities for effective mixing of dust particles and local pollution aerosols near the surface. Consequently, the large extinction and AOD observed in the winter are caused by the combined contributions of dust particles and other aerosol types near the surface. In contrast, the mean base height of 1.59 km in the spring reduces the probability of the mixing of dust and near-ground aerosols.The meanδ of 0.10 in winter is slightly smaller than that in spring (0.12). Particularly, the value of δ
p in winter (0.18) is much smaller than that in spring (0.27). This finding is potentially attributed to the mixing of dust particles and water vapor in the lower troposphere. Figure 11 shows the correlation between δ
p and RH in the spring and winter. The δ
p and RH are weakly linearly correlated in the spring (correlation coefficient of −0.19); however, the δ
p and RH exhibit a strong, negative linear correlation in the winter (correlation coefficient of −0.71), which indicates the strong effect of RH on the δ
p. In the winter, plenty of dust aerosols transported at a relatively low horizontal velocity and low height are determined by regional weather systems and topography. Therefore, dust particles have abundant time to mix with the local water vapor and pollution aerosols. Dust particles may interact with acidic pollutants (such as SO2 and NO
x) under high RH conditions. The neutralization process (such as the reaction between CaCO3 and H2SO4) promotes the hygroscopic growth of dust particles [16, 25]. Under high relative humidity conditions, activated dust particle would experience hygroscopic growth to produce enhanced scattering and extinction [25, 26, 49], as well as reducing the depolarization ratio [16, 50, 51]. Murayama et al. [50] observed the δ
p decreased with the increasing of RH. Ikegami et al. [51] found that dust particles could be covered with sulfuric acid solution in the troposphere. The internal mixture would lead to the change of depolarization ratio and its dependence to RH.Scatter diagrams of the correlation between particle depolarization ratio and relative humidity measured during winter (a) and spring (b) dust events over Wuhan from October 2010 to September 2012. The relative humidity and particle depolarization ratio were simultaneously measured using the radiosondes and ground-based depolarization lidar. Here,R represents the linear correlation coefficient.
(a)
(b)With the prevailing westerlies, most Asian dust events occur in the spring [14–17]; few events occur in the winter [52, 53]. Kim et al. [53] and Kim [54] observed that more than 80% of dust emissions occur in the spring, based on long-term observations made in Korea, which indicates that dust affects Japan and Korea mostly in the spring. In contrast, approximately half of the dust events over Wuhan are observed in the winter, indicating that winter is another season when plenty of intrusive dust plumes occur. Two years of CALIPAO level 1B data (the same period of the polarization lidar observations) were employed to validate the frequency of seasonal dust events (Table 4). Two regions were chosen for the comparison: southeast China (28°–33°N, 109°–119°E) and the area of Korea and Japan (33°–38°N, 127°–137°E). Wuhan is located in the middle of the selected southeast China region. Cloud screening was first performed using CALIPSO level 2 333 m and 5 km cloud layer products following the method of Liu et al. [23]. After screening for clouds, there remained 6230 and 8537 5 km profiles for southeast China and the Japan-Korea area, respectively. In southeast China, the dust occurrence is highest in the spring and slightly lower in winter; the frequency is lowest in the summer. The Korea-Japan region exhibits similar seasonal variations; however, for every season, the occurrence of dust is relatively less there than in southeast China. There is another difference between the two regions. The dust occurrence probability over southeast China in the winter (0.68) is nearly the same as that as in the spring (0.73). However, the dust occurrence probability in the Korea-Japan region in the winter (0.26) is slightly less than half of that in the spring (0.56). These findings based on two years of CALIPSO data agree with other reports from Korea and Japan and the results from our ground-based depolarization lidar. In the winter, dust from Mongolia and northern China is always entrained by more southeastward cold frontal systems.Table 4
Comparison of seasonal dust occurrences in southeast China (25–35°N, 110–118°E) and Korea-Japan (33°–38°N, 127°–137°E), based on CALIPSO data from October 2010 to September 2012.
Dust occurrence
Spring
Summer
Fall
Winter
Southeast China
0.73
0.37
0.56
0.68
Korea-Japan
0.56
0.13
0.25
0.26


## 4.2. Comparisons of Optical Properties with Other Observations

The volume depolarization ratios listed in Table2 are smaller than those of other Asian dust events reported in previous studies. The dust aerosol volume depolarization ratios are 0.21 on the Tibet Plateau [20], 0.17–0.25 in Japan [15], 0.18 in Korea [55], 0.35 in Beijing [56], and 0.19 in Taiwan [25]. The dust plumes are always transported at a low horizontal velocity as they leave the westerly jet and move southward to Wuhan. Enough time is provided for the dust particles to externally and internally mix with other urban aerosols. Cottle et al. [57] observed the maximum volume depolarization ratio within the dust layers reduced substantially from 0.27 over Vancouver to 0.15 over Egbert by the time the dust layers travelled across North America during April 2010. They explained this finding that was the result of continued loss of dust, external mixing with the background aerosol layer or more likely a combination of both. Large quantity of pollutant (such as soot, sulfate, and nitrate) is released due to the heating from November to next March over northern China. The atmospheric diffusion condition in winter is usually unfavorable, which causes the persistence of aerosols aggregation [58]. As a result, a plenty of anthropogenic aerosols intrude into the transported dust plume and weaken its overall nonspherical shape property. Local emissions from vehicles and coal-fired power plants, industrial waste gas, and city construction plant dust yield abundant local pollution aerosols over Wuhan and nearby areas. Thus the lidar-observed dust plumes over Wuhan were actually mixed with those aerosols that intruded into during transport and were produced locally. Chemical processes may occur between the dust particles and the acidic pollution under suitable RH conditions. Then, the hydrophobic dust particles may become hydrophilic [25]. These mixing and chemical processes may make the dust particles more spherical (reducing the particle depolarization ratio).The extinctions measured over Wuhan are larger than those measured at other sites. The mean extinction measured on 29–31 May 2008 in Seoul, Korea, was 0.08–0.12 km−1. Further statistical dust characteristics were examined based on 14 dust events occurring between April 2006 and March 2009. The extinctions ranged from 0.10 km−1 to 0.20 km−1; the mean was 0.14 km−1 [22]. A similar extinction range (0.02–0.23 km−1) was measured using Raman lidar in Thessaloniki, Greece, based on 12 cases observed during the period 2001-2002 [59]. The mean AOD (0.32) in the dust layer over Wuhan is relatively large compared to other downstream regions. Liu et al. [39] observed an extensive dust storm originating on 17 August 2006 in North Africa using CALIPSO lidar. Dust layer AODs (0.08-0.09) were measured on 28 August 2006 in the Gulf of Mexico using the NASA high spectral resolution lidar. Papayannis et al. [60] measured dust layer AODs (0.10–0.25) using 8 lidar stations in Europe from May 2000 to December 2002. The extinctions and AODs over Wuhan were larger than those reported in other studies, which was caused by high local aerosol loadings. Pure dust particles are removed over time with continual sedimentation during transport. Hence, dust particles contribute increasingly less to the total extinction and AOD, while other aerosols inside the dust layer contribute increasingly more. Hänel et al. [61] generally observed an elevated aerosol layer at 2.5–5 km height caused by the complex mixture of aged desert dust, biomass burning smoke, and industrial pollution over eastern Asia at the Global Atmospheric Watch station of Shangdianzi near Beijing, China. In conclusion, the mixture of dust aerosol and locally emitted aerosols cause the large extinction and large AOD in Wuhan.


## 5. Summary and Conclusions

In this study, two years of Asian dust transported over long distances were simultaneously observed using a ground-based depolarization lidar and CALIPSO at Wuhan (30.5°N, 114.4°E), China. The vertical distribution, horizontal range, and optical properties of the dust in 13 dust events were examined from October 2010 to September 2012. The results are summarized as follows.The first case, which was a long-duration, intense dust event, occurred on 22–29 April 2011. The dust layer primarily occurred at a height of 1.52–5.02 km and had a mean centroid height of 3.04 km. The mean volume depolarization ratio, particle depolarization ratio, extinction, and optical depth were 0.12, 0.24, 0.11 km−1, and 0.38, respectively. The double-layer structure observed using the depolarization lidar during the period 22–25 April was also seen in the CALIPSO data at heights of 1–4 km and 3–9 km. The horizontal ranges of the dust layers were ~950 km and ~750 km on 23 April. The second case, which occurred on 7–9 December 2010, contained a single layer extending from the ground surface to a height of 2.22 km and had a mean centroid height of 1.02 km. The mean volume depolarization ratio, particle depolarization ratio, extinction, and optical depth were 0.15, 0.27, 0.13 km−1, and 0.26, respectively. The dust plume was observed at a height of 0–2 km and had a horizontal range of ~1300 km, based on the CALIPSO data.The statistical characteristics of the dust were also examined. The dust primarily occurred at a height of approximately 1.45–3.51 km and had a mean centroid height of 2.34 km. A meanδ of 0.12, as well as the mean δ
p of 0.22, was observed, indicating the mixing of other spherical aerosols and the hygroscopic growth of the dust particles. The extinction and AOD exhibited large variations; the mean extinction and AOD were 0.19 km−1 and 0.32, respectively. The CALIPSO data indicate that the dust plumes had horizontal ranges of approximately 750–2400 km in 9 cases.The dust occurrence probability over Wuhan in the winter (0.68) is slightly smaller than that in spring (0.73). Moreover, there are a few dust features that differ between the spring and winter. The springtime dust layers tend to be higher due to more intense Asian dust outbreaks, which are caused by the higher frequency of strong winds. The mean extinction and AOD are 140% and 31% larger in the winter, respectively. Moreover, theδ and δ
p are 20% and 50% smaller in the winter, respectively. These differences may be due to the temperature inversion in the lower troposphere, the mixing of local aerosols, and hygroscopic growth under suitable RH conditions. Compared with other downstream dust regions, the dust aerosols over Wuhan exhibit a larger extinction and AOD and a smaller δ.Dust largely affects Wuhan during the spring and winter every year, and different spring and winter dust properties were observed in this study. However, the precise explanations for these different characteristics require more extensive study and further analysis. The optical properties of dust obtained only by polarization lidar are limited. The lidar ratio and Ångström exponent are two other important parameters to depict the microphysics properties of dust aerosols, which may help us to further understand the particle size and absorbing property. Therefore, multiwavelength Raman/polarization lidar observations are needed at out site. We will perform such studies in the future.



---

*Source: 536762-2015-01-26.xml*