Investigating the quality of modeled aerosol profiles based on combined lidar and sunphotometer data

In this study we present an evaluation of the Comprehensive Air Quality Model with extensions (CAMx) for Thessaloniki using radiometric and lidar data. The aerosol mass concentration profiles of CAMx are compared against the PM2.5 and PM2. 5−10 concentration profiles retrieved by the Lidar-Radiometer Inversion Code (LIRIC). The CAMx model and the LIRIC algorithm results were compared in terms of mean mass concentration profiles, center of mass and integrated mass concentration in the boundary layer and the free troposphere. The mean mass concentration comparison resulted in profiles within the same order of magnitude and similar vertical structure for the PM2. 5 particles. The mean centers of mass values are also close, with a mean bias of 0.57 km. On the opposite side, there are larger differences for the PM2. 5−10 mode, both in the boundary layer and in the free troposphere. In order to grasp the reasons behind the discrepancies, we investigate the effect of aerosol sources that are not properly included in the model's emission inventory and in the boundary conditions such as the wildfires and the desert dust component. The identification of the cases that are affected by wildfires is performed using wind backward trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model in conjunction with satellite fire pixel data from MODerate-resolution Imaging Spectroradiometer (MODIS) Terra and Aqua global monthly fire location product MCD14ML. By removing those cases the correlation coefficient improves from 0.69 to 0.87 for the PM2. 5 integrated mass in the boundary layer and from 0.72 to 0.89 in the free troposphere. The PM2.5 center of mass fractional bias also decreases to 0.38 km. Concerning the analysis of the desert dust component, the simulations from the Dust Regional Atmospheric Model (BSC-DREAM8b) were deployed. When only the Saharan dust cases are taken into account, BSC-DREAM8b generally outperforms CAMx when compared with LIRIC, achieving a correlation of 0.91 and a mean bias of −29.1 % for the integrated mass in the free troposphere and a correlation of 0.57 for the center of mass. CAMx, on the other hand, underestimates the integrated mass in the free troposphere. Consequently, the accuracy of CAMx is limited concerning the transported Saharan dust cases. We conclude that the performance of CAMx appears to be best for the PM2.5 particles, both in the boundary layer and in the free troposphere. Sources of particles not properly taken into account by the model are confirmed to negatively affect its performance, especially for the PM2. 5−10 particles.The authors would like to acknowledge the EU projects MACC-III (Monitoring Atmospheric Composition and Climate – III, grant agreement no. 633080) and MACC-II project (Monitoring Atmospheric Composition and Climate – Interim Implementation, grant agreement no. 283576). The simulated results presented in this research paper have been produced using the EGI and HellasGrid infrastructures. The ACTRIS-2 project from the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 654109 is gratefully acknowledged. The authors would also like to acknowledge the support provided by the Scientific Computing Center at Aristotle University of Thessaloniki throughout the progress of the work on air quality forecasting. BSC-DREAM8b simulations were performed on the Mare Nostrum supercomputer hosted by Barcelona Supercomputing Center-Centro Nacional de Supercomputacion (BSC-CNS). S. Basart wants to acknowledge the CICYT project (CGL2013-46736). Elina Giannakaki acknowledges the support of the Academy of Finland (project no. 270108).Peer ReviewedPostprint (published version

Detection and characterization of birch pollen in the atmosphere using a multiwavelength Raman polarization lidar and Hirst-type pollen sampler in Finland

We present the results of birch pollen characterization using lidar and in situ measurements based on a 11 d pollination period from 5 to 15 May 2016 at the European Aerosol Research Lidar Network (EARLINET) station in Vehmasmäki (Kuopio; 62∘44′ N, 27∘33′ E), Finland. The ground-based multiwavelength Raman polarization lidar PollyXT performed continuous measurements at this rural forest site and has been combined with a Hirst-type volumetric air sampler, which measured the pollen type and concentration at roof level (4 m). The period was separated into two parts due to different atmospheric conditions and detected pollen types. During the first period, high concentrations of birch pollen were measured with a maximum 2 h average pollen concentration of 3700 grains m−3. Other pollen types represented less than 3 % of the total pollen count. In observed pollen layers, the mean particle depolarization ratio at 532 nm was 10±6 % during the intense birch pollination period. Mean lidar ratios were found to be 45±7 and 55±16 sr at 355 and 532 nm, respectively. During the second period, birch pollen was still dominant, but a significant contribution of spruce pollen was observed as well. Spruce pollen grains are highly nonspherical, leading to a larger mean depolarization ratio of 26±7 % for the birch–spruce pollen mixture. Furthermore, higher lidar ratios were observed during this period with mean values of 60±3 and 62±10 sr at 355 and 532 nm, respectively. The presented study shows the potential of the particle depolarization ratio to track pollen grains in the atmosphere.</p

Profiling water vapor mixing ratios in Finland by means of a Raman lidar, a satellite and a model

We present tropospheric water vapor profiles measured with a Raman lidar during three field campaigns held in Finland. Co-located radio soundings are available throughout the period for the calibration of the lidar signals. We investigate the possibility of calibrating the lidar water vapor profiles in the absence of co-existing on-site soundings using water vapor profiles from the combined Advanced InfraRed Sounder (AIRS) and the Advanced Microwave Sounding Unit (AMSU) satellite product; the Aire Limitee Adaptation dynamique Developpement INternational and High Resolution Limited Area Model (ALADIN/HIRLAM) numerical weather prediction (NWP) system, and the nearest radio sounding station located 100 km away from the lidar site (only for the permanent location of the lidar). The uncertainties of the calibration factor derived from the soundings, the satellite and the model data are <2.8, 7.4 and 3.9 %, respectively. We also include water vapor mixing ratio intercomparisons between the radio soundings and the various instruments/model for the period of the campaigns. A good agreement is observed for all comparisons with relative errors that do not exceed 50% up to 8 km altitude in most cases. A 4-year seasonal analysis of vertical water vapor is also presented for the Kuopio site in Finland. During winter months, the air in Kuopio is dry (1.15 +/- 0.40 g kg(-1)); during summer it is wet (5.54 +/- 1.02 g kg(-1)); and at other times, the air is in an intermediate state. These are averaged values over the lowest 2 km in the atmosphere. Above that height a quick decrease in water vapor mixing ratios is observed, except during summer months where favorable atmospheric conditions enable higher mixing ratio values at higher altitudes. Lastly, the seasonal change in disagreement between the lidar and the model has been studied. The analysis showed that, on average, the model underestimates water vapor mixing ratios at high altitudes during spring and summer.Peer reviewe

Optical characterization of pure pollen types using a multi-wavelength Raman polarization lidar

We present a novel algorithm for characterizing the optical properties of pure pollen particles, based on the depolarization ratio values obtained in lidar measurements. The algorithm was first tested and validated through a sim-ulator and then applied to the lidar observations during a 4-month pollen campaign from May to August 2016 at the Eu-ropean Aerosol Research Lidar Network (EARLINET) station in Kuopio (62 • 44 N, 27 • 33 E), in Eastern Finland. With a Burkard sampler, 20 types of pollen were observed and identified from concurrent measurements, with birch (Be-tula), pine (Pinus), spruce (Picea), and nettle (Urtica) pollen being the most abundant, contributing more than 90 % of the total pollen load, regarding number concentrations. Mean values of lidar-derived optical properties in the pollen layer were retrieved for four intense pollination periods (IPPs). Li-dar ratios at both 355 and 532 nm ranged from 55 to 70 sr for all pollen types, without significant wavelength dependence. An enhanced depolarization ratio was found when there were pollen grains in the atmosphere, and an even higher depo-larization ratio (with mean values of 0.25 or 0.14) was observed with the presence of the more non-spherical spruce or pine pollen. Under the assumption that the backscatter-related Ångström exponent between 355 and 532 nm should be zero for pure pollen, the depolarization ratio of pure pollen particles at 532 nm was assessed, resulting in 0.24±0.01 and 0.36 ± 0.01 for birch and pine pollen, respectively. Pollen optical properties at 1064 and 355 nm were also estimated. The backscatter-related Ångström exponent between 532 and 1064 nm was assessed to be ∼ 0.8 (∼ 0.5) for pure birch (pine) pollen; thus the longer wavelength would be a better choice to trace pollen in the air. Pollen depolarization ratios of 0.17 and 0.30 at 355 nm were found for birch and pine pollen, respectively. The depolarization values show a wavelength dependence for pollen. This can be the key parameter for pollen detection and characterization.</p

Spectral dependence of birch and pine pollen optical properties using a synergy of lidar instruments

Active remote sensors equipped with the capability to detect polarization, a shape-relevant parameter, are essential to aerosol particle identification in the vertical domain. Most commonly, the linear particle depolarization ratio has been available at the shorter wavelengths of 355 and/or 532 nm. Recently, linear particle depolarization ratios at longer wavelengths (910, 1064, and 1565 nm) have emerged in lidar aerosol research. In this study, a synergy of three lidars, namely a PollyXT lidar, a Vaisala CL61 ceilometer, and a HALO Photonics StreamLine Pro Doppler lidar, as well as in situ aerosol and pollen observations have been utilized to investigate the spectral dependence of birch and pine pollen particles. We found that, regardless of the pollen type, the linear particle depolarization ratio was subject to the amount of pollen and its relative contribution to the aerosol mixture in the air. More specifically, during birch pollination, characteristic linear particle depolarization ratios of 5 ± 2 % (355 nm), 28 ± 6 % (532 nm), 23 ± 6 % (910 nm), and 33 ± 4 % (1565 nm) were retrieved at the pollen layer. Regarding the pine-dominant period, characteristic linear particle depolarization ratios of 6 ± 2 %, 43 ± 11 %, 22 ± 6 %, and 26 ± 3 % were determined at wavelengths of 355, 532, 910, and 1565 nm, respectively. For birch, the linear particle depolarization ratio at 1565 nm was the highest, followed by the 532 and 910 nm wavelengths, respectively. A sharp decrease at 355 nm was evident for birch pollen. For pine pollen, a maximum at the 532 nm wavelength was observed. There was no significant change in the linear particle depolarization ratio at 910 nm for the pollen types considered in this study. Given the low concentration of pollen in the air, the inclusion of the longer wavelengths (910 and 1565 nm) for the detection of birch and pine can be beneficial due to their sensitivity to trace large aerosol particles.</p

A methodology for investigating dust model performance using synergistic EARLINET/AERONET dust concentration retrievals

Systematic measurements of dust concentration profiles at a continental scale were recently made possible by the development of synergistic retrieval algorithms using combined lidar and sun photometer data and the establishment of robust remote-sensing networks in the framework of Aerosols, Clouds, and Trace gases Research InfraStructure Network (ACTRIS)/European Aerosol Research Lidar Network (EARLINET). We present a methodology for using these capabilities as a tool for examining the performance of dust transport models. The methodology includes considerations for the selection of a suitable data set and appropriate metrics for the exploration of the results. The approach is demonstrated for four regional dust transport models (BSC-DREAM8b v2, NMMB/BSC-DUST, DREAMABOL, DREAM8-NMME-MACC) using dust observations performed at 10 ACTRIS/EARLINET stations. The observations, which include coincident multi-wavelength lidar and sun photometer measurements, were processed with the Lidar-Radiometer Inversion Code (LIRIC) to retrieve aerosol concentration profiles. The methodology proposed here shows advantages when compared to traditional evaluation techniques that utilize separately the available measurements such as separating the contribution of dust from other aerosol types on the lidar profiles and avoiding model assumptions related to the conversion of concentration fields to aerosol extinction values. When compared to LIRIC retrievals, the simulated dust vertical structures were found to be in good agreement for all models with correlation values between 0.5 and 0.7 in the 1–6 km range, where most dust is typically observed. The absolute dust concentration was typically underestimated with mean bias values of -40 to -20 μg m−3 at 2 km, the altitude of maximum mean concentration. The reported differences among the models found in this comparison indicate the benefit of the systematic use of the proposed approach in future dust model evaluation studies

A methodology for investigating dust model performance using synergistic EARLINET/AERONET dust concentration retrievals

Systematic measurements of dust concentration profiles at a continental scale were recently made possible by the development of synergistic retrieval algorithms using combined lidar and sun photometer data and the establishment of robust remote-sensing networks in the framework of Aerosols, Clouds, and Trace gases Research Infra-Structure Network (ACTRIS)/European Aerosol Research Lidar Network (EARLINET). We present a methodology for using these capabilities as a tool for examining the performance of dust transport models. The methodology includes considerations for the selection of a suitable data set and appropriate metrics for the exploration of the results. The approach is demonstrated for four regional dust transport models (BSC-DREAM8b v2, NMMB/BSC-DUST, DREAM-ABOL, DREAM8-NMME-MACC) using dust observations performed at 10 ACTRIS/EARLINET stations. The observations, which include coincident multi-wavelength lidar and sun photometer measurements, were processed with the Lidar-Radiometer Inversion Code (LIRIC) to retrieve aerosol concentration profiles. The methodology proposed here shows advantages when compared to traditional evaluation techniques that utilize separately the available measurements such as separating the contribution of dust from other aerosol types on the lidar profiles and avoiding model assumptions related to the conversion of concentration fields to aerosol extinction values. When compared to LIRIC retrievals, the simulated dust vertical structures were found to be in good agreement for all models with correlation values between 0.5 and 0.7 in the 1-6 km range, where most dust is typically observed. The absolute dust concentration was typically underestimated with mean bias values of -40 to -20 mu g m(-3) at 2 km, the altitude of maximum mean concentration. The reported differences among the models found in this comparison indicate the benefit of the systematic use of the proposed approach in future dust model evaluation studies

Aerosol Characterization by Means of Polly(XT) Raman Lidar at South Africa, India and Finland

A multi-wavelength Raman lidar PollyXT supplied by the Finnish Meteorological Institute (FMI) is operated remotely at Finland since November 2012. In addition the lidar performed measurements in two long-term aerosol experimental campaigns, one close to New Delhi in India (March 2008-March 2009) and one at Elandsfontein about 150 km from Johannesburg in South Africa (December 2009-January 2011). The selected Raman lidar data obtained over these periods have been analyzed. We study the seasonal behavior of vertical profiles of intensive and extensive optical properties at the different measurement sites. The range of the retrieved aerosol profiles were studied in terms of different aerosol sources as well as in differences in the mixing state of aerosols

Lidar Observations of Birch and Spruce Pollen in Finland

Pollen has various effects on human health and the environment. To understand phenomena behind atmospheric pollen transport and hence improve pollen forecasts, vertically resolved optical properties and geometrical characteristics of the pollen distribution need to be studied. Lidar measurements and especially the particle depolarization ratio have been found to be an excellent tool to track pollen grains. In this study we present first results of atmospheric pollen characterization based on a 11 days period of birch and spruce pollination events

Are EARLINET and AERONET climatologies consistent? the case of Thessaloniki, Greece

In this study we investigate the climatological behavior of the aerosol optical properties over Thessaloniki during the years 2003–2017. For this purpose, measurements of two independent instruments, a lidar and a sunphotometer, were used. These two instruments represent two individual networks, the European Lidar Aerosol Network (EARLINET) and the Aerosol Robotic Network (AERONET). They include different measurement schedules. Fourteen years of lidar and sunphotometer measurements were analyzed, independently of each other, in order to obtain the annual cycles and trends of various optical and geometrical aerosol properties in the boundary layer, in the free troposphere, and for the whole atmospheric column. The analysis resulted in consistent statistically significant and decreasing trends of aerosol optical depth (AOD) at 355 nm of ĝ'23.2 and ĝ'22.3 % per decade in the study period over Thessaloniki for the EARLINET and the AERONET datasets, respectively. Therefore, the analysis indicates that the EARLINET sampling schedule can be quite effective in producing data that can be applied to long-term climatological studies. It is also shown that the observed decreasing trend is mainly attributed to changes in the aerosol load inside the boundary layer. Seasonal profiles of the most dominant aerosol mixture types observed over Thessaloniki have been generated from the lidar data. The higher values of the vertically resolved extinction coefficient at 355 nm appear in summer, while the lower ones appear in winter. The dust component is more dominant in the free troposphere than in the boundary layer during summer. The biomass burning layers tend to arrive in the free troposphere during spring and summer. This kind of information can be quite useful for applications that require a priori aerosol profiles. For instance, they can be utilized in models that require aerosol climatological data as input, in the development of algorithms for satellite products, and also in passive remote-sensing techniques that require knowledge of the aerosol vertical distribution