Search NASA⌕ Search

Engineering topics

Ilya Slutsker

Publications and source records attributed to Ilya Slutsker.

The Extreme Forest Fires in California/Oregon in 2020: Aerosol Optical and Physical Properties and Comparisons of Aged Versus Fresh Smoke

Wildfire activity in the western United States during August to October 2020 was exceptional in terms of the fire severity and area burned. Extremely dry biomass fuels from near historic multi-year drought conditions were further exacerbated with very hot and dry conditions in 2020. These conditions when coupled with strong offshore flow allowed many ignitions to grow into extremely large and severe wildfires. Long-term monitoring at a few AERONET sites in California showed that the number of days with high Aerosol Optical Depth at 440 nm (AOD 440 >1) in 2020 was greater than any other year going back to the beginning of the data records in 2002. A wide range of fine mode particle volume median radii were retrieved from AERONET data over the course of these fires suggesting significant variability in combustion conditions and aging processes. Additionally, the fine mode radii in some of these smoke plumes in 2020 were very large especially at high AOD (∼0.22–0.32 μm volume median radius), likely due to both coagulation and condensation occurring during aging at very high particulate concentrations. The largest fine mode particle radii combined with narrow distributions resulted in some very rare AOD spectra showing peak AOD at 500 nm and decreasing to lower AOD at both shorter and longer wavelengths. The most extreme retrieved size distributions and associated measured AOD spectra were principally observed in long-distance transported smoke plumes from these western United States fires at sites in Colorado, Maryland and Virginia, possibly due to further aging during transport. Additionally, strong absorption was sometimes observed at short wavelengths with much lower single scattering albedo at 440 nm compared to 675 nm in some plumes consistent with significant brown carbon (BrC) and/or coated black carbon (BC) absorption in biomass burning particles. This strong spectral absorption signature observed at some California sites and dates remained similarly strong in some smoke plumes observed at some east coast sites in Maryland and Virginia, thereby suggesting that the lifetime of these particular BrC and/or coated BC absorbing species was greater than 5 days.

Aerosol↗

Employing Relaxed Smoothness Constraints on Imaginary Part of Refractive Index in AERONET Aerosol Retrieval Algorithm

In the Aerosol Robotic Network (AERONET) retrieval algorithm, smoothness constraints on the imaginary part of the refractive index provide control of retrieved spectral dependence of aerosol absorption by preventing the inversion code from fitting the noise in optical measurements and thus avoiding unrealistic oscillations of retrievals with wavelength. The history of implementation of the smoothness constraints in the AERONET retrieval algorithm is discussed. It is shown that the latest version of the smoothness constraints on the imaginary part of refractive index, termed standard and employed by Version 3 of the retrieval algorithm, should be modified to account for strong variability of light absorption by brown-carbon-containing aerosols in UV through mid-visible parts of the solar spectrum. In Version 3 strong spectral constraints were imposed at high values of the Ångström exponent (440–870 nm) since black carbon was assumed to be the primary absorber, while the constraints became increasingly relaxed as aerosol exponent deceased to allow for wavelength dependence of absorption for dust aerosols. The new version of the smoothness constraints on the imaginary part of the refractive index assigns different weights to different pairs of wavelengths, which are the same for all values of the Ångström exponent. For example, in the case of four-wavelength input, the weights assigned to short-wavelength pairs (440–675, 675–870 nm) are small so that smoothness constraints do not suppress natural spectral variability of the imaginary part of the refractive index. At longer wavelengths (870–1020 nm), however, the weight is 10 times higher to provide additional constraints on the imaginary part of refractive index retrievals of aerosols with a high Ångström exponent due to low sensitivity to aerosol absorption for longer channels at relatively low aerosol optical depths. The effect of applying the new version of smoothness constraints, termed relaxed, on retrievals of single-scattering albedo is analyzed for case studies of different aerosol types: black- and brown-carbon-containing fine mode aerosols, mineral dust coarse mode aerosols, and urban industrial fine mode aerosol. It is shown that for brown-carbon-containing aerosols employing the relaxed smoothness constraints resulted in significant reduction in retrieved single-scattering albedo and spectral residual errors (compared to standard) at the short wavelengths. For example, biomass burning smoke cases showed a reduction in single-scattering albedo and spectral residual error at 380 nm of ∼ 0.033 and ∼ 17 %, respectively, for the Rexburg site and ∼ 0.04 and ∼ 12.7 % for the Rimrock site, both AERONET sites in Idaho, USA. For a site with very high levels of black-carbon-containing aerosols (Mongu, Zambia), the effect of modification in the smoothness constraints was minor. For mineral dust aerosols at small Ångström exponent values (Mezaira site, UAE), the spectral constraint on the imaginary part of the refractive index was already relaxed in Version 3; therefore the new relaxed constraint results in minimal change. In the case of weakly absorbing urban industrial aerosols at the GSFC site, there are significant changes in retrieved single-scattering albedo using relaxed assumption, especially reductions at longer wavelengths: ∼ 0.016 and ∼ 0.02 at 875 and 1020 nm, respectively, for 440 nm aerosol optical depth (AOD) ∼ 0.3. The modification of smoothness constraints on the imaginary part of the refractive index has a minor effect on retrievals of other aerosol parameters such as the real part of the refractive index and parameters of the aerosol size distribution. The implementation of the relaxed smoothness constraints on the imaginary part of the refractive index in the next version of the AERONET inversion algorithm will produce significant impacts at some sites in short wavelength channels (380 and 440 nm) for some biomass burning smoke cases with significant brown carbon content and possibly in mid-visible channels (500 and 675 nm) to near-infrared channels (870 to 1020 nm) for some urban industrial aerosol types. However, most differences in single-scattering albedo retrievals between those applying the new relaxed constraint and the standard constraint will be within the uncertainty of the single-scattering albedo retrievals, depending on the level of aerosol optical depth, Ångström exponent, brown carbon content and wavelength.

Aliaksandr Siniuk↗

Analysis of Aeronet Extended Wavelength Retrievals of Aerosol Absorption Parameters Including 380 NM and 500 NM for Detection of Brown Carbon in Biomass Burning and Iron Oxides in Desert Dust

To investigate near UV aerosol absorbers, such as iron in mineral dust and brown carbon from biomass burning, AERONET spectral imaginary refractive index and size distributions retrievals were performed with the input of additional wavelength data. This includes adding both measured spectral AOD and sky radiances at 380 nm and 500 nm to the four channels (440, 675, 870 and 1020 nm) that are used in the standard retrieval data product. The vector radiative transfer code in the AERONET Version 3 retrievals enabled more accurate computations in the ultraviolet region (380 nm). Calibration of the 380 nm channel radiances was done by the vicarious method since the integrating sphere radiance source that is used for all other wavelengths does not produce sufficient energy at 380 nm. Spectral single scattering albedo, imaginary refractive indices and aerosol size distributions from these six wavelength retrievals are compared to the standard four wavelength retrieval values.

Aerosol↗

Characterization and Calibration of an AERONET-OC Radiometer

The Ocean Color component of the global Aerosol Robotic Network (AERONET-OC) utilizes CE-318 sun photometers modified for above-water radiometry from fixed structures such as oil rigs, lighthouses and service platforms. Primarily, AERONET-OC measurements allow determination of the water-leaving radiance required for the validation of ocean color satellite data products. In this study, we characterized and calibrated one instrument from the AERONET-OC network, identified as AERONET #080. A laser-illuminated integrating sphere of known radiance enabled determination of the absolute radiance responsivity for 7 bands. We compared the results to calibrations from the AERONET facility at the Goddard Space Flight Center of the National Aeronautics and Space Administration and from the Joint Research Centre of the European Commission. These results agree within the estimated mean comparison uncertainty of 1.88 % (k = 2). We also assessed these results using calibrated lamp-illuminated integrating spheres and observed a spectral dependence to the comparison results that is unexplained.

B. C. Johnson↗

Advances in the Ocean Color Component of the Aerosol Robotic Network (AERONET-OC)

The Ocean Color Component of the Aerosol Robotic Network(AERONET-OC) supports activities related to ocean color such as validation of satellite data products, assessment of atmospheric correction schemes, and evaluation of bio-optical models through globally distributed standardized measurements of water-leaving radiance and aerosol optical depth. In view of duly assisting the AERONET-OC data user community, this work (i) summarizes the latest investigations on a number of scientific issues related to above-water radiometry, (ii) emphasizes the network expansion that from 2002 until the end of 2020 integrated 31 effective measurement sites, (iii) shows the equivalence of data product accuracy across sites and time for measurements performed with different instrument series, (iv) illustrates the variety of water types represented by the network sites ensuring validation activities across a diversity of observation conditions, and (v) documents the availability of water-leaving radiance data corrected for bidirectional effects by applying a method specifically developed for chlorophyll-a-dominated waters and an alternative one that is likely suitable for any water type.

In situ oceanic observations↗

The AERONET Version 3 Aerosol Retrieval Algorithm, Associated Uncertainties and Comparisons to Version 2

The Aerosol Robotic Network (AERONET) Version 3 (V3) aerosol retrieval algorithm is described, which is based on the Version 2 (V2) algorithm with numerous updates. Comparisons of V3 aerosol retrievals to those of V2 are presented, along with a new approach to estimate uncertainties in many of the retrieved aerosol parameters. Changes in the V3 aerosol retrieval algorithm include (1) a new polarized radiative transfer code (RTC), which replaced the scalar RTC of V2, (2) detailed characterization of gas absorption by adding NO2 and H2O to specify total gas absorption in the atmospheric column, specification of vertical profiles of all the atmospheric species, (3) new bidirectional reflectance distribution function (BRDF) parameters for land sites adopted from the MODIS BRDF/Albedo product, (4) a new version of the extraterrestrial solar flux spectrum, and (5) a new temperature correction procedure of both direct Sun and sky radiance measurements. The potential effect of each change in V3 on single scattering albedo (SSA) retrievals was analyzed. The operational almucantar retrievals of V2 versus V3 were compared for four AERONET sites: GSFC, Mezaira, Mongu, and Kanpur. Analysis showed very good agreement in retrieved parameters of the size distributions. Comparisons of SSA retrievals for dust aerosols (Mezaira) showed a good agreement in 440 nm SSA, while for longer wavelengths V3 SSAs are systematically higher than those of V2, with the largest mean difference at 675 nm due to cumulative effects of both extraterrestrial solar flux and BRDF changes. For non-dust aerosols, the largest SSA deviation is at 675 nm due to differences in extraterrestrial solar flux spectrums used in each version. Further, the SSA 675 nm mean differences are very different for weakly (GSFC) and strongly (Mongu) absorbing aerosols, which is explained by the lower sensitivity to a bias in aerosol scattering optical depth by less absorbing aerosols. A new hybrid (HYB) sky radiance measurement scan is introduced and discussed. The HYB combines features of scans in two different planes to maximize the range of scattering angles and achieve scan symmetry, thereby allowing for cloud screening and spatial averaging, which is an advantage over the principal plane scan that lacks robust symmetry. We show that due to an extended range of scattering angles, HYB SSA retrievals for dust aerosols exhibit smaller variability with solar zenith angles (SZAs) than those of almucantar (ALM), which allows extension of HYB SSA retrievals to SZAs less than 50∘ to as small as 25∘. The comparison of SSA retrievals from closely time-matched HYB and ALM scans in the 50 to 75∘ SZA range showed good agreement with the differences below ∼0.005. We also present an approach to estimate retrieval uncertainties which utilizes the variability in retrieved parameters generated by perturbing both measurements and auxiliary input parameters as a proxy for retrieval uncertainty. The perturbations in measurements and auxiliary inputs are assumed as estimated biases in aerosol optical depth (AOD), radiometric calibration of sky radiances combined with solar spectral irradiance, and surface reflectance. For each set of Level 2 Sun/sky radiometer observations, 27 inputs corresponding to 27 combinations of biases were produced and separately inverted to generate the following statistics of the inversion results: average, standard deviation, minimum and maximum values. From these statistics, standard deviation (labeled U27) is used as a proxy for estimated uncertainty, and a lookup table (LUT) approach was implemented to reduce the computational time. The U27 climatological LUT was generated from the entire AERONET almucantar (1993–2018) and hybrid (2014–2018) scan databases by binning U27s in AOD (440 nm), Angström exponent (AE, 440–870 nm), and SSA (440, 675, 870, 1020 nm). Using this LUT approach, the uncertainty estimates U27 for each individual V3 Level 2 retrieval can be obtained by interpolation using the corresponding measured and inverted combination of AOD, AE, and SSA.

Aliaksandr Siniuk↗

Advancements in the Aerosol Robotic Network (AERONET) Version 3 database – automated near-real-time quality control algorithm with improved cloud screening for Sun photometer aerosol optical depth (AOD) measurements

The Aerosol Robotic Network (AERONET) has provided highly accurate, ground-truth measurements of the aerosol optical depth (AOD) using Cimel Electronique Sun–sky radiometers for more than 25 years. In Version 2 (V2) of the AERONET database, the near-real-time AOD was semiautomatically quality controlled utilizing mainly cloud-screening methodology, while additional AOD data contaminated by clouds or affected by instrument anomalies were removed manually before attaining quality-assured status (Level 2.0). The large growth in the number of AERONET sites over the past 25 years resulted in significant burden to the manual quality control of millions of measurements in a consistent manner. The AERONET Version 3 (V3) algorithm provides fully automatic cloud screening and instrument anomaly quality controls. All of these new algorithm updates apply to near-real-time data as well as post-field-deployment processed data, and AERONET reprocessed the database in 2018. A full algorithm redevelopment provided the opportunity to improve data inputs and corrections such as unique filter-specific temperature characterizations for all visible and near-infrared wavelengths, updated gaseous and water vapor absorption coefficients, and ancillary data sets. The Level 2.0 AOD quality-assured data set is now available within a month after post-field calibration, reducing the lag time from up to several months. Near-real-time estimated uncertainty is determined using data qualified as V3 Level 2.0 AOD and considering the difference between the AOD computed with the pre-field calibration and AOD computed with pre-field and post-field calibration. This assessment provides a near-real-time uncertainty estimate for which average differences of AOD suggest a +0.02 bias and one sigma uncertainty of 0.02, spectrally, but the bias and uncertainty can be significantly larger for specific instrument deployments. Long-term monthly averages analyzed for the entire V3 and V2 databases produced average differences (V3–V2) of +0.002 with a ±0.02 SD (standard deviation), yet monthly averages calculated using time-matched observations in both databases were analyzed to compute an average difference of −0.002 with a ±0.004 SD. The high statistical agreement in multiyear monthly averaged AOD validates the advanced automatic data quality control algorithms and suggests that migrating research to the V3 database will corroborate most V2 research conclusions and likely lead to more accurate results in some cases.

David M. Giles↗

The NASA Micro Pulse Lidar Network (MPLNET): Early Results from Development of Diurnal Climatologies

The NASA Micro Pulse Lidar Network (MPLNET) is a federated network of micro pulse lidar sites, mostly co-located with the NASA Aerosol Robotic Network (AERONET), providing information on the vertical properties of both aerosols and clouds. MPLNET began in 2000, and has grown to include ~80 sites worldwide. MPLNET data is collected continuously, offering the ability to examine diurnal changes in aerosol and cloud properties. To date, twelve sites have over 10 years of data and more will hit this milestone soon. Here we present early results from development of Level 3 monthly diurnal climatologies from MPLNET.

Ellsworth J Welton↗