Articles | Volume 12, issue 6
Atmos. Meas. Tech., 12, 3269–3288, 2019

Special issue: Holistic Analysis of Aerosol in Littoral Environments - A...

Atmos. Meas. Tech., 12, 3269–3288, 2019

Research article 20 Jun 2019

Research article | 20 Jun 2019

Detecting layer height of smoke aerosols over vegetated land and water surfaces via oxygen absorption bands: hourly results from EPIC/DSCOVR in deep space

Detecting layer height of smoke aerosols over vegetated land and water surfaces via oxygen absorption bands: hourly results from EPIC/DSCOVR in deep space
Xiaoguang Xu1,2, Jun Wang1, Yi Wang1, Jing Zeng1, Omar Torres3, Jeffrey S. Reid4, Steven D. Miller5, J. Vanderlei Martins2, and Lorraine A. Remer2 Xiaoguang Xu et al.
  • 1Department of Chemical and Biochemical Engineering, Center for Global and Regional Environmental Research, and Informatics Initiative, The University of Iowa, Iowa City, Iowa 52241, USA
  • 2Joint Center for Earth Systems Technology and Department of Physics, University of Maryland Baltimore County, Baltimore, Maryland 21250, USA
  • 3Atmospheric Chemistry and Dynamics Laboratory, NASA Goddard Space Flight Center, Greenbelt, Maryland 20770, USA
  • 4Marine Meteorology Division, Naval Research Laboratory, Monterey, California 93943, USA
  • 5Cooperative Institute for Research in the Atmosphere, Colorado State University, Fort Collins, Colorado 80523, USA

Correspondence: Jun Wang ( and Xiaoguang Xu (


We present an algorithm for retrieving aerosol layer height (ALH) and aerosol optical depth (AOD) for smoke over vegetated land and water surfaces from measurements of the Earth Polychromatic Imaging Camera (EPIC) onboard the Deep Space Climate Observatory (DSCOVR). The algorithm uses Earth-reflected radiances in six EPIC bands in the visible and near-infrared and incorporates flexible spectral fitting that accounts for the specifics of land and water surface reflectivity. The fitting procedure first determines AOD using EPIC atmospheric window bands (443, 551, 680, and 780 nm), then uses oxygen (O2) A and B bands (688 and 764 nm) to derive ALH, which represents an optical centroid altitude. ALH retrieval over vegetated surface primarily takes advantage of measurements in the O2 B band. We applied the algorithm to EPIC observations of several biomass burning events over the United States and Canada in August 2017. We found that the algorithm can be used to obtain AOD and ALH multiple times daily over water and vegetated land surface. Validation is performed against aerosol extinction profiles detected by the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) and against AOD observed at nine Aerosol Robotic Network (AERONET) sites, showing, on average, an error of 0.58 km and a bias of −0.13km in retrieved ALH and an error of 0.05 and a bias of 0.03 in retrieved AOD. Additionally, we show that the aerosol height information retrieved by the present algorithm can potentially benefit the retrieval of aerosol properties from EPIC's ultraviolet (UV) bands.

1 Introduction

Aerosol vertical distribution is an important but poorly constrained variable that strongly influences how aerosol particles affect Earth's energy budget. In particular, absorption of solar radiation by smoke and dust particles can result in diabatic heating, alter atmospheric stability, and affect cloud formation and life cycle. These effects depend critically on the altitude of aerosol layers (Babu et al., 2011; Ge et al., 2014; Koch and Del Genio, 2010; Satheesh et al., 2008; Wendisch et al., 2008). Consequently, the aerosol profile factors into the magnitude and even the sign of aerosol direct and indirect effects. An accurate representation of aerosol altitude is thus essential for model prediction of weather and climate (Choi and Chung, 2014; Samset et al., 2013). The thermal signature of dust, in particular, can likewise influence the Earth's longwave budget and, through the interference of retrievals of water vapor and temperature, influence the measurement of atmospheric state (Maddy et al., 2012). Additionally, knowledge of ALH is essential for retrieving aerosol absorption properties in the ultraviolet (UV) channels (Torres et al., 1998), aerosol microphysical properties from multi-angular photopolarimetric measurements (Chowdhary et al., 2005; Waquet et al., 2009), and for atmospheric correction for ocean color remote sensing (Duforêt et al., 2007).

Despite the importance of aerosol vertical distribution, the simulation of aerosol layer height (ALH) in current climate models is subject to large inter-model variation and uncertainty (Kipling et al., 2016; Koffi et al., 2012). Although the assimilation of spaceborne lidar (e.g., CALIOP, or the Cloud–Aerosol Lidar with Orthogonal Polarization) observations can improve the prediction of vertical allocation (Zhang et al., 2011), a significant challenge remains due to the sparsity of lidar observations over space and time. Therefore, frequent satellite observations of global aerosol vertical distribution based on more abundant observation sources are critically needed.

One such source, the recently launched Deep Space Climate Observatory (DSCOVR) mission, has introduced an unprecedented opportunity to acquire ALH information multiple times daily. Launched on 11 February 2015, the DSCOVR spacecraft flies a Lissajous orbit around the Sun–Earth L1 Lagrangian point about 1.5 million kilometers from Earth in the direction of the Sun. While its primary payload is the plasma magnetometer that measures solar wind, DSCOVR also carries the Earth Polychromatic Imaging Camera (EPIC) to image the sunlit disk of Earth every 60–100 min from its stable L1 vantage point. EPIC captures its imagery at a resolution of 2048 by 2018 pixels, with about 8 km of pixel resolution at the disk's center. EPIC was equipped with a double-wheel filter to sequentially measure Earth-reflected solar radiance in 10 narrow bands spanning from the UV to the near-infrared (NIR) spectrum (Fig. 1a), with central wavelengths at 318, 325, 340, 388, 443, 551, 680, 688, 764, and 780 nm (Marshak et al., 2018).

Figure 1EPIC instrument filter response function (blue) and atmospheric spectral transmission (orange). Panel (a) includes all 10 EPIC bands, whereas panels (b) and (c) show zoomed-in views for the 688 nm channel in the O2 B band and the 764 nm channel in the O2 A band, respectively. Here the atmospheric transmission is simulated by the UNL-VRTM, with a spectral step size and a spectral full-width at half maximum of 0.02 nm.


As shown in Fig. 1b–c, two of EPIC's bands are located within the oxygen (O2) “A” and “B” bands (764 and 688 nm), each associated with a reference continuum band at 780 and 680 nm, respectively. These four bands, offering spectral contrasts between absorption bands and continuum bands (known as differential optical absorption spectroscopic, or DOAS, ratios), were originally designed for determining cloud height (Yang et al., 2013). Recently, Xu et al. (2017) presented an algorithm to simultaneously retrieve aerosol optical depth (AOD) and ALH using the EPIC measurements via these four bands and, for the first time, demonstrated EPIC's promising application for determining dust plume height over ocean surfaces during daytime hours.

The present study, building upon the development of Xu et al. (2017) for determining dust ALH over ocean, extends the algorithm to retrieve ALH from EPIC measurements over land surfaces as well. The new developments in this study include implementing smoke aerosol optical properties, land surface characterization, and more robust strategies for the procedures of pixel selection and spectral fitting. The augmentation of the Xu et al. (2017) algorithm takes an important additional step towards our goal of providing more frequent global ALH and AOD information for multi-species global aerosol.

The paper is outlined as follows. Section 2 briefly reviews the physical principle for remote sensing of ALH in the O2 absorption bands and explains the challenges for retrieving ALH over land surfaces, limiting our algorithm development to water and vegetated land surface. Section 3 describes the key elements, procedures, and assumptions of the updated ALH retrieval algorithm. ALH retrievals of smoke events over Canada and the United States in August 2017 are demonstrated in Sect. 4. Section 5 evaluates the retrieved smoke ALH and AOD from EPIC against aerosol extinction profiles measured by CALIOP and AOD measured at AERONET sites. Section 6 discusses the potential application of ALH retrievals for determining UV aerosol properties. Section 7 summarizes our findings and discusses potential broader applications of this new ALH information.

2 Remote sensing principle and challenges

The ALH retrieval from EPIC takes advantage of the absorption features of molecular O2 in the A band at 759–771 nm and the B band at 686–695 nm. The electronic transitions from the ground state to two excited states of an O2 molecule, coupled with vibrational–rotational transitions, produce O2 absorption lines in the A and B bands with distinct spectral variability (Fig. 1b–c). As a result, the large variability of atmospheric opacity within the O2 A and B bands leads to a wide range of penetration depths of solar radiation (Ding et al., 2016). The spectroscopic characteristics in the O2 bands are related to how scattered light from aerosol particles interacts with O2 absorption at different altitudes. Furthermore, since O2 molecules are well mixed in the atmospheric column, the altitude dependence of O2 absorption can provide information on the path length of light scattered by aerosol particles and is thus related to the height of aerosol layers. The premise of this retrieval dates back to Hanel (1961) and Yamamoto and Wark (1961), who estimated cloud-top pressure based on the amount of absorption by molecular CO2 and O2 above the cloud layer. More recently, a number of satellite sensors have been designed to capture the O2 absorption feature to provide cloud-top height and ALH; among those are GOME, SCIAMACHY, MERIS, POLDER, DSCOVR, and most recently the TROPOspheric Monitoring Instrument (TROPOMI; Xu et al., 2018a, and references therein).

Figure 2Schematic diagram of the physical principle of the remote sensing of aerosol layer height in the O2 absorption bands. Shown is the same aerosol layer located at two different altitudes in the atmosphere. Due to the scattering by aerosol particles, light scattered by a higher aerosol layer (b) would travel a shorter path length to reach the satellite sensor than those scattered by the lower-altitude aerosol (a), resulting in less absorption by O2 and thus a higher radiance value detected by the satellite (i.e., IH>IL).


Figure 2 illustrates the physical principle for sensing ALH using the O2 absorption spectroscopic approach, which relies on the fact that a scattering aerosol layer can scatter sunlight back to space, shortening the path length of light traveling in the atmosphere and reducing the chance of that light being absorbed by O2 molecules. As a result, an elevated scattering layer enhances the top-of-atmosphere (TOA) reflectance within the O2 absorption bands as detected by a satellite. In other words, for a given aerosol layer of fixed AOD placed at different altitudes, the higher the altitude, the larger the TOA reflectance, i.e., IH>IL in Fig. 2. This relationship contrasts with that of the reference continuum band, wherein the TOA reflectance is not sensitive to the vertical location of the aerosol layer but depends only on the column AOD, i.e., IHIL for a given AOD. Based on this principle, the ratio of TOA reflectance between in-band and the continuum band, or the DOAS ratio, provides a practical way to infer the ALH (e.g., Dubuisson et al., 2009b; Duforêt et al., 2007; Xu et al., 2017).

Figure 3Sensitivity of DOAS ratios of TOA reflectance around the O2 A and B bands to smoke ALH as simulated by UNL-VRTM for different AOD and surface reflectance (As) values. Simulations were performed for a typical biomass burning aerosol model as observed by EPIC for the geometry of θ0,θ,Δϕ=[42, 37, 165], where θ0 and θ are solar and viewing zenith angles, and Δϕ the relative azimuth angle.


To build the links between DOAS ratios and ALH, we simulate TOA reflectance as observed by EPIC measurements with the state-of-the-art Unified Linearized Vector Radiative Transfer Model (UNL-VRTM; described in Sect. 3.5). Figure 3 shows simulations of the relationship between smoke ALH (y axis) and DOAS ratios (x axis) for EPIC channels in O2 A and B bands and for various AOD values and surface reflectance (As, indicated by different colors). Those simulations were performed for a typical biomass burning aerosol model as observed by EPIC for the geometry of [θ0,θ,Δϕ]=[42, 37, 165], where θ0 and θ are solar and viewing zenith angles, and Δϕ the relative azimuth angle between the Sun and the viewer. As seen from Fig. 3, the DOAS ratios in general increase with the rise of ALH. Meanwhile, the relationship is strongly dependent on the aerosol loading and surface reflectivity.

The sensitivity of DOAS to ALH is enhanced for lower surface reflectance and larger AOD. Conversely, the sensitivity decreases rapidly for smaller AOD (e.g., AOD=0.1 in Fig. 3a, d) or as As increases, a condition under which it is difficult to discriminate between the aerosol scattering contribution to TOA reflectance and the contribution from the surface. Therefore, it is challenging to obtain ALH information for low aerosol loading conditions and over bright surfaces. This is consistent with the findings of many previous information content studies for determining ALH from O2 A and B band measurements (Ding et al., 2016; Dubuisson et al., 2009b). Duforêt et al. (2007) showed that the centroid altitude of a single aerosol layer over a dark surface can be retrieved from Polarization and Directionality of the Earth's Reflectances (POLDER) measurements with an error of less than 1 km when AOD is over 0.2. Similarly, Dubuisson et al. (2009b) suggested AOD>0.3 and As<0.06 in order to achieve a retrieval accuracy of 0.5 km from POLDER.

With the above retrieval principles and challenges in mind, for our algorithm development we need to determine over what underlying Earth surfaces EPIC measurements have sufficient information to allow for a practical retrieval of ALH. First and foremost, the surface reflectance must be low enough in one of the O2 A and B bands to provide an unambiguous signal from the aerosol. According to Fig. 3, the DOAS–ALH sensitivity becomes substantial when As<0.1 in moderate aerosol loading conditions (i.e., τ680=0.4). Hence, As=0.1 is set as the upper threshold, above which an ALH retrieval will not be attempted, despite the sensitivity increases in high aerosol loading (e.g., τ680=1.0) for higher As.

Figure 4Surface reflectance spectra (from the ASTER spectral library; Baldridge et al., 2009) for various surface types in the visible to NIR range, with selected EPIC spectral bands within this range shown in gray.


Figure 4 shows the magnitude of surface reflectance from the ASTER spectral library for typical Earth surface types (Baldridge et al., 2009). The water surface is dark, with reflectance lower than 0.03 in both the O2 A and B bands. Over land, green vegetation is the only surface type that offers a reflectance below 0.1 in the O2 B band (As=0.1 is indicated by the red dotted line). In fact, As in the O2 A band is considerably above 0.1 for all land types considered in Fig. 4. These findings provide important guidance on the ALH algorithm design. Specifically, EPIC can retrieve ALH only over water and vegetated surfaces. In particular, ALH retrieval over vegetated surface would enlist DOAS ratios in the O2 B band. Additionally, retrieval of ALH requires a sufficiently high aerosol loading. We set AOD thresholds for retrieving ALH over both the water and vegetation per the algorithm overview to follow.

Figure 5Flowchart of the AOD–ALH retrieval algorithm.


3 EPIC aerosol layer height retrieval algorithm

3.1 Algorithm overview

Figure 5 illustrates the processing of EPIC data and the ALH retrieval procedure. Briefly, the retrieval algorithm entails the following steps.

  1. Calculate TOA reflectance in six EPIC visible and NIR bands (443, 551, 680, 688, 764, and 780 nm) from the calibrated EPIC level 1B digital data.

  2. Identify EPIC pixels that are suitable for aerosol height retrieval. Through various tests, this step screens out pixels having clouds, over-water sun glints, and bright land surfaces, which are performed separately for water and land pixels. Surface pressure comes from MERRA-2 reanalysis data and we determine surface reflectance in EPIC bands using GOME-2 and Moderate-resolution Imaging Spectroradiometer (MODIS) surface products.

  3. Aggregate the original EPIC pixels into a box of 3×3 individual pixels, an area with a size of about 24 km at nadir. In many cases, not all pixels within a box are suitable for retrieval (i.e., cloud, glint, and bright land). If the number of available pixels within a box is not less than four (of the total of nine), calculate mean values of TOA reflectance, satellite geometries, and surface reflectance for the available pixels. Otherwise, do not conduct an aerosol retrieval for the box.

  4. Invert the aggregated EPIC observations using pre-calculated lookup tables to obtain smoke ALH and AOD. The inversion uses a flexible spectral fitting strategy that considers the specific surface type.

While the retrieval procedure is based on our algorithm for retrieving dust ALH over ocean from EPIC measurements (Xu et al., 2017), it was upgraded in several ways. First, the algorithm is extended to retrieve ALH over vegetated land surface. The capability of inferring ALH over vegetation is predicated on the O2 B band, where the surface reflectance is sufficiently low. Accordingly, the new algorithm uses methods separated for land and water surfaces in determining surface reflectance and screening cloud. Second, the algorithm incorporates a smoke aerosol optical model in order to retrieve biomass burning smoke ALH. Third, rather than aggregating EPIC data from satellite pixels to regular latitude–longitude grids, the new algorithm retrieves over aggregated boxes, each consisting of an array of 3×3 EPIC pixels. This change is based on the consideration that the geolocation of EPIC data has made significant improvements in the new version of level 1B data (Geogdzhayev and Marshak, 2018). Lastly, to obtain AOD and ALH the new algorithm employs a flexible spectral fitting strategy by considering the underlying surface reflectivity. In the following procedural subsections, we describe these changes in full detail.

3.2 Obtaining EPIC TOA reflectance

For this algorithm development, we used the EPIC level 1B (L1B) version 02 imagery data, available from NASA's Atmospheric Science Data Center (ASDC) at (last access: 26 November 2018). Various prelaunch and postlaunch calibrations were applied to the L1B EPIC data. EPIC visible and near-IR (NIR) channels at 443, 551, 680, and 780 nm were cross-calibrated with independent LEO satellite instruments, including the MODIS onboard the Terra and Aqua satellites (Geogdzhayev and Marshak, 2018). The two O2 absorption channels (688 and 764 nm) were calibrated using lunar surface reflectivity with EPIC lunar observations at the time of full moon as seen from Earth (Ohtake et al., 2010). The image data in EPIC L1B products are in digital units of counts per second that are converted into reflectance for each visible and NIR channel using calibration factors provided at ASDC online: (last access: 26 November 2018). The TOA reflectance values are calculated as

(1) R λ = K ( λ ) C ( λ ) μ 0 ,

where C(λ) is the EPIC measured signal in units of counts per second at the wavelength of λ, K(λ) is a calibration factor, and μ0 is the cosine of solar zenith angle, θ0. Additionally, pixels that are far from nadir and strongly distorted by the Earth ellipsoid (e.g., view zenith angle θ>70 or solar zenith angle θ0>70) are excluded.

3.3 Determining surface albedo and pressure

Past studies show that an accurate ALH retrieval critically depends on the appropriate assumption of surface reflectance (Corradini and Cervino, 2006; Ding et al., 2016; Dubuisson et al., 2009a). As discussed in the preceding section, our retrievals are confined to water and dark vegetated surfaces. Land and water mask information for each EPIC pixel was determined using the GSHHG coastline database (Wessel and Smith, 1996). Following Xu et al. (2017), surface reflectance over water surface was obtained from the GOME-2 surface Lambert-equivalent reflectivity (LER) database (Koelemeijer et al., 2003; Tilstra et al., 2017). GOME-2 LER products contain spectral LER albedo at 21 channels that are 1 nm wide, available globally for each month of the year at 0.25× 0.25 grid spacing (25×25km cells). From these data, we selected LER albedos at 440, 555, 670, 758, and 780 nm to represent surface reflectance for the nearest EPIC bands.

The Moderate-resolution Imaging Spectroradiometer (MODIS) BRDF/Albedo product (MCD43) provides parameters that can be used to determine surface reflectance over land in the first seven MODIS channels (Lucht et al., 2000; Schaaf et al., 2002). Those parameters, including surface bidirectional reflectance distribution function (BRDF) and albedos, were inverted from the atmospherically corrected (i.e., Rayleigh and aerosol components removed) MODIS reflectance observations from both the Terra and Aqua satellites over a 16 d period. We used black-sky albedo and white-sky albedo compiled in the level 3 daily Climate Modeling Grid (CMG) Albedo Product (MCD43C3) at a spatial resolution of 0.1× 0.1. Here, black-sky albedo is the directional hemispherical reflectance if the illumination comes only from the Sun at solar noon. White-sky albedo is the bi-hemispherical reflectance under isotropic skylight illumination. The actual bi-hemispherical reflectance can be computed from white-sky and black-sky albedos via a weighting coefficient that depends primarily on solar angle and columnar optical depth (Kokhanovsky et al., 2005; Lewis and Barnsley, 1994; Schaaf et al., 2002). Following Kokhanovsky et al. (2005), we calculated the Lambertian surface albedo in the 469, 555, 645, and 858 nm MODIS bands by weighting the contributions of white-sky and black-sky albedos. It should be noted that the effect of non-Lambertian surface reflection may bias the ALH retrieval because uncertainty in surface reflectance can substantially affect the ALH retrieval accuracy (see Appendix A). Nevertheless, this type of impact could be limited as EPIC's Earth observations are confined within an almost constant viewing geometry with scattering angles between 165 and 178. Further studies are needed to examine the detailed impacts, which will be one of our future efforts.

The differences in the spectral position and width of corresponding EPIC and MODIS channels may result in discrepancies. To compensate, we adjust MODIS reflectance values to equivalent EPIC bands. The adjustment factors, in the form of linear regression coefficients, were determined from analyzing USGS (United States Geological Survey) hyperspectral data for green vegetation. Figure 6a shows the selected reflectance spectra of 47 vegetation samples from the USGS Spectral Library version 7 (Kokaly et al., 2017) and the spectral locations of EPIC (blue) and MODIS (red dotted) bands. The linear regression of reflectance values in each EPIC band and corresponding MODIS band is illustrated as scatterplots in Fig. 6c–h. The regression slope and offset, as well as coefficient of determination (R2), are listed in Fig. 6b. As shown, the reflectances between EPIC bands and their closest MODIS bands are highly correlated (R2 from 0.96 to 1.0), and therefore the best-fitting equations in Fig. 6b can be used to empirically derive EPIC surface reflectance from MODIS surface reflectance products.

Figure 6Determination of surface reflectance from MODIS surface albedo products. (a) USGS reflectance spectra for different vegetation samples and (b) statistics for regression (i.e., slope, offset, and coefficient of determination, R2) of reflectance between different EPIC bands, marked as red dotted lines in (a), and the respective closest MODIS bands, marked as blue lines in (a), according to the spectral library. (c–f) Scatterplots of reflectance in each EPIC band versus reflectance in the corresponding MODIS bands. Linear regression lines are shown in red.


The retrieval also requires auxiliary meteorological information to realistically characterize O2 absorption properties. Surface pressure strongly affects the reflectance in O2 absorption bands as observed from space because it is a direct proxy for the columnar loading of air (and thus oxygen) molecules. Here, we enlisted surface pressure information from the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) datasets (Gelaro et al., 2017). MERRA-2's 1-hourly surface pressure at 0.5 by 0.675 grids were interpolated to the location and scan time of each EPIC pixel. In addition, the atmospheric temperature–pressure profile also impacts the width and strength of O2 absorption lines. However, such an influence on the radiative transfer is negligible for EPIC's 1 to 2 nm wide bands. In this study, our algorithm employs a standard temperature–pressure profile representing the midlatitude summer atmosphere (McClatchey et al., 1972).

3.4 Masking clouds, sun glints, and bright land surface

After acquiring EPIC TOA reflectance and surface reflectance, the algorithm conducts a masking exercise to select clear-sky pixels suitable for aerosol retrieval. This includes the screenings of clouds, sun glint over water surfaces, and highly reflective land surfaces. Over both land and water, the cloud mask combines a set of brightness and homogeneity tests following Martins et al. (2002). The brightness tests screen out cloud pixels for which EPIC TOA reflectance in two or three visible bands (443 and 680 nm over land; 443, 680, and 780 nm over water) exceeds prescribed thresholds. The homogeneity test, on the other hand, identifies cloud pixels by evaluating the standard deviation of TOA reflectance at 443 and 551 nm within 3×3 neighboring pixels. Still, thin and small subpixel clouds can be missed in the relatively coarse resolution 8 km EPIC pixels, leading to an overestimation of AOD and an underestimation of ALH. Therefore, contamination by small-scale clouds is one of the known issues of retrieval quality, and coupling higher-resolution cloud mask information from geostationary sensors may help to overcome this issue.

In addition, EPIC pixels affected by sun glint and highly reflective land surface are also removed from consideration by the mask. We identify the over-water glint area as those pixels having a glint angle smaller than 30 (Levy et al., 2013). Highly reflective land surface pixels are identified using MODIS land surface products. Any EPIC pixel having a normalized difference vegetation index (NDVI) value below 0.2 or having a 680 nm surface reflectance larger than 0.1 is flagged as being a bright surface and excluded from the retrievals.

3.5 Lookup tables

Following Xu et al. (2017), the revised algorithm retrieves smoke ALH and AOD from EPIC measurements via a set of lookup tables (LUTs) constructed using the UNL-VRTM. The UNL-VRTM (, last access: 15 June 2019) is a radiative transfer test bed developed specifically for atmospheric remote sensing (Wang et al., 2014). By integrating the VLIDORT vector radiative transfer model (Spurr, 2006) and particulate scattering codes (Spurr et al., 2012), UNL-VRTM can perform simulations for two or more sets of aerosol microphysical properties, typically with aerosols in one fine mode and one coarse mode. It also incorporates HITRAN spectroscopic gaseous absorption with up to 22 trace gases (Rothman et al., 2013), allowing for accurate hyperspectral simulations of remote sensing observations (Xu et al., 2018b).

The LUTs for the current retrieval consist of simulated EPIC TOA reflectance at selected spectral bands for a set of AOD and ALH values under various atmospheric and observation scenarios (e.g., Sun–Earth–sensor geometry, surface reflectance, and surface pressure) as shown in Table 1. To build the LUTs, hyperspectral monochromatic radiances were simulated using UNL-VRTM and convolved with EPIC instrumental filter response functions for the selected six bands. Table 2 provides the spectral range, interval, and resolution (in terms of full-with at half maximum; FWHM) of the UNL-VRTM simulations.

Table 1Adopted parameters for generating LUTs.

a θ0 indicates the solar zenith angle, θ the satellite viewing zenith angle, and Δφ the relative azimuth angle.

Download Print Version | Download XLSX

Table 2Spectral settings of UNL-VRTM for constructing the LUTs.

Download Print Version | Download XLSX

The UNL-VRTM simulations assume a Gaussian-like aerosol profile characterized by a centroid altitude and a half-width parameter (Spurr and Christi, 2014; Xu et al., 2017). The centroid altitude is taken to represent ALH. A half-width of 1 km is assumed following Xu et al. (2017). This value was also used to derive AOD from UV observations by the Total Ozone Mapping Spectrometer (TOMS) and the Ozone Monitoring Instrument (OMI) (Torres et al., 1998). For this study, we implemented smoke optical properties calculated based on the Lorenz–Mie scattering theory. The smoke refractive index was assumed to be 1.5–0.012i, following Dubovik et al. (2002), and its spectral dependence was neglected. We assumed smoke particles followed a bi-lognormal size distribution, as adopted for the MODIS dark-target aerosol algorithm (Levy et al., 2013; Remer et al., 2005). Specifically, fine-mode volume median radius and standard deviation were prescribed as 0.14+0.01τ680µm and 0.44, respectively. The coarse-mode counterparts are 0.14+0.01τ680µm and 0.80, respectively. The volume ratio between fine and coarse modes is (0.01+0.3τ680)/(0.01+0.09τ680). Here, τ680 is the total AOD at 680 nm.

3.6 Strategy of flexible spectral fitting

The retrieval procedure involves two steps over both water and land surfaces and, at the same time, incorporates flexible spectral fitting that accounts for the specifics of surface reflectivity. First, the TOA reflectance in the EPIC “atmospheric window” channels are matched with LUTs to determine AOD, since the TOA reflectance does not depend on ALH in these channels. Second, based on this estimated AOD, the DOAS ratios around the O2 bands are fitted to estimate ALH. In each step, least-squared fitting is applied, and fitting residuals are reported.

In general, retrieving aerosol information from reflected solar radiation in visible to NIR wavelengths over land is more challenging than over ocean, since the satellite signal tends to be dominated by surface contributions over land. It is thus often the case that, for the same satellite instrument, separate over-land and over-ocean algorithms are developed with different strategies for characterizing surface reflectivity and band selection for fitting, e.g., MODIS aerosol algorithms (Remer et al., 2005). Similarly, our retrieval of smoke ALH from EPIC uses different band combinations over land versus over ocean, adjusted by surface type and the spectral signature of smoke aerosol in the TOA reflectance.

In general, reflectance of water surfaces does not exceed 0.03 across the entire visible to NIR spectrum (Fig. 4). Therefore, AOD can be determined by fitting TOA reflectance in all four atmospheric window channels at 443, 551, 680, and 780 nm. In contrast, the 780 nm band is excluded for spectral fitting over vegetated land surface because of the high chlorophyll reflectance. Similarly, in the fitting of DOAS ratios, different weights are given for the ratios in the O2 A and B bands for different surfaces, adjusted by the sensitivity of the DOAS ratio to ALH (e.g., Fig. 3). The weighting coefficients are 0.4 for R688∕R680 and 0.6 for R764∕R780 over water surface, whereas these values are 0.9 and 0.1 over vegetated surface.

Figure 7Demonstration of ALH and AOD retrievals for three EPIC scenes acquired on 25 August 2017 around UTC times of 16:10, 17:26, and 18:32. (a) Red–green–blue (RGB) composites of EPIC 443, 551, and 680 nm data. The gold line indicates the CALIOP sub-orbital track with an overpass time of 19:05 UTC. (b) Retrieved smoke ALH when 680 nm AOD is larger than 0.2. (c) Retrieved 680 nm AOD. Small circles on AOD maps represent AOD values observed at corresponding AERONET sites within 2 h of the EPIC scan time.


Figure 8Same as Fig. 7 but for three EPIC scenes captured on 26 August 2017 at UTC times of 16:11, 17:17, and 18:22. The CALIOP overpass was at 18:15 UTC.


4 Retrieval demonstration

We apply our algorithm to six EPIC scenes over the Hudson Bay–Great Lakes area obtained during 25–26 August 2017, with three consecutive scenes considered on each day. We chose these cases primarily because smoke plumes prevailed over both water and vegetated surface, and there were available aerosol vertical profile measurements from CALIOP that could be used to validate the ALH retrievals. As shown in the EPIC RGB images (Figs. 7a and 8a), smoke aerosols appeared as diffuse plumes emitted from wildfires in western Canada and transported over the western and southern portions of Hudson Bay (Peterson et al., 2018). The retrieved smoke ALH is shown in Figs. 7b and 8b, and retrieved 680 nm AOD in Figs. 7c and 8c. It is noted that data gaps in the AOD maps represent screened-out bright pixels due to either cloud or bright land surface; ALH retrievals were only available when AOD was larger than 0.2.

Obvious spatial variations are noted in retrieved smoke ALH and AOD. On 25 August, smoke plumes had AOD values ranging from 0.1 to 0.45, with higher loading found at downwind regions in the south (Fig. 7c). An ALH of 4–5 km was found over Hudson Bay, whereas the smoke altitudes decreased to 2–4 km over land off the bay's western and southern shores (Fig. 7b). Southward, the ALH increased rapidly to 4–6 km towards the Great Lakes. By 26 August, the smoke plumes had traveled southeast (Fig. 8c). The smoke altitudes remained at 3–5 km over the eastern part of Hudson Bay and 2–4 km over the bay's south side (Fig. 8b). Altitudes of smoke plumes over the coast of the northeastern US were higher than 5 km. Aside from spatial variations, the retrievals also revealed diurnal changes in ALH and the evolution of the smoke plumes. For instance, the ALH of smoke plumes over the Hudson Bay and the north side of the Great Lakes rose by about 0.5 km within 2 h from local morning to afternoon on 25 August (Fig. 7b).

5 Retrieval validation

For the validation of the EPIC-retrieved ALH, we used CALIOP level 2 version 4.10 aerosol extinction profiles at 5 km spatial resolution, retrieved from CALIOP observations of attenuated backscatter at 532 nm (Young and Vaughan, 2009). To facilitate quantitative comparison of aerosol altitude, we used a mean extinction height calculated from the CALIOP extinction profile, following Koffi et al. (2012):

(2) ALH CALIOP = i = 1 n β ext , i Z i i = 1 n β ext , i .

Here, βext,i is the aerosol extinction coefficient (km−1) at 532 nm at level i, and Zi is the altitude (km) of level i. Thus, ALHCALIOP represents an effective ALH weighted by the aerosol extinction signal at each level and is consistent with ALH as defined in our EPIC algorithm.

In the CALIOP level 2 products, aerosol extinction is only retrieved for the layers in which aerosols are detected, as permitted by the instrument's signal-to-noise ratio (SNR). Atmospheric layers with no aerosols detected are classified as “clear air” and assigned an aerosol extinction coefficient of 0 km−1. The detection limits are defined in terms of backscattering ratio, which depends on an aerosol lidar ratio (Vaughan et al., 2009). As such, aerosol layers that span a large altitude range frequently remain undetected (Toth et al., 2018), particularly for absorbing aerosols in the daytime. As indicated by Winker et al. (2013), the aerosol extinction threshold in a daytime CALIOP scan is 0.01–0.03 km−1 for 80 km horizontal averaging resolution and increases to 0.07 km−1 for 5 km horizontal averaging resolution. In reality, aerosols are ubiquitous throughout the troposphere, though the concentration can be very low in the free troposphere. However, excluding the clear-air layers would lead to a biased estimate of ALHCALIOP, particularly for cases of predominantly clear-air layers present below an elevated aerosol layer. To compensate for this bias, we use an exponentially decayed background aerosol extinction profile to provide a proxy for the aerosol extinction coefficients of these undetected aerosol layers within the troposphere. The background profile has a columnar AOD of about 0.07 at 532 nm according to Tomasi and Petkov (2014), who use various lidar measurements to characterize the summertime Arctic atmosphere. Though aerosol extinction within the undetected aerosol layers by no means follows the background extinction profile identically, adding such a background aerosol extinction profile could, to the first order, correct the potential bias in ALHCALIOP. Therefore, we consider ALHCALIOP with and without the imposed background aerosol in the following validation analysis.

Figure 9Comparison of ALH retrieved from EPIC and CALIOP level 2 aerosol extinction profile in both “curtain” view (a, c) and “scatterplot” view (b, d) for 2 d overpasses. The CALIOP orbital tracks are marked on EPIC RGB images in Figs. 7–8. Different EPIC scan times are marked with different line and symbol colors in the comparison. The error bar for EPIC ALH represents the standard deviation for an array of 3×3 24 km retrieval pixels, while that for CALIOP ALH represents the standard deviation of over five adjacent CALIOP 5 km footprints. Also shown in panels (a, c) is the CALIOP* ALH (thin black curve) calculated with a background aerosol profile imposed for undetected aerosol layers (see text for detail), whereas the bold black curve represents the CALIOP ALH without considering background aerosol.


Figure 9 presents the evaluation of our ALH retrieval against CALIOP observations, showing that the ALH retrieval captures the overall spatial variability of ALH as seen by CALIOP. Quantitatively, 67 % and 59 % of the retrievals are within 0.5 km of difference from the counterparts of ALHCALIOP on 25 and 26 August, respectively. Considering all EPIC–CALIOP ALH pairs, over 65 % of ALH retrievals are within an uncertainty envelope of 0.5 km (Fig. 10a–b), and EPIC ALH has a mean bias of −0.13km and a root mean square error (RMSE) of 0.58 km, capturing 52 % of the variation of ALHCALIOP (Fig. 10a). If background aerosol is imposed in the CALIOP ALH calculation, our retrieved EPIC ALH is found to have a bias of 0.23 km and an RMSE of 0.57 km (Fig. 10b).

Figure 10Comparison of EPIC ALH and AOD with the corresponding CALIOP and AERONET measurements. (a) Scatterplot of EPIC ALH versus CALIOP ALH by including all EPIC–CALIOP pairs shown in Fig. 9. The color of each scatter point represents the EPIC 680 nm AOD value for the same EPIC pixels. (b) Same as panel (a) but for CALIOP ALH calculated with a background aerosol profile imposed for undetected aerosol layers. (c) Scatterplot of EPIC 680 nm AOD versus AERONET 675 nm AOD, collocated at the nine AERONET sites listed in Table 3. The dotted lines in both the scatterplots represent error envelops, i.e., ±0.5km for ALH and ±(0.05+10%) for AOD. Also annotated are the one-to-one line (solid black line), linear regression fit (red line), number of scatter points (N), coefficient of determination (R2), significance level (P), linear regression equation, bias, percentage of scatter points within the error envelop (EE %), and root mean square error (RMSE).


Table 3AERONET sites selected for AOD validation.

a NAOD is the number of collocated AERONET AOD values.

Download Print Version | Download XLSX

Figure 11EPIC UV aerosol index (UVAI) for the three scenes on 25 August 2017 shown in Fig. 7. Here the UVAI was obtained from EPIC level 2 UV aerosol products.


The EPIC 680 nm AOD retrievals during the two days were compared against 675 nm AOD observations at nine AERONET sites (Table 3). The collocation method follows Ichoku et al. (2002) but was modified to compare EPIC AOD retrievals over 3×3 pixels centered at the AERONET sites with 1 h AERONET AOD observations around the EPIC scan time. Collocated AERONET AOD values at each site are shown as circles in Figs. 7c and 8c. A comparison of EPIC and AERONET AODs is shown in Fig. 10c. The collocated AOD pairs, though with limited data samplings, have over 77 % falling in an uncertainty envelope of ±(0.05+0.1AOD) with a coefficient of determination (R2) of 0.54. The EPIC AOD shows a positive bias of 0.03 and an RMSE of 0.05, which was dominated by four subsets of EPIC AODs and was likely caused by cloud contamination.

Figure 12Scatterplots of EPIC UVAI (Fig. 11) versus current retrievals of ALH (Fig. 7) on 25 August 2017 (17:26 UTC) under different AOD values as indicated in the title of each panel. Also annotated are the linear regression fit, number of scatter points (N), coefficient of determination (R2), and significance level (P).


6 Implication to the retrieval of UV absorbing aerosol properties

One important implication of this study is that the retrieved ALH can provide complementary height information for determining absorbing aerosol properties from EPIC's UV bands. EPIC also measures backscattered UV radiances at 340 and 388 nm, which was designed to detect and retrieve the optical properties of UV-absorbing aerosols like mineral dust and smoke by using the UV aerosol index (UVAI) (Marshak et al., 2018). UVAI quantifies the difference in spectral dependence between measured and calculated near-UV radiances assuming a purely molecular atmosphere (Torres et al., 1998). Physically, inferring aerosol properties from those UV bands requires the characterization of aerosol height, since UV radiance is sensitive to aerosol vertical distribution. For example, Jeong and Hsu (2008) retrieved aerosol single-scattering albedo (SSA) from OMI radiance via synergic use of AOD from MODIS and ALH from CALIOP. Currently, EPIC's UV aerosol algorithm utilizes a climatological aerosol height dataset derived from CALIOP observations (Torres et al., 2013). However, these static climatological data may fail to capture the dynamic variation of aerosol height and thus induce uncertainties in the UV-retrieved aerosol properties. With the aerosol height and loading available from EPIC's O2 A and B bands, these closures are now possible with measurements from a single instrument. Such synergy can be also applied to the TROPOspheric Monitoring Instrument (TROPOMI) and the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) satellite, both of which obtain hyperspectral measurements from the UV to the NIR covering the O2 A and B bands (Omar et al., 2018; Veefkind et al., 2012).

Here, we compare our retrieved ALHs with EPIC UVAI to illustrate the importance, as well the potential benefit, of the ALH retrievals to EPIC's UV aerosol products. The EPIC UVAI data are publicly available at (last access: 26 November 2018). Figure 11 shows the UVAI of the same EPIC scenes on 25 August. These smoke plumes are highlighted by large values of UVAI, which are in contrast to clouds that typically exhibit a UVAI close to zero (Torres et al., 1998). Since UVAI is a function of ALH, AOD, and SSA, its correlation with ALH varies with AOD. As shown in Fig. 12, the sensitivity of UVAI to ALH, as well as the correlation between them, increases as AOD increases. In particular, high AOD values (e.g., over 0.4) may result in UVAI values from less than 1 to about 4, depending on the ALH (Fig. 12, bottom right panel). Therefore, the use of ALH derived from the EPIC O2 bands is expected to improve EPIC UV aerosol retrievals.

7 Summary and discussion

We extend our retrieval algorithm for retrieving over-water dust ALH from EPIC measurements (Xu et al., 2017) to biomass burning smoke aerosols over both water and vegetated land surfaces. The new algorithm uses Earth-reflected radiances in EPIC's visible and NIR bands and incorporates flexible spectral fitting that accounts for the specifics of vegetation and water surface reflectivity. The fitting procedure first determines AOD using four EPIC atmospheric window bands (443, 551, 680, and 780 nm) and then uses the DOAS ratios formulated in the O2 A and B bands (688 and 764 nm) to derive the ALH that represents an optical centroid altitude. ALH retrieval over vegetated surface primarily takes advantage of the measurements in the O2 B band, where surface reflectance is sufficiently low to yield aerosol height information. Surface reflectance values are specified using MODIS surface products over land and GOME-2 LER products over ocean.

We applied the algorithm to six EPIC images, with three images on each day, having biomass burning plumes over the United States and Canada acquired on 25–26 August 2017. The algorithm is able to retrieve AOD and ALH multiple times daily over both water and vegetated land surfaces. The retrieved ALHs were validated against CALIOP extinction-weighted aerosol height (ALHCALIOP), showing that EPIC-retrieved ALH has an RMSE of 0.58 km and captures 52 % of the variation of ALHCALIOP, and 65 % of EPIC and CALIOP ALH pairs are within an uncertainty envelope of 0.5 km. EPIC-retrieved AODs are validated against AERONET AOD observed at nine sites, indicating an RMSE of 0.05, and over 77 % of EPIC AOD retrievals fall within the error envelope of ±(0.05+0.1AOD). In addition, by comparing the retrieved ALH and EPIC UVAI, we show that the aerosol height information retrieved by the present algorithm can potentially benefit the retrieval of aerosol properties from the EPIC UV bands.

Over 3 years of data recorded thus far by EPIC (July 2015–present) offer an opportunity to characterize aerosol height from diurnal to seasonal scales. Such datasets with well-characterized uncertainties are valuable for evaluating aerosol vertical distribution simulated by climate models. A follow-on study, based on the current work and Xu et al. (2017), will examine the long-term global retrieval of ALH for dust and smoke aerosols. Challenges for such a global retrieval include the heterogeneity of aerosol species, whereas our efforts have been mainly focused on single-species (dust or smoke) aerosol assumptions. In fact, mixtures of aerosol are more commonly found in nature. Furthermore, the currently available EPIC UVAI products can help identify elevated absorbing aerosols (mainly dust and smoke), which can then be used to determine the retrieval targets on any EPIC image.

Data availability

The data presented in this paper can be obtained by emailing the corresponding authors (Xiaoguang Xu and Jun Wang).

Appendix A: Sensitivity and error analysis

In this Appendix, we investigate the sensitivity of EPIC measurements in the O2 A and B bands to the aerosol vertical profile shape, aerosol optical properties, and surface reflectance assumed in the retrieval algorithm. Then, we estimate potential ALH retrieval errors due to the uncertainties in these parameters. For this case, ALH is derived from two DOAS ratios (ρ) in the O2 A and B bands. Mathematically, the retrieval error (ϵ^) of ALH can be estimated using the optimal estimation approach (Xu and Wang, 2015):

(A1) ϵ ^ - 2 = K T S ϵ - 1 K ,

where K is the Jacobian matrix of ρ with respect to ALH, and Sϵ is the observation error covariance matrix for ρ. In this study, we assume no error correlation between the two DOAS ratios. That is, Sϵ is a diagonal matrix with its elements equal to the error variance for ρ. Observation error variance matrix Sϵ can be characterized as a sum of two terms:

(A2) S ϵ = S ρ + S m .

Here, Sρ is the error covariance matrix describing the uncertainty for EPIC measurements. Sm is the covariance matrix for forward model errors caused by inaccurate model assumptions and uncertainties in model parameterizations, and Sm can be calculated by

(A3) S m = K b T S b K b ,

where Sb is the error covariance matrix describing uncertainties of the vector of model parameters b, and Kb is Jacobian matrix of ρ with respect to b. In the following analysis, we consider four parameters for the vector b, i.e., the half-width parameter (σH) defining the quasi-Gaussian aerosol vertical profile, the 680 nm AOD (τ680), aerosol single-scattering albedo (SSA), and surface reflectance As.

Figure A1UNL-VRTM-simulated DOAS ratios ρ in EPIC O2 A and B band (a) and their Jacobian gradients with respect to ALH (b), half-width parameter σH for the quasi-Gaussian aerosol vertical profile (c), 680 nm AOD τ680 (d), SSA (e), and surface reflectance As (f). The y axis represents aerosols being present at various ALH values. Three colors indicate ρ and its Jacobians for the O2 A band with As=0.05 (blue), the O2 A band with As=0.30 (green), and the O2 B band with As=0.05 (red). Two AOD loadings are indicated by the solid (τ680=0.4) and dotted (τ680=1.0) lines, respectively. Simulations are performed with the same aerosol model and observation geometry in Fig. 3.


Figure A2Estimated ALH retrieval errors owing to EPIC measurement errors and uncertainties in aerosol and surface assumptions. Panels (a) and (c) are for a surface with As=0.05 reflectance in both the O2 A and B bands (close to a water surface), and panels (b) and (d) are for As=0.30 in the A band and As=0.05 in the B band (close to a vegetation surface). (a–b) ALH retrieval errors due to 2 % EPIC measurement error alone (black curves) and with model uncertainty from all four parameters added (red curves). (c–d) The error in retrieved ALH due to uncertainty from each of the parameters: 0.5 km in σH (blue), 0.05 in τ680 (orange), 0.03 in SSA (green), and 0.01 in As (purple). In all panels, the solid and dotted curves represent 0.4 and 1.0 τ680, respectively.


Figure A1 presents the UNL-VRTM-simulated DOAS ratios ρ in the EPIC O2 A and B band (panel a) and their Jacobian gradients with respect to ALH (panel b) and with respect to the abovementioned four model parameters (panels c–f) for the vector b for various surface and aerosol loading conditions. Figure A2 shows the estimated ALH retrieval errors owing to EPIC observation error and uncertainties in the assumed model parameters for water (with As=0.05 in both O2 A and B bands) and vegetation (with As=0.30 in the A band and 0.05 in the B band) surfaces. The findings from these Jacobian and error analysis results can be summarized as follows.

  • The sensitivity of ρ with respect to ALH (Fig. A1b) is weak for aerosols confined in the lower atmosphere (ALH below 1.5 km). The sensitivity increases rapidly with the increase in ALH, peaks for ALH about 4 km, and then deceases slightly for higher elevated aerosols. The sensitivity in the O2 A band is stronger than that in the O2 B band, and it is stronger for low surface reflectance and high aerosol loading. By considering 2 % uncertainty for the EPIC DOAS ratios (Geogdzhayev and Marshak, 2018), the ALH retrieval error (ϵ^0) is less than 0.5 km for the water surface and is less than 0.75 km for the vegetated surface if ALH is over 1.5 km (black curves in Fig. A2a–b). This error characterization is consistent with the retrieval error (0.57–0.58 km) as estimated through the validation with CALIOP data.

  • The DOAS ratios ρ have a negative sensitivity to σH for elevated aerosol, and the sensitivity turns to positive for ALH below 1.5 km (Fig. A1c). An error of 0.5 km σH may cause up to 0.3 km of retrieval error for ALH (blue curves in Fig. A2c–d).

  • The DOAS ratios ρ have a positive sensitivity to τ680, and the sensitivity increases as ALH increases (Fig. A1d). An error of 0.05 in τ680 can lead to a retrieval error from 0 to 0.6 km for ALH, depending on the aerosol loading and surface type (orange curves in Fig. A2c–d).

  • DOAS ratios are sensitive to SSA to some degree, especially for large AOD values (Fig. A1e). However, the sensitivity to SSA is much less significant than the sensitivity to AOD and surface reflectance because the reflectance at TOA depends more on surface reflectance and AOD (than SSA in a relative sense). As a result, SSA only has a marginal impact on the ALH retrieval error (green curves in Fig. A2c–d), which is consistent with findings by Sanders et al. (2015).

  • The DOAS ratio ρ has a negative sensitivity to As, and the sensitivity increases as ALH increases (Fig. A1f). An error of 0.01 in As can result in a retrieval error of 0.1–0.6 km for ALH over water surface and 0.1–0.4 km over vegetation surface, depending on the aerosol loading (purple curves in Fig. A2c–d).

  • If including possible uncertainties from these four parameters in vector b and EPIC measurement error, the estimated ALH retrieval errors (ϵ^All) are shown as red curves in Fig. A2a–b. For elevated smoke aerosols with τ680=0.4, ϵ^All ranges from 0.4 to 1.0 km for water surface and from 0.7 to 1.25 km for vegetation surface; the ϵ^All range for τ680=1.0 is reduced to 0.3–0.5 km for water and 0.4–0.6 km for vegetation surface.

Author contributions

JW and XX designed the research; XX and JW conducted the research and wrote the paper; YW contributed data analysis for validation; JZ, OT, JSR, SDM, JVM, and LAR contributed to the writing.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Holistic Analysis of Aerosol in Littoral Environments – A Multidisciplinary University Research Initiative (ACP/AMT inter-journal SI)”. It is not associated with a conference.


This research is supported in part by the NASA DSCOVR Earth Science Algorithms Program (grant no. NNX17AB05G, managed by Richard S. Eckman) and in part by the Office of Naval Research (ONR) Multidisciplinary University Research Initiatives (MURI) program under award no. N00014-16-1-2040. We acknowledge the computational support from the High-Performance Computing group at the University of Iowa and the UMBC High-Performance Computing Facility (HPCF). EPIC level 1B (L1B) version 02 imagery data are obtained from the NASA Atmospheric Science Data Center (ASDC). We also acknowledge the AERONET program for providing AOD data retrieved from sun photometers. We thank Feng Xu and two other anonymous reviewers for their insightful comments on the paper.

Financial support

This research has been supported by NASA (grant no. NNX17AB05G) and the Office of Naval Research (grant no. N00014-16-1-2040).

Review statement

This paper was edited by Sebastian Schmidt and reviewed by Feng Xu and two anonymous referees.


Babu, S. S., Moorthy, K. K., Manchanda, R. K., Sinha, P. R., Satheesh, S. K., Vajja, D. P., Srinivasan, S., and Kumar, V. H. A.: Free tropospheric black carbon aerosol measurements using high altitude balloon: Do BC layers build “their own homes” up in the atmosphere?, Geophys. Res. Lett., 38, L08803,, 2011. 

Baldridge, A. M., Hook, S. J., Grove, C. I., and Rivera, G.: The ASTER spectral library version 2.0, Remote Sens. Environ., 113, 711–715,, 2009. 

Choi, J.-O. and Chung, C. E.: Sensitivity of aerosol direct radiative forcing to aerosol vertical profile, Tellus B, 66, 24376,, 2014. 

Chowdhary, J., Cairns, B., Mishchenko, M. I., Hobbs, P. V., Cota, G. F., Redemann, J., Rutledge, K., Holben, B. N., and Russell, E.: Retrieval of Aerosol Scattering and Absorption Properties from Photopolarimetric Observations over the Ocean during the CLAMS Experiment, J. Atmos. Sci., 62, 1093–1117,, 2005. 

Corradini, S. and Cervino, M.: Aerosol extinction coefficient profile retrieval in the oxygen A-band considering multiple scattering atmosphere. Test case: SCIAMACHY nadir simulated measurements, J. Quant. Spectrosc. Ra., 97, 354–380,, 2006. 

Ding, S., Wang, J., and Xu, X.: Polarimetric remote sensing in oxygen A and B bands: sensitivity study and information content analysis for vertical profile of aerosols, Atmos. Meas. Tech., 9, 2077–2092,, 2016. 

Dubovik, O., Holben, B., Eck, T. F., Smirnov, A., Kaufman, Y. J., King, M. D., Tanré, D., and Slutsker, I.: Variability of Absorption and Optical Properties of Key Aerosol Types Observed in Worldwide Locations, J. Atmos. Sci., 59, 590–608,<0590:voaaop>;2, 2002. 

Dubuisson, P., Frouin, R., Dessailly, D., Duforêt, L., Léon, J.-F., Voss, K., and Antoine, D.: Estimating the altitude of aerosol plumes over the ocean from reflectance ratio measurements in the O2 A-band, Remote Sens. Environ., 113, 1899–1911,, 2009a. 

Dubuisson, P., Frouin, R., Dessailly, D., Duforêt, L., Léon, J.-F., Voss, K., and Antoine, D.: Estimating the altitude of aerosol plumes over the ocean from reflectance ratio measurements in the O2 A-band, Remote Sens. Environ., 113, 1899–1911,, 2009b. 

Duforêt, L., Frouin, R., and Dubuisson, P.: Importance and estimation of aerosol vertical structure in satellite ocean-color remote sensing, Appl. Optics, 46, 1107–1119,, 2007. 

Ge, C., Wang, J., and Reid, J. S.: Mesoscale modeling of smoke transport over the Southeast Asian Maritime Continent: coupling of smoke direct radiative effect below and above the low-level clouds, Atmos. Chem. Phys., 14, 159–174,, 2014. 

Gelaro, R., McCarty, W., Suárez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., Wargan, K., Coy, L., Cullather, R., Draper, C., Akella, S., Buchard, V., Conaty, A., Silva, A. M. D., Gu, W., Kim, G.-K., Koster, R., Lucchesi, R., Merkova, D., Nielsen, J. E., Partyka, G., Pawson, S., Putman, W., Rienecker, M., Schubert, S. D., Sienkiewicz, M., and Zhao, B.: The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), J. Climate, 30, 5419–5454,, 2017. 

Geogdzhayev, I. V. and Marshak, A.: Calibration of the DSCOVR EPIC visible and NIR channels using MODIS Terra and Aqua data and EPIC lunar observations, Atmos. Meas. Tech., 11, 359–368,, 2018. 

Hanel, R. A.: Determination of cloud altitude from a satellite, J. Geophys. Res., 66, 1300–1300,, 1961. 

Ichoku, C., Chu, D. A., Mattoo, S., Kaufman, Y. J., Remer, L. A., Tanré, D., Slutsker, I., and Holben, B. N.: A spatio-temporal approach for global validation and analysis of MODIS aerosol products, Geophys. Res. Lett., 29, 1616,, 2002. 

Jeong, M.-J. and Hsu, N. C.: Retrievals of aerosol single-scattering albedo and effective aerosol layer height for biomass-burning smoke: Synergy derived from “A-Train” sensors, Geophys. Res. Lett., 35, L24801,, 2008. 

Kipling, Z., Stier, P., Johnson, C. E., Mann, G. W., Bellouin, N., Bauer, S. E., Bergman, T., Chin, M., Diehl, T., Ghan, S. J., Iversen, T., Kirkevåg, A., Kokkola, H., Liu, X., Luo, G., van Noije, T., Pringle, K. J., von Salzen, K., Schulz, M., Seland, Ø., Skeie, R. B., Takemura, T., Tsigaridis, K., and Zhang, K.: What controls the vertical distribution of aerosol? Relationships between process sensitivity in HadGEM3–UKCA and inter-model variation from AeroCom Phase II, Atmos. Chem. Phys., 16, 2221–2241,, 2016. 

Koch, D. and Del Genio, A. D.: Black carbon semi-direct effects on cloud cover: review and synthesis, Atmos. Chem. Phys., 10, 7685–7696,, 2010. 

Koelemeijer, R. B. A., de Haan, J., and Stammes, P.: A database of spectral surface reflectivity in the range 335–772 nm derived from 5.5 years of GOME observations, J. Geophys. Res., 108, 4070,, 2003. 

Koffi, B., Schulz, M., Bréon, F.-M., Griesfeller, J., Winker, D., Balkanski, Y., Bauer, S., Berntsen, T., Chin, M., Collins, W. D., Dentener, F., Diehl, T., Easter, R., Ghan, S., Ginoux, P., Gong, S., Horowitz, L. W., Iversen, T., Kirkevåg, A., Koch, D., Krol, M., Myhre, G., Stier, P., and Takemura, T.: Application of the CALIOP layer product to evaluate the vertical distribution of aerosols estimated by global models: AeroCom phase I results, J. Geophys. Res.-Atmos., 117, D10201,, 2012. 

Kokaly, R. F., Clark, R. N., Swayze, G. A., Livo, K. E., Hoefen, T. M., Pearson, N. C., Wise, R. A., Benzel, W. M., Lowers, H. A., Driscoll, R. L., and Klein, A. J.: USGS Spectral Library Version 7, Reston, VA, Report 1035, 68, 2017. 

Kokhanovsky, A. A., Mayer, B., and Rozanov, V. V.: A parameterization of the diffuse transmittance and reflectance for aerosol remote sensing problems, Atmos. Res., 73, 37–43,, 2005. 

Levy, R. C., Mattoo, S., Munchak, L. A., Remer, L. A., Sayer, A. M., Patadia, F., and Hsu, N. C.: The Collection 6 MODIS aerosol products over land and ocean, Atmos. Meas. Tech., 6, 2989–3034,, 2013. 

Lewis, P. and Barnsley, M. J.: Influence of the sky radiance distribution on various formulations of the earth surface albedo, in: Proc. Conf. Phys. Meas. Sign. Remote Sens, Val d'Isere, France, 707–715, 1994. 

Lucht, W., Schaaf, C. B., and Strahler, A. H.: An algorithm for the retrieval of albedo from space using semiempirical BRDF models, Geoscience and Remote Sensing, IEEE Transactions on, 38, 977–998,, 2000. 

Maddy, E. S., DeSouza-Machado, S. G., Nalli, N. R., Barnet, C. D., Strow, L., Wolf, W. W., Xie, H., Gambacorta, A., King, T. S., Joseph, E., Morris, V., Hannon, S. E., and Schou, P.: On the effect of dust aerosols on AIRS and IASI operational level 2 products, Geophys. Res. Lett., 39, L10809,, 2012. 

Marshak, A., Herman, J., Szabo, A., Blank, K., Carn, S., Cede, A., Geogdzhayev, I., Huang, D., Huang, L.-K., Knyazikhin, Y., Kowalewski, M., Krotkov, N., Lyapustin, A., McPeters, R., Meyer, K., Torres, O., and Yang, Y.: Earth Observations from DSCOVR/EPIC Instrument, B. Am. Meteorol. Soc., 99, 1829–1850,, 2018. 

Martins, J. V., Tanré, D., Remer, L., Kaufman, Y., Mattoo, S., and Levy, R.: MODIS Cloud screening for remote sensing of aerosols over oceans using spatial variability, Geophys. Res. Lett., 29, 1619,, 2002. 

McClatchey, R. A., Fenn, R. W., Selby, J. E. A., Volz, F. E., and Garing, J. S.: Optical Properties of the Atmosphere (Third Edition), Air Force Cambridge Research Labs, Hanscom AFB, MA, 110 pp., 1972. 

Ohtake, M., Matsunaga, T., Yokota, Y., Yamamoto, S., Ogawa, Y., Morota, T., Honda, C., Haruyama, J., Kitazato, K., Takeda, H., Iwasaki, A., Nakamura, R., Hiroi, T., Kodama, S., and Otake, H.: Deriving the Absolute Reflectance of Lunar Surface Using SELENE (Kaguya) Multiband Imager Data, Space Sci. Rev., 154, 57–77,, 2010. 

Omar, A. H., Tzortziou, M., Coddington, O., and Remer, L. A.: Plankton Aerosol, Cloud, ocean Ecosystem mission: atmosphere measurements for air quality applications, J. Appl. Remote Sens., 12, 042608,, 2018. 

Peterson, D. A., Campbell, J. R., Hyer, E. J., Fromm, M. D., Kablick, G. P., Cossuth, J. H., and DeLand, M. T.: Wildfire-driven thunderstorms cause a volcano-like stratospheric injection of smoke, Climate and Atmospheric Science, 1, 30,, 2018. 

Remer, L. A., Kaufman, Y. J., Tanre, D., Mattoo, S., Chu, D. A., Martins, J. V., Li, R.-R., Ichoku, C., Levy, R. C., Kleidman, R. G., Eck, T. F., Vermote, E., and Holben, B. N.: The MODIS Aerosol Algorithm, Products, and Validation, J. Atmos. Sci., 62, 947–973,, 2005. 

Rothman, L. S., Gordon, I. E., Babikov, Y., Barbe, A., Chris Benner, D., Bernath, P. F., Birk, M., Bizzocchi, L., Boudon, V., Brown, L. R., Campargue, A., Chance, K., Cohen, E. A., Coudert, L. H., Devi, V. M., Drouin, B. J., Fayt, A., Flaud, J. M., Gamache, R. R., Harrison, J. J., Hartmann, J. M., Hill, C., Hodges, J. T., Jacquemart, D., Jolly, A., Lamouroux, J., Le Roy, R. J., Li, G., Long, D. A., Lyulin, O. M., Mackie, C. J., Massie, S. T., Mikhailenko, S., Müller, H. S. P., Naumenko, O. V., Nikitin, A. V., Orphal, J., Perevalov, V., Perrin, A., Polovtseva, E. R., Richard, C., Smith, M. A. H., Starikova, E., Sung, K., Tashkun, S., Tennyson, J., Toon, G. C., Tyuterev, V. G., and Wagner, G.: The HITRAN2012 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 130, 4–50,, 2013. 

Samset, B. H., Myhre, G., Schulz, M., Balkanski, Y., Bauer, S., Berntsen, T. K., Bian, H., Bellouin, N., Diehl, T., Easter, R. C., Ghan, S. J., Iversen, T., Kinne, S., Kirkevåg, A., Lamarque, J.-F., Lin, G., Liu, X., Penner, J. E., Seland, Ø., Skeie, R. B., Stier, P., Takemura, T., Tsigaridis, K., and Zhang, K.: Black carbon vertical profiles strongly affect its radiative forcing uncertainty, Atmos. Chem. Phys., 13, 2423–2434,, 2013. 

Sanders, A. F. J., de Haan, J. F., Sneep, M., Apituley, A., Stammes, P., Vieitez, M. O., Tilstra, L. G., Tuinder, O. N. E., Koning, C. E., and Veefkind, J. P.: Evaluation of the operational Aerosol Layer Height retrieval algorithm for Sentinel-5 Precursor: application to O2 A band observations from GOME-2A, Atmos. Meas. Tech., 8, 4947–4977,, 2015. 

Satheesh, S. K., Moorthy, K. K., Babu, S. S., Vinoj, V., and Dutt, C. B. S.: Climate implications of large warming by elevated aerosol over India, Geophys. Res. Lett., 35, L19809,, 2008. 

Schaaf, C. B., Gao, F., Strahler, A. H., Lucht, W., Li, X., Tsang, T., Strugnell, N. C., Zhang, X., Jin, Y., Muller, J.-P., Lewis, P., Barnsley, M., Hobson, P., Disney, M., Roberts, G., Dunderdale, M., Doll, C., d'Entremont, R. P., Hu, B., Liang, S., Privette, J. L., and Roy, D.: First operational BRDF, albedo nadir reflectance products from MODIS, Remote Sens. Environ., 83, 135–148,, 2002. 

Spurr, R.: VLIDORT: A linearized pseudo-spherical vector discrete ordinate radiative transfer code for forward model and retrieval studies in multilayer multiple scattering media, J. Quant. Spectrosc. Ra., 102, 316–342, 2006. 

Spurr, R. and Christi, M.: On the generation of atmospheric property Jacobians from the (V)LIDORT linearized radiative transfer models, J. Quant. Spectrosc. Ra., 142, 109–115,, 2014. 

Spurr, R., Wang, J., Zeng, J., and Mishchenko, M. I.: Linearized T-matrix and Mie scattering computations, J. Quant. Spectrosc. Ra., 113, 425–439,, 2012. 

Tilstra, L. G., Tuinder, O. N. E., Wang, P., and Stammes, P.: Surface reflectivity climatologies from UV to NIR determined from Earth observations by GOME-2 and SCIAMACHY, J. Geophys. Res.-Atmos., 122, 4084–4111,, 2017. 

Tomasi, C. and Petkov, B. H.: Calculations of relative optical air masses for various aerosol types and minor gases in Arctic and Antarctic atmospheres, J. Geophys. Res.-Atmos., 119, 1363–1385,, 2014. 

Torres, O., Bhartia, P. K., Herman, J. R., Ahmad, Z., and Gleason, J.: Derivation of aerosol properties from satellite measurements of backscattered ultraviolet radiation: Theoretical basis, J. Geophys. Res., 103, 17099–17110,, 1998. 

Torres, O., Ahn, C., and Chen, Z.: Improvements to the OMI near-UV aerosol algorithm using A-train CALIOP and AIRS observations, Atmos. Meas. Tech., 6, 3257–3270,, 2013. 

Toth, T. D., Campbell, J. R., Reid, J. S., Tackett, J. L., Vaughan, M. A., Zhang, J., and Marquis, J. W.: Minimum aerosol layer detection sensitivities and their subsequent impacts on aerosol optical thickness retrievals in CALIPSO level 2 data products, Atmos. Meas. Tech., 11, 499–514,, 2018. 

Vaughan, M. A., Powell, K. A., Winker, D. M., Hostetler, C. A., Kuehn, R. E., Hunt, W. H., Getzewich, B. J., Young, S. A., Liu, Z., and McGill, M. J.: Fully Automated Detection of Cloud and Aerosol Layers in the CALIPSO Lidar Measurements, J. Atmos. Ocean. Tech., 26, 2034–2050,, 2009. 

Veefkind, J. P., Aben, I., McMullan, K., Förster, H., de Vries, J., Otter, G., Claas, J., Eskes, H. J., de Haan, J. F., Kleipool, Q., van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P., Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P. F.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83,, 2012. 

Wang, J., Xu, X., Ding, S., Zeng, J., Spurr, R., Liu, X., Chance, K., and Mishchenko, M.: A numerical testbed for remote sensing of aerosols, and its demonstration for evaluating retrieval synergy from a geostationary satellite constellation of GEO-CAPE and GOES-R, J. Quant. Spectrosc. Ra., 146, 510–528,, 2014. 

Waquet, F., Cairns, B., Knobelspiesse, K., Chowdhary, J., Travis, L. D., Schmid, B., and Mishchenko, M. I.: Polarimetric remote sensing of aerosols over land, J. Geophys. Res., 114, D01206,, 2009. 

Wendisch, M., Hellmuth, O., Ansmann, A., Heintzenberg, J., Engelmann, R., Althausen, D., Eichler, H., Müller, D., Hu, M., Zhang, Y., and Mao, J.: Radiative and dynamic effects of absorbing aerosol particles over the Pearl River Delta, China, Atmos. Environ., 42, 6405–6416,, 2008. 

Wessel, P. and Smith, W. H. F.: A global, self-consistent, hierarchical, high-resolution shoreline database, J. Geophys. Res.-Sol. Ea., 101, 8741–8743,, 1996. 

Winker, D. M., Tackett, J. L., Getzewich, B. J., Liu, Z., Vaughan, M. A., and Rogers, R. R.: The global 3-D distribution of tropospheric aerosols as characterized by CALIOP, Atmos. Chem. Phys., 13, 3345–3361,, 2013. 

Xu, X. and Wang J.: Retrieval of aerosol microphysical properties from AERONET photopolarimetric measurements: 1. Information content analysis, J. Geophys. Res.-Atmos., 120, 7059–7078, 2015. 

Xu, X., Wang, J., Wang, Y., Zeng, J., Torres, O., Yang, Y., Marshak, A., Reid, J., and Miller, S.: Passive remote sensing of altitude and optical depth of dust plumes using the oxygen A and B bands: First results from EPIC/DSCOVR at Lagrange-1 point, Geophys. Res. Lett., 44, 7544–7554,, 2017. 

Xu, X., Wang, J., Wang, Y., and Kokhanovsky, A.: Chapter 1 – Passive Remote Sensing of Aerosol Height, in: Remote Sensing of Aerosols, Clouds, and Precipitation, Elsevier, Cambridge, MA, 1–22, 2018a. 

Xu, X., Wang, J., Zeng, J., Hou, W., Meyer, K. G., Platnick, S. E., and Wilcox, E. M.: A pilot study of shortwave spectral fingerprints of smoke aerosols above liquid clouds, J. Quant. Spectrosc. Ra., 221, 38–50,, 2018b. 

Yamamoto, G. and Wark, D. Q.: Discussion of the letter by R. A. Hanel, “Determination of cloud altitude from a satellite”, J. Geophys. Res., 66, 3596–3596,, 1961. 

Yang, Y., Marshak, A., Mao, J., Lyapustin, A., and Herman, J.: A method of retrieving cloud top height and cloud geometrical thickness with oxygen A and B bands for the Deep Space Climate Observatory (DSCOVR) mission: Radiative transfer simulations, J. Quant. Spectrosc. Ra., 122, 141–149,, 2013. 

Young, S. A. and Vaughan, M. A.: The Retrieval of Profiles of Particulate Extinction from Cloud-Aerosol Lidar Infrared Pathfinder Satellite Observations (CALIPSO) Data: Algorithm Description, J. Atmos. Ocean. Tech., 26, 1105–1119,, 2009. 

Zhang, J., Campbell, J. R., Reid, J. S., Westphal, D. L., Baker, N. L., Campbell, W. F., and Hyer, E. J.: Evaluating the impact of assimilating CALIOP-derived aerosol extinction profiles on a global mass transport model, Geophys. Res. Lett., 38, L14801,, 2011. 

Short summary
Detecting aerosol layer height from space is challenging. The traditional method relies on active sensors such as lidar that provide the detailed vertical structure of the aerosol profile but is costly with limited spatial coverage (more than 1 year is needed for global coverage). Here we developed a passive remote sensing technique that uses backscattered sunlight to retrieve smoke aerosol layer height over both water and vegetated surfaces from a sensor 1.5 million kilometers from the Earth.