Search

In remote sensing applications, clouds are generally characterized by two properties: cloud optical thickness (COT) and effective radius of water–ice particles (Reff), as well as additionally by geometric properties when specific information is available. Most of the current operational passive remote sensing algorithms use a mono-angular bispectral method to retrieve COT and Reff. They are based on pre-computed lookup tables while assuming a homogeneous plane-parallel cloud layer. In this work, we use the formalism of the optimal estimation method, applied to airborne near-infrared high-resolution multi-angular measurements, to retrieve COT and Reff as well as the corresponding uncertainties related to the measurement errors, the non-retrieved parameters, and the cloud model assumptions. The measurements used were acquired by the airborne radiometer OSIRIS (Observing System Including PolaRization in the Solar Infrared Spectrum), developed by the Laboratoire d'Optique Atmosphérique. It provides multi-angular measurements at a resolution of tens of meters, which is very suitable for refining our knowledge of cloud properties and their high spatial variability. OSIRIS is based on the POLDER (POlarization and Directionality of the Earth's Reflectances) concept as a prototype of the future 3MI (Multi-viewing Multi-channel Multi-polarization Imager) planned to be launched on the EUMETSAT-ESA MetOp-SG platform in 2024. The approach used allows the exploitation of all the angular information available for each pixel to overcome the radiance angular effects. More consistent cloud properties with lower uncertainty compared to operational mono-directional retrieval methods (traditional bispectral method) are then obtained. The framework of the optimal estimation method also provides the possibility to estimate uncertainties of different sources. Three types of errors were evaluated: (1) errors related to measurement uncertainties, which reach 6 % and 12 % for COT and Reff, respectively, (2) errors related to an incorrect estimation of the ancillary data that remain below 0.5 %, and (3) errors related to the simplified cloud physical model assuming independent pixel approximation. We show that not considering the in-cloud heterogeneous vertical profiles and the 3D radiative transfer effects leads to an average uncertainty of 5 % and 4 % for COT and 13 % and 9 % for Reff.

The features detected in monolayer atmospheric columns sounded by the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and classified as cloud or aerosol layers by the CALIOP version 4 (V4) cloud and aerosol discrimination (CAD) algorithm are reassessed using perfectly collocated brightness temperatures measured by the Imaging Infrared Radiometer (IIR) aboard the same satellite. Using the IIR's three wavelength measurements of layers that are confidently classified by the CALIOP CAD algorithm, we calculate two-dimensional (2-D) probability distribution functions (PDFs) of IIR brightness temperature differences (BTDs) for different cloud and aerosol types. We then compare these PDFs with 1-D radiative transfer simulations for ice and water clouds and dust and marine aerosols. Using these IIR 2-D BTD signature PDFs, we develop and deploy a new IIR-based CAD algorithm and compare the classifications obtained to the results reported by the CALIOP-only V4 CAD algorithm. IIR observations are shown to be able to identify clouds with a good accuracy. The IIR cloud identifications agree very well with layers classified as confident clouds by the V4 CAD algorithm (88 %). More importantly, simultaneous use of IIR information reduces the ambiguity in a notable fraction of “not confident” V4 cloud classifications. 28 % and 14 % of the ambiguous V4 cloud classifications are reclassified more appropriately as confident cloud layers through the use of the IIR observations in the tropics and in the midlatitudes, respectively. IIR observations are of relatively little help in deriving high-confidence classifications for most aerosols, as the low altitudes and small optical depths of aerosol layers yield IIR signatures that are similar to those of clear skies. However, misclassifications of aerosol layers, such as dense dust or elevated smoke layers, by the V4 CAD algorithm can be corrected to cloud layer classification by including IIR information. 10 %, 16 %, and 6 % of the ambiguous V4 dust, polluted dust, and tropospheric elevated smoke, respectively, are found to be misclassified cloud layers by the IIR measurements.

The algorithm presented in this paper was developed to retrieve ice water content (IWC) profiles in cirrus clouds. It is based on optimal estimation theory and combines ground-based visible lidar and thermal infrared (TIR) radiometer measurements in a common retrieval framework in order to retrieve profiles of IWC together with a correction factor for the backscatter intensity of cirrus cloud particles. As a first step, we introduce a method to retrieve extinction and IWC profiles in cirrus clouds from the lidar measurements alone and demonstrate the shortcomings of this approach due to the backscatter-to-extinction ambiguity. As a second step, we show that TIR radiances constrain the backscattering of the ice crystals at the visible lidar wavelength by constraining the ice water path (IWP) and hence the IWC, which is linked to the optical properties of the ice crystals via a realistic bulk ice microphysical model. The scattering phase function obtained from the microphysical model is flat around the backscatter direction (i.e., there is no backscatter peak). We show that using this flat backscattering phase function to define the backscatter-to-extinction ratio of the ice crystals in the retrievals with the lidar-only algorithm results in an overestimation of the IWC, which is inconsistent with the TIR radiometer measurements. Hence, a synergy algorithm was developed that combines the attenuated backscatter profiles measured by the lidar and the measurements of TIR radiances in a common optimal estimation framework to retrieve the IWC profile together with a correction factor for the phase function of the bulk ice crystals in the backscattering direction. We show that this approach yields consistent lidar and TIR results. The resulting lidar ratios for cirrus clouds are found to be consistent with previous independent studies.

Accurate calculations of shortwave reflectances in clear-sky aerosol-laden atmospheres are necessary for various applications in atmospheric sciences. However, computational cost becomes increasingly important for some applications such as data assimilation of top-of-atmosphere reflectances in models of atmospheric composition. This study aims to provide a benchmark that can help in assessing these two requirements in combination. We describe a protocol and input data for 44 080 cases involving various solar and viewing geometries, four different surfaces (one oceanic bidirectional reflectance function and three albedo values for a Lambertian surface), eight aerosol optical depths, five wavelengths, and four aerosol types. We first consider two models relying on the discrete ordinate method: VLIDORT (in vector and scalar configurations) and DISORT (scalar configuration only). We use VLIDORT in its vector configuration as a reference model and quantify the loss of accuracy due to (i) neglecting the effect of polarization in DISORT and VLIDORT (scalar) models and (ii) decreasing the number of streams in DISORT. We further test two other models: the 6SV2 model, relying on the successive orders of scattering method, and Forward-Lobe Two-Stream Radiance Model (FLOTSAM), a new model under development by two of the authors. Typical mean fractional errors of 2.8 % and 2.4 % for 6SV2 and FLOTSAM are found, respectively. Computational cost depends on the input parameters but also on the code implementation and application as some models solve the radiative transfer equations for a range of geometries while others do not. All necessary input and output data are provided as a Supplement as a potential resource for interested developers and users of radiative transfer models.

Simulations of total and polarized cloud reflectance angular signatures such as the ones measured by the multi-angular and polarized radiometer POLDER3/PARASOL are used to evaluate cloud heterogeneity effects on cloud parameter retrievals. Effects on optical thickness, albedo, effective radius and variance of the cloud droplet size distribution and aerosol parameters above cloud are analyzed. Three different clouds that have the same mean optical thicknesses were generated: the first with a flat top, the second with a bumpy top and the last with a fractional cloud cover. At small scale (50 m), for oblique solar incidence, the illumination effects lead to higher total but also polarized reflectances. The polarized reflectances even reach values that cannot be predicted by the 1-D homogeneous cloud assumption. At the POLDER scale (7 km × 7 km), the angular signature is modified by a combination of the plane–parallel bias and the shadowing and illumination effects. In order to quantify effects of cloud heterogeneity on operational products, we ran the POLDER operational algorithms on the simulated reflectances to retrieve the cloud optical thickness and albedo. Results show that the cloud optical thickness is greatly affected: biases can reach up to −70, −50 or +40 % for backward, nadir and forward viewing directions, respectively. Concerning the albedo of the cloudy scenes, the errors are smaller, between −4.7 % for solar incidence angle of 20° and up to about +8 % for solar incidence angle of 60°. We also tested the heterogeneity effects on new algorithms that allow retrieving cloud droplet size distribution and cloud top pressures and also aerosol above clouds. Contrary to the bi-spectral method, the retrieved cloud droplet size parameters are not significantly affected by the cloud heterogeneity, which proves to be a great advantage of using polarized measurements. However, the cloud top pressure obtained from molecular scattering in the forward direction can be biased up to about 60 hPa (around 550 m). Concerning the aerosol optical thickness (AOT) above cloud, the results are different depending on the available angular information. Above the fractional cloud, when only side scattering angles between 100 and 130° are available, the AOT is underestimated because of the plane–parallel bias. However, for solar zenith angle of 60° it is overestimated because the polarized reflectances are increased in forward directions.

Information content analyses on cloud top altitude (CTOP) and geometrical thickness (CGT) from multi-angular A-band measurements in the case of monolayer homogeneous clouds are conducted. In the framework of future multi-angular radiometer development, we compared the potential performances of the 3MI (Multi-viewing, Multi-channel and Multi-polarization Imaging) instrument developed by EUMETSAT, which is an extension of POLDER/PARASOL instrument and MSPI (Multiangle SpectroPolarimetric Imager) developed by NASA's Jet Propulsion Laboratory. Quantitative information content estimates were realized for thin, moderately opaque and opaque clouds for different surface albedo and viewing geometry configurations. Analyses show that retrieval of CTOP is possible with a high accuracy in most of the cases investigated. Retrieval of CGT is also possible for optically thick clouds above a black surface, at least when CGT > 1–2 km and for thin clouds for CGT > 2–3 km. However, for intermediate optical thicknesses (COT ≃ 4), we show that the retrieval of CGT is not simultaneously possible with CTOP. A comparison between 3MI and MSPI shows a higher information content for MSPI's measurements, traceable to a thinner filter inside the oxygen A-band, yielding higher signal-to-noise ratio for absorption estimation. Cases of cloud scenes above bright surfaces are more complex but it is shown that the retrieval of CTOP remains possible in almost all situations while the information content on CGT appears to be insufficient in many cases, particularly for COT

In this study, various ice particle habits are investigated in conjunction with inferring the optical properties of ice clouds for use in the Global Change Observation Mission-Climate (GCOM-C) satellite programme. We develop a database of the single-scattering properties of five ice habit models: plates, columns, droxtals, bullet rosettes, and Voronoi. The database is based on the specification of the Second Generation Global Imager (SGLI) sensor on board the GCOM-C satellite, which is scheduled to be launched in 2017 by the Japan Aerospace Exploration Agency. A combination of the finite-difference time-domain method, the geometric optics integral equation technique, and the geometric optics method is applied to compute the single-scattering properties of the selected ice particle habits at 36 wavelengths, from the visible to the infrared spectral regions. This covers the SGLI channels for the size parameter, which is defined as a single-particle radius of an equivalent volume sphere, ranging between 6 and 9000 µm. The database includes the extinction efficiency, absorption efficiency, average geometrical cross section, single-scattering albedo, asymmetry factor, size parameter of a volume-equivalent sphere, maximum distance from the centre of mass, particle volume, and six nonzero elements of the scattering phase matrix. The characteristics of calculated extinction efficiency, single-scattering albedo, and asymmetry factor of the five ice particle habits are compared. Furthermore, size-integrated bulk scattering properties for the five ice particle habit models are calculated from the single-scattering database and microphysical data. Using the five ice particle habit models, the optical thickness and spherical albedo of ice clouds are retrieved from the Polarization and Directionality of the Earth's Reflectances-3 (POLDER-3) measurements, recorded on board the Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar (PARASOL) satellite. The optimal ice particle habit for retrieving the SGLI ice cloud properties is investigated by adopting the spherical albedo difference (SAD) method. It is found that the SAD is distributed stably due to the scattering angle increases for bullet rosettes with an effective diameter (Deff) of 10 µm and Voronoi particles with Deff values of 10, 60, and 100 µm. It is confirmed that the SAD of small bullet-rosette particles and all sizes of Voronoi particles has a low angular dependence, indicating that a combination of the bullet-rosette and Voronoi models is sufficient for retrieval of the ice cloud's spherical albedo and optical thickness as effective habit models for the SGLI sensor. Finally, SAD analysis based on the Voronoi habit model with moderate particle size (Deff = 60 µm) is compared with the conventional general habit mixture model, inhomogeneous hexagonal monocrystal model, five-plate aggregate model, and ensemble ice particle model. The Voronoi habit model is found to have an effect similar to that found in some conventional models for the retrieval of ice cloud properties from space-borne radiometric observations.

Moderate Resolution Imaging Spectroradiometer (MODIS) retrieves cloud droplet effective radius ( r e ) and optical thickness ( τ ) by projecting observed cloud reflectances onto a precomputed look-up table (LUT). When observations fall outside of the LUT, the retrieval is considered "failed" because no combination of τ and r e within the LUT can explain the observed cloud reflectances. In this study, the frequency and potential causes of failed MODIS retrievals for marine liquid phase (MLP) clouds are analyzed based on 1 year of Aqua MODIS Collection 6 products and collocated CALIOP and CloudSat observations. The retrieval based on the 0.86 µm and 2.1 µm MODIS channel combination has an overall failure rate of about 16% (10% for the 0.86 µm and 3.7 µm combination). The failure rates are lower over stratocumulus regimes and higher over the broken trade wind cumulus regimes. The leading type of failure is the " r e too large" failure accounting for 60%-85% of all failed retrievals. The rest is mostly due to the " r e too small" or τ retrieval failures. Enhanced retrieval failure rates are found when MLP cloud pixels are partially cloudy or have high subpixel inhomogeneity, are located at special Sun-satellite viewing geometries such as sunglint, large viewing or solar zenith angles, or cloudbow and glory angles, or are subject to cloud masking, cloud overlapping, and/or cloud phase retrieval issues. The majority (more than 84%) of failed retrievals along the CALIPSO track can be attributed to at least one or more of these potential reasons. The collocated CloudSat radar reflectivity observations reveal that the remaining failed retrievals are often precipitating. It remains an open question whether the extremely large r e values observed in these clouds are the consequence of true cloud microphysics or still due to artifacts not included in this study.