Abstract. This paper presents the operational cloud retrieval algorithms for the TROPOspheric Monitoring Instrument (TROPOMI) on board the European Space Agency Sentinel-5 Precursor (S5P) mission scheduled for launch in 2017.Two algorithms working in tandem are used for retrieving cloud properties: OCRA (Optical Cloud Recognition Algorithm) and ROCINN (Retrieval of Cloud Information using Neural Networks). OCRA retrieves the cloud fraction using TROPOMI measurements in the ultraviolet (UV) and visible (VIS) spectral regions, and ROCINN retrieves the cloud top height (pressure) and optical thickness (albedo) using TROPOMI measurements in and around the oxygen A-band in the near infrared (NIR).Cloud parameters from TROPOMI/S5P will be used not only for enhancing the accuracy of trace gas retrievals but also for extending the satellite data record of cloud information derived from oxygen A-band measurements, a record initiated with the Global Ozone Monitoring Experiment (GOME) on board the second European Remote-Sensing Satellite (ERS-2) over 20 years ago.The OCRA and ROCINN algorithms are integrated in the S5P operational processor UPAS (Universal Processor for UV/VIS/NIR Atmospheric Spectrometers), and we present here UPAS cloud results using the Ozone Monitoring Instrument (OMI) and GOME-2 measurements. In addition, we examine anticipated challenges for the TROPOMI/S5P cloud retrieval algorithms, and we discuss the future validation needs for OCRA and ROCINN.

In this paper we analyze the accuracy and efficiency of several radiative transfer models for inferring cloud parameters from radiances measured by the Earth Polychromatic Imaging Camera (EPIC) on board the Deep Space Climate Observatory (DSCOVR). The radiative transfer models are the exact discrete ordinate and matrix operator methods with matrix exponential, and the approximate asymptotic and equivalent Lambertian cloud models. To deal with the computationally expensive radiative transfer calculations, several acceleration techniques such as, for example, the telescoping technique, the method of false discrete ordinate, the correlated k-distribution method and the principal component analysis (PCA) are used. We found that, for the EPIC oxygen A-band absorption channel at 764 nm, the exact models using the correlated k-distribution in conjunction with PCA yield an accuracy better than 1.5 % and a computation time of 18 s for radiance calculations at 5 viewing zenith angles.

The new generation of atmospheric composition sensors such as TROPOMI is capable of providing spectra of high spatial and spectral resolution. To process this vast amount of spectral information, fast radiative transfer models (RTMs) are required. In this regard, we analyzed the efficiency of two acceleration techniques based on the principal component analysis (PCA) for simulating the Hartley-Huggins band spectra. In the first one, the PCA is used to map the data set of optical properties of the atmosphere to a lower-dimensional subspace, in which the correction function for an approximate but fast RTM is derived. The second technique is based on the dimensionality reduction of the data set of spectral radiances. Once the empirical orthogonal functions are found, the whole spectrum can be reconstructed by performing radiative transfer computations only for a specific subset of spectral points. We considered a clear-sky atmosphere where the optical properties are defined by Rayleigh scattering and trace gas absorption. Clouds can be integrated into the model as Lambertian reflectors. High computational performance is achieved by combining both techniques without losing accuracy. We found that for the Hartley-Huggins band, the combined use of these techniques yields an accuracy better than 0.05% while the speedup factor is about 20. This innovative combination of both PCA-based techniques can be applied in future works as an efficient approach for simulating the spectral radiances in other spectral regions.

This paper outlines the matrix exponential description of radiative transfer. The eigendecomposition method which serves as a basis for computing the matrix exponential and for representing the solution in a discrete ordinate setting is considered. The mathematical equivalence of the discrete ordinate method, the matrix operator method, and the matrix Riccati equations method is proved rigorously by means of the matrix exponential formalism. For optically thin layers approximate solution methods relying on the Padé and Taylor series approximations to the matrix exponential, as well as on the matrix Riccati equations are presented. For optically thick layers, the asymptotic theory with higher-order corrections is derived, and parametrizations of the asymptotic functions and constants for a water-cloud model with a Gamma size distribution are obtained.

Current atmospheric composition sensors provide a large amount of high spectral resolution data. The accurate processing of this data employs time-consuming line-by-line (LBL) radiative transfer models (RTMs). In this paper, we describe a method to accelerate hyperspectral radiative transfer models based on the clustering of the spectral radiances computed with a low-stream RTM and the regression analysis performed for the low-stream and multi-stream RTMs within each cluster. This approach, which we refer to as the Cluster Low-Streams Regression (CLSR) method, is applied for computing the radiance spectra in the O 2 A-band at 760 nm and the CO 2 band at 1610 nm for five atmospheric scenarios. The CLSR method is also compared with the principal component analysis (PCA)-based RTM, showing an improvement in terms of accuracy and computational performance over PCA-based RTMs. As low-stream models, the two-stream and the single-scattering RTMs are considered. We show that the error of this approach is modulated by the optical thickness of the atmosphere. Nevertheless, the CLSR method provides a performance enhancement of almost two orders of magnitude compared to the LBL model, while the error of the technique is below 0.1% for both bands. simulations. In [7,8], the transmission function for a given spectral interval is fitted by a sum of exponentials, while the corresponding fitting coefficients are computed from a reduced number of monochromatic computations. A similar approach is described in Moncet et al. [9], where the fitting weights and the most representative wavelengths are chosen appropriately.The state-of-the-art hyperspectral RTMs employ dimensionality reduction techniques such as principal component analysis (PCA). In [10,11], PCA is applied to the spectral radiance data to establish a set of empirical orthogonal functions (EOFs), so that an arbitrary spectrum at full spectral resolution can be reconstructed as a weighted sum of EOFs. The weights are found by performing monochromatic simulations at a reduced number of wavelengths. To accelerate the computations in the O 2 A-band, Natraj et al. [12] proposed a fundamentally different PCA-based radiative transfer model, in which the dimensionality of the optical properties data is reduced. A two-stream radiative transfer model was used as an approximate model, and the dependency of the corresponding correction factor on the optical parameters was modeled by a second-order Taylor expansion about the mean value of the optical parameters in the reduced optical data space. This approach was extended to other dimensionality reduction techniques [13] and spectral ranges [14-16]; moreover, it was implemented in conjunction with PCA for spectral radiances [17] and with the k-distribution method [18]. The errors of these approaches are usually below 0.1% for the spectral radiances, while the performance enhancement may reach several orders of magnitude depending on the spectral region and the required level of accuracy.In Efremenko et al. [19] it was shown that, aft...