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.

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

This article is the second in a series of studies investigating the benefits of multispectral measurements to improve the atmospheric parameter retrievals. In the first paper, we presented an information content (IC) analysis from the thermal infrared (TIR) and shortwave infrared (SWIR) bands of Thermal And Near infrared Sensor for carbon Observations–Fourier Transform Spectrometer (TANSO-FTS) instrument dedicated to greenhouse gas retrieval in clear sky conditions. This second paper presents the potential of the spectral synergy from TIR to visible for aerosol characterization, and their impact on the retrieved CO2 and CH4 column concentrations. The IC is then used to determine the most informative spectral channels for the simultaneous retrieval of greenhouse gas total columns and aerosol parameters. The results show that a channel selection spanning the four bands can improve the computation time and retrieval accuracy. Therefore, the spectral synergy allows obtaining up to almost seven different aerosol parameters, which is comparable to the most informative dedicated instruments. Moreover, a channel selection from the TIR to visible bands allows retrieving CO2 and CH4 total columns simultaneously in the presence of one aerosol layer with a similar accuracy to using all channels together to retrieve each gas separately in clear sky conditions.

The Greenhouse gases Observing SATellite (GOSAT) mission, and in particular the Thermal And Near infrared Sensor for carbon Observations–Fourier Transform Spectrometer (TANSO-FTS) instrument, has the advantage of being able to measure simultaneously the same field of view in different spectral ranges with a high spectral resolution. These features allow studying the benefits of using multispectral measurements to improve the CO2 and CH4 retrievals. In order to quantify the impact of the spectral synergy on the retrieval accuracy, we performed an information content (IC) analysis from simulated spectra corresponding to the three infrared bands of TANSO-FTS. The advantages and limitations of using thermal and shortwave infrared simultaneously are discussed according to surface type and state vector composition. The IC is then used to determine the most informative spectral channels for the simultaneous retrieval of CO2 and CH4. The results show that a channel selection spanning the three infrared bands can improve the computation time and retrieval accuracy. Therefore, a selection of less than 700 channels from the thermal infrared (TIR) and shortwave infrared (SWIR) bands allows retrieving CO2 and CH4 simultaneously with a similar accuracy to using all channels together to retrieve each gas separately.

Most of the current aerosol retrievals from passive sensors are restricted to cloud-free scenes, which strongly reduces our ability to monitor the aerosol properties at a global scale and to estimate their radiative forcing. The presence of aerosol above clouds (AAC) affects the polarized light reflected by the cloud layer, as shown by the spaceborne measurements provided by the POlarization and Directionality of Earth Reflectances (POLDER) instrument on the PARASOL satellite. In a previous work, a first retrieval method was developed for AAC scenes and evaluated for biomass-burning aerosols transported over stratocumulus clouds. The method was restricted to the use of observations acquired at forward scattering angles (90–120°) where polarized measurements are highly sensitive to fine-mode particle scattering. Non-spherical particles in the coarse mode, such as mineral dust particles, do not much polarize light and cannot be handled with this method. In this paper, we present new developments that allow retrieving also the properties of mineral dust particles above clouds. These particles do not much polarize light but strongly reduce the polarized cloud bow generated by the liquid cloud layer beneath and observed for scattering angles around 140°. The spectral attenuation can be used to qualitatively identify the nature of the particles (i.e. accumulation mode versus coarse mode, i.e. mineral dust particles versus biomass-burning aerosols), whereas the magnitude of the attenuation is related to the optical thickness of the aerosol layer. We also use the polarized measurements acquired in the cloud bow to improve the retrieval of both the biomass-burning aerosol properties and the cloud microphysical properties. We provide accurate polarized radiance calculations for AAC scenes and evaluate the contribution of the POLDER polarization measurements for the simultaneous retrieval of the aerosol and cloud properties. We investigate various scenes with mineral dust particles and biomass-burning aerosols above clouds. For clouds, our results confirm that the droplet size distribution is narrow in high-latitude ocean regions and that the droplet effective radii retrieved from both polarization measurements and from total radiance measurements are generally close for AAC scenes (departures smaller than 2 μm). We found that the magnitude of the primary cloud bow cannot be accurately estimated with a plane parallel transfer radiative code. The errors for the modeling of the polarized cloud bow are between 4 and 8% for homogenous cloudy scenes, as shown by a 3-D radiative transfer code. These effects only weakly impact the retrieval of the Aerosol Optical Thickness (AOT) performed with a mineral dust particle model for which the microphysical properties are entirely known (relative error smaller than 6%). We show that the POLDER polarization measurements allow retrieving the AOT, the fine-mode particle size, the Ångström exponent and the fraction of spherical particles. However, the complex refractive index and the coarse-mode particle size cannot be accurately retrieved with the present polarization measurements. Our complete and accurate algorithm cannot be applied to process large amounts of data, so a simpler algorithm was developed to retrieve the AOT and the Ångström exponent above clouds in an operational way. Illustrations are provided for July–August 2008 near the African coast. Large mean AOTs above clouds at 0.865 μm (>0.3) are retrieved for oceanic regions near the coasts of South Africa that correspond to biomass-burning aerosols, whereas even larger mean AOTs above clouds for mineral dust particles (>0.6) are also retrieved near the coasts of Senegal. For these regions and time period, the direct AAC radiative forcing is likely to be significant. The final aim of this work is the global monitoring of the AAC properties and the estimation of the direct aerosol radiative forcing in cloudy scenes.