Seasonal Variation of Dust Aerosol Vertical Distribution in Arctic Based on Polarized Micropulse Lidar Measurement

: This study investigates the seasonal variation of dust aerosol vertical distribution using polarized Micropulse lidar (MPL) measurements at the Atmospheric Radiation Measurement (ARM) North Slope of Alaska (NSA) observatory from January 2013 to September 2017. For the ﬁrst time, multi-year aerosol backscatter coefﬁcients are retrieved at the ARM NSA site from MPL measurements and are consistent with co-located high spectral resolution lidar (HSRL) measurements. The high-quality aerosol backscatter coefﬁcient retrievals are used to derive the particle depolarization ratio (PDR) at the wavelength of 532 nm, which is used to identify the presence of dust aerosols. The annual cycles of the vertical distributions of dust backscatter coefﬁcient and PDR and dust aerosol optical depth (DAOD) show that aerosol loading has a maximum in late winter and early spring but a minimum in late summer and early autumn. Vertically, dust aerosol occurs in the entire troposphere in spring and winter and in the low and middle troposphere in summer and autumn. Because dust aerosols are effective ice nuclei, the seasonality of dust aerosol vertical distribution has important implications for the Arctic climate through aerosol–cloud–radiation interactions, primarily through impacting mixed-phase cloud processes.


Introduction
The Arctic region has been undergoing tremendous changes, including rising temperature and melting sea ice over the past several decades [1].However, the potentially critical processes have not yet been fully characterized in some model schemes, especially for Arctic aerosols and their effect on the radiance balance [2][3][4][5], which induces significant uncertainties in predicting future Arctic climate changes [6,7].
Dust aerosol, as one of important aerosol types, affects the radiative budget directly, and modifies cloud properties and lifetimes indirectly by acting as ice nuclei particles (INPs) [8,9] and cloud condensation nuclei (CCN) [10].For example, dust aerosols can significantly enhance ice concentrations in mixed-phased clouds and further impact the properties of mixed-phased clouds [11,12], which are the dominated Arctic clouds, particularly distinct in spring and fall [11,[13][14][15][16] and also contribute large equilibrium climate sensitivity in the general circulation model (GCM) [17][18][19].The vertical distribution of aerosols is one of the key factors to disentangle their effects on earth due to their importance for radiative effects [20], surface air quality studies [21] and model validation [22,23].And there is strong seasonal cycle in aerosol optical properties over Arctic [6,24,25].Considering the high variability of aerosol properties in time and space, it is necessary to evaluate their impacts in climate models [26].Thus, the information on detailed seasonal cycle of tropospheric dust aerosols including its vertical distribution is critical to improve our understanding of dust aerosol impacts on Arctic climate and environment [24,27].
Arctic aerosols were reported when Mitchell et al. [28] investigated Arctic haze sixty years ago.The long-term continuous surface observations began in the mid-1970s at Arctic sites [29], while the most of lidar measurements are limited to short-term campaigns [30][31][32][33][34][35].Surface measurements showed that the scattering and absorption of Arctic aerosols are the largest in late winter and early spring and smallest in summer [29].Treffeisen et al. [36] used Stratospheric Gas Experiment (SAGE) to provide the first climatology of Arctic upper tropospheric aerosol extinction with the maximum through the spring and a base level value in mid-summer.Nevertheless, SAGE measurements are not available in the low and middle troposphere, which can be provided with lidar measurements.The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) lidar measurements are used widely for Arctic aerosol characterizations.Devasthale et al. [35] characterized the macro-physical properties of aerosol layers over the Arctic based on the 4-year CALIPSO data and their results showed that the most part of aerosol are confined in the boundary layer in winter and summer.Di Pierro et al. [37] presented the first pan-Arctic view of Arctic haze seasonality in the low troposphere using 6-year CLIPSO data: maximum during winter and early spring and minimum in summer and early autumn.Yang et al. [38] analyzed the aerosol optical depth (AOD) with 13-year CALIPSO data, and showed the largest from late autumn to early spring over the Arctic.In addition, 4-year ground-based Mie polarized lidar was used to show that the highest aerosol concentration is in late spring and summer and the lowest is in late summer and fall, while the largest aerosol particle depolarization ratio (PDR) is in winter and spring and the smallest is in summer at Ny-Ålesund site [39].Zhang et al. [40] used the 6-year high spectra resolution lidar (HSRL) and found that Arctic haze plays a governing role in later winter and spring at the North Slope of Alaska (NSA) site.Overall, the seasonality of aerosol properties at Arctic sites depends on the geographical locations, which are driven by long-range transportation patterns in winter and spring in combination with frequent precipitation removal mechanisms in summer [26,41,42].However, dust aerosols over Arctic region attract little attention.
It is known that the considerable dust events occur over Arctic [43,44].They are from either the long-range transportation of low-latitude dust aerosols [44,45], or the cold and high latitudes [26,43].Most of the vertical dust aerosol observation studies in the Arctic are limited to short-term campaigns [46], and the studies on seasonality of Arctic dust aerosol vertical distribution are still missing.Although micro-pulse lidar (MPL) has been operated continuously at Atmospheric Radiation Measurement (ARM) NSA (71.3 • N, 156.6 • W) observatory since 1998 and polarization measurement was added in 2006, the polarization measurement has large uncertainties, especially in middle and upper troposphere.Our recent work on the MPL afterpulse correction greatly improved not only the MPL middle and upper tropospheric aerosol backscatter observations, but also the MPL polarization measurements [47].The latter is critically important to detect nonspherical particles such as dust aerosols and ice particles using the PDR.Based on these facts, this study presented the long-term and reliable MPL retrieval results at the ARM NSA observatory for the first time, and the seasonal variations of dust aerosol vertical distributions are documented.
In Section 2, the technical aspects of the lidar measurements, retrieval method and the dust aerosol identification method are described.The validation of MPL retrievals with collocated HSRL observations and the seasonal variations of dust aerosols are given in Section 3. In Section 4, the summary and conclusions are presented.

Instrumentation and Datasets
The ARM NSA site, which is an ideal location for understanding Arctic aerosol processes [48], is located to the northeast of Utqiagvik (formerly Barrow), Alaska (71.3 • N, 156.6 • W, 8 m ASL).There has been a suite of atmospheric measurements to observe the changing Arctic since 1997.In this study, the main datasets used are the measurements from MPL (532 nm), HSRL (532 nm) and Ka-band ARM Zenith Radar (KAZR, 35 GHz) at the NSA observatory.
MPL is an eye-safe lidar system that provides highly sensitive continuous measurements of vertical backscattering profiles from optically thin clouds, aerosols and molecules.A polarized MPL system had been operating since August 2006 at the NSA observatory.Starting in September 2010, the MPL was upgraded to a "fast-switching" system to improve the temporal matching between co-pol and cross-pol polarization channel measurements [49].The upgraded MPL has been used to avoid the high depolarization values at the NSA site since October 2017 [50], which were caused by the afterpulse noise from the MPL avalanche photodiode (APD) detector and have been corrected with the new afterpulse method developed by Xie et al. [47].More detailed extensive descriptions of the MPL system are given by Campbell et al. [51].
HSRL technique takes advantage of the spectral distribution of the lidar return signal to separate aerosol and molecular signals [52].As a result, the calibrated particulate backscatter, extinction coefficient and depolarization ratio are derived using the retrieval scheme for the instrument with a vertical/temporal resolution of 30 m/30 s [53].HSRL measurements are used as the reference for validating the retrieval result from the elastic lidars [54].The HSRL at the NSA site worked from March 2011 to July 2019.There is much missing data from 2011 to 2012.Thus, HSRL measurements after January 2013 are analyzed and only the HSRL backscatter coefficient observations are used to validate the MPL backscatter retrievals due to the noisy HSRL aerosol extinction and depolarization ratio measurements.
The KAZR measures the spectra of backscattered power (reflectivity) as a function of Doppler velocity of the cloud and precipitation particles in the atmosphere, and started in November 2011 at the NSA site.The best-estimated radar reflectivity factor in the Active Remote Sensing of Cloud Layers (ARSCL) algorithm [55] was used in this study to better separate optically thin ice clouds and precipitation from aerosols because KAZR is more sensitive to large particles.Here, the four-year and nine-month data from January 2013 to September 2017 were selected to document the seasonal vertical distribution of dust aerosols over the Arctic after considering the overlap period from multi-instrument observations.In addition, the air pressure, air temperature, and relative humidity from radiosondes at the NSA site were used to calculate lidar molecular backscatter signals [56].

Automated Retrievals of Aerosol Optical Properties with MPL Measurements
The first step to obtain aerosol optical products with MPL measurements is the preprocessing of the data, i.e., the raw signal corrections, including background noise, afterpulse, overlap and emitted laser energy normalization corrections [47,51].Then, a one-hour averaging is performed to increase the signal noise ratio (SNR) for middle and upper tropospheric aerosol detection, and the atmospheric feature is detected [57].Here, the signals are used only when SNR > 3 during the nighttime and SNR > 1.5 during the daytime.The primary classification of cloud and aerosol is undertaken to select the proper lidar ratio using the Fernald method [58], based on attenuated backscattering ratio (threshold is set as 3).Then, particular backscatter coefficients at 532 nm are determined using the Fernald method [58] (with a fixed lidar ratio of 55 sr as a good representative for the most aerosol types [59]).After retrieval, the cloud and aerosol are separated based on height-dependent thresholds determined from the height-dependent probability density functions (PDFs) of particular backscatter coefficients, because clouds are optically thicker than aerosol layers at the same altitude [60].To further remove potential misclassification of optically thin clouds or precipitation as aerosols, the collocated KAZR reflectivity measurements are used and pixels are removed when reflectivity is more than −50 dBZ, which could indicate low concentration but large hydrometers.
The critical technique in the Fernald method for cloud-free profiles is to find the proper reference height [61].It is assumed that there is only molecular at the location of the reference height, where the range-corrected signal in height range X(r) should be completely proportional to the attenuated molecular backscatter coefficient β att m (r).
where E is the output energy monitor pulse, C is the calibration constant of the instrument, β p (r) and β m (r) are particulate and molecular backscatter coefficients, respectively, and α p (r) and α m (r) are particulate and molecular extinction coefficients, respectively.
This processing is called the Rayleigh fits [62].The Rayleigh fits can find all possible reference height intervals.The quality tests are needed to find the optimum reference height in each profile.The more details about the quality tests are described in the appendix of Baars et al. [61].

Dust Aerosol Identification
The PDR (δ p ) [63] is calculated as: where δ m is the molecular volume depolarization ratio (VDR), which is a constant with height [64].It is set to be 0.05 in this study.δ v is the total VDR, which is calculated based on the afterpulse-corrected (δ v_AC ) depolarization ratio [47,65]: where G dep is the electro-optical gain ratio between the copol and crosspol channels and is set to be 0.67 empirically in this study.The R is the scattering ratio.
According to the Equation ( 3), the systematic uncertainties of PDR are from R and δ v and δ m .In general, the uncertainty from δ m is always small, so it can be neglected.The propagation factor is typically between 1 and 2 for the δ v uncertainty.It is noticeable that the propagation factor for uncertainties in R varies significantly with the R and can be significantly larger when there is fairly small R (R < 1.1) [65,66].In addition, the default range of PDR is 0-0.6 [67].The PDR is set as NAN once R < 1.1 or PDR > 0.6 in this study.
High PDR is an indicator of non-spherical particles [68], thus, PDRs are widely used to identify dust aerosols [65,[69][70][71].Figure 1 shows the PDF of PDR for all aerosols.The PDF for all aerosols is a left-skewed distribution from 0.0-0.6 with a peak location at 0.03, which includes the hypothetical contributions from dust and non-dust aerosols as indicated with the black and green lines, respectively.Observations show that the PDR of the non-dust aerosol (e.g., marine, smoke, urban) is always 0.0-0.1 [67,[72][73][74][75][76][77].Liu et al. [78] showed that pure Africa dust aerosols have a mean PDR about 0.3 based on CALIPSO lidar measurements.Amiridis et al. [79] showed the peak of PDR distribution for pure dust aerosols is 0.28 with a range mainly between 0.15 and 0.4, while the peak of polluted dust distribution is 0.09 with a range mainly between 0.05 and 0.25.Dust aerosols in polar regions are from local and long-range transportation sources [43,80].Thus, Arctic dust aerosols have different levels of mixing states including pure dust, polluted dust and dusty marine, which results in a wide range of observed dust aerosol PDR.It is obvious that there are overlaps between PDR distributions of non-dust aerosols and dust aerosols.A threshold of 0.09 was selected to empirically identify dust aerosol according to the established literature in this work [67,74].Uncertainties associated with the threshold selection are considered in the results.cated with the black and green lines, respectively.Observations show that the PDR of the non-dust aerosol (e.g., marine, smoke, urban) is always 0.0-0.1 [67,[72][73][74][75][76][77].Liu et al. [78] showed that pure Africa dust aerosols have a mean PDR about 0.3 based on CALIPSO lidar measurements.Amiridis et al. [79] showed the peak of PDR distribution for pure dust aerosols is 0.28 with a range mainly between 0.15 and 0.4, while the peak of polluted dust distribution is 0.09 with a range mainly between 0.05 and 0.25.Dust aerosols in polar regions are from local and long-range transportation sources [43,80].Thus, Arctic dust aerosols have different levels of mixing states including pure dust, polluted dust and dusty marine, which results in a wide range of observed dust aerosol PDR.It is obvious that there are overlaps between PDR distributions of non-dust aerosols and dust aerosols.A threshold of 0.09 was selected to empirically identify dust aerosol according to the established literature in this work [67,74].Uncertainties associated with the threshold selection are considered in the results.

Figure 1.
The probability density functions (PDF) of particle depolarization ratio (PDR) for all aerosols observed by MPL from January 2013 to September 2017 at the NSA site.The black line is for all the aerosols derived from MPL observations, while the red and green lines indicate non-dust and dust aerosol modes fitted to the shape of the PDF of all aerosol PDR.A PDR of 0.10 is the cross point of the dust and non-dust aerosol modes.

The Retrieval Case and Validation
To demonstrate the performance of the retrieval results, HSRL observations and MPL retrievals during 11-14 February 2015 are shown in Figure 2. The MPL (Figure 2b) can well capture the aerosol structures detected by HSRL (Figure 2a) although the peak values of strong aerosol layers are slightly high, and weak aerosol regions are slightly small.This reflects the larger uncertainties using MPL measurements than HSRL measurements due to assumptions related to boundary conditions and lidar ratio as discussed earlier.VDR shows the consistent aerosol layers with MPL backscatter coefficients.There are two strong aerosol layers within 0.0-4.0km on 11 February, followed by the two strong and

The Retrieval Case and Validation
To demonstrate the performance of the retrieval results, HSRL observations and MPL retrievals during 11-14 February 2015 are shown in Figure 2. The MPL (Figure 2b) can well capture the aerosol structures detected by HSRL (Figure 2a) although the peak values of strong aerosol layers are slightly high, and weak aerosol regions are slightly small.This reflects the larger uncertainties using MPL measurements than HSRL measurements due to assumptions related to boundary conditions and lidar ratio as discussed earlier.VDR shows the consistent aerosol layers with MPL backscatter coefficients.There are two strong aerosol layers within 0.0-4.0km on 11 February, followed by the two strong and three weak aerosol layers in 0.0-4.0km on 12-14 February, and a strong aerosol layer in 4.0-8.0km at about 11:00 UTC of 14 February, which may be from long-range transportation.PDR (Figure 2d) is calculated with particular backscatter coefficients (Figure 2b) and VDR (Figure 2c).PDR shows the aerosol layers more clearly.Dust aerosols (Figure 2e) are detected when PDR is more than 0.09.The strong aerosol layers in 0.0-4.0km on 11-12 and 15 February and in 4.0-8.0km at 11:00 UTC of 14 February are identified as dust aerosols using the selected threshold.
The statistical comparison of MPL retrievals with co-located HSRL observations is shown in Figure 3a based on all measurements from January 2013 to September 2017.The median backscatter coefficient profiles and the numeric range of main backscatter coefficient distributions with 25th-75th percentiles from MPL (red line and blue box) and HSRL (cyan blue and magenta box) are generally consistent, which indicates the high reliability of our MPL backscatter retrievals.However, MPL is larger in 0.0-3.0km and 7.2-9.0km, while HSRL is slightly larger in 3.0-7.2km.To show the relative difference between MPL and HSRL clearly, the ratios of MPL and HSRL median backscatter coefficients are shown in Figure 3b.The absolute differences are mainly within 20%, with the maximum below 40%.The differences may be due to the assumption of a constant lidar ratio for all aerosols in this study while it is a fact that different types of aerosols have different lidar ratios.However, lidar ratios with different aerosol types are difficult to implement with MPL measurements.On the other hand, issues with the instrument itself may also contribute to these differences.For example, both lidar systems use small fields of view, which leads to a significant overlap function below 5 km and introduces uncertainties in backscatter coefficient measurements.Considering these factors, the differences between MPL and HSRL aerosol measurement uncertainties at the NSA site are acceptable in this study.
Remote Sens. 2022, 14, x FOR PEER REVIEW 6 of 16 three weak aerosol layers in 0.0-4.0km on 12-14 February, and a strong aerosol layer in 4.0-8.0km at about 11:00 UTC of 14 February, which may be from long-range transportation.PDR (Figure 2d) is calculated with particular backscatter coefficients (Figure 2b) and VDR (Figure 2c).PDR shows the aerosol layers more clearly.Dust aerosols (Figure 2e) are detected when PDR is more than 0.09.The strong aerosol layers in 0.0-4.0km on 11-12 and 15 February and in 4.0-8.0km at 11:00 UTC of 14 February are identified as dust aerosols using the selected threshold.The statistical comparison of MPL retrievals with co-located HSRL observations is shown in Figure 3a based on all measurements from January 2013 to September 2017.The median backscatter coefficient profiles and the numeric range of main backscatter coefficient distributions with 25th-75th percentiles from MPL (red line and blue box) and HSRL (cyan blue and magenta box) are generally consistent, which indicates the high reliability of our MPL backscatter retrievals.However, MPL is larger in 0.0-3.0km and 7.2-9.0km, while HSRL is slightly larger in 3.0-7.2km.To show the relative difference between MPL and HSRL clearly, the ratios of MPL and HSRL median backscatter coefficients are shown in Figure 3b.The absolute differences are mainly within 20%, with the maximum below 40%.The differences may be due to the assumption of a constant lidar ratio for all aerosols in this study while it is a fact that different types of aerosols have different lidar ratios.However, lidar ratios with different aerosol types are difficult to implement with MPL measurements.On the other hand, issues with the instrument itself may also contribute to these differences.For example, both lidar systems use small fields of view, which leads to a significant overlap function below 5 km and introduces uncertainties in backscatter The 5th and 95th percentiles of MPL backscatter coefficients indicate broader distributions of MPL retrieved backscatter coefficients in the entire troposphere, which implies the larger random noises of MPL backscatter coefficient retrievals over HSRL backscatter coefficient observations.This is illustrated in Figure 2 and is consistent with similar differences between HSRL and elastic lidars (e.g., MPL and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP)) reported in other studies [54,81,82].The lower limit values (5th percentile) may be caused by the selections of reference heights or the constant selection of lidar ratio values under strong aerosol layers.The differences of noise-level and photon counting detection dynamical range between MPL and HSRL measurements could also contribute to these differences.
To our knowledge, the high quality MPL aerosol retrievals and their consistency with the co-located HSRL from the long-term observations are shown here for the first time.Most of the previous research has used range-corrected signals to show the retrieval results with mean seasonal profiles, or validated the retrieval results with column integral variables such as AOD [83,84].This highlights the value of long-term ARM MPL measurements at the NSA site for understanding Arctic climate changes through aerosol-cloud-radiation interaction.
Remote Sens. 2022, 14, x FOR PEER REVIEW 7 of 16 coefficient measurements.Considering these factors, the differences between MPL and HSRL aerosol measurement uncertainties at the NSA site are acceptable in this study.The 5th and 95th percentiles of MPL backscatter coefficients indicate broader distributions of MPL retrieved backscatter coefficients in the entire troposphere, which implies the larger random noises of MPL backscatter coefficient retrievals over HSRL backscatter coefficient observations.This is illustrated in Figure 2 and is consistent with similar differences between HSRL and elastic lidars (e.g., MPL and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP)) reported in other studies [54,81,82].The lower limit values (5th percentile) may be caused by the selections of reference heights or the constant selection of lidar ratio values under strong aerosol layers.The differences of noise-level and photon counting detection dynamical range between MPL and HSRL measurements could also contribute to these differences.
To our knowledge, the high quality MPL aerosol retrievals and their consistency with the co-located HSRL from the long-term observations are shown here for the first time.Most of the previous research has used range-corrected signals to show the retrieval results with mean seasonal profiles, or validated the retrieval results with column integral variables such as AOD [83,84].This highlights the value of long-term ARM MPL measurements at the NSA site for understanding Arctic climate changes through aerosolcloud-radiation interaction.

Seasonal Variations of Dust Aerosols
Figure 4 shows the seasonal cycle of MPL-derived dust aerosol backscatter coefficient and PDR profiles from January 2013 to September 2017 at the NSA site.This analysis is bases on cloud-free profiles.Dust aerosol is an important component of the Arctic haze.Dust aerosols occur in winter and spring in the entire troposphere, and in summer and autumn in the low (0-2 km) and middle (2-5 km) troposphere.These results are consistent

Seasonal Variations of Dust Aerosols
Figure 4 shows the seasonal cycle of MPL-derived dust aerosol backscatter coefficient and PDR profiles from January 2013 to September 2017 at the NSA site.This analysis is bases on cloud-free profiles.Dust aerosol is an important component of the Arctic haze.Dust aerosols occur in winter and spring in the entire troposphere, and in summer and autumn in the low (0-2 km) and middle (2-5 km) troposphere.These results are consistent with the seasonality of Arctic haze [29,37].In the upper troposphere (5-10 km) dust aerosols occurring during the late winter and early spring (January-April) are likely associated with the long-range transportation of natural and anthropogenic aerosols from low-latitude areas including Eurasia, Siberia, North America, and Southeast Asia [6,7,25,26].Further, Yang et al. [44] showed that polar jets can carry dust aerosols from source regions to the Arctic with CALIOP measurement.There are few low-dust aerosol occurrences and small aerosol samples (red line of Figure 4c) during August-October.A high occurrence of thin and opaque clouds during summer and autumn is responsible for the low samples during this period [36,85].Long-range dust transportation from middle latitudes is weak [44], but wet deposition is obvious in this period [86].Thus, dust aerosols during August-October are primarily from local and regional emissions [42].
Figure 5 shows the vertical distributions (every 150 m) of seasonal particulate backscatter coefficient (Figure 5a-d) and PDR PDFs (Figure 5e-h) based on MPL measurements at the NSA site from January 2013 to September 2017.The corresponding seasonal averaged profiles (black lines) of backscatter coefficient and PDR are over-plotted, separately.In summer and autumn, the backscatter coefficients of dust aerosols are relatively large, but PDR is small, implying that they may be the local dusty marine aerosols [87].In winter and spring, dust aerosols are detected in the entire troposphere.The mean profiles of aerosol backscatter coefficients and PDRs of dust aerosols show that the strongest are in middle and upper troposphere in spring.The vertical PDF distribution in spring is broader than in winter.This indicates that spring dust aerosols have the more complex mixture states of dust and non-dust aerosols due to the long-range transportation from the mid-latitudes [37,41,78,[88][89][90][91].
sols occurring during the late winter and early spring (January-April) are likely associated with the long-range transportation of natural and anthropogenic aerosols from low-latitude areas including Eurasia, Siberia, North America, and Southeast Asia [6,7,25,26].Further, Yang et al. [44] showed that polar jets can carry dust aerosols from source regions to the Arctic with CALIOP measurement.There are few low-dust aerosol occurrences and small aerosol samples (red line of Figure 4c) during August-October.A high occurrence of thin and opaque clouds during summer and autumn is responsible for the low samples during this period [36,85].Long-range dust transportation from middle latitudes is weak [44], but wet deposition is obvious in this period [86].Thus, dust aerosols during August-October are primarily from local and regional emissions [42].Figure 5 shows the vertical distributions (every 150 m) of seasonal particulate backscatter coefficient (Figure 5a-d) and PDR PDFs (Figure 5e-h) based on MPL measurements at the NSA site from January 2013 to September 2017.The corresponding seasonal averaged profiles (black lines) of backscatter coefficient and PDR are over-plotted, separately.In summer and autumn, the backscatter coefficients of dust aerosols are relatively large, but PDR is small, implying that they may be the local dusty marine aerosols [87].In winter and spring, dust aerosols are detected in the entire troposphere.The mean profiles of aerosol backscatter coefficients and PDRs of dust aerosols show that the strongest are in middle and upper troposphere in spring.The vertical PDF distribution in spring is broader than in winter.This indicates that spring dust aerosols have the more complex mixture states of dust and non-dust aerosols due to the long-range transportation from the mid-latitudes [37,41,78,[88][89][90][91]. Figure 6 shows the seasonal aerosol samples and dust occurrences.The number of the aerosol samples (black line in Figure 6) decreases with height for all seasons.The maximum amount of aerosol samples is in winter, and the minimum amount of aerosol samples is in autumn due to the seasonal cycle of clouds in the Arctic [7].For the dust aerosol occurrences using PDR threshold 0.09 (red line in Figure 6), the maximum occurs in spring and the minimum occurs in autumn [92].The maximum of occurrence profiles indicate that Arctic dust aerosols occur in three dust aerosol layers: 0-4 km, 4-8 km and 8-10 km during spring and winter, which is consistent with the example given in Figure 2. The dust aerosol occurrence is less than 0.1 in autumn and summer.
To quantify the uncertainty of the selected threshold on dust aerosol occurrences, the results within 30% of selected thresholds (0.09 ± 0.09 * 30%) (green, red, and blue lines in Figure 6) were compared.The PDR values can be interpreted as the different ratios of pure dust in all aerosols [93], which means that higher PDR values reflect the more dust aerosols mixed.For example, the dust occurrences with 0.06 are higher, especially below 4 km, which infers that there are more non-dust aerosols included.Below 4 km, the difference between dust aerosol occurrences with 0.12 and 0.09 are much smaller than the difference between 0.06 and 0.09.Above 4 km, dust occurrences have small difference among the three thresholds, which implies that dust aerosols are the dominant aerosol types in the Arctic upper troposphere.To balance between including as many dust aerosols as possible and avoiding non-dust aerosols, a PDR threshold of 0.09 was selected.Figure 6 shows the seasonal aerosol samples and dust occurrences.The number of the aerosol samples (black line in Figure 6) decreases with height for all seasons.The maximum amount of aerosol samples is in winter, and the minimum amount of aerosol samples is in autumn due to the seasonal cycle of clouds in the Arctic [7].For the dust aerosol occurrences using PDR threshold 0.09 (red line in Figure 6), the maximum occurs in spring and the minimum occurs in autumn [92].The maximum of occurrence profiles indicate that Arctic dust aerosols occur in three dust aerosol layers: 0-4 km, 4-8 km and 8-10 km during spring and winter, which is consistent with the example given in Figure 2. The dust aerosol occurrence is less than 0.1 in autumn and summer.
To quantify the uncertainty of the selected threshold on dust aerosol occurrences, the results within 30% of selected thresholds (0.09  0.09 * 30%) (green, red, and blue lines in Figure 6) were compared.The PDR values can be interpreted as the different ratios of pure dust in all aerosols [93], which means that higher PDR values reflect the more dust aerosols mixed.For example, the dust occurrences with 0.06 are higher, especially below 4 km, which infers that there are more non-dust aerosols included.Below 4 km, the difference between dust aerosol occurrences with 0.12 and 0.09 are much smaller than the To derive the dust aerosol optical depth (DAOD), the dust aerosol extinction profile is needed, which is calculated by the multiplication of the retrieved particulate backscatter coefficients and assumed lidar ratio.The DAOD is the vertical integration of extinction profiles in the entire height range.For the seasonal variation of DAOD (Figure 7), the maximum is in later winter and early spring and the minimum is in late summer and early autumn, which is consistent with the trend and values of DAOD in the previous studies [94,95].
Here, the different lidar ratios are used to show the effect of lidar ratios on monthly DAOD retrievals.Figure 7a shows that the seasonality of DAODs with different lidar ratios is almost consistent except for the absolute values.To prove the small effect of lidar ratio on the monthly DAOD retrievals, the ratio of DAOD with different lidar ratios (black and magenta lines in Figure 7b) is calculated.The ratio is almost 5% except for the values in July-October (5-10%), which may result from small samples during these months.difference between 0.06 and 0.09.Above 4 km, dust occurrences have small difference among the three thresholds, which implies that dust aerosols are the dominant aerosol types in the Arctic upper troposphere.To balance between including as many dust aerosols as possible and avoiding non-dust aerosols, a PDR threshold of 0.09 was selected.To derive the dust aerosol optical depth (DAOD), the dust aerosol extinction profile is needed, which is calculated by the multiplication of the retrieved particulate backscatter coefficients and assumed lidar ratio.The DAOD is the vertical integration of extinction profiles in the entire height range.For the seasonal variation of DAOD (Figure 7), the maximum is in later winter and early spring and the minimum is in late summer and early autumn, which is consistent with the trend and values of DAOD in the previous studies [94,95].
Here, the different lidar ratios are used to show the effect of lidar ratios on monthly DAOD retrievals.Figure 7a shows that the seasonality of DAODs with different lidar ratios is almost consistent except for the absolute values.To prove the small effect of lidar ratio on the monthly DAOD retrievals, the ratio of DAOD with different lidar ratios (black and magenta lines in Figure 7b) is calculated.The ratio is almost 5% except for the values in July-October (5-10%), which may result from small samples during these months.

Conclusions
A four-year and nine-month altitude-resolved distribution of MPL retrieved backscatter coefficient and PDR of Arctic dust aerosols at NSA site from January 2013 to September 2017 is presented.The long-term and stable MPL retrieved results and the comparison of MPL and co-located HSRL are shown for the first time, indicating the high quality of MPL retrievals.The dust aerosols are identified with PDR.The seasonality of Arctic dust aerosols shows the maximum is in late winter and early spring extending in the entire troposphere, which is from the long-range transportation of the mid-latitude dust aerosols, and the minimum is in the early summer and late autumn occurring in the low and middle troposphere, which originates from the local and regional sources.
Our understanding of the processes that control the Arctic aerosol properties remains highly uncertain.Research efforts are needed to constrain the seasonality and vertical profiles of Arctic aerosols to improve the model skills in simulation [7,22,96,97] and the ability of understanding aerosol-cloud-climate impacts [26].The high occurrence of dust aerosols in Arctic highlighted the importance to simulate dust aerosols and their impacts on mixed-phase and ice clouds reliably.Our multi-year temporal distribution of MPL backscatter coefficient and PDR over the Arctic will provide a new dataset to validate these processes in global models.The limitation of ground-based MPL measurements is just one-point observations.In the future, the satellite such as CALIOP and SAGE can be combined to provide more information for the Arctic dust aerosol distributions.

Figure 1 .
Figure 1.The probability density functions (PDF) of particle depolarization ratio (PDR) for all aerosols observed by MPL from January 2013 to September 2017 at the NSA site.The black line is for all the aerosols derived from MPL observations, while the red and green lines indicate non-dust and dust aerosol modes fitted to the shape of the PDF of all aerosol PDR.A PDR of 0.10 is the cross point of the dust and non-dust aerosol modes.

Figure 2 .
Figure 2. The time-height section of (a) HSRL particular backscatter coefficients, (b) MPL retrieved particular backscatter coefficients, (c) MPL total volume depolarization ratio (VDR), and (d) MPL particle depolarization ratio (PDR) from February 11 to 14, 2015 at the NSA site.(a-d) are for all aerosols.(e) MPL backscatter coefficient for dust aerosols that are identified using the PDR threshold 0.09.

Figure 2 .
Figure 2. The time-height section of (a) HSRL particular backscatter coefficients, (b) MPL retrieved particular backscatter coefficients, (c) MPL total volume depolarization ratio (VDR), and (d) MPL particle depolarization ratio (PDR) from February 11 to 14, 2015 at the NSA site.(a-d) are for all aerosols.(e) MPL backscatter coefficient for dust aerosols that are identified using the PDR threshold 0.09.

Figure 3 .
Figure 3. (a) The comparison of MPL (blue box) and HSRL (magenta box) backscatter coefficient vertical quartile distributions at the NSA site from January 2013 to September 2017.Box and whisker plots include the median, 25th and 75th percentiles (left and right of the boxes), and 5th and 95th percentiles (left and right of the whiskers).The red and cyan lines are for the medians of HSRL and MPL, respectively.(b) The ratio of MPL and HSRL backscatter medians.

Figure 3 .
Figure 3. (a) The comparison of MPL (blue box) and HSRL (magenta box) backscatter coefficient vertical quartile distributions at the NSA site from January 2013 to September 2017.Box and whisker plots include the median, 25th and 75th percentiles (left and right of the boxes), and 5th and 95th percentiles (left and right of the whiskers).The red and cyan lines are for the medians of HSRL and MPL, respectively.(b) The ratio of MPL and HSRL backscatter medians.

Figure 4 .
Figure 4.The climatological month-altitude sections of dust aerosols (a) particular backscatter coefficients and (b) PDR, and (c) aerosol samples (red line) and dust occurrence (blue line) for cloudfree profiles from MPL at the NSA site from January 2013 to September 2017.The dust occurrence is defined as the ratio of dust aerosol samples and all aerosol samples.The samples are based on the data with 1 h and 30 m bin.

Figure 4 .
Figure 4.The climatological month-altitude sections of dust aerosols (a) particular backscatter coefficients and (b) PDR, and (c) aerosol samples (red line) and dust occurrence (blue line) for cloudfree profiles from MPL at the NSA site from January 2013 to September 2017.The dust occurrence is defined as the ratio of dust aerosol samples and all aerosol samples.The samples are based on the data with 1 h and 30 m bin.

Figure 5 .
Figure 5.The vertical PDF distribution of seasonal particulate backscatter coefficients (a-d) and PDR for MPL (e-h) at the NSA site from January 2013 to September 2017.Rows from top to bottom are for spring (a,e), summer (b,f), autumn (c,g) and winter (d,h), separately.The seasonal average profiles (black lines) are plotted in each sub-figure.The color bar is the log scale of PDF distribution values of particular backscatter coefficients and PDR in each 150 m.

Figure 5 .
Figure 5.The vertical PDF distribution of seasonal particulate backscatter coefficients (a-d) and PDR for MPL (e-h) at the NSA site from January 2013 to September 2017.Rows from top to bottom are for spring (a,e), summer (b,f), autumn (c,g) and winter (d,h), separately.The seasonal average profiles (black lines) are plotted in each sub-figure.The color bar is the log scale of PDF distribution values of particular backscatter coefficients and PDR in each 150 m.

Figure 6 .
Figure 6.The aerosol samples (black line) in (a) spring, (b) summer, (c) autumn and (d) winter, and dust aerosol occurrences using the PDR thresholds of 0.06 (green line), 0.09 (red line) and 0.12 (blue line) for MPL measurements at the NSA site from January 2013 to September 2017.The samples are based on the data with 1 h and 30 m bin.

Figure 6 . 16 Figure 7 .
Figure 6.The aerosol samples (black line) in (a) spring, (b) summer, (c) autumn and (d) winter, and dust aerosol occurrences using the PDR thresholds of 0.06 (green line), 0.09 (red line) and 0.12 (blue line) for MPL measurements at the NSA site from January 2013 to September 2017.The samples are based on the data with 1 h and 30 m bin.Remote Sens. 2022, 14, x FOR PEER REVIEW 11 of 16 year and nine-month altitude-resolved distribution of MPL retrieved backscatter coefficient and PDR of Arctic dust aerosols at NSA site from January 2013 to September 2017 is presented.The long-term and stable MPL retrieved results and the comparison of MPL and co-located HSRL are shown for the first time, indicating the high

Figure 7 .
Figure 7.The seasonal cycle of (a) dust aerosol optical depth (DAOD) with lidar ratios of 40 sr (red line), 50 sr (green line) and 55 sr (blue line) for MPL measurements at the NSA site from January 2013 to September 2017.(b) The ratio of DAODs with different lidar ratios (cyan line is the reference line, black line is the ratio of 40 and 55, magenta line is the ratio of 50 and 55).