Sun-sky Photometer Research Articles

Evaluation of on-site calibration procedures for SKYNET Prede POM sun–sky photometers

Abstract. To retrieve columnar intensive aerosol properties from sun–sky photometers, both irradiance and radiance calibration factors are needed. For the irradiance the solar calibration constant, V0, which denotes the instrument counts for a direct normal solar flux extrapolated to the top of the atmosphere, must be determined. The solid view angle, SVA, is a measure of the field of view of the instrument, and it is important for obtaining the radiance from sky diffuse irradiance measurements. Each of the three sun-photometer networks considered in the present study (SKYNET, AERONET, WMO GAW) adopts different protocols of calibration, and we evaluate the performance of the on-site calibration procedures, applicable to every kind of sun–sky photometer but tested in this analysis only on SKYNET Prede POM01 instruments, during intercomparison campaigns and laboratory calibrations held in the framework of the Metrology for Aerosol Optical Properties (MAPP) European Metrology Programme for Innovation and Research (EMPIR) project. The on-site calibration, performed as frequently as possible (ideally monthly) to monitor changes in the device conditions, allows operators to track and evaluate the calibration status on a continuous basis, considerably reducing the data gaps incurred by the periodic shipments for performing centralized calibrations. The performance of the on-site calibration procedures for V0 was very good at sites with low turbidity, showing agreement with a reference calibration between 0.5 % and 1.5 % depending on wavelengths. In the urban area, the agreement decreases between 1.7 % and 2.5 %. For the SVA the difference varied from a minimum of 0.03 % to a maximum of 3.46 %.

Aerosol optical and radiative characteristics under contrasting pollution conditions in a typical urban valley in Northwestern China

Aerosol optical and radiative properties play a key role in climate change. Precisely pinpointing the optical and radiative properties of various types of aerosols during diverse pollution events in urban settings remains a challenge. In this study, aerosol optical and radiative properties under contrasting pollution conditions were investigated, including dusty, haze, fireworks, and clean days, in a typical urban valley in Northwestern China, based on observations from the Sun-Sky photometers and simulations from libRadtran. The results show that on dusty days aerosols have high absorption and low backscattering, while on haze days they are characterized by fine, absorptive organic aerosol particles with pronounced forward scattering in the ultraviolet and visible spectra. On fireworks days during Chinese New Year (CNY), fine-mode aerosols from fireworks dominate, with the highest scattering and the lowest absorption under the four pollution conditions, and the particle peak radius growth responds rapidly to changes in relative humidity. Aerosols generally cause Earth's surface cooling and atmospheric warming across various pollution conditions. Notably, dusty days, clean days, and haze days all exhibit a lower positive aerosol radiative forcing (ARF) in the ultraviolet (UV) spectrum at the top of the atmosphere (TOA) due to lower single-scattering albedo (SSA). In contrast, CNY uniquely displays negative UV_ARF at TOA, attributed to high-SSA non-absorptive particles. Aerosols on dusty days have higher positive radiative forcing than on haze days, emphasizing the absorptive impact of mixed dust aerosols. These findings provide valuable insights into the behavior of aerosols under various pollution conditions in a typical urban valley, contributing to a better understanding of the environmental effects of various aerosols in arid urban regions.

Evaluation of the vertically-resolved aerosol radiative effect on shortwave and longwave ranges using sun-sky photometer and ceilometer measurements

The aerosol radiative effect (ARE) is one of the atmospheric components still affected by large uncertainty. One of the causes is related to the fact that the longwave (LW) component is usually neglected, even though it is necessary for an accurate quantification of the ARE together with the shortwave component (SW). In this study we have developed a methodology based on the GAME (Global Atmospheric Model) radiative transfer model (RTM) that allows to obtain the radiative effect of the atmospheric aerosol for both spectral ranges in an automated way. The microphysical and optical properties necessary to feed the RTM have been obtained through the GRASP (Generalized Retrieval of Aerosol and Surface Properties) algorithm, with the combination of ceilometer and sun-sky photometer data. Data measured in Granada (Spain) during 2017 have been used for the evaluation and implementation of this methodology. According to the results, the ARE in the SW spectral range (ARESW) varies between 0 and − 50 Wm−2 for most of the data, whereas the ARE in the LW range (ARELW) varies between 0 and 5 Wm−2, at heights near the surface. In general, the obtained results agree with those found in the literature, with negative values in the SW range (cooling effect) and positive values in the LW (heating effect). The seasonal analysis shows that, for both components, the ARE is more important during the spring and summer seasons, when the aerosol load is greater, as expected. The analysis of the aerosol heating rate (AHR) shows positive values in the SW and negative values in the LW range. The majority of the AHRSW data varies between 0 and 1 Kd−1 during the year whereas the AHRLW does it between 0 and − 0.15 Kd−1. The seasonal analysis of the AHR shows that the greater monthly average values are found during spring, however there is not much variability along the year, with the exception of February, under the effects of an extreme dust intrusion. The mineral dust particles in this event cause an ARESW of −130 Wm−2 and an ARELW of 23 Wm−2 (ARELW/ARESW = 17%), thus pointing out that the LW component should not be neglected for coarse mode particles. Additionally, it is observed that the vertical distribution of the aerosol layers strongly influences the ARE and the AHR obtained profiles, affecting the way the atmospheric cooling/heating occurs in the vertical coordinate.

Evaluation of the Vertically-Resolved Aerosol Radiative Effect on Shortwave and Longwave Ranges Using Sun-Sky Photometer and Ceilometer Measurements

Evaluation of the Vertically-Resolved Aerosol Radiative Effect on Shortwave and Longwave Ranges Using Sun-Sky Photometer and Ceilometer Measurements

Overview of the SLOPE I and II campaigns: aerosol properties retrieved with lidar and sun–sky photometer measurements

Abstract. The Sierra Nevada Lidar aerOsol Profiling Experiment I and II (SLOPE I and II) campaigns were intended to determine the vertical structure of aerosols by remote sensing instruments and test the various retrieval schemes for obtaining aerosol microphysical and optical properties with in situ measurements. The SLOPE I and II campaigns were developed during the summers of 2016 and 2017, respectively, combining active and passive remote sensing with in situ measurements at stations belonging to the AGORA observatory (Andalusian Global ObseRvatory of the Atmosphere) in the Granada area (Spain). In this work, we use the in situ measurements of these campaigns to evaluate aerosol properties retrieved by the GRASP code (Generalized Retrieval of Atmosphere and Surface Properties) combining lidar and sun–sky photometer measurements. We show an overview of aerosol properties retrieved by GRASP during the SLOPE I and II campaigns. In addition, 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 the absorption coefficient. The statistical analysis of 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 North African intrusions. The vertically resolved analysis denotes a decay of the atmospheric aerosols with an altitude up to 5 km a.s.l. Finally, desert dust and biomass burning events were chosen to show the high potential of GRASP to retrieve vertical profiles of aerosol properties (e.g. absorption coefficient and single scattering albedo) for different aerosol types. The aerosol properties retrieved by GRASP show good agreement with simultaneous in situ measurements (nephelometer, aethalometer, scanning mobility particle sizer, and aerodynamic particle sizer) performed at the Sierra Nevada Station (SNS) in Granada. In general, GRASP overestimates the in situ data at the SNS with a mean difference lower than 6 µm3 cm−3 for volume concentration, and 11 and 2 Mm−1 for the scattering and absorption coefficients. On the other hand, the comparison of GRASP with airborne measurements also shows an overestimation with mean absolute differences of 14 ± 10 and 1.2 ± 1.2 Mm−1 for the scattering and absorption coefficients, showing a better agreement for the absorption (scattering) coefficient with higher (lower) aerosol optical depth. The potential of GRASP shown in this study will contribute to enhancing the representativeness of the aerosol vertical distribution and provide information for satellite and global model evaluation.

Improving the sectional Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) aerosols of the Weather Research and Forecasting-Chemistry (WRF-Chem) model with the revised Gridpoint Statistical Interpolation system and multi-wavelength aerosol optical measurements: the dust aerosol observation campaign at Kashi, near the Taklimakan Desert, northwestern China

Abstract. The Gridpoint Statistical Interpolation data assimilation (DA) system was developed for the four size bin sectional Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) aerosol mechanism in the Weather Research and Forecasting-Chemistry (WRF-Chem) model. The forward and tangent linear operators for the aerosol optical depth (AOD) analysis were derived from WRF-Chem aerosol optical code. We applied three-dimensional variational DA to assimilate the multi-wavelength AOD, ambient aerosol scattering coefficient, and aerosol absorption coefficient, measured by the sun–sky photometer, nephelometer, and aethalometer, respectively. These measurements were undertaken during a dust observation field campaign at Kashi in northwestern China in April 2019. The results showed that the DA analyses decreased the model aerosols' low biases; however, it had some deficiencies. Assimilating the surface particle concentration increased the coarse particles in the dust episodes, but AOD and the coefficients for aerosol scattering and absorption were still lower than those observed. Assimilating aerosol scattering coefficient separately from AOD improved the two optical quantities. However, it caused an overestimation of the particle concentrations at the surface. Assimilating the aerosol absorption coefficient yielded the highest positive bias in the surface particle concentration, aerosol scattering coefficient, and AOD. The positive biases in the DA analysis were caused by the forward operator underestimating aerosol mass scattering and absorption efficiency. As compensation, the DA system increased particle concentrations excessively to fit the observed optical values. The best overall improvements were obtained from the simultaneous assimilation of the surface particle concentration and AOD. The assimilation did not substantially change the aerosol chemical fractions. After DA, the clear-sky aerosol radiative forcing at Kashi was −10.4 W m−2 at the top of the atmosphere, which was 55 % higher than the radiative forcing value before DA.

The characterization of Taklamakan dust properties using a multiwavelength Raman polarization lidar in Kashi, China

Abstract. The Taklamakan desert is an important dust source for the global atmospheric dust budget and a cause of the dust weather in East Asia. The characterization of Taklamakan dust in the source region is still very limited. To fill this gap, the DAO (dust aerosol observation) was conducted in April 2019 in Kashi, China. The Kashi site is about 150 km from the western rim of the Taklamakan desert and is strongly impacted by desert dust aerosols, especially in spring time, i.e., April and May. According to sun–sky photometer measurements, the aerosol optical depth (at 500 nm) varied in the range of 0.07–4.70, and the Ångström exponent (between 440 and 870 nm) in the range of 0.0–0.8 in April 2019. In this study, we provide the first profiling of the 2α+3β+3δ parameters of Taklamakan dust based on a multiwavelength Mie–Raman polarization lidar. For Taklamakan dust, the Ångström exponent related to the extinction coefficient (EAE, between 355 and 532 nm) is about 0.01 ± 0.30, and the lidar ratio is found to be 45 ± 7 sr (51 ± 8–56 ± 8 sr) at 532 (355) nm. The particle linear depolarization ratios (PLDRs) are about 0.28–0.32 ± 0.07 at 355 nm, 0.36 ± 0.05 at 532 nm and 0.31 ± 0.05 at 1064 nm. Both lidar ratios and depolarization ratios are higher than the typical values of Central Asian dust in the literature. The difference is probably linked to the fact that observations in the DAO campaign were collected close to the dust source; therefore, there is a large fraction of coarse-mode and giant particles (radius >20 µm) in the Taklamakan dust. Apart from dust, fine particles coming from local anthropogenic emissions and long-range transported aerosols are also non-negligible aerosol components. The signatures of pollution emerge when dust concentration decreases. The polluted dust (defined by PLDR532≤0.30 and EAE355-532≥0.20) is featured with reduced PLDRs and enhanced EAE355−532 compared to Taklamakan dust. The mean PLDRs of polluted dust generally distributed in the range of 0.20–0.30. Due to the complexity of the nature of the involved pollutants and their mixing state with dust, the lidar ratios exhibit larger variabilities compared to those of dust. The study provides the first reference of novel characteristics of Taklamakan dust measured by Mie–Raman polarization lidar. The data could contribute to complementing the dust model and improving the accuracy of climate modeling.

Ground-based retrievals of aerosol column absorption in the UV spectral region and their implications for GEMS measurements

Quantifying the spectral variation of column aerosol absorption in the ultraviolet (UV) and visible (Vis) wavelengths is required for accurate satellite-based aerosol and trace-gas retrievals. Retrievals of the column-averaged imaginary part of refractive index and single scattering albedo (SSA) in the UV–Vis range have been performed at Yonsei University, Seoul, Korea, since 2016 by combining co-located measurements from the NASA Aerosol Robotic Network (AERONET) Cimel sun-sky photometer, the Ultraviolet Multifilter Rotating Shadowband Radiometer (UV-MFRSR), the SKYNET Prede sky radiometer, and the NASA Pandora sun spectrometer. We investigated the spectral variation of column-averaged imaginary part of refractive index for UV–Vis wavelengths to refine models used in our aerosol retrieval algorithm to process measurements from the upcoming Geostationary Environment Monitoring Satellite (GEMS). The retrieved imaginary part of refractive index for highly absorbing fine pollution particles (BC), dust (DS), and non-absorbing (NA) particles in the selected UV–Vis range (380–440 nm) showed 0–20%, 30%, and 0–40% of spectral dependence, respectively. Retrievals of Ozone Monitoring Instrument (OMI) measurement data using the improved aerosol model showed improved correlation with AERONET data compared to the old algorithm that did not properly account for aerosol absorption effects. These results corroborate the advantage of using local climatology derived from ground-based UV–Vis spectral aerosol absorption measurements for satellite GEMS aerosol retrievals over East Asia. Moreover, this study reveals that spectral variations in the UV column aerosol absorption in East Asia differ from those in other regions.

Evaluation of LIRIC Algorithm Performance Using Independent Sun-Sky Photometer Data at Two Altitude Levels

This work evaluates the Lidar-Radiometer Inversion Code (LIRIC) using sun-sky photometers located at different altitudes in the same atmospheric column. Measurements were acquired during an intensive observational period in summer 2012 at Aerosols, Clouds, and Trace gases Research InfraStructure Network (ACTRIS)/Aerosol Robotic Network (AERONET) Granada (GRA; 37.16°N, 3.61°W, 680 m above sea level (a.s.l.)) and Cerro Poyos (CP; 37.11°N, 3.49°W, 1820 m a.s.l.) sites. Both stations operated AERONET sun-photometry, with an additional lidar system operating at Granada station. The extended database of simultaneous lidar and sun-photometry measurements from this study allowed the statistical analysis of vertically resolved microphysical properties retrieved with LIRIC, with 70% of the analyzed cases corresponding to mineral dust. Consequently, volume concentration values were 46 μm3/cm3 on average, with a value of ~30 μm3/cm3 corresponding to the coarse spheroid mode and concentrations below 10 μm3/cm3 for the fine and coarse spherical modes. According to the microphysical properties’ profiles, aerosol particles reached altitudes up to 6000 m a.s.l., as observed in previous studies over the same region. Results obtained from comparing the LIRIC retrievals from GRA and from CP revealed good agreement between both stations with differences within the expected uncertainties associated with LIRIC (15%). However, larger discrepancies were found for 10% of the cases, mostly due to the incomplete overlap of the lidar signal and/or to the influence of different aerosol layers advected from the local origin located between both stations, which is particularly important in cases of low aerosol loads. Nevertheless, the results presented here demonstrate the robustness and self-consistency of LIRIC and consequently its applicability to large databases such as those derived from ACTRIS-European Aerosol Research Lidar Network (EARLINET) observations.

Different strategies to retrieve aerosol properties at night-time with the GRASP algorithm

Abstract. This study evaluates the potential of the 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 ground-based lidar 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 the 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 aerosol optical depth values.

A 10-year characterization of the Saharan Air Layer lidar ratio in the subtropical North Atlantic

Abstract. Particle extinction-to-backscatter ratio (lidar ratio) is a key parameter for a correct interpretation of elastic lidar measurements. Of particular importance is the determination of the lidar ratio of the Saharan Air Layer mineral dust transported into the free troposphere over the North Atlantic region. The location of the two sun photometer stations managed by the Izaña Atmospheric Research Centre (IARC) on the island of Tenerife and a decade of available micropulse lidar (MPL) data allow us to determine the lidar ratio under almost pure-dust conditions. This result can be considered representative of the Saharan dust transported westward over the North Atlantic in the subtropical belt. Three different methods have been used to calculate the lidar ratio in this work: (1) using the inversion of sky radiance measurements from a sun–sky photometer installed at the Izaña Observatory (2373 m a.s.l.) under free-troposphere conditions; (2) the one-layer method, a joint determination using a micropulse lidar sited at the Santa Cruz de Tenerife sea-level station and photometric information considering one layer of aerosol characterized by a single lidar ratio; and (3) the two-layer method, a joint determination using the micropulse lidar and photometric information considering two layers of aerosol with two different lidar ratios. The one-layer method only uses data from a co-located photometer at Santa Cruz de Tenerife, while the two-layer conceptual approach incorporates photometric information at two heights from the observatories of Izaña and Santa Cruz de Tenerife. The almost pure-dust lidar ratio retrieval from the sun–sky photometer and from the two-layer method give similar results, with lidar ratios at 523 nm of 49 ± 6 and 50 ± 11 sr. These values obtained from a decade of data records are coincident with other studies in the literature reporting campaigns in the subtropical North Atlantic region. This result shows that the two-layer method is an improved conceptual approach compared to the single-layer approach, which matches the real lower-troposphere structure well. The two-layer method is able to retrieve reliable lidar ratios and therefore aerosol extinction profiles despite the inherent limitations of the elastic lidar technique. We found a lack of correlation between lidar ratio and Ångström exponent (α), which indicates that the dust lidar ratio can be considered independent of dust size distribution in this region. This finding suggests that dust is, under most atmospheric conditions, the predominant aerosol in the North Atlantic free troposphere, which is in agreement with previous studies conducted at the Izaña Observatory.

Long-Term Ground-Based Measurements of Aerosol Optical Depth over Kuwait City

We analyze ten years (2008–2017) of ground-based observations of the Aerosol Optical Depth (AOD) in the atmosphere of Kuwait City, in Middle East. The measurements were conducted with a CIMEL sun-sky photometer, at various wavelengths. The daily average AOD at 500 nm (AOD500) is 0.45, while the mean Ångström coefficient (AE), calculated from the pair of wavelengths 440 and 870 nm, is 0.61. The observed high AOD500 values (0.75–2.91), are due to regional sand and dust storm events, which are affecting Kuwait with a mean annual frequency of almost 20 days/year. The long-term record analysis of AOD500 and AE, shows a downward and upward tendency respectively, something which could be attributed to the continuous expansion and industrialization of the main city of Kuwait, in combination with the simultaneous increase of soil moisture over the area. By utilizing back trajectories of air masses for up to 4 days, we assessed the influence of various regions to the aerosol load over Kuwait. The high aerosol loads during spring, are attributed to the dominance of coarse particles from Saudi Arabia (AOD500 0.56–0.74), a source area that contributes the 56% to the mean annual AOD500. Other dust sources affecting significantly Kuwait originated from the regions of Iraq and Iran with contribution of 21%.

A multiwavelength numerical model in support of quantitative retrievals of aerosol properties from automated lidar ceilometers and test applications for AOT and PM<sub>10</sub> estimation

Abstract. The use of automated lidar ceilometer (ALC) systems for the aerosol vertically resolved characterization has increased in recent years thanks to their low construction and operation costs and their capability of providing continuous unattended measurements. At the same time there is a need to convert the ALC signals into usable geophysical quantities. In fact, the quantitative assessment of the aerosol properties from ALC measurements and the relevant assimilation in meteorological forecast models is amongst the main objectives of the EU COST Action TOPROF (“Towards operational ground-based profiling with ALCs, Doppler lidars and microwave radiometers for improving weather forecasts”). Concurrently, the E-PROFILE program of the European Meteorological Services Network (EUMETNET) focuses on the harmonization of ALC measurements and data provision across Europe. Within these frameworks, we implemented a model-assisted methodology to retrieve key aerosol properties (extinction coefficient, surface area, and volume) from elastic lidar and/or ALC measurements. The method is based on results from a large set of aerosol scattering simulations (Mie theory) performed at UV, visible, and near-IR wavelengths using a Monte Carlo approach to select the input aerosol microphysical properties. An average “continental aerosol type” (i.e., clean to moderately polluted continental aerosol conditions) is addressed in this study. Based on the simulation results, we derive mean functional relationships linking the aerosol backscatter coefficients to the abovementioned variables. Applied in the data inversion of single-wavelength lidars and/or ALCs, these relationships allow quantitative determination of the vertically resolved aerosol backscatter, extinction, volume, and surface area and, in turn, of the extinction-to-backscatter ratios (i.e., the lidar ratios, LRs) and extinction-to-volume conversion factor (cv) at 355, 532, and 1064 nm. These variables provide valuable information for visibility, radiative transfer, and air quality applications. This study also includes (1) validation of the model simulations with real measurements and (2) test applications of the proposed model-based ALC inversion methodology. In particular, our model simulations were compared to backscatter and extinction coefficients independently retrieved by Raman lidar systems operating at different continental sites within the European Aerosol Research Lidar Network (EARLINET). This comparison shows good model–measurement agreement, with LR discrepancies below 20 %. The model-assisted quantitative retrieval of both aerosol extinction and volume was then tested using raw data from three different ALCs systems (CHM 15k Nimbus), operating within the Italian Automated LIdar-CEilometer network (ALICEnet). For this purpose, a 1-year record of the ALC-derived aerosol optical thickness (AOT) at each site was compared to direct AOT measurements performed by colocated sun–sky photometers. This comparison shows an overall AOT agreement within 30 % at all sites. At one site, the model-assisted ALC estimation of the aerosol volume and mass (i.e., PM10) in the lowermost levels was compared to values measured at the surface level by colocated in situ instrumentation. Within this exercise, the ALC-derived daily-mean mass concentration was found to reproduce the corresponding (EU regulated) PM10 values measured by the local air quality agency well in terms of both temporal variability and absolute values. Although limited in space and time, the good performances of the proposed approach suggest it could possibly represent a valid option to extend the capabilities of ALCs to provide quantitative information for operational air quality and meteorological monitoring.

Comparative assessment of GRASP algorithm for a dust event over Granada (Spain) during ChArMEx-ADRIMED 2013 campaign

Abstract. In this study, vertical profiles and column-integrated aerosol properties retrieved by the GRASP (Generalized Retrieval of Atmosphere and Surface Properties) algorithm are evaluated with in situ airborne measurements made during the ChArMEx-ADRIMED field campaign in summer 2013. In the framework of this campaign, two different flights took place over Granada (Spain) during a desert dust episode on 16 and 17 June. The GRASP algorithm, which combines lidar and sun–sky photometer data measured at Granada, was used to retrieve aerosol properties. Two sun-photometer datasets are used: one co-located with the lidar system and the other in the Cerro Poyos station, approximately 1200 m higher than the lidar system but at a short horizontal distance. Column-integrated aerosol microphysical properties retrieved by GRASP are compared with AERONET products showing a good agreement. Differences between GRASP retrievals and airborne extinction profiles are in the range of 15 to 30 %, depending on the instrument on board the aircraft used as reference. On 16 June, a case where the dust layer was coupled to the aerosol layer close to surface, the total volume concentration differences between in situ data and GRASP retrieval are 15 and 36 % for Granada and Cerro Poyos retrievals, respectively. In contrast, on 17 June the dust layer was decoupled from the aerosol layer close to the surface, and the differences are around 17 % for both retrievals. In general, all the discrepancies found are within the uncertainly limits, showing the robustness and reliability of the GRASP algorithm. However, the better agreement found for the Cerro Poyos retrieval with the aircraft data and the vertical homogeneity of certain properties retrieved with GRASP, such as the scattering Ångström exponent, for cases with aerosol layers characterized by different aerosol types, shows that uncertainties in the vertical distribution of the aerosol properties have to be considered. The comparison presented here between GRASP and other algorithms (i.e. AERONET and LIRIC) and with airborne in situ measurements shows the potential to retrieve the optical and microphysical profiles of the atmospheric aerosol properties. Also, the advantage of GRASP versus LIRIC is that GRASP does not assume the results of the AERONET inversion as a starting point.

Satellite and ground-based monitoring of smoke in the atmosphere during the summer wildfires in European Russia in 2010 and Siberia in 2012

Using Moderate Resolution Imaging Spectroradiometer (MODIS) (Aqua and Terra satellites) and in situ observations, a comparative analysis of two large-scale smoke events caused by the summer wildfires in European Russia (ER) in 2010 and Western Siberia (WS) in 2012 was carried out. In the 5-day periods of the extreme smoke pollution (5–9 August 2010 in ER and 27–31 July 2012 in WS), the number of active fires in the equal territories, confined by the coordinates 47°–65° N, 25°–55° E and 51°–70° N, 71°–104° E, was found to be 4754 for ER and 3823 for WS. With this, the regional mean aerosol optical depths (AODs) were found to be (1.02 ± 0.02) and (1.00 ± 0.04), not much differing for both the events. The regional mean aerosol radiative forcing effects at the top (R1) and the bottom (R2) of the atmosphere over ER/WS according to MODIS observations were estimated to be (−61 ± 1) and (−54 ± 2) W m−2, and (−107 ± 2) and (−96 ± 3) W m−2, respectively. At the same time, the local values of AOD and the local absolute values of R1 and R2 over WS were considerably higher than those over ER. MODIS AOD (L3) data during the wildfires of 2010 were validated by AOD data obtained by the sun-sky photometer CIMEL, operating at the AERONET station Zvenigorod. The rates of radiative heating of the smoky atmosphere over ER and WS were also estimated and compared with the existed temperature anomalies, obtained using National Centers for Environmental Prediction National Center for Atmospheric Research reanalysis data. Optical and microphysical properties of smoke aerosols during the wildfires in ER and WS also revealed some similar characteristics. The aerosols were mostly found in the submicron-size fraction and characterized by very high single-scattering albedos (0.95–0.98). In the dense smoke conditions, the degree of linear polarization at the scattering angle 90° during both the events decreased to negative values ranging between −0.1 and −0.15. The optical properties of smoke aerosols were mainly conditioned by unusually narrow particle size distribution.

Aerosol optical and microphysical properties observed by the lidar technique from a forest-fire smoke event over Madrid

On 27 August 2012, a wildfire occurred in the western zone of the Madrid region. Consequently, a significant release of smoke and aerosols was injected into the free troposphere and advected by the synoptic circulation in trajectories that passed over the capital. This event was detected by the Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT; Research Centre for Energy, Environment and Technology), Madrid lidar station at about 5 km height. While the main aerosol plume did not pass over the lidar station, the eastern portion of the dispersed emissions did. No aerosol layers were detected by lidar in the free-troposphere for the next few days. Since few observations of fresh biomass burning aerosol have been carried out by lidar so far, this study contributes to previous studies on biomass burning aerosols, in particular regarding the first stages of the smoke transport.This article studies the optical and microphysical properties of the smoke aerosols retrieved by the lidar technique along with the recently developed Lidar/Radiometer Inversion Code (LiRIC), which combines both lidar and sun–sky photometer data. The sun–sky photometer data are provided by the Agencia Estatal de Meteorología (AEMET; Spanish Meteorological Agency) station close to the CIEMAT site. We suggest that the aerosol properties retrieved might be linked to the water uptake phenomenon in view of the fact that simultaneous water vapour lidar measurements exhibited a close relationship with backscatter coefficient profiles. This article is concerned not only with showing the capability of combining different remote-sensing techniques and the LiRIC code to provide aerosol vertical distribution for fine and coarse modes, but also with providing signs about the plausible interaction between aerosols and humidity that leads fresh biomass burning aerosols to act as cloud condensation nuclei. This key role, known as the indirect effect, remains the major source of uncertainty when the global radiation budget is assessed.

Development of a new data-processing method for SKYNET sky radiometer observations

Abstract. In order to reduce uncertainty in the estimation of Direct Aerosol Radiative Forcing (DARF), it is important to improve the estimation of the single scattering albedo (SSA). In this study, we propose a new data processing method to improve SSA retrievals for the SKYNET sky radiometer network, which is one of the growing number of networks of sun-sky photometers, such as NASA AERONET and others. There are several reports that SSA values from SKYNET have a bias compared to those from AERONET, which is regarded to be the most accurate due to its rigorous calibration routines and data quality and cloud screening algorithms. We investigated possible causes of errors in SSA that might explain the known biases through sensitivity experiments using a numerical model, and also using real data at the SKYNET sites at Pune (18.616° N/73.800° E) in India and Beijing (39.586° N/116.229° E) in China. Sensitivity experiments showed that an uncertainty of the order of ±0.03 in the SSA value can be caused by a possible error in the ground surface albedo or solid view angle assumed for each observation site. Another candidate for possible error in the SSA was found in cirrus contamination generated by imperfect cloud screening in the SKYNET data processing. Therefore, we developed a new data quality control method to get rid of low quality or cloud contamination data, and we applied this method to the real observation data at the Pune site in SKYNET. After applying this method to the observation data, we were able to screen out a large amount of cirrus-contaminated data and to reduce the deviation in the SSA value from that of AERONET. We then estimated DARF using data screened by our new method. The result showed that the method significantly reduced the difference of 5 W m−2 that existed between the SKYNET and AERONET values of DARF before screening. The present study also suggests the necessity of preparing suitable a priori information on the distribution of coarse particles ranging in radius between 10 μm and 30 μm for the analysis of heavily dust-laden atmospheric cases.

Columnar aerosol properties in a Northeastern Atlantic site (Plymouth, United Kingdom) by means of ground based skyradiometer data during years 2000–2008

Between 2000 and 2008, columnar optical and radiative properties were measured at the Plymouth Marine Laboratory (PML), UK (50° 21.95′N, 4° 8.85′W) using an automatic Prede POM01L sun–sky photometer. The database was analyzed for aerosol optical properties using the SKYRAD radiative inversion algorithm and calibrated using the in situ SKYIL calibration method. Retrievals include aerosol optical depth, Ångström wavelength exponent, aerosol volume distribution, refractive index and single scattering albedo. The results show that the Plymouth site is characterized by low values of aerosol optical depth with low variability (0.18 ± 0.08 at 500 nm) and a mean annual Ångström exponent of 1.03 ± 0.21. The annual mean of the single scattering albedo is 0.97, indicative of non-absorbing aerosols. The aerosol properties were classified in terms of air mass back trajectories: the area is mainly affected by Atlantic air masses and the dominant aerosol type is a mixture of maritime particles, present in low burdens with variable size. The maritime air masses were defined by annual mean values for the AOD (at 500 nm) of 0.13–0.14 and a wavelength exponent of 0.96–1.03. Episodic anthropogenic and mineral dust intrusions occasionally occur, but they are sporadic and dilute (AOD at 500 nm about 0.20). Tropical continental air masses were characterized by the highest AOD at 500 nm (0.34) and the lowest wavelength exponent (0.83), although they were the least represented in the analysis.

Dust layer effects on the atmospheric radiative budget and heating rate profiles

The effect of mineral aerosol optical properties and vertical distribution on clear-sky, instantaneous and daily-average aerosol direct radiative effects (DREs) and heating rates (HRs) is analyzed in the solar (S, 0.3–4 μm) and terrestrial (T, 4–80 μm) spectral domain, respectively. The used radiative transfer model is based on lidar, sun-sky photometer, and radiosonde measurements. The study focuses on the Sahara dust outbreak of July 16, 2009 which advected dust particles from north-western Africa over south-eastern Italy. Clear-sky, instantaneous aerosol DREs and HRs undergo large changes within few hours, for the variability of the dust aerosol properties and vertical distribution. The daily-average, clear-sky aerosol S-DRE is near −5 Wm−2 and −12 Wm−2 at the top of the atmosphere (ToA) and surface (sfc), respectively. The daily-average aerosol T-DRE offsets the S-DRE by about one third at the ToA and by about one half at the surface. The daily average aerosol HR integrated over the whole aerosol column is 0.5 and −0.3 K day−1 in the S and T domain, respectively. Thus, the all-wave integrated HR is 0.2 K day−1. These results highlight the importance of accounting for the interaction of dust particles with T and S radiation. Sensitivity tests indicate that the uncertainties of the aerosol refractive index, size distribution, and vertical distribution have on average a large impact on aerosol HRs in the S and T domain, respectively. Refractive index and aerosol size distribution uncertainties also have a large impact on S- and T-DREs. The aerosol vertical distribution that has a negligible impact on aerosol S-DREs, is important for aerosol T-DREs. It is also shown that aerosol HRs and DREs in the terrestrial domain are affected by the water vapour vertical distribution.

Direct radiative forcing during Sahara dust intrusions at a site in the Central Mediterranean: Anthropogenic particle contribution

The clear-sky, instantaneous Direct Radiative Effect (DRE) by all and anthropogenic particles is calculated during selected dust intrusion events at a site in the Central Mediterranean, to evaluate the role of anthropogenic particle's radiative effects and obtain a better estimate of the DRE by desert dust. The Mediterranean basin is continuously influenced by long-range transported emissions from continental Europe and northern Africa. The clear-sky aerosol DRE is calculated by a two stream radiative transfer model in the solar (0.3–4 μm) and infrared (4–200 μm) spectral range, at the top of the atmosphere (ToA) and at the Earth's surface (sfc), during typical weak and strong dust intrusion events. Aerosol optical properties by AERONET sun-sky photometer measurements and aerosol vertical profiles by EARLINET lidar measurements, both performed at Lecce (40.33°N, 18.10° E) in the Central Mediterranean, are used to perform radiative transfer simulations. Instantaneous solar DREs are negative as a consequence of the cooling effect by aerosol particles and span the – (25–11) W m − 2 and – (53–25) W m − 2 range at the ToA and surface, respectively during the investigated dusty days. Infrared DREs offset the solar DRE from 8% up to 14% at the ToA and from 9% up to 22% at the surface. A methodology has been implemented to estimate the contribution of anthropogenic particles during dusty days since their contribution can be significant over the Central Mediterranean. We have found that the aerosol optical depth (AOD) by anthropogenic particles is from 35% up to 65% of the AOD by all particles. The solar DRE by anthropogenic particles represents from 44% up to 83% and from 28% up to 60% of the DRE by all particles, at the ToA and surface, respectively. If we assume that the DRE by natural and anthropogenic particles can be linearly added, the contribution of natural particles to the whole aerosol DRE is lower than 56% and 72% at the ToA and surface, respectively during the tested dusty days.