Search

The eddy covariance technique (EC) is used worldwide to measure surface fluxes of greenhousegases and energy balance components. Nevertheless, in the scientific community it is well-knownthat EC presents an imbalance of the energy components. In this regard, the atmospheric boundarylayer (ABL) directly influences the mass and energy transfer between surface and atmosphere.Thus, knowledge and characterization of the ABL might be essential to disentangling the driverscausing the imbalance of the energy components measured by EC. This work aims to relate ABLcharacteristics to the accuracy of the surface energy balance closure (SEBC) obtained by EC.The study was carried out in an irrigated olive orchard (Olea europaea var. europaea L.) in theSoutheastern Iberian Peninsula (37.9427º N, 3.3002º W and 370 m asl) during the intensiveBLOOM (turBulence and oLea pOllen prOperties experiMent) campaign, from May 19th to June20th 2022. In order to characterise ABL dynamics, remote sensing techniques are commonly used.One of them is via Doppler lidar, which provides measurements of wind components andturbulence-related products at high spatial and temporal resolutions. In particular, a Doppler lidarStream Line (HALO photonics) with a temporal resolution of 2 seconds and 10 minutes (for verticaland scanning measurements, respectively) and 30 m of vertical spatial resolution was used toretrieve the turbulent kinetic energy dissipation rate (ε) and wind shear (sh) with a commonresolution of 3 minutes as indicators of convective and mechanical sources of turbulence,respectively. To assess the SEBC, we used (1) a three-axis sonic anemometer (CSAT-3, CampbellScientific, Logan, UT, USA) and an enclosed path infrared gas analyser (IRGA, Li-Cor 7200;Lincoln, NE, USA) on a tower 9 m tall to measure latent heat (λE) and sensible heat (H) fluxes; (2)an incoming and outgoing short- and long-wave 4-component radiometer (CNR-4, Kipp & Zonen,Delft, Netherlands) to measure net radiation (Rn); and (3) two each: soil moisture probes (CS616,CSI), thermocouples (TCAV, CSI) and heat flux plates (HFP01, Hukseflux, Delft, the Netherlands)at 0.10 m, 0.04 m and 0.08 m depth, respectively to calculate the soil heat flux (G).Preliminary results show that SEBC enhances under turbulent conditions (slope and R2 from 0.49and 0.92 to 0.81 and 0.94, respectively). However, when the source of turbulence is mechanical theSEBC is less accurate (slope and R2 from 0.82 and 0.94 to 0.77 and 0.9, respectively). A moredetailed study based on principal component analysis of the ABL height, skewness of the verticalwind velocity and vertical profile of horizontal wind, together with ε and sh, among other Dopplerlidar-derived products is expected to offer reliable information that is highly relevant regarding theinfluence of the ABL on the SEBC.ACKNOWLEDGMENTS: This work was supported by theSpanish Ministry of Science and Innovation through project PID2020-117825GB-C21(INTEGRATYON3), the Andalusian regional Development through the projects B-RNM-60-UGR20 (OLEAGEIs) and P18-RT-3629 (ICAERSA), including European Union ERDF funds.

A statistical validation campaign of Aeolus L2B wind products has been performed with a ground-based Doppler lidar system and radiosondes at the ACTRIS/AGORA facility in Granada (Spain). The validation activities with the automatic ground-based lidar system lasted from the release of the reprocessed Aeolus data to the orbit shift of June 2021. This validation was performed using 109 B10 and B11 Aeolus products (within 100 km around the station) and 30-min averages of coincident ground-based lidar measurements (mean bin distance of ~50 km). The comparison yielded an approximately equal overestimation and underestimation of Aeolus HLOS wind speed during that period for the Rayleigh clear and Mie cloudy configurations. However, the reliability of the results was constrained to Aeolus’ lowermost bins (roughly up to 3.5 km asl), due to the limited vertical coverage of the ground-based lidar measurements. Several spin-off analyses were performed varying the maximum distance to consider an Aeolus bin into the comparison and the average of the ground-based lidar measurements, in order to confirm the reliability of the criteria considered. An additional study was performed with Aeolus products after the orbit shift (B12 with a mean bin distance of ~75 km). A set of 7 radiosondes were launched between June 2021 and July 2022 (B12, mean bin distance of ~75 km) aiming to increase their coincidence in space and time with the satellite overpass (~30 minutes before). The radiosondes could provide full vertical coverage of Aeolus profiles (from surface up to ~20 km asl) and the comparison yielded that the Rayleigh clear HLOS wind speed presented an approximately equal overestimation and underestimation, while the Mie cloudy HLOS wind speed was significantly overestimated. Spin-off analyses were performed in order to test how the spatiotemporal colocation of the radiosonde affected the results.

We use Doppler lidar wind profiles from six locations around the globe to evaluate the wind profile forecasts in the boundary layer generated by the operational global Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts (ECMWF). The six locations selected cover a variety of surfaces with different characteristics (rural, marine, mountainous urban, and coastal urban). We first validated the Doppler lidar observations at four locations by comparison with co-located radiosonde profiles to ensure that the Doppler lidar observations were of sufficient quality. The two observation types agree well, with the mean absolute error (MAE) in wind speed almost always less than 1 m s−1. Large deviations in the wind direction were usually only seen for low wind speeds and are due to the wind direction uncertainty increasing rapidly as the wind speed tends to zero. Time–height composites of the wind evaluation with 1 h resolution were generated, and evaluation of the model winds showed that the IFS model performs best over marine and coastal locations, where the mean absolute wind vector error was usually less than 3 m s−1 at all heights within the boundary layer. Larger errors were seen in locations where the surface was more complex, especially in the wind direction. For example, in Granada, which is near a high mountain range, the IFS model failed to capture a commonly occurring mountain breeze, which is highly dependent on the sub-grid-size terrain features that are not resolved by the model. The uncertainty in the wind forecasts increased with forecast lead time, but no increase in the bias was seen. At one location, we conditionally performed the wind evaluation based on the presence or absence of a low-level jet diagnosed from the Doppler lidar observations. The model was able to reproduce the presence of the low-level jet, but the wind speed maximum was about 2 m s−1 lower than observed. This is attributed to the effective vertical resolution of the model being too coarse to create the strong gradients in wind speed observed. Our results show that Doppler lidar is a suitable instrument for evaluating the boundary layer wind profiles in atmospheric models.

The Atmospheric Boundary Layer (ABL) over two middle-latitude rural sites was characterized in terms of mean horizontal wind and turbulence sources using a standard classification methodology based on Doppler lidar. The first location was an irrigated olive orchard in Úbeda (Southern Spain), representing one of the most important crops in the Mediterranean basin and a typical site with Mediterranean climate. The second location was PolWET peatland site in Rzecin (Northwestern Poland), representing one of the largest natural terrestrial carbon storages that have a strong interaction with the climate system. The results showed typical situations for non cloud-topped ABL cases, where ABL is fully developed during daytime due to convection, with high turbulent activity and strong positive skewness indicating frequent and powerful updrafts. The cloud-topped cases showed the strong influence that clouds can have on ABL development, preventing it to reach the same maximum height and introducing top-down movements as an important contribution to mixing. The statistical analysis of turbulent sources allowed for finding a common diurnal cycle for convective mixing at both sites, but nocturnal wind shear driven turbulence with marked differences in its vertical distribution. This analysis demonstrates the Doppler lidar measurements and the classification algorithm strong potential to characterize the dynamics of ABL in its full extent and with high temporal resolution. Moreover, some recommendations for future improvement of the classification algorithm were provided on the basis of the experience gained., Fundacion Ramon Areces, European Space Agency 4000119961/16/NL/FF/mg, Polish National Science Centre (NCN) 2021/40/C/ST10/00023, Spanish Government CGL2015-73250-JIN CGL201681092-R CGL2017-83538-C3-1-R CGL2017-90884-REDT PID2020117825GB-C21 PID2020-120015RB-100, Andalusian Regional Government P18-RT-3820, FEDER-UGR program ARNM-430-UGR20, University of Granada, ACTRIS-2 Research Infrastructure Project of the European Union's Horizon 2020 research and innovation program 654109, European Cooperation in Science and Technology (COST) CA18235, Universidad de Granada/CBUA

This work focuses on the analysis of the seasonal cycle of temperature and relative humidity (RH) profiles and integrated water vapor (IWV) obtained from microwave radiometer (MWR) measurements over the mid-latitude city of Granada, southern Spain. For completeness the study, the maximum atmospheric boundary layer height (ABLHmax) is also included. To this end, we have firstly characterized the HATPRO-RPG MWR errors using 55 co-located radiosondes (RS) by means of the mean-bias (biasbar) profile and the standard deviation (SDbias) profile classified under all-weather conditions and cloud-free conditions. This characterization pointed out that temperature from HATPRO-MWR presents a very low biasbar respects RS mostly below 2.0 km agl, ranging from positive to negative values under all-weather conditions (from 1.7 to -0.4 K with SDbias up to 3.0 K). Under cloud-free conditions, the bias was very similar to that found under all-weather conditions (1.8 to -0.4 K) but with smaller SDbias (up to 1.1 K). The same behavior is also seen in this lower part (ground to 2.0 km agl) for RH. Under all-weather conditions, the mean RH bias ranged from 3.0 to -4.0% with SDbias between 10 and 16.3% while under cloud-free conditions the bias ranged from 2.0 to -0.4% with SDbias from 0.5 to 13.3%. Above 2.0 km agl, the SDbias error increases considerably up to 4 km agl (up to -20%), and then decreases slightly above 7.0 km agl (up to -5%). In addition, IWV values from MWR were also compared with the values obtained from the integration of RS profiles, showing a better linear fit under cloud-free conditions (R2 = 0.96) than under all-weather conditions (R2 = 0.82). The mean bias under cloud-free conditions was -0.80 kg/m2 while for all-weather conditions it was -1.25 kg/m2. Thus, the SDbiasfor all the statistics (temperature, RH and IWV) of the comparison between MWR and RS presented higher values for all-weather conditions than for cloud-free conditions ones. It points out that the presence of clouds is a key factor to take into account when MWR products are used. The second part of this work is devoted to a seasonal variability analysis over five years, leading us to characterize thermodynamically the troposphere over our site. This city atmosphere presents a clear seasonal cycle where temperature, ABLHmax and IWV increase from winter to summer and decrease in autumn, meanwhile RH decreases along the warmer seasons. This city presents cold winters (mean daily maximum temperature: 10.6 ± 1.1 °C) and dry/hot summers (mean daily maximum temperature of 28.8 ± 0.9 °C and mean daily maximum of surface RH up to 55.0 ± 6.0%) at surface (680 m asl). Moreover, considering temporal trends, our study pointed out that only temperature and RH showed a linear increase in winters with a mean-rate of (0.5 ± 0.1) °C/year and (3.4 ± 1.7) %/year, respectively, from ground to 2.0 km agl, meanwhile IWV presented a linear increase of 1.0 kg·m-2/year in winters, 0.78 kg·m-2/year in summers and a linear decrease in autumns of -0.75 kg·m-2/year., Andalusia Regional Government through project P12-RNM-2409, Spanish Ministry of Economy and Competitiveness through projects CGL2013-45410-R, CGL2015-73250-JIN and CGL2016-81092-R, Juan de la Cierva grant IJCI-2016-30007

The Sierra Nevada Lidar aerOsol Profiling Experiment I and II (SLOPE I and II) campaigns were intended to determine the vertical structure of the aerosol by remote sensing instruments and test the various retrieval schemes for obtaining aerosol microphysical and optical properties with in-situ measurements. These campaigns deployed a set of in-situ and remote sensing instruments at the stations include in AGORA observatory (Andalusian Global ObseRvatory of the Atmosphere) in the Granada area (Spain) along summer in 2016 and 2017. In this work, using the in-situ measurements performed at a high-altitude station, Sierra Nevada station, and airborne flights, we evaluate the retrievals of aerosol properties by GRASP code (Generalized Retrieval of Atmosphere and Surface Properties) combining lidar and sun-sky photometer measurements. Besides, we show an overview of aerosol properties retrieved by GRASP during SLOPE I and II campaigns. We evaluate the GRASP retrievals of total aerosol volume concentration (discerning between fine and coarse modes), extinction and scattering coefficients, and for the first time we present an evaluation of absorption coefficient. The statistical analysis of the aerosol optical and microphysical properties, both column-integrated and vertically-resolved, from May to July 2016 and 2017 shows a large variability in aerosol load and types. The results show a strong predominance of desert dust particles due to the North African intrusions. The vertically-resolved analysis denotes a decay of the atmospheric aerosols with altitude up to 5 km a.s.l. Finally, two events of desert dust and biomass burning were used to show the high potential of GRASP to retrieve and study the aerosol properties profiles such as absorption coefficient and single scattering albedo for different aerosol types. The aerosol properties retrieved by GRASP show good agreement with simultaneous in situ measurements performed at Sierra Nevada Station (SNS) in Granada. In general, GRASP overestimates the in situ data at SNS with a mean difference lower than 6 µm3/cm3 for volume concentration, 11 Mm−1 and 2 Mm−1 for scattering and absorption coefficient. On the other hand, the comparison of GRASP with airborne measurements also shows an overestimation with mean absolute differences of 14 ± 10 Mm−1 and 1.2 ± 1.2 Mm−1 for scattering and absorption coefficients, showing a better agreement for absorption (scattering) coefficient with higher (lower) aerosol optical depth. The potentiality of GRASP showed in this study will contribute to enhancing the representativeness of the aerosol vertical distribution and provide information for satellite and global model evaluation.

Remote sensing measurements of aerosols using depolarization Raman Lidar systems from 4 EARLINET (European Aerosol research Lidar Network) stations are used for a comprehensive analysis of Saharan dust events over the Mediterranean basin in the period 2014–2017. In this period, we selected to study 51 dust events regarding the geometrical, optical and microphysical properties of dust particles, classifying them and assessing their radiative forcing effect on the atmosphere. From West to East, the stations of Granada, Potenza, Athens and Limassol were selected as representative Mediterranean cities regularly affected by Saharan dust intrusions. Emphasis was given on lidar measurements in the visible (532 nm) and specifically on the consistency of the particle linear depolarization ratio (δp532), the extinction-to-backscatter lidar ratio (LR532) and the Aerosol Optical Thickness (AOT532) within the observed dust layers. We found mean δp532 values of 0.24 ± 0.05, 0.26 ± 0.06, 0.28 ± 0.05 and 0.28 ± 0.04, mean LR532 values of 52 ± 8 sr, 51 ± 9 sr, 52 ± 9 sr and 49 ± 6 sr, and mean AOT532 values of 0.40 ± 0.31, 0.11 ± 0.07, 0.12 ± 0.10 and 0.32 ± 0.17, for Granada, Potenza, Athens and Limassol, respectively. The mean layer thickness values were found to range from ~1700 to ~3400 m. Additionally, based also on a previous aerosol type classification scheme provided by airborne High Spectral Resolution Lidar (HSRL) observations and on air mass backward trajectory analysis, a clustering analysis was performed in order to identify the major mixing aerosol types over the studied area. Furthermore, a synergy of lidar measurements and modeling was used to deeply analyze the solar and thermal radiative forcing of airborne dust. In total, a cooling behavior in the solar range and a significantly lower heating behavior in the thermal range was estimated. Depending on the dust optical and geometrical properties, the load intensity and the solar zenith angle (SZA), the estimated solar radiative forcing values range from −59 to −22 W m−2 at the surface and from −24 to −1 W m−2 at the top of the atmosphere (TOA). Similarly, in the thermal spectral range these values range from +2 to +4 W m−2 for the surface and from +1 to +3 W m−2 for the TOA. Finally, the radiative forcing seems to be inversely proportional to the dust mixing ratio, since higher absolute values are estimated for less mixed dust layers.

Vibrational Raman effect is widely used in atmospheric lidar systems, but rotational Raman present several advantages. We have implemented a new setup in the ultraviolet branch of an existing multiwavelength lidar system to collect signal from rotational Raman lines of Oxygen and Nitrogen. We showed that, with an appropriate filter wavelength selection, the systematic error introduced in the particle optical properties due to temperature dependence was less than 4%. With this new setup, we have been able to retrieve aerosol extinction and backscatter coefficients profiles at 355 nm with 1-h time resolution during daytime and up to 1-min time resolution during nighttime., Grupo de Física de la Atmósfera (RNM119)

The financial support by the European Union's Horizon 2020 research and innovation program (ACTRIS-2, grant agreement no. 654109) is gratefully acknowledged. The investigation was supported by the Belarusian Republican Foundation for Fundamental, Research Agreement No. F18EA-011, by the Russian Foundation for Basic Research (grant No. 18-55-81001), and also by Vietnam Academy of Science and Technology (project code QTRU05.01/18-20). The Portuguese team acknowledges the support from the Portuguese Science Foundation, in the frame of the European Regional Development Fund -COMPETE 2020, under the projects UID/GEO/04683/2013 (POCI-01-0145-FEDER-007690), The paper presents the preliminary results of the lidar&radiometer measurement campaign (LRMC2017), estimation of statistical relations between aerosol mode concentrations retrieved from CALIOP and ground-based lidar stations and case study of fire smoke events in the Eurasian regions using combined ground-based and space lidar and radiometer observations., European Union's Horizon 2020 research and innovation program (ACTRIS-2) 654109, Belarusian Republican Foundation for Fundamental, Research F18EA-011, Russian Foundation for Basic Research (RFBR) 18-55-81001, Vietnam Academy of Science and Technology QTRU05.01/18-20, Portuguese Science Foundation, in the frame of the European Regional Development Fund -COMPETE 2020 UID/GEO/04683/2013 (POCI-01-0145-FEDER-007690)

The Planetary Boundary Layer (PBL) is a relevant part of the atmosphere with a variable extension that clearly plays an important role in fields like air quality or weather forecasting. Passive and active remote sensing systems have been widely applied to analyze PBL characteristics. The combination of different remote sensing techniques allows obtaining a complete picture on the PBL dynamic. In this study, we analyze the PBL using microwave radiometer, elastic lidar and Doppler lidar data. We use co-located data simultaneously gathered in the framework of SLOPE-I (Sierra Nevada Lidar aerOsol Profiling Experiment) campaign at Granada (Spain) during a 90- day period in summer 2016. Firstly, the PBL height (PBLH) obtained from microwave radiometer data is validated against PBLH provided by analyzing co-located radiosondes, showing a good agreement. In a second stage, active remote sensing systems are used for deriving the PBLH. Thus, an extended Kalman filter method is applied to data obtained by the elastic lidar while the vertical wind speed variance method is applied to the Doppler lidar. PBLH′s derived by these approaches are compared to PBLH retrieved by the microwave radiometer. The results show a good agreement among these retrievals based on active remote sensing in most of the cases, although some discrepancies appear in instances of intense PBL changes (either growth and/or decrease).

This study evaluates the potential of GRASP algorithm (Generalized Retrieval of Aerosol and Surface Properties) to retrieve continuous day-to-night aerosol properties, both column-integrated and vertically-resolved. The study is focused on the evaluation of GRASP retrievals during an intense Saharan dust event that occurred during the Sierra Nevada Lidar aerOsol Profiling Experiment I (SLOPE I) field campaign. For daytime aerosol retrievals, we combined the measurements of the lidar ground-based from EARLINET (European Aerosol Research Lidar Network) station and sun/sky photometer from AERONET (Aerosol Robotic Network), both instruments co-located in Granada (Spain). However, for night-time retrievals three different combinations of active and passive remote sensing measurements are proposed. The first scheme (N0) uses lidar night-time measurements in combination with the interpolation of sun/sky daytime measurements. The other two schemes combine lidar night-time measurements with night-time aerosol optical depth obtained by lunar photometry either using intensive properties of the aerosol retrieved during sun/sky daytime measurements (N1) or using the moon aureole radiance obtained by sky camera images (N2). Evaluations of the columnar aerosol properties retrieved by GRASP are done versus standard AERONET retrievals. The coherence of day-to-night evolutions of the different aerosol properties retrieved by GRASP is also studied. The extinction coefficient vertical profiles retrieved by GRASP are compared with the profiles calculated by Raman technique at night-time with differences below 30% for all schemes at 355, 532 and 1064 nm. Finally, the volume concentration and scattering coefficient retrieved by GRASP at 2500 m a.s.l. are evaluated by in-situ measurements at this height at Sierra Nevada Station. The differences between GRASP and in-situ measurements are similar for the different schemes, with differences below 30 % for both volume concentration and scattering coefficient. In general, for the scattering coefficient, the GRASP N0 and N1 show better results than the GRASP N2 schemes, while for volume concentration, GRASP N2 shows the lowest differences against in-situ measurements, around 10 %, for high the aerosol optical depth values.

In this paper we extract the aerosol microphysical properties for a collection of mineral dust cases measured by multi-wavelength depolarization Raman lidar systems located at the National Technical University of Athens (NTUA, Athens, Greece) and the Andalusian Institute for Earth System Research (IISTA-CEAMA, Granada, Spain). The lidar-based retrievals were carried out with the Spheroidal Inversion eXperiments software tool (SphInX) developed at the University of Potsdam (Germany). The software uses regularized inversion of a two-dimensional enhancement of the Mie model based on the spheroid-particle approximation with the aspect ratio determining the particle shape. The selection of the cases was based on the transport time from the source regions to the measuring sites. The aerosol optical depth as measured by AERONET ranged from 0.27 to 0.54 (at 500 nm) depending on the intensity of each event. Our analysis showed the hourly mean particle linear depolarization ratio and particle lidar ratio values at 532 nm ranging from 11 to 34% and from 42 to 79 sr respectively, depending on the mixing status, the corresponding air mass pathways and their transport time. Cases with shorter transport time showed good agreement in terms of the optical and SphInX-retrieved microphysical properties between Athens and Granada providing a complex refractive index value equal to 1.4 + 0.004i. On the other hand, the results for cases with higher transport time deviated from the aforementioned ones as well as from each other, providing, in particular, an imaginary part of the refractive index ranging from 0.002 to 0.005. Reconstructions of two-dimensional shape-size distributions for each selected layer showed that the dominant effective particle shape was prolate with diverse spherical contributions. The retrieved volume concentrations reflect overall the intensity of the episodes., Spanish Ministry of Sciences, Innovation and Universities through project CGL2016-81092, Spanish Ministry of Education, Culture and Sports through grant FPU14/03684

The research leading to these results has received funding from the H2020 program from the European Union (grant agreement no. 654109, 778349) and also from the Spanish Ministry of Industry, Economy and Competitiviness (MINECO, ref. CGL2013-45410-R, CGL2016-81092-R, CGL2017-85344-R, TEC2015-63832-P), the Spanish Ministry of Science, Innovation and Universities (ref. CGL2017-90884-REDT); the CommSensLab "Maria de Maeztu" Unity of Excellence (ref. MDM-2016-0600) financed by the Spanish Agencia Estatal de Investigación. Co-funding was also provided by the European Union through the European Regional Development Fund (ref. POCI-01-0145-FEDER-007690, ALT20-03-0145-FEDER-000004, ALT20-03-0145-FEDER-000011); by the Andalusia Regional Government (ref. P12-RNM-2409); by the Madrid Regional Government (projects TIGAS-CM, ref. Y2018/EMT-5177 and AIRTEC-CM, ref. P2018/EMT4329); by the University of Granada through “Plan Propio. Programa 9 Convocatoria 2013” and by the Portuguese Foundation for Science and Technology and national funding (ref. SFRH/BSAB/143164/2019). The BSC-DREAM8b and NNMB/BSC-Dust (now NMMB-MONARCH) model simulations were performed by the Mare Nostrum supercomputer hosted by the Barcelona Supercomputer Center (BSC). S. Basart acknowledges the AXA Research Fund for supporting aerosol research at the BSC through the AXA Chair on Sand and Dust Storms Fund, as well as the InDust project (COST Action CA16202). The authors gratefully acknowledge the NOAA Air Resources Laboratory (ARL) for the provision of the HYSPLIT transport and dispersion model and/or READY website (http://www.ready.noaa.gov) used in this publication., An unprecedented extreme Saharan dust event was registered in winter time from 20 to 23 February 2017 over the Iberian Peninsula (IP). We report on aerosol optical properties observed under this extreme dust intrusion through passive and active remote sensing techniques. For that, AERONET (AErosol RObotic NETwork) and EARLINET (European Aerosol Research LIdar NETwork) databases are used. The sites considered are: Barcelona (41.38°N, 2.17°E), Burjassot (39.51°N, 0.42°W), Cabo da Roca (38.78°N, 9.50°W), Évora (38.57°N, 7.91°W), Granada (37.16°N, 3.61°W) and Madrid (40.45°N, 3.72°W). Large aerosol optical depths (AOD) and low Ångström exponents (AE) are observed. An AOD of 2.0 at 675 nm is reached in several stations. A maximum peak of 2.5 is registered in Évora. During and around the peak of AOD, AEs close to 0 and even slightly negative are measured. With regard to vertically-resolved aerosol optical properties, particle backscatter coefficients as high as 15 Mm−1 sr−1 at 355 nm are recorded at the lidar stations. Layer-mean lidar ratios are found in the range 40–55 sr at 355 nm and 34–61 sr at 532 nm during the event. The particle depolarization ratios are found to be constant inside the dust layer, and consistent from one site to another. Layer-mean values vary in the range 0.19–0.31. Another remarkable aspect of the event is the limited vertical distribution of the dust plume which never exceeds 5 km. The extreme aspect of the event also presented a nice case for testing the ability of two dust forecast models, BSC-DREAM8b and NMMB/BSC-Dust, to reproduce the arrival, the vertical distribution and the intensity of the dust plume over a long-range transport region. In the particular case of the February 2017 dust event, we found a large underestimation in the forecast of the extinction coefficient provided by BSC-DREAM8b at all heights independently of the site. In contrast NMMB/BSC-Dust forecasts presented a better agreement with the observations, especially in southwestern part of the IP. With regard to the forecast skill as a function of lead time, no clear degradation of the prognostic is appreciated at 24, 48 and 72 h for Évora and Granada stations (South). However the prognostic does degrade (bias increases and/or correlation decreases) for Barcelona (North), which is attributed to the fact that Barcelona is at a greater distance from the source region and to the singularity of the event., Funding from the H2020 program from the European Union (grant agreement no. 654109, 778349), Spanish Ministry of Industry, Economy and Competitiviness (MINECO, ref. CGL2013-45410-R, CGL2016-81092-R, CGL2017-85344-R, TEC2015-63832-P), Spanish Ministry of Science, Innovation and Universities (ref. CGL2017-90884-REDT), CommSensLab "Maria de Maeztu" Unity of Excellence (ref. MDM-2016-0600) financed by the Spanish Agencia Estatal de Investigación, Co-funding was also provided by the European Union through the European Regional Development Fund (ref. POCI-01-0145-FEDER-007690, ALT20-03-0145-FEDER-000004, ALT20-03-0145-FEDER-000011), Andalusia Regional Government (ref. P12-RNM-2409); by the Madrid Regional Government (projects TIGAS-CM, ref. Y2018/EMT-5177 and AIRTEC-CM, ref. P2018/EMT4329), Portuguese Foundation for Science and Technology and national funding (ref. SFRH/BSAB/143164/2019)

The AEAOD– ΔAEAOD grid proposed by Gobbi et al. (2007) is a graphical method used to visually represent the spectral characterization of aerosol optical depth (AOD), i.e. Angström exponent (AE) and its curvature, in order to infer the fine mode contribution (η) to the total AOD and the size of the fine mode aerosol particles. Perrone et al. (2014) applied this method for the wavelengths widely used in lidar measurements. However, in neither case does the method allow for a direct relationship between η and the fine mode fraction contribution to the total aerosol population. Some discussions are made regarding the effect of shape and composition to the classical AE-ΔAE plot. The potential use of particle backscatter measurements, widely used in aerosol characterization methods together with extinction measurements, is also discussed in the AE-ΔAE grid context. A modification is proposed that yields the submicron contribution to the total volume concentration by using particle extinction data, and a comparison to experimental measurements is made. Our results indicate that the use of a modified AE-ΔAE grid plot to directly obtain submicrometric and micrometric mode fraction to the total aerosol population is feasible if a volume-based bimodal particle size distribution is used instead of a number-based one., Andalusia Regional Government through project P12-RNM-2409, Spanish Ministry of Sciences, Innovation and Universities (CGL2016-81092 and CGL2017 -90884 - REDT)

Four years (2014-2017) of observations of depolarization Raman Lidar systems of four EARLINET (European Aerosol research Lidar Network) stations [from West to East: Granada (Spain), Potenza (Italy), Athens (Greece) and Limassol (Cyprus)] were collected and used to a statistical analysis of Saharan dust events over Mediterranean basin. In this study, emphasis is given to the consistency of the particle linear depolarization ratio (δp532), the extinction-to-backscatter ratio mentioned as Lidar Ratio (LR532) and the Aerosol Optical Thickness (AOT532) within the observed Saharan dust layers, corresponding to the visible range (532 nm). Geometrical properties and clusters of aerosol mixtures are also presented. Our clustering was based on previous classification by airborne High Spectral Resolution Lidar (HSRL) observations and was further supported by backward trajectory analysis. We found mean δp532 values of 0.24±0.05, 0.26±0.06, 0.28±0.05 and 0.28±0.04, mean LR532 values of 52±8 sr, 51±9 sr, 52±9 sr and 49±6 sr, mean AOT532 values of 0.40±0.31, 0.11±0.07, 0.12±0.10 and 0.32±0.17 and mean layer thicknesses of 3392±1458 m, 2150±1082 m, 1872±816 m and 1716±567 m for Granada, Potenza, Athens and Limassol respectively. This work could assist in bridging the existing gaps related to the extensive and intensive dust aerosol properties over the Mediterranean and enriching the bibliography about mixed aerosol layers from different sources (e.g. dust and biomass burning (BB) aerosols, dust and urban/ industrial aerosols).

An unprecedented extreme Saharan dust event was registered in winter time from 20 to 23 February 2017 over the Iberian Peninsula (IP). We report on aerosol optical properties observed under this extreme dust outbreak through remote sensing (active and passive) techniques. For that, EARLINET (European Aerosol Research LIdar NETwork) lidar and AERONET (AErosol RObotic NETwork) Sun-photometer Cimel CE 318 measurements are used. The sites considered are: Barcelona (41.38º N, 2.17º E), Burjassot (39.51º N, 0.42º W), Cabo da Roca (38.78º N, 9.50º W), Évora (38.57º N, 7.91º W), Granada (37.16º N, 3.61º W) and Madrid (40.45º N, 3.72º W). In general, large aerosol optical depths (AOD) and low Ångström exponents (AE) are observed. An AOD of 2.0 at 675 nm is reached in several stations. Maximum values of AOD675 of 2.5 are registered in Évora. During and around the peak of AOD675, AEs close to 0 are measured. With regard to vertically-resolved aerosol optical properties, particle backscatter coefficients as high as 1.5∙10−5 m−1 sr−1 at 355 nm are recorded at every lidar stations. Mean lidar ratios are found in the range 40–55 sr at 355 nm and 34–61 sr at 532 nm during the event inside the dust layer. Mean particle and volume depolarization ratios are found to be very consistent between lidar stations. They range 0.19–0.31 and 0.12–0.26 respectively. The optical properties are also found very stable with height in the dust layer. Another remarkable aspect of the event is the limited height of the dust transport which is found between the ground and 5 km. Our vertically-resolved aerosol properties are also used to estimate the performances of two dust models, namely BSC-DREAM8b and NMMB/BSC-Dust, in order to evaluate their forecast skills in such intense dust outbreaks. We found that forecasts provided by the NMMB/BSC-Dust show a better agreement with observations than the ones from BSC-DREAM8b. The BSC-DREAM8b forecasts (24 h) present a large underestimation during the event. No clear degradation of the prognostics is appreciated in 24, 48, 72 h except for the Barcelona station.

The Planetary Boundary Layer (PBL) is the lowermost region of troposphere and endowed with turbulent characteristics, which can have mechanical or thermodynamic origins. Such behavior gives to this layer great importance, mainly in studies about pollutant dispersion and weather forecasting. However, the instruments usually applied in studies about turbulence in the PBL have limitations in spatial resolution (anemometer towers) or temporal resolution (aircrafts). In this study we propose the synergetic use of remote sensing systems (microwave radiometer [MWR], Doppler lidar [DL] and elastic lidar [EL]) to analyze the PBL behavior. Furthermore, we show how some meteorological variables such as air temperature, aerosol number density, vertical wind, relative humidity and net radiation might influence the PBL dynamic. The statistical moments of the high frequency distributions of the vertical velocity, derived from DL and of the backscattered coefficient derived from EL, are corrected by two methodologies, namely first lag and −2/3 correction. The corrected profiles present small differences when compare against the uncorrected profiles, showing low influence of noise and the viability of the proposed methodology. Two case studies were analyzed in detail, one corresponding to a well-defined PBL and another one corresponding to a situation with presence of a Saharan dust lofted aerosol layer and clouds. In both cases the results provided by the different instruments are complementary, thus the synergistic use of the different systems allow us performing a detailed monitoring of the PBL.

Two cases from SLOPE campaign at Granada are analyzed in terms of particle microphysical properties using novel software developed at Potsdam University. Multiwavelength Raman lidar measurements of particle extinction and backscatter coefficients as well as linear particle depolarization ratios are used as input for the software. The result of the retrieval is a 2-dimensional particle volume distribution as a function of radius and aspect ratio, from which the particle microphysical properties are obtained.

This work focuses on the study of aerosol hygroscopic growth during Sierra Nevada Lidar AerOsol Profiling Experiment (SLOPE I) campaign by using the synergy of active and passive remote sensors at Granada valley station (IISTA-CEAMA station) and in-situ instrumentation at a mountain station (Sierra Nevada station, SNS), located ~ 20 km away from IISTA-CEAMA and about 1 km above this station. To this end, a methodology based on the combination of calibrated water vapour mixing ratio (r) profiles, retrieved from an EARLINET multiwavelength Raman lidar (RL), and continuous temperature profiles from a microwave radiometer (MWR) for obtaining relative humidity (RH) profiles with high temporal resolution is used. This methodology is validated against an approach using radiosounding (RS) data at the EARLINET IISTA-CEAMA station, achieving differences in hygroscopic growth parameter (γ) lower than 5 % between the methodology based on RS and that based on remote sensing. During SLOPE I the remote sensing methodology used for aerosol hygroscopic growth studies has been checked using Mie calculations with in-situ measurements at SNS. For SLOPE I case with remote sensing instrumentation, an increase in particle backscatter coefficient at 355 and 532 nm is observed from 1.5 to 2.4 km a.s.l. in the relative humidity range of 78–98 %, but also a decrease on Ångström exponent (AE) and particle linear depolarization ratio (PLDR). This fact indicates that particles become larger and more spherical as relative humidity increases, what is strongly related to aerosol hygroscopic growth. In addition, atmospheric stability is checked ensuring well-mixed layers. Vertical and horizontal wind analysis is performed by means of a co-located Doppler lidar system at IISTA-CEAMA station, in order to evaluate the horizontal and vertical dynamics of the air masses. Finally, the Hänel parameterization is applied to experimental data for both stations, obtaining relative differences between RL and Mie simulations up to 13 % and 10 % for 532 nm and 355 nm, respectively.

Strong events of long-range transported biomass burning aerosol were detected during July 2013 at three EARLINET (European Aerosol Research Lidar Network) stations, namely Granada (Spain), Leipzig (Germany) and Warsaw (Poland). Satellite observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and CALIOP (Cloud-Aerosol Lidar with Orthogonal Polarization) instruments, as well as modeling tools such as HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) and NAAPS (Navy Aerosol Analysis and Prediction System), have been used to estimate the sources and transport paths of those North American forest fire smoke particles. A multiwavelength Raman lidar technique was applied to obtain vertically resolved particle optical properties, and further inversion of those properties with a regularization algorithm allowed for retrieving microphysical information on the studied particles. The results highlight the presence of smoke layers of 1–2 km thickness, located at about 5 km a.s.l. altitude over Granada and Leipzig and around 2.5 km a.s.l. at Warsaw. These layers were intense, as they accounted for more than 30 % of the total AOD (aerosol optical depth) in all cases, and presented optical and microphysical features typical for different aging degrees: color ratio of lidar ratios (LR532 ∕ LR355) around 2, α-related ångström exponents of less than 1, effective radii of 0.3 µm and large values of single scattering albedos (SSA), nearly spectrally independent. The intensive microphysical properties were compared with columnar retrievals form co-located AERONET (Aerosol Robotic Network) stations. The intensity of the layers was also characterized in terms of particle volume concentration, and then an experimental relationship between this magnitude and the particle extinction coefficient was established.

espanolEste trabajo presenta la caracterizacion, durante un evento especial de transporte a gran distancia, de particulas de combustion de biomasa, procedentes de incendios forestales en Canada. El evento se detecto entre el 11 y el 17 de julio de 2013 en la estacion de EARLINET en Granada (37.16° N, 3.61° W, 680 snm), asi como en otras estaciones de EARLINET. Para su caracterizacion se han usado diferentes tecnicas, incluyendo activas/pasivas y terrestres/por satelite. En particular, el estudio implico el procesado de senales lidar tanto elasticas como Raman, y datos de fotometro solar de la red AERONET. Ademas se han usado modelos como HYSPLIT, NAAPS y BSC-DREAM8b para investigar el comportamiento de las plumas de aerosol a escala sinoptica. Los resultados destacan la presencia de una estructura con multiples capas que muestran un exponente de Angstrom de coeficiente de retrodispersion alto (aproximadamente 2) y una despolarizacion lineal de particulas muy baja (menos de 0.04) en las capas a 5-6 km (snm) de altura, lo que indica que son particulas considerablemente pequenas y esfericas. Durante el periodo mas intenso, el 14 de julio, se detectaron valores del coeficiente de retrodispersion de hasta 2.9×10ˉ⁶ (m sr) ˉ¹ para 532 nm. EnglishThis work presents the characterization of aerosols during a special episode of long-range transported biomass burning particles from Canadian forest fires. This event was detected at the EARLINET Granada station (37.16° N, 3.61° W, 680 m asl) from 11ᵗʰ to 17ᵗʰ July 2013, and also over other EARLINET stations. Different techniques, covering active/passive and ground based/satellite ones, have been used. In particular, the study has involved processing of both Raman and elastic lidar signals and sun-photometer AERONET data. Furthermore, modeling tools as HYSPLIT, NAAPS and BSC-DREAM8b have been used to understand the synoptic behavior of the aerosol plumes. The results highlight the presence of a multilayered structure showing large backscatter-related Angstrom Exponent (around 2) and very low particle linear depolarization (less than 0.04) at layers around 5-6 km (asl) which means considerably small and spherical particles. During the strongest period, on 14ᵗʰ July, aerosol backscatter values up to 2.9×10ˉ⁶ (m sr) ˉ¹ at 532 nm were detected in this layer.

The Planetary Boundary Layer (PBL) is the lowermost part of the troposphere. In this work, we analysed some high order moments and PBL height detected continuously by three remote sensing systems: an elastic lidar, a Doppler lidar and a passive Microwave Radiometer, during the SLOPE-2016 campaign, which was held in Granada from May to August 2016. This study confirms the feasibility of these systems for the characterization of the PBL, helping us to justify and understand its behaviour along the day.

An intense biomass-burning (BB) event from North America in July 2015 was observed over Ny-Ålesund (Spitsbergen, European Arctic). An extreme air pollution took place and aerosol optical depth (AOD) of more than 1 at 500nm occurs in middle and lower troposphere. We analyse data from the multi-wavelength Raman-lidar KARL of Alfred Wegener Institute to derive microphysical properties of the aerosol of one interesting layer from 3186 to 3306 m via regularization. We found credible and confidential microphysical parameters.