**Research article**
06 Oct 2020

**Research article** | 06 Oct 2020

# Cloud-top pressure retrieval with DSCOVR EPIC oxygen A- and B-band observations

Bangsheng Yin Qilong Min Emily Morgan Yuekui Yang Alexander Marshak and Anthony B. Davis

^{1},

^{1},

^{1},

^{2},

^{2},

^{3}

**Bangsheng Yin et al.**Bangsheng Yin Qilong Min Emily Morgan Yuekui Yang Alexander Marshak and Anthony B. Davis

^{1},

^{1},

^{1},

^{2},

^{2},

^{3}

^{1}Atmospheric Sciences Research Center, University at Albany, Albany, NY, USA^{2}NASA Goddard Space Flight Center, Climate and Radiation Laboratory, Greenbelt, MD, USA^{3}Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA

^{1}Atmospheric Sciences Research Center, University at Albany, Albany, NY, USA^{2}NASA Goddard Space Flight Center, Climate and Radiation Laboratory, Greenbelt, MD, USA^{3}Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA

**Correspondence**: Qilong Min (qmin@albany.edu)

**Correspondence**: Qilong Min (qmin@albany.edu)

Received: 02 Oct 2019 – Discussion started: 04 Feb 2020 – Revised: 31 Jul 2020 – Accepted: 06 Aug 2020 – Published: 06 Oct 2020

An analytic transfer inverse model for Earth Polychromatic Imaging Camera
(EPIC) observations is proposed to retrieve the cloud-top pressure (CTP) with the
consideration of in-cloud photon penetration. In this model, an analytic equation
was developed to represent the reflection at the top of the atmosphere from above cloud, in cloud, and below cloud. The coefficients of this analytic equation
can be derived from a series of EPIC simulations under different atmospheric
conditions using a nonlinear regression algorithm. With estimated cloud
pressure thickness, the CTP can be retrieved from EPIC observation data by
solving the analytic equation. To simulate the EPIC measurements, a program
package using the double-*k* approach was developed. Compared to line-by-line
calculation, this approach can calculate high-accuracy results with a
100-fold computation time reduction. During the retrieval processes,
two kinds of retrieval results, i.e., baseline CTP and retrieved CTP, are
provided. The baseline CTP is derived without considering in-cloud photon
penetration, and the retrieved CTP is derived by solving the analytic
equation, taking into consideration in-cloud and below-cloud
interactions. The retrieved CTPs for the oxygen A and B bands are smaller
than their related baseline CTP. At the same time, both baseline CTP and
retrieved CTP at the oxygen B band are larger than those at the oxygen
A band. Compared to the difference in baseline CTP between the B band and
A band, the difference in retrieved CTP between these two bands is generally
reduced. Out of around 10 000 cases, in retrieved CTP between the A and B bands
we found an average bias of 93 mb with a standard deviation of 81 mb. The
cloud layer top pressure from Cloud–Aerosol Lidar and Infrared Pathfinder
Satellite Observations (CALIPSO) measurements is used for validation.
Under single-layer cloud situations, the retrieved CTPs for the oxygen
A band agree well with the CTPs from CALIPSO, the mean difference of which
within 5 mb in the case study. Under multiple-layer cloud situations, the
CTPs derived from EPIC measurements may be larger than the CTPs of high-level thin clouds due to the effect of photon penetration.

The Deep-Space Climate Observatory (DSCOVR) satellite is an observation platform orbiting within the first Sun–Earth Lagrange point (L1), 1.5 million km from the Earth, carrying a suite of instruments oriented both
earthward and sunward. One of the earthward instruments is the Earth
Polychromatic Imaging Camera (EPIC), which can take images of the
Earth with a spatial resolution of 10 km at nadir. EPIC continuously
monitors the entire sunlit Earth for backscatter, with a nearly constant
scattering angle between 168.5 and 175.5^{∘} from sunrise to sunset with 10 narrowband filters: 317, 325, 340, 388, 443, 552, 680, 688, 764, and 779 nm (Marshak et al., 2018). Of the 10 narrowband channels, there are two oxygen absorption and reference pairs, 764 nm versus 779.5 nm and 680 nm versus 687.75 nm, for the oxygen A and B bands. The cloud-top pressure (CTP) or cloud-top height (CTH) is an important cloud property for climate and weather studies. Based on differential oxygen absorption, both EPIC oxygen A-band and B-band pairs can be used to retrieve CTP. It is worth noting that although CTP and CTH reference the same characteristic of clouds, the conversion between the two depends on the related atmospheric profile.

Although the theory of using oxygen absorption bands to retrieve CTP was proposed decades ago (Yamamoto and Wark, 1961), it is still very challenging to perform retrievals accurately due to the complicated in-cloud penetration effect (Yang et al., 2019, 2013; Davis et al., 2018a, b; Richardson and Stephens, 2018; Loyola et al., 2018; Lelli et al., 2014, 2012; Schuessler et al., 2013; Ferlay et al., 2010; Rozanov and Kokhanovsky, 2004; Kokhanovsky and Rozanov, 2004; Daniel et al., 2003; Koelemeijer et al., 2001; Kuze and Chance, 1994; O'Brien and Mitchell, 1992; Fischer and Grassl, 1991). To estimate the CTP from satellite measurements, many approaches have been designed to retrieve cloud effective top pressures without considering in-cloud photon penetration. These approaches do not consider light penetrating cloud, and therefore the derived CTH is lower than the cloud top, and the effective top pressure is higher than CTP. In the meantime, to improve the retrieval accuracy of CTP, various techniques have been applied to retrieval methods with in-cloud photon penetration. For example, Kokhanovsky and Rozanov (2004) proposed a simple semi-analytical model for the calculation of the top-of-atmosphere (TOA) reflectance of an underlying surface–atmosphere system accounting for both aerosol and cloud scattering. Based on the work of Kokhanovsky and Rozanov (2004), Rozanov and Kokhanovsky (2004) developed an asymptotic algorithm for the CTH and the geometrical thickness determination using measurements of the cloud reflection function. This retrieval method was applied by Lelli et al. (2012, 2014) to derive CTH using measurements from the GOME instrument onboard the ESA ERS-2 space platform.

Currently, based on the measurements of DSCOVR EPIC, the Atmospheric
Science Data Center (ASDC) at the National Aeronautics and Space Administration
(NASA) Langley Research Center archives both calibrated EPIC reflectance ratio data and processed Level 2 cloud retrieval products, including cloud
cover, cloud optical depth (COD), and cloud effective top pressure at oxygen A and B bands (Yang et al., 2019; Holdaway and Yang, 2016; Meyer et al., 2016). By using EPIC reflectance ratio data at oxygen A-band and B-band absorption to reference channels, Yang et al. (2013) developed a method to retrieve CTH and cloud geometrical thickness simultaneously for a fully cloudy scene over ocean surface. First, their method calculates cloud centroid heights for both A- and B-band channels using the ratios between the reflectance of the absorption and reference channels, then derives the CTH and the cloud geometrical thickness from the two-dimensional lookup tables that relate the sum and the difference between the retrieved centroid heights for the A and B bands to the CTH and the cloud geometrical thickness. The difference in the O_{2} A- and B-band cloud centroid heights results from the different penetration depths of the two bands. Compared to the cloud height variability, the penetration depth differences are much smaller, and the retrieval accuracy from this method can be affected by instrument noise (Davis et al., 2018a, b).

In this paper, to address the issue of in-cloud penetration, we propose an analytic method to retrieve the CTP by using DSCOVR EPIC oxygen A- and B-band observations. This analytical method adopts ideas of the semi-analytical model (Kokhanovsky and Rozanov, 2004; Rozanov and Kokhanovsky, 2004) with a quadratic EPIC analytic radiative transfer equation to analyze the radiative transfer in oxygen A- and B-band channels. The structure of this paper is as follows: Sect. 2 describes the theory and methods, which includes several subsections, i.e., the introduction of DSCOVR EPIC oxygen A- and B-band filters, the theory of CTP retrieval based on EPIC oxygen A- and B-band observations, and the detailed retrieval algorithm. Section 3 describes the application and validation of the CTP retrieval method, which also includes several subsections, i.e., case studies of CTP retrieval, validation of the retrieval method, and retrieval of global observation. Section 4 states the conclusions of this study.

## 2.1 DSCOVR EPIC oxygen A- and B-band filters

EPIC filters at 764 and 779 nm cover the oxygen A-band absorption and
reference bands, respectively (Fig. 1a). The high-resolution absorption
optical depth spectrum at the oxygen A band and B band is calculated by the
Line-By-Line Radiative Transfer Model (LBLRTM; Clough et al., 2005) with
the HITRAN 2016 database (Gordon et al., 2017) for the US standard atmosphere.
In this wavelength range, the O_{3} absorption is very weak (O_{3} optical depth <0.003) and there are no other gas absorptions. The background aerosol and Rayleigh scattering optical depth vary smoothly within the A-band range; the differences between in-band and the reference band are
negligible at nominal EPIC response functions. EPIC filters at 688 and
680 nm cover the oxygen B-band absorption and reference band, respectively
(Fig. 1b). Compared to the oxygen A band, O_{3} absorption is slightly stronger in the oxygen B-band range, with an O_{3} optical depth around 0.01. Any water vapor absorption in the B-band range is negligible. In the standard atmospheric model, from the oxygen B-band reference band to the absorption band, the O_{3} absorption and Rayleigh scattering optical depth decreased by approximately 0.002 and 0.002, respectively. This may have some impacts on the CTP retrieval from the oxygen B band (more discussion in later sections). It is worth noting that for EPIC measurements at both oxygen A and B bands, the surface influence cannot be ignored. For example, in snow- or ice-covered area the surface albedo is high; in plant-covered area, the surface albedo changes substantially between the oxygen A band and B band due to the impact of the spectral red edge (Seager et al., 2005).

In general, if we use the pair of oxygen A and B absorption and reference
bands together, the impact of other absorption lines, background Rayleigh
scattering, and aerosol optical depth is very limited. At the same time, as
a well-mixed major atmospheric component, the vertical distribution of
oxygen in the atmosphere, is very stable under varying atmospheric
conditions. Thus, we can use the ratio of reflected radiance (or
reflectance) at the TOA of oxygen absorption to reference bands (i.e.,
*R*_{764} and *R*_{779}, *R*_{688}, and *R*_{680}) to study the photon path length distribution and derive cloud information. Also, compared to any specific EPIC oxygen absorption bands (i.e., *R*_{764} and *R*_{688}), the ratios of absorption to reference channels (i.e., *R*_{764}∕*R*_{779} and *R*_{688}∕*R*_{680}) are less impacted by instrument calibration and other measurement error. This can be explained by the following reasons: first, the EPIC measurements at the oxygen A and B absorption and reference bands share the same sensor and optical system, so when calculating the ratios, some preprocessing calibration errors can be reduced. Second, to calculate *R*_{764} and *R*_{688}, the ratio of lunar reflectance at neighboring channels (i.e., *F*(764, 779) and *F*(688, 680)) to the calibration factors of the oxygen A and B reference bands (i.e., *K*_{779} and *K*_{680}) are used (Geogdzhayev and Marshak, 2018; Marshak et al., 2018). Therefore, the accuracy of *R*_{764} and *R*_{688} is determined by the stability of *F*(764, 779) and *F*(688, 680) and the accuracy of *K*_{779} and *K*_{680} together. But the accuracy of the absorption to reference ratios is
only determined by the stability of *F*(764, 779) and *F*(688, 680).

## 2.2 Theory of CTP retrieval based on EPIC oxygen A- and B-band observation

In our study, we tried two methods to retrieve the CTP based on EPIC oxygen A-band and B-band measurements: (1) build a lookup table (LUT) for various atmospheric conditions and perform the retrieval by searching the LUT; (2) develop an analytic transfer inverse model for EPIC observations and calculate the related coefficients based on a series of simulated values, then use this analytic transfer inverse model to retrieve the CTP. In this paper, we mainly focus on the second method.

### 2.2.1 Method 1: LUT-based approach

One commonly used method of retrieval for satellite observation is through the building and usage of LUTs (Loyola et al., 2018, Gastellu-Etchegorry and Esteve, 2003). The LUT-based approach can be fast because the most computationally expensive part of the inversion procedure is completed before the retrieval itself. For DSCOVR EPIC observations, we can build LUTs by simulating EPIC measurements under various atmospheric conditions, such as different surface albedo, solar zenith and viewing angles, COD, CTP, and cloud pressure thickness. Comparing the related simulated reflectance at the oxygen absorption and reference bands, we can obtain two LUTs for the reflectance ratios of absorption to reference at the EPIC oxygen A band and B band, respectively, which can be used for the CTP retrieval. Detailed information on the simulated reflectance ratio of absorption to reference is stated in Sect. 2.3.3.

During the retrieval process, the EPIC measurements (e.g., reflectance at
oxygen A and B bands) with related solar zenith and viewing angles can be
obtained from the EPIC Level 1B data; COD information (retrieved from other
EPIC channels) can be obtained from EPIC Level 2 data. At the same time, we
can get surface albedo from Global Ozone Monitoring Experiment 2 (GOME-2)
surface Lambertian-equivalent reflectivity (LER) data (Tilstra et al., 2017). At this point the CTP and cloud pressure thickness are the only
unknown variables. The cloud pressure thickness or the cloud vertical
distribution has a substantial impact on the accuracy of the CTP retrievals
(Carbajal Henken et al., 2015; Fischer and Grassl, 1991; Rozanov and
Kokhanovsky, 2004; Preusker and Lindstrot, 2009). In this study, the cloud
pressure thickness is used as an input parameter to retrieve the CTP.
However, no related accurate cloud pressure thickness is currently provided by other
satellite sensors. To constrain the error from the estimation of cloud
pressure thickness, we related it to the cloud optical thickness. It is
reasonable because clouds with higher optical thickness normally have higher
values of pressure thickness. To explore the correlation between cloud
pressure thickness and cloud optical thickness, we use the related cloud
data from the Modern-Era Retrospective analysis for Research and Applications
version 2 (MERRA-2; Gelaro et al., 2017), which is a NASA atmospheric
reanalysis for the satellite era using the Goddard Earth Observing System
Model version 5 (GEOS-5) with Atmospheric Data Assimilation System (ADAS).
Based on statistical analysis of 1-year single-layer liquid water clouds
over an oceanic region (23.20^{∘} S, 170.86^{∘} W;
2.11^{∘} S, 144.14^{∘} W) in 2017, we can get an
equation for cloud pressure thickness approximation, i.e., cloud pressure
thickness (mb) $=\mathrm{2.5}\times \text{COD}+\mathrm{23}$. The derived correlation coefficients are dependent on the case region and time selections. Due to the complexity of cloud vertical distribution, whatever the accuracy of the correlation coefficients is, the estimation will certainly bring in error.

With an estimated cloud pressure thickness, a multivariable LUT searching method can then be used to interpolate and obtain the CTP. It is worth noting that the reflectance ratio of absorption to reference can be seen as a function of surface albedo, solar zenith and viewing angles, COD, CTP, and cloud pressure thickness. Some atmospheric variables will have a nonlinear effect on the reflectance ratio. For example, the reflectance ratio is more sensitive to the variation of COD when COD is small. Overall, the reflectance ratio varies monotonically and smoothly with these variables (shown in Fig. 3). With a relatively high-resolution simulated table, we can use a localized linear interpolation method to estimate the proper values. Multiple interpolations are needed for this method to decrease the number of LUT dimensions, which will cost more time than the analytic transfer inverse model method. The retrieval error of this method is determined by the resolution of the LUT; i.e., the higher the resolution, the higher the retrieval accuracy. However, for multiple dimensional LUTs, the increase in resolution will increase the table size exponentially, which will increase computational cost substantially for the table building and inverse searching. Another possible method to increase the retrieval accuracy is using different interpolation methods. For example, if the value of LUT varies nonlinearly with a variable, using a high-order interpolation method may be better than using a linear interpolation method (Dannenberg, 1998).

### 2.2.2 Method 2: analytic transfer inverse model

For a long time, various efforts have been devoted to the study of radiative transfer in the atmosphere, including scattering, absorption, and emission (Chandrasekhar, 1960; Irvine, 1964; Ivanov and Gutshabash, 1974; van de Hulst, 1980, 2012; Ishimaru, 1999; Thomas and Stamnes, 2002; Davis and Marshak, 2002; Kokhanovsky et al., 2003; Marshak and Davis, 2005; Pandey et al., 2012). In this study, we develop an analytic radiative transfer equation to analyze the radiative transfer at the oxygen A and B bands. Through solving the analytic equation, we can retrieve the CTP information directly. The theory of CTP retrieval is similar for EPIC oxygen A-band and B-band observations. Here we use the oxygen A band as an example to study the radiative transfer model. For the oxygen A band, the photon path length distribution is capable of describing vital information related to a variety of cloud and atmospheric characteristics:

where *p*(*l*) is the photon path length distribution, *κ*_{ν} is the
gaseous absorption coefficient at wavenumber *v*, $\mathit{\mu}=\mathrm{cos}\left(\mathit{\theta}\right),{\mathit{\mu}}_{\mathrm{0}}=\mathrm{cos}\left({\mathit{\theta}}_{\mathrm{0}}\right)$, (*θ*,*ϕ*) and (*θ*_{0},*ϕ*_{0}) are zenith and azimuth angles for the solar
and sensor view, respectively, and *I*_{0} and *I*_{ν} are incident solar
radiation and sensor-measured solar radiation, respectively.

When clouds exist, the incident solar radiation is reflected to TOA in three
primary ways. First, incident solar radiation is reflected by the cloud-top
layer directly as a result of single scattering. Second, the incident solar
radiation will penetrate into the cloud and be reflected back to TOA through
the cloud top via multiple scattering. Third, the incident solar radiation will
pass through the cloud and arrive at the surface; after that, it is reflected
back into the cloud and finally scattered back to TOA through the cloud
top. Due to the position of the EPIC instrument and the long distance between
EPIC and Earth, we can consider the solar zenith angle and sensor view
angle to be nearly reversed. At the oxygen A band, the reflected solar radiation
will be reduced due to oxygen absorption depending on photon path length
distributions. Absorption is negligible in the oxygen A-band reference band.
The oxygen A band and its reference band are also attenuated by air masses and
aerosol through Rayleigh scattering and aerosol extinction. In the standard
atmospheric model, the optical depth of Rayleigh scattering (*τ*_{Ray}) at the oxygen A band (B band) and its reference band is 0.026 (0.040) and 0.024 (0.042), respectively (Bodhaine et al., 1999). The absolute
difference in Rayleigh scattering optical depth ($\mathrm{\Delta}{\mathit{\tau}}_{\mathrm{Ray}}={\mathit{\tau}}_{\mathrm{Ray}}^{\mathrm{In}-\mathrm{band}}-{\mathit{\tau}}_{\mathrm{Ray}}^{\mathrm{Ref}}$) between them is within 0.002. Compared to Rayleigh scattering, the difference in background aerosol optical depth (Δ*τ*_{Aer}) between absorbing and reference bands is smaller, within 0.0005. Therefore, the attenuations from Rayleigh scattering and aerosol extinction at EPIC oxygen absorption and its reference band are close to each other. Thus, when we use the ratio of EPIC-measured reflectance at the oxygen A band and its reference band to derive the photon path length distribution and retrieve cloud information such as CTP, the impact of Rayleigh scattering and aerosol extinction can be simplified in the analytic transfer inverse model.

To simplify the analytic transfer inverse model for EPIC observations, we
made a series of assumptions, e.g., isotropic component, a plane-parallel
homogenous cloud assumption with quasi-Lambertian reflecting surfaces. These
assumptions have been widely used in radiative transfer calculation for
cloud studies. In this model, *μ* and *μ*_{0} are the same as in Eq. (1); *ϕ* is the relative azimuth angle between the Sun and satellite sensors; *A*_{surf} is the surface albedo; ${\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}$, ${\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Base}}$, and ${\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Surface}}$ are oxygen A-band absorption optical depth from TOA to the cloud-top layer, cloud-bottom layer, and surface, respectively; $\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Above}-\mathrm{Cld}}$, $\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{In}-\mathrm{Cld}}$, and $\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Below}-\mathrm{Cld}}$ are layered oxygen A-band
absorption optical depth above cloud, in cloud, and below cloud,
respectively; and the functions *f* indicate their contribution to the ratio of measured
reflectance at the oxygen A band (*R*_{A}) and reference band (*R*_{f}). A detailed analysis of the EPIC analytic transfer inverse model is shown as follows.

- 1.
*Above cloud.*The reflected solar radiation is determined by the oxygen absorption optical depth above the cloud and air mass directly.$$\begin{array}{}\text{(2)}& \begin{array}{rl}& f\left(\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Above}-\mathrm{Cld}},{\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=f\left(\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Above}-\mathrm{Cld}}\right)f\left({\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi},\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}={a}_{\mathrm{0}}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}\left({\displaystyle \frac{\mathrm{1}}{\mathit{\mu}}}+{\displaystyle \frac{\mathrm{1}}{{\mathit{\mu}}_{\mathrm{0}}}}\right)\end{array}\end{array}$$Here,

*a*_{0}is a weight coefficient. - 2.
*Within cloud.*The reflected solar radiation is not only determined by oxygen absorption optical depth above cloud and in cloud, but also by penetration-related factors, e.g., COD. Due to photon penetration, the oxygen parameter ${\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}$ influences the enhanced path length absorption.$$\begin{array}{}\text{(3)}& \mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{In}-\mathrm{Cld}}={\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Base}}-{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}\end{array}$$Equivalence theorem (Irvine, 1964; Ivanov and Gutshabash, 1974; van de Hulst, 1980) is used to separate absorption from scattering.

$$\begin{array}{}\text{(4)}& \begin{array}{rl}& f\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}},\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{In}-\mathrm{Cld}},{\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=f\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}},\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{In}-\mathrm{Cld}}\right)f\left({\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=f\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}\right){f}_{\mathrm{1}}\left({\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+f\left(\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{In}-\mathrm{Cloud}}\right){f}_{\mathrm{2}}\left({\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\end{array}\end{array}$$$f\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}\right)$ is determined by two absorption dependences: strong ($\sim \sqrt{{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}}$) and weak ($\sim {\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}$).

$$\begin{array}{}\text{(5)}& f\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}\right)={a}_{\mathrm{1}}\sqrt{{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}}+{b}_{\mathrm{1}}\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}\right)\end{array}$$Based on asymptotic approximation (Kokhanovsky et al., 2003; Pandey et al., 2012), the reflection of a cloud without considering below-cloud interactions is given by Eq. (6).

$$\begin{array}{}\text{(6)}& \begin{array}{rl}R\left(\mathit{\tau},\mathit{\mu},{\mathit{\mu}}_{\mathrm{0}},T\right)& ={R}_{\mathrm{0}}^{\mathrm{\infty}}\left(\mathit{\tau},\mathit{\mu},{\mathit{\mu}}_{\mathrm{0}}\right)-TK\left(\mathit{\mu}\right)K\left({\mathit{\mu}}_{\mathrm{0}}\right)\\ & ={R}_{\mathrm{0}}^{\mathrm{\infty}}\left(\mathit{\tau},{f}_{\mathrm{1}}\left(\mathit{\mu},{\mathit{\mu}}_{\mathrm{0}}\right)\right)-T{f}_{\mathrm{2}}\left(\mathit{\mu},{\mathit{\mu}}_{\mathrm{0}}\right)\end{array}\end{array}$$Here, ${R}_{\mathrm{0}}^{\mathrm{\infty}}$ is the reflectance of a semi-infinite cloud,

*K*(*μ*) is the escape function of*μ*, and*T*is the global transmittance of a cloud.*T*can be estimated by Eq. (7), with the cloud optical thickness*τ*_{cld}, the asymmetry parameter*g*, and a numerical constant*α*=1.07.$$\begin{array}{}\text{(7)}& T={\displaystyle \frac{\mathrm{1}}{\mathrm{0.75}{\mathit{\tau}}_{\mathrm{cld}}\left(\mathrm{1}-g\right)+\mathit{\alpha}}}\end{array}$$The

*f*_{1}and*f*_{2}functions have a quadratic form as follows.$$\begin{array}{}\text{(8)}& \begin{array}{rl}{f}_{i-\mathrm{1}}& \phantom{\rule{0.25em}{0ex}}={a}_{i}T+{b}_{i}\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)\\ & \phantom{\rule{0.25em}{0ex}}+{c}_{i}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{d}_{i}\mathit{\mu},{\mathit{\mu}}_{\mathrm{0}},i=\mathrm{2},\mathrm{3}\end{array}\end{array}$$Combining Eqs. (4), (5), and (8), we can get Eq. (9).

$$\begin{array}{}\text{(9)}& \begin{array}{rl}& f\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}},\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Cld}},{\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=\left({a}_{\mathrm{1}}\sqrt{{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}}+{b}_{\mathrm{1}}\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}\right)\right)({a}_{\mathrm{2}}T+{b}_{\mathrm{2}}\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{c}_{\mathrm{2}}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{d}_{\mathrm{2}}\mathit{\mu}{\mathit{\mu}}_{\mathrm{0}})\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{In}-\mathrm{Cloud}}({a}_{\mathrm{3}}T+{b}_{\mathrm{3}}\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{c}_{\mathrm{3}}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{d}_{\mathrm{3}}\mathit{\mu}{\mathit{\mu}}_{\mathrm{0}})\end{array}\end{array}$$ - 3.
*Below cloud.*The equivalence theorem used for below cloud is similar to within cloud (Kokhanovsky et al., 2003; Pandey et al., 2012).$$\begin{array}{}\text{(10)}& \begin{array}{rl}& f\left(\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Below}-\mathrm{Cld}},{\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}=T{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Surface}}{\displaystyle \frac{{A}_{\mathrm{Surf}}}{\mathrm{1}+\left({e}_{\mathrm{4}}\cdot T+{f}_{\mathrm{4}}\right)\cdot {A}_{\mathrm{Surf}}}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\cdot \left({a}_{\mathrm{4}}T+{b}_{\mathrm{4}}\left(\mathit{\mu}+\mathit{\mu}\mathrm{0}\right)+{c}_{\mathrm{4}}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{d}_{\mathrm{4}}\mathit{\mu}{\mathit{\mu}}_{\mathrm{0}}\right)\end{array}\end{array}$$Combining Eqs. (2), (9), and (10), we can get the total EPIC analytic transfer equation as follows.

$$\begin{array}{}\text{(11)}& \begin{array}{rl}& -\mathrm{log}\left({\displaystyle \frac{{R}_{A}}{{R}_{f}}}\right)=f\left(\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Above}-\mathrm{Cld}},{\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+f\left({\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}},\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Cld}},{\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+f\left(\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Below}-\mathrm{Cld}},{\mathit{\mu}}_{\mathrm{0}},\mathit{\mu},\mathit{\varphi}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+\mathrm{\Delta}{\mathit{\tau}}_{\mathrm{BG}}\left({\displaystyle \frac{\mathrm{1}}{\mathit{\mu}}}+{\displaystyle \frac{\mathrm{1}}{{\mathit{\mu}}_{\mathrm{0}}}}\right)\end{array}\end{array}$$In Eq. (11), Δ

*τ*_{BG}represents the sum of the optical depth difference in background extinction (i.e., Rayleigh scattering Δ*τ*_{Ray}, aerosol extinction Δ*τ*_{Aer}, and O_{3}$\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{3}}}$) between oxygen in-band and the reference band, as shown in Eq. (12).$$\begin{array}{}\text{(12)}& \mathrm{\Delta}{\mathit{\tau}}_{\mathrm{BG}}=\mathrm{\Delta}{\mathit{\tau}}_{\mathrm{Ray}}+\mathrm{\Delta}{\mathit{\tau}}_{\mathrm{Aer}}+\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{3}}}\end{array}$$As stated in the previous subsection, in the standard atmospheric model with background aerosol loading, $\left(\mathrm{\Delta}{\mathit{\tau}}_{\mathrm{Ray}},\mathrm{\Delta}{\mathit{\tau}}_{\mathrm{Aer}},\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{3}}}\right)$ is approximately (0.002, 0.0005, −0.0005) and (−0.002, −0.0005, −0.002), respectively, at the oxygen A and B bands, and thus Δ

*τ*_{BG}is approximately 0.002 and −0.0045, respectively, at these two bands.In this total analytic equation, there are 17 coefficients (${a}_{\mathrm{0}},{a}_{\mathrm{1}},{b}_{\mathrm{1}},{a}_{\mathrm{2}},\mathrm{\dots}{d}_{\mathrm{4}},{e}_{\mathrm{4}},{f}_{\mathrm{4}}$), which can be calculated through a nonlinear regression algorithm according to a series of simulated values for different atmospheric conditions. Based on Eq. (11), we can finally obtain a quadratic equation,

$$A{\sqrt{{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}}}^{\mathrm{2}}+B\sqrt{{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Top}}}+C=\mathrm{0},$$where the parameters

*A*,*B*, and*C*can be derived from Eq. (11) directly, as shown in Eq. (13).$$\begin{array}{}\text{(13a)}& \begin{array}{rl}A& \phantom{\rule{0.25em}{0ex}}={a}_{\mathrm{0}}\left({\displaystyle \frac{\mathrm{1}}{\mathit{\mu}}}+{\displaystyle \frac{\mathrm{1}}{{\mathit{\mu}}_{\mathrm{0}}}}\right)+{b}_{\mathrm{1}}({a}_{\mathrm{2}}T+{b}_{\mathrm{2}}\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{c}_{\mathrm{2}}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{d}_{\mathrm{2}}\mathit{\mu}{\mathit{\mu}}_{\mathrm{0}})\end{array}\end{array}$$$$\begin{array}{}\text{(13b)}& \begin{array}{rl}B& \phantom{\rule{0.25em}{0ex}}={a}_{\mathrm{1}}({a}_{\mathrm{2}}T+{b}_{\mathrm{2}}\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{c}_{\mathrm{2}}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{d}_{\mathrm{2}}\mathit{\mu}{\mathit{\mu}}_{\mathrm{0}})\end{array}\end{array}$$$$\begin{array}{}\text{(13c)}& \begin{array}{rl}C& \phantom{\rule{0.25em}{0ex}}=-\mathrm{log}\left({\displaystyle \frac{{R}_{A}}{{R}_{f}}}\right)-\mathrm{\Delta}{\mathit{\tau}}_{\mathrm{BG}}\left({\displaystyle \frac{\mathrm{1}}{\mathit{\mu}}}+{\displaystyle \frac{\mathrm{1}}{{\mathit{\mu}}_{\mathrm{0}}}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}-\mathrm{\Delta}{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{In}-\mathrm{Cloud}}({a}_{\mathrm{3}}T+{b}_{\mathrm{3}}\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+{c}_{\mathrm{3}}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{d}_{\mathrm{3}}\mathit{\mu}{\mathit{\mu}}_{\mathrm{0}})\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}-T{\mathit{\tau}}_{{\mathrm{O}}_{\mathrm{2}}}^{\mathrm{Surface}}{\displaystyle \frac{{A}_{\mathrm{Surf}}}{\mathrm{1}+\left({e}_{\mathrm{4}}\cdot T+{f}_{\mathrm{4}}\right)\cdot {A}_{\mathrm{Surf}}}}\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}({a}_{\mathrm{4}}T+{b}_{\mathrm{4}}\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{c}_{\mathrm{4}}T\left(\mathit{\mu}+{\mathit{\mu}}_{\mathrm{0}}\right)+{d}_{\mathrm{4}}\mathit{\mu}{\mathit{\mu}}_{\mathrm{0}})\end{array}\end{array}$$When these parameters (i.e.,

*A*,*B*, and*C*) are obtained from EPIC observation data and other data sources, we can easily solve the quadratic equation to retrieve the cloud-top O_{2}absorption depth and then derive CTP.

## 2.3 Detailed retrieval algorithm

As previously stated, in method 2, the analytic EPIC equation (i.e., Eq. 11) is key for the CTP retrieval. To derive the coefficients of Eq. (11), a series of model simulations for various atmospheric conditions is needed. Thus, developing a radiative transfer model to simulate the EPIC measurements at A and B bands and their reference bands is the first thing we need to complete.

### 2.3.1 Oxygen A- and B-band absorption coefficient calculation

To simulate the EPIC measurements, one of the most important steps is calculating oxygen absorption coefficients at the oxygen A band and B band. In this step, the HITRAN 2016 database is used to provide the absorption parameters, and the LBLRTM package is used to calculate oxygen absorption coefficients layer by layer. In our algorithm, the whole Earth atmosphere is divided by 63 layers.

Since oxygen absorption coefficients are pressure- (or pressure-squared) and
temperature-dependent, and the line shapes (*k*_{i}) of oxygen A and
B bands are well fitted as Lorentzian in the lower atmosphere, the
relationship can be written as follows:

where *S*_{i} is the line intensity, *v*_{i} and *α*_{i} are the
line center wavenumber and half-width, respectively, and *P*_{0} and *T*_{0}
are standard atmospheric pressure and temperature, respectively.

In the simulation of EPIC measurements, the atmospheric layer at a given layer-average pressure can have drastically different temperature depending on the atmospheric profile in use. To ensure the accuracy of simulation, we need to use the LBLRTM package to calculate oxygen absorption coefficients for each pressure–temperature profile, which is a time-consuming process. Our goal has been to find a simple and fast method to calculate oxygen absorption coefficients for different atmospheric profiles. Based on the study of Chou and Kouvaris (1986), Min et al. (2014) proposed a fast method to calculate oxygen absorption optical depth for any given atmosphere by using a polynomial fitting function, as shown in Eq. (16):

where *A*_{vLM} is optical depth for layer *L*, spectral point *v*, and atmosphere model *M*; ${\mathit{\rho}}_{{\mathrm{O}}_{\mathrm{2}}}$ is molecular column density
($\frac{\text{molecules}}{{\mathrm{cm}}^{\mathrm{2}}}\times {\mathrm{10}}^{-\mathrm{23}}$); *T*_{LM} is the average temperature for layer *L* for a given atmosphere; and *T*_{mL} is average temperature over all six typical geographic–seasonal model atmospheres (M1 to M6, i.e., tropical model, midlatitude summer model, midlatitude winter model, subarctic summer model, subarctic winter model, and the US standard 1976 model) for layer *L*. To derive the coefficients *a*_{0}, *a*_{1}, and *a*_{2}, we first calculated oxygen optical depth coefficients for all typical atmospheres (M1 to M6) by using the LBLRTM package and then selected three of them (e.g., M1, M5, and M6) to calculate the polynomial fitting coefficients. This method has been successfully used by Min et al. (2014) to simulate high-resolution oxygen A-band measurements.

### 2.3.2 Fast radiative transfer model for simulating high-resolution oxygen A and B bands

At oxygen A and B absorption bands, there are lots of absorption lines,
therefore we cannot simply calculate narrowband mean optical depth and then
calculate the radiation for various atmospheric conditions when simulating
EPIC narrowband measurements. The correct way is described as follows:
firstly, simulate the solar radiation spectrum *S*(*k*(*λ*))
under specific atmospheric conditions, then integrate the spectrum with EPIC narrowband filter *R*(*k*(*λ*)) to obtain simulated
narrowband measurements (Eq. 17).

With the high-spectral-resolution oxygen absorption coefficient data, we can simulate the high-resolution upward diffuse oxygen A-band or B-band spectrum through DISORT code (Stamnes et al., 1988) for any given atmospheric condition, which has various surface albedo, solar zenith angle (SZA), COD, CTH (CTP), and cloud geometric (pressure) thickness. However, due to the high spectral resolution, it is very time-consuming when performing line by line (LBL) calculations. Thus, developing a fast radiative transfer model for simulating the high-resolution oxygen A-band and B-band spectrum is necessary.

In this project, the double-*k* approach is used to develop a fast radiative transfer model for the oxygen A band and B band, respectively. Min and Harrison (2004) and Duan et al. (2005) proposed a fast radiative transfer model. In their approach, the radiation from absorption and scattering processes of cloud and aerosol are split into the single- and multiple-scattering components: the single scattering component is computed line by line (LBL), while multiple scattering (second order and higher) radiance is approximated.

Equation (18) is from Eq. (1) in Duan et al. (2005): ss and ms indicate single and multiple scattering, respectively. *Z* is the optical properties of the
atmosphere as a function of pressure *p* and temperature *T*, with *P* being the phase function of that layer. h and l represent higher and lower number of layers and streams, respectively. *F* is the transform function between wavenumber space and *k* space, defined from a finite set of *k*(*λ*_{i}).

The application of the double-*k* approach in the oxygen A band has been presented in detail in Duan et al. (2005). Here we take the oxygen B band as an example. The detailed fast radiative transfer model for simulating the high-resolution oxygen B band is as follows: the first-order scattering radiance is calculated accurately by using a higher number of layers and streams for all required wavenumber grid points. The multiple-scattering component is extrapolated and/or interpolated from a finite set of calculations in the space of two integrated gaseous absorption optical depths to the wavenumber grids: a double-*k* approach. The double-*k* approach substantially reduces the error due to
the uncorrelated nature of overlapping absorption lines. More importantly,
these finite multiple-scattering radiances at specific *k* values are computed
with a reduced number of layers and/or streams in the forward radiative
transfer model. To simulate an oxygen B-band spectrum with high accuracy, 33 *k* values and 99 calculations of radiative transfer are chosen in our program.
This results in around a 100-fold time reduction with respect to the
standard forward radiative transfer calculation.

As shown in Fig. 2, under clear-sky and thin liquid water cloud situations,
the simulated high-resolution upward diffuse oxygen B-band spectra from LBL
calculation and the double-*k* approach are compared. The spectrum difference
between LBL calculation and the double-*k* approach is very small (Fig. 2a). Under
both situations, most of the relative difference between these two methods
is under 0.5 %. The obvious relative difference (>1 %) occurs only in the wavelength range with high absorption optical depth,
which has little contribution to the integrated solar radiation. Therefore,
for the simulated narrowband measurements at the EPIC oxygen B band, the
relative difference between LBL and the double-*k* approach is much smaller than that of the high-resolution spectrum, which is less than 0.1 % for a clear day. Compared to a clear-sky situation, the relative difference for cloud situations can be bigger. As shown in Table 1, the relative difference is −0.06 % and −0.32 % for typical high-level optical thin cloud and low-level thick cloud situations, respectively. The comparison of simulated narrowband measurements at the EPIC oxygen A-band channel (764 nm) is also shown in Table 1; the relative differences between LBL and the double-*k* approach are −0.06 %, 0.21 %, and 0.23 % for a clear day, high-level thin cloud, and
low-level thick cloud cases, respectively. In general, the accuracy of
double-*k* approach for both oxygen A and B absorption bands is high.

### 2.3.3 Simulation of oxygen A and B bands for different atmospheric conditions

Using the EPIC measurement simulation package, we made a series of simulations with different settings for surface albedo, solar zenith angle, COD, CTH (CTP), and cloud geometric (pressure) thickness (or cloud-bottom height). The results of these simulations consist of a data table, which can be used not only to calculate the coefficients for the analytic equation, but also to study the sensitivity of every variable.

According to the previous theory analysis, the ratio of reflectance radiance (i.e., absorption to the reference) at TOA is determined by the photon path length distribution at oxygen A and B bands: the larger the mean photon path length, the stronger the absorption and the smaller the reflectance ratio. To make the figures easy to view and understand, we use the cloud-top and cloud-bottom geometric height to represent CTP and thickness information in Fig. 3. As shown in Fig. 3a, the ratio of upward diffuse radiance at the oxygen B band and its reference band is sensitive to the cloud-top height (pressure). The higher the CTH, the larger the ratio. At the same time, this ratio is affected by the cloud-bottom height (or cloud geometric thickness) when the other cloud parameters are fixed; the lower the cloud bottom (or the larger the cloud geometric thickness), the smaller the ratio. It is consistent with the theory analysis: (1) the higher the CTH, the shorter the mean photon path length and the weaker the absorption. (2) When the COD is given, larger cloud geometric thickness means smaller cloud density, and sunlight can penetrate deeper into the cloud, which results in a longer mean photon path length. In Fig. 3b, for clouds with given CTH, COD, and geometric thickness, the ratio decreases with the solar and view angles. This can be understood to mean that the larger the solar and view angles, the longer the mean photon pathlength and the stronger the absorption. In Fig. 3c, for clouds with given CTH and geometric thickness, when the COD is small (e.g., COD <5), the reflectance ratio increases with COD. However, when COD is larger than 16, the effect of COD is small. This is because the larger the COD, the shallower the sunlight penetration, and the shorter the mean photon pathlength. In Fig. 3d, for clouds with given COD, CTP, and geometric thickness, the ratio decreases with surface albedo. The smaller the COD, the stronger the impact of the surface albedo. This is because thick cloud prevents incident sunlight from passing through it to reach the surface and also prevents reflected light from going back to the TOA.

For the oxygen A band, the ratio of upward diffuse radiance at absorption and reference bands shows similar characteristics as the oxygen B band. Compared to the oxygen B band, under the same atmospheric conditions, the oxygen absorption at the A band is stronger, and the ratio of the A band to its reference band has smaller values (shown in Fig. 4). As stated previously, for land area covered with plants, the surface albedo may change substantially from the oxygen B band to A band due to the presence of the red edge. Therefore, accurate spectral data on surface albedo for CTP retrieval are vitally important, especially for optically thin clouds.

## 3.1 Case studies of CTP retrieval

The dataset of DSCOVR EPIC measurements at 00:17:51 GMT on 25 July 2016 is
used for the case studies. The reflectance at oxygen A and B bands with
related solar zenith and viewing angles is obtained from the EPIC Level 1B
data; COD information (retrieved from other EPIC channels) is obtained from
EPIC Level 2 data. The surface albedo data are obtained from Global Ozone
Monitoring Experiment 2 (GOME-2) surface Lambertian-equivalent reflectivity
(LER) data. Detailed information on the datasets is shown in the
“Data availability” section. To reduce the impact of the Earth's surface, we
selected the region located in the spatial range of 75^{∘} S to
85^{∘} N and 177 to 175^{∘} W for case studies, which is mainly covered by ocean. To constrain the influence of surface albedo and broken clouds, only pixels with total cloud cover (i.e., EPIC cloud mask 4), surface albedo less than 0.05, and liquid assumed COD
larger than 3 are considered. In the selected region, around 10 000 pixels
are finally chosen for case studies.

In our retrieval algorithm, we have two kinds of retrieval results: baseline CTP and retrieved CTP. The baseline CTP is used as a reference for the retrieved CTP. It is similar to the effective CTP in Yang et al. (2019), which does not consider cloud penetration. The retrieved CTP is calculated by the analytic equation, which considers in-cloud and below-cloud interactions.

During the baseline CTP calculation, the impact of in-cloud penetration is
ignored, and the incident light that reached cloud top is assumed to be reflected
back directly. As shown in Eq. (19), the baseline absorption optical
depth *τ*_{base} is derived from the ratio of upward diffuse radiance
at absorption bands and their reference bands directly. According to the
model-calculated oxygen A- and B-band absorption optical depth profile at
the specific solar zenith angle, the baseline CTP can be derived directly.

As shown in Fig. 5, the baseline CTP value at the A band is slightly higher than the effective CTP from NASA ASDC L2 data. But the baseline CTP value at the B band is substantially higher than the effective CTP from NASA ASDC L2 data. For both the A band and B band, the difference between baseline CTP and effective CTP increases with the CTP. For low-level clouds, the mean differences are up to 60 and 100 mb at the A band and B band, respectively. The difference may be mainly from the calculation of oxygen A- and B-band absorption coefficients or the absorption optical depth profile.

Based on the simulated reflectance ratio under different atmospheric conditions, we can calculate the coefficients for the analytic radiative transfer equations by using a nonlinear fitting algorithm. The coefficients for different SZAs are calculated individually to reduce the fitting error. Based on the calculated coefficients, we can retrieve the CTP with DSCOVR EPIC observation data at the oxygen A and B bands.

During the CTP retrieval, with the exception of the previously mentioned
analytic equation coefficients, we can get the surface albedo data from
GOME and obtain reflectance data, solar zenith angle, view angles, and COD from
the NASA ASDC data file. Another very important step in the retrieval
processing is the acquisition of cloud pressure thickness data, which have a
substantial impact on the retrieval results. We currently use a statistical
approach (i.e., cloud pressure thickness (mb) $=\mathrm{2.5}\times \text{COD}+\mathrm{23}$) to
estimate the cloud pressure thickness based on COD. As shown in Fig. 6a–d,
the retrieved CTP when considering cloud penetration is smaller than
baseline CTP. For this case, the mean differences between baseline CTP and
retrieved CTP for the oxygen A band and B band are around 57 and 85 mb,
respectively, which is consistent with theoretical expectations. For clouds
with a given CTP, the mean photon path length will increase substantially
when considering cloud penetration. A decrease in retrieved CTP will result
in order to match the measurement ratio of absorption to reference. Compared
to the O_{2} A band, both baseline CTP and retrieved CTP for the O_{2} B band are larger (Fig. 6e–h). This is because the absorption of solar radiation in the O_{2} B band is weaker than that of the O_{2} A band, and the incident light at the oxygen B band can penetrate deeper into the cloud, allowing more light to pass through. The difference in retrieved CTP between the B band and A band (approx. 93 mb with a standard deviation of 83 mb) is generally reduced in comparison to baseline B band and A band (approx. 114 mb with a standard deviation of 73 mb). This indicates, as expected, more photon penetration correction for the B band than the A band.

We also used the LUT-based method to perform the retrieval for the same observation data because both methods share the same EPIC simulation package and the same simulated data table, the results of which are similar.

## 3.2 Validation of the retrieval method

To validate the analytic transfer inverse model method for CTP retrieval, we
used another independent measurement of CTP, i.e., cloud layer top pressure
from Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations
(CALIPSO; Vaughan et al., 2004), as a reference. For the previously stated
case, i.e., DSCOVR EPIC measurements at 00:17:51 GMT on 25 July 2016, we
used the cloud layer data from the CALIPSO IIR version 4.2 Level 2 product with
5 km resolution at 00:01:47 GMT on 25 July 2016 as its reference for
validation. To constrain the error from spatial differences between
different satellite measurements, we only chose the pixels of EPIC and
CALIPSO measurements with a spatial distance of within 0.1^{∘} (degree of latitude or longitude) to make comparisons. For the EPIC measurements, the same as previously stated, only pixels with total cloud cover (i.e., EPIC cloud mask 4), surface albedo less than 0.05, and liquid assumed COD larger than 3 are considered. As shown in Fig. 7a, there is a series of pixels (around 400 cases) from EPIC and CALIPSO measurements that can be used for the validation analysis. For reader convenience, we perform the analyses by using the case number as the *x* axis. Figure 7b shows the
comparisons of cloud layer top pressure from CALIPSO and different CTPs
(i.e., effective CTP, baseline CTP, and retrieved CTP) from EPIC
measurements. Figure 7c shows the cloud layer number measured by CALIPSO.
According to Fig. 7b and c, we can get some results: under single-layer
cloud situations, the CTPs derived from EPIC measurements are close to the
CTP from CALIPSO; under multilayer cloud situations, the CTPs derived from
EPIC measurements are larger than the CTP from CALIPSO. Figure 7d shows the
expanded view of Fig. 7b for some cases under single-layer cloud
situations. For these single-layer cloud cases (with case numbers 46–156), the mean values of CTP for CALIPSO, EPIC effective,
EPIC baseline, and EPIC retrieval are 846, 834, 866, and 850 mb, respectively.
Compared to the CTP from CALIPSO measurements, the EPIC effective and
baseline CTPs are 12 mb smaller and 20 mb larger, respectively; the EPIC
retrieval with consideration of photon penetration is only 4 mb larger. This
shows that our method for the CTP retrieval is valid and accurate under
single-layer cloud situations with COD >3 and low surface albedo. Under multilevel cloud situations, high-level clouds are often thin clouds, which can be detected by CALIPSO but are hard to derive by our retrieval method. This is because the EPIC-retrieved CTP mainly shows the pressure of cloud layer that reflects the major part of incident sunlight.

## 3.3 Retrieval of global observation

We applied our retrieval algorithm to global DSCOVR EPIC measurement data at oxygen A and B bands. During the retrieval, only pixels with total cloud cover (i.e., cloud mask index of 4), surface albedo <0.25, and COD ≥3 are considered. To make the pictures easy to visualize and analyze, all invalid values are plotted as white (or blank) pixels.

Figure 8a shows the synthesized RGB (red, green, blue) picture of EPIC measurements at 00:17:51 GMT on 25 July 2016. At this point in time sunlight covers most of the Pacific Ocean. In this figure, the white pixels represent cloud cover. Figure 8b shows the global COD (NASA ASDC L2 data), in which the white areas and colorful areas indicate clear-sky areas and cloudy areas, respectively. On the whole, the cloudy areas are consistent with the RGB image. The highlighted (red) areas indicate that the cloud systems there contain optically heavy clouds. Figure 8c shows the A-band effective CTP (NASA ASDC L2 data); the white areas indicate clear sky or no valid values, and warm (brown) and cold (blue) color areas indicate high-level and low-level clouds, respectively. According to the A-band effective CTP, high-level clouds are dominant in the equatorial area, and low-level clouds play a major role in the cloud systems in the North Pacific area. Figure 8d and e show the baseline and retrieved CTP at the A band, respectively, in which cloudy areas are consistent with the A-band effective CTP image on the whole. Due to the filtering setting in the CTP retrieval algorithm, there are more white pixels (invalid values) in these two figures. The difference in the A-band retrieved CTP and A-band effective CTP is shown in Fig. 8d. The A-band retrieved CTP is overall smaller than A-band effective CTP, with a difference within 100 mb. The highlighted (brown or red) areas are located in high-level cloud areas or large COD areas. This indicates that the complexity of a cloud system has a significant impact on the CTP retrieval. Figure 8g and h show the baseline and retrieved CTP in the B band, respectively, which are similar to but greater than the A band. As shown in Fig. 8i, the retrieved CTP at the EPIC B band is overall significantly larger than the retrieved CTP at the EPIC A band, for which the mean difference is up to 200 mb.

As previously stated in Sect. 3.2, under single-layer cloud situations, the CTPs derived from EPIC A-band measurements have good agreement with the CTP from CALIPSO measurements; under multiple-layer cloud situations, the CTPs derived from EPIC measurements may be larger than the CTPs of high-level thin clouds due to the effect of photon penetration. Therefore, in the global range for large-scale low-level stratus clouds, the retrieved CTPs from EPIC A-band measurements should agree well with the actual value of CTPs, but for a complex cloud system with multiple-layer clouds, the CTPs derived from EPIC A-band measurements may be larger than those of high-level thin clouds.

In-cloud photon penetration has significant impacts on the CTP retrieval when using DSCOVR EPIC oxygen A- and B-band measurements. To address this issue, we proposed two methods: (1) the LUT-based method and (2) the analytic transfer inverse model method for CTP retrieval with consideration of in-cloud photon penetration. In the analytic transfer inverse model method, we build an analytic equation that represents the reflection at TOA from above cloud, in cloud, and below cloud. The coefficients of this analytic equation can be derived from a series of EPIC simulations under different atmospheric conditions using a nonlinear regression algorithm. With EPIC observation data, the related solar zenith and sensor view angle, surface albedo data, COD, and estimated cloud pressure thickness, we can retrieve the CTP by solving the analytic equation.

We developed a package for the DSCOVR EPIC measurement simulation. The high-resolution radiation spectrum must be simulated first and then integrated
with the EPIC filter function in order to accurately simulate EPIC
measurements. Because this process is highly time-consuming, a polynomial
fitting function is used when calculating the oxygen absorption coefficients
under different atmospheric conditions. At the same time, the double-*k*
approach is applied to the high-resolution spectrum simulation to further
reduce time costs, which can obtain high-accuracy results with a 100-fold
time reduction. The results of the EPIC simulation measurements are
consistent with theoretical analysis.

Based on the EPIC simulation measurements, we derived a series of coefficients from various solar zenith angles for the analytic EPIC equations. Using these coefficients, we performed CTP retrieval for real EPIC observation data. We have two kinds of retrieval results: baseline CTP and retrieved CTP. The baseline CTP is similar to the effective CTP in Yang et al. (2019), which does not consider cloud penetration. The retrieved CTP is derived by solving the analytic equation with the consideration of in-cloud and below-cloud interactions. Compared to the effective CTP provided by NASA ASDC L2 data, the baseline CTP value at the A band is slightly higher, but the baseline CTP value at the B band is substantially higher. The retrieved CTP for both the oxygen A and B bands is smaller than the related baseline CTP. At the same time, compared to the oxygen A band, both baseline CTP and retrieved CTP at the oxygen B band are larger. The cloud layer top pressure from CALIPSO measurements is used to validate the CTP derived from EPIC measurements. Under single-layer cloud situations, the retrieved CTPs for the oxygen A band agree well with the CTPs from CALIPSO, for which the mean difference is within 5 mb in the case study. Under multiple-layer cloud situations, the CTPs derived from EPIC measurements may be larger than the CTPs of high-level thin clouds due to the effect of photon penetration.

Currently, this analytical transfer model method can only retrieve CTP, and it still requires cloud pressure thickness as an input parameter. However, in the satellite observations, both CTP and cloud pressure thickness are unknown. The estimation or assumption of cloud pressure thickness will bring extra error into CTP retrieval. In the near future, we plan to address this issue.

The dataset of DSCOVR EPIC Level 1B is available by visiting the website of ASDC of NASA (https://doi.org/10.5067/EPIC/DSCOVR/L1B.002, NASA LARC ASDC DAAC, 2018a; https://doi.org/10.1175/BAMS-D-17-0223.1, Marshak et al., 2018). The dataset of EPIC Level 2 is available by visiting the website of ASDC of NASA (https://doi.org/10.5067/EPIC/DSCOVR/L2_Cloud_01; NASA LARC ASDC DAAC, 2018b; https://doi.org/10.5194/amt-12-2019-2019, Yang et al., 2019). The dataset of surface albedo from GOME is available by visiting the website of Tropospheric Emission Monitoring Internet Service (TEMIS) (http://temis.nl/surface/gome2_ler/databases/, last access: 13 September 2017, TEMIS, 2017; https://doi.org/10.1002/2016JD025940, Tilstra et al., 2017). The dataset of cloud layer data from CALIPSO is available by visiting the website of NASA by using EOSIDS Earthdata account (https://subset.larc.nasa.gov/calipso/index.php, last access: 29 September 2020, NASA LARC, 2017; https://doi.org/10.1117/12.572024, Vaughan et al., 2004).

All authors contributed to planning and writing the paper. QM initiated and led the EPIC CTP retrieval project and designed the analytic transfer inverse model for EPIC observation. BY developed a fast radiative transfer model for simulating the high-resolution oxygen B band, implemented the CTP retrieval algorithms by using the analytic transfer inverse model and lookup table method, and drafted the paper. EM calculated the high-resolution oxygen absorption coefficients by using the LBLRTM model and HITRAN database. YY provided access to the EPIC Level 1B and Level 2 products and provided guidance for the evaluation of EPIC CTP retrieval. AM and ABD conducted the EPIC-related studies and provided guidance for the design of the EPIC CTP retrieval algorithm with the consideration of in-cloud photon penetration.

The authors declare that they have no conflict of interest.

This work was partially supported by NASA's Research Opportunities in Space and Earth Science (ROSES) program element for DSCOVR Earth Science Algorithms managed by Richard Eckman, the National Science Foundation (NSF), and the National Oceanic and Atmospheric Administration (NOAA) Educational Partnership Program with Minority Serving Institutions. This work was partially performed at the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration (contract no. 80NM0018D0004).

This research has been supported by the NSF (grant no. AGS-1608735) and the NOAA Educational Partnership Program with Minority Serving Institutions (grant no. NA11SEC4810003).

This paper was edited by Alexander Kokhanovsky and reviewed by three anonymous referees.

Bodhaine, B. A., Wood, N. B., Dutton, E. G., and Slusser, J. R.: On Rayleigh optical depth calculations, J. Atmos. Ocean. Tech., 16, 1854–1861, https://doi.org/10.1175/1520-0426(1999)016<1854:ORODC>2.0.CO;2, 1999.

Carbajal Henken, C. K., Doppler, L., Lindstrot, R., Preusker, R., and Fischer, J.: Exploiting the sensitivity of two satellite cloud height retrievals to cloud vertical distribution, Atmos. Meas. Tech., 8, 3419–3431, https://doi.org/10.5194/amt-8-3419-2015, 2015.

Chandrasekhar, S.: Radiative transfer, Dover, New York, 1960.

Chou, M. D. and Kouvaris, L.: Monochromatic calculations of atmospheric radiative transfer due to molecular line absorption, J. Geophys. Res., 91, 4047–4055, 1986.

Clough, S. A., Shephard, M. W., Mlawer, E. J., Delamere, J. S., Iacono, M. J., Cady-Pereira, K., Boukabara, S., and Brown, P. D.: Atmospheric radiative transfer modeling: a summary of the AER codes, J. Quant. Spectrosc. Ra., 91, 233–244, 2005.

Daniel, J. S., Solomon, S., Miller, H. L., Langford, A. O., Portmann, R. W., and Eubank, C. S.: Retrieving cloud information from passive measurements of
solar radiation absorbed by molecular oxygen and O_{2}−O_{2}, J. Geophys.
Res., 108, 4515, https://doi.org/10.1029/2002JD002994, 2003.

Dannenberg, R. B.: Interpolation error in waveform table lookup, in: Proceedings of the International Computer Music Conference, International Computer Music Association, San Francisco, 1998.

Davis, A. B. and Marshak, A.: Space–time characteristics of light transmitted through dense clouds: A Green's function analysis, J. Atmos. Sci., 59, 2713–2727, 2002.

Davis, A. B., Merlin, G., Cornet, C., Labonnote, L. C., Riédi, J., Ferlay, N., Dubuisson, P., Min, Q., Yang, Y., and Marshak, A.: Cloud information content in EPIC/DSCOVR's oxygen A-and B-band channels: An optimal estimation approach, J. Quant. Spectrosc. Ra., 216, 6–16, 2018a.

Davis, A. B., Ferlay, N., Libois, Q., Marshak, A., Yang, Y., and Min, Q.: Cloud information content in EPIC/DSCOVR's oxygen A-and B-band channels: A physics-based approach, J. Quant. Spectrosc. Ra., 220, 84–96, 2018b.

Duan, M., Min, Q., and Li, J.: A fast radiative transfer model for simulating high-resolution absorption bands, J. Geophys. Res., 110, D15201, https://doi.org/10.1029/2004JD005590, 2005.

Ferlay, N., Thieuleux, F., Cornet, C., Davis, A. B., Dubuisson, P., Ducos, F., Parol, F., Riédi, J., and Vanbauce, C.: Toward new inferences about cloud structures from multidirectional measurements in the oxygen A band: middle-of-cloud pressure and cloud geometrical thickness from POLDER-3/PARASOL, J. Appl. Meteorol. Clim., 49, 2492–2507, 2010.

Fischer, J. and Grassl, H.: Detection of cloud-top height from backscattered radiances within the oxygen A band. Part 1: Theoretical study, J. Appl. Meteorol., 30, 1245–1259, 1991.

Gastellu-Etchegorry, J. P., Gascon, F., and Esteve, P.: An interpolation procedure for generalizing a look-up table inversion method, Remote Sens. Environ., 87, 55–71, 2003.

Gelaro, R., McCarty, W., Suárez, M. J., Todling, R., Molod, A., Takacs, L., Randles, C. A., Darmenov, A., Bosilovich, M. G., Reichle, R., and Wargan, K.: 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, https://doi.org/10.5194/amt-11-359-2018, 2018.

Gordon, I. E., Rothman, L. S., Hill, C., Kochanov, R. V., Tan, Y., Bernath, P. F., Birk, M., Boudon, V., Campargue, A., Chance, K. V., and Drouin, B. J.: The HITRAN2016 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 203, 3–69, 2017.

Holdaway, D. and Yang, Y.: Study of the effect of temporal sampling frequency on DSCOVR observations using the GEOS-5 nature run results (Part II): Cloud Coverage, Remote Sens., 8, 431, https://doi.org/10.3390/rs8050431, 2016.

Ishimaru, A.: Wave propagation and scattering in random media, Wiley-IEEE-Press, New York, 1999.

Irvine, W. M.: The formation of absorption bands and the distribution of photon optical paths in a scattering atmosphere, B. Astron. I. Neth., 17, 266–279, 1964.

Ivanov, V. V. and Gutshabash, S. D.: Propagation of brightness wave in an optically thick atmosphere, Phys. Atmos. Okeana, 10, 851–863, 1974.

Koelemeijer, R. B. A., Stammes, P., Hovenier, J. W., and Haan, J. D.: A fast method for retrieval of cloud parameters using oxygen A band measurements from the Global Ozone Monitoring Experiment, J. Geophys. Res., 106, 3475–3490, 2001.

Kokhanovsky, A. A. and Rozanov, V. V.: The physical parameterization of the top-of-atmosphere reflection function for a cloudy atmosphere–underlying surface system: the oxygen A-band case study, J. Quant. Spectrosc. Ra., 85, 35–55, https://doi.org/10.1016/S0022-4073(03)00193-6, 2004.

Kokhanovsky, A. A., Rozanov, V. V., Zege, E. P., Bovesmann, H., and Burrows, J. P.: A semi analytical cloud retrieval algorithm usinfg backscattered radiation in 0.4–2.4 µm spectral region, J. Geophys. Res., 108, 4008, https://doi.org/10.1029/2001JD001543, 2003.

Kuze, A. and Chance, K. V.: Analysis of cloud top height and cloud coverage from satellites using the O_{2}*A* and *B* bands, J. Geophys. Res., 99, 14481–14491, 1994.

Lelli, L., Kokhanovsky, A. A., Rozanov, V. V., Vountas, M., and Burrows, J. P.: Linear trends in cloud top height from passive observations in the oxygen A-band, Atmos. Chem. Phys., 14, 5679–5692, https://doi.org/10.5194/acp-14-5679-2014, 2014.

Lelli, L., Kokhanovsky, A. A., Rozanov, V. V., Vountas, M., Sayer, A. M., and Burrows, J. P.: Seven years of global retrieval of cloud properties using space-borne data of GOME, Atmos. Meas. Tech., 5, 1551–1570, https://doi.org/10.5194/amt-5-1551-2012, 2012.

Loyola, D. G., Gimeno García, S., Lutz, R., Argyrouli, A., Romahn, F., Spurr, R. J. D., Pedergnana, M., Doicu, A., Molina García, V., and Schüssler, O.: The operational cloud retrieval algorithms from TROPOMI on board Sentinel-5 Precursor, Atmos. Meas. Tech., 11, 409–427, https://doi.org/10.5194/amt-11-409-2018, 2018.

Marshak, A. and Davis, A. (Eds.).: 3D radiative transfer in cloudy atmospheres, Springer Science & Business Media, Springer, Berlin, Heidelberg, New York, 2005.

Marshak, A., Herman, J., Adam, S., Carn, S., Cede, A., Geogdzhayev, I., Huang, D., Huang, L. K., Knyazikhin, Y., Kowalewski, M., and Krotkov, N.: Earth observations from DSCOVR EPIC instrument, B. Am. Meteorol. Soc., 99, 1829–1850, https://doi.org/10.1175/BAMS-D-17-0223.1, 2018.

Meyer, K., Yang, Y., and Platnick, S.: Uncertainties in cloud phase and optical thickness retrievals from the Earth Polychromatic Imaging Camera (EPIC), Atmos. Meas. Tech., 9, 1785–1797, https://doi.org/10.5194/amt-9-1785-2016, 2016.

Min, Q. and Harrison, L. C.: Retrieval of atmospheric optical depth profiles
from downward-looking high-resolution O_{2} A-band measurements: Optically thin conditions, J. Atmos. Sci., 61, 2469–2477, 2004.

Min, Q., Yin, B., Li, S., Berndt, J., Harrison, L., Joseph, E., Duan, M., and Kiedron, P.: A high-resolution oxygen A-band spectrometer (HABS) and its radiation closure, Atmos. Meas. Tech., 7, 1711–1722, https://doi.org/10.5194/amt-7-1711-2014, 2014.

NASA ASDC: CAL_LID_L2_05kmCLay-Standard-V4-20, available at: https://www-calipso.larc.nasa.gov/products/inventory.php, last access: 14 July 2017.

NASA LARC ASDC DAAC (Langley Research Center's Atmospheric Science Data Center Distributed Active Archive Centers): DSCOVR EPIC Level 1B Version 2, https://doi.org/10.5067/EPIC/DSCOVR/L1B.002, 2018a.

NASA LARC ASDC DAAC (Langley Research Center's Atmospheric Science Data Center Distributed Active Archive Centers): DSCOVR EPIC Cloud Products, https://doi.org/10.5067/EPIC/DSCOVR/L2_Cloud_01, 2018b.

NASA LARC: Data product: LIDAR Level2 Version 4.20 5-km Cloud Layer, available at: https://subset.larc.nasa.gov/calipso/index.php, last access: 29 September 2020.

O'Brien, D. M. and Mitchell, R. M.: Error estimates for retrieval of cloud-top pressure using absorption in the A band of oxygen, J. Appl. Meteorol., 31, 1179–1192, 1992.

Pandey, P., De Ridder, K., Gillotay, D., and van Lipzig, N. P. M.: Estimating cloud optical thickness and associated surface UV irradiance from SEVIRI by implementing a semi-analytical cloud retrieval algorithm, Atmos. Chem. Phys., 12, 7961–7975, https://doi.org/10.5194/acp-12-7961-2012, 2012.

Preusker, R. and Lindstrot, R.: Remote Sensing of Cloud-Top Pressure Using Moderately Resolved Measurements within the Oxygen A Band-A Sensitivity Study, J. Appl. Meteorol. Clim., 48, 1562–1574, 2009.

Richardson, M. and Stephens, G. L.: Information content of OCO-2 oxygen A-band channels for retrieving marine liquid cloud properties, Atmos. Meas. Tech., 11, 1515–1528, https://doi.org/10.5194/amt-11-1515-2018, 2018.

Rozanov, V. V. and Kokhanovsky, A. A.: Semianalytical cloud retrieval algorithm as applied to the cloud top altitude and the cloud geometrical thickness determination from top-of-atmosphere reflectance measurements in the oxygen A band, J. Geophys. Res., 109, 4070, https://doi.org/10.1029/2003JD004104, 2004.

Schuessler, O., Rodriguez, D. G. L., Doicu, A., and Spurr, R.: Information Content in the Oxygen A-Band for the Retrieval of Macrophysical Cloud Parameters, IEEE T. Geosci. Remote, 52, 3246–3255, 2013.

Seager, S., Turner, E. L., Schafer, J., and Ford, E. B.: Vegetation's red edge: a possible spectroscopic biosignature of extraterrestrial plants, Astrobiology, 5, 372–390, 2005.

Stamnes, K., Tsay, S. C., Wiscombe, W., and Jayaweera, K.: Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media, Appl. Optics, 27, 2502–2509, 1988.

Thomas, G. E. and Stamnes, K.: Radiative transfer in the atmosphere and ocean, Cambridge University Press, Cambridge, 2002.

Tilstra, L. G., Wang, P., and Stammes, P.: Surface reflectivity climatologies from UV to NIR determined from Earth observations by GOME-2 and SCIAMACHY, J. Geophys. Res., 122, 4084–4111, https://doi.org/10.1002/2016JD025940, 2017.

Tropospheric Emission Monitoring Internet Service (TEMIS): GOME-2 surface LER, available at: http://temis.nl/surface/gome2_ler/databases/, last access: 13 September 2017.

Van de Hulst, H. C.: Multiple Light Scattering: Tables, Formulas, and Applications, Academic Press, 299 pp., 1980.

Van de Hulst, H. C.: Multiple light scattering: tables, formulas, and applications, Elsevier, Burlington Elsevier Science, 2012.

Vaughan, M. A., Young, S. A., Winker, D. M., Powell, K. A., Omar, A. H., Liu, Z., Hu, Y., and Hostetler, C. A.: Fully automated analysis of space-based lidar data: An overview of the CALIPSO retrieval algorithms and data products, in: Laser radar techniques for atmospheric sensing, International Society for Optics and Photonics, 5575, 16–30, https://doi.org/10.1117/12.572024, 2004.

Yamamoto, G. and Wark, D. Q.: Discussion of letter by 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.

Yang, Y., Meyer, K., Wind, G., Zhou, Y., Marshak, A., Platnick, S., Min, Q., Davis, A. B., Joiner, J., Vasilkov, A., Duda, D., and Su, W.: Cloud products from the Earth Polychromatic Imaging Camera (EPIC): algorithms and initial evaluation, Atmos. Meas. Tech., 12, 2019–2031, https://doi.org/10.5194/amt-12-2019-2019, 2019.