The cloud droplet number concentration (Nd) is of central interest to improve the understanding of cloud physics and for quantifying the effective radiative forcing by aerosol‐cloud interactions. ...Current standard satellite retrievals do not operationally provide Nd, but it can be inferred from retrievals of cloud optical depth (τc) cloud droplet effective radius (re) and cloud top temperature. This review summarizes issues with this approach and quantifies uncertainties. A total relative uncertainty of 78% is inferred for pixel‐level retrievals for relatively homogeneous, optically thick and unobscured stratiform clouds with favorable viewing geometry. The uncertainty is even greater if these conditions are not met. For averages over 1° ×1° regions the uncertainty is reduced to 54% assuming random errors for instrument uncertainties. In contrast, the few evaluation studies against reference in situ observations suggest much better accuracy with little variability in the bias. More such studies are required for a better error characterization. Nd uncertainty is dominated by errors in re, and therefore, improvements in re retrievals would greatly improve the quality of the Nd retrievals. Recommendations are made for how this might be achieved. Some existing Nd data sets are compared and discussed, and best practices for the use of Nd data from current passive instruments (e.g., filtering criteria) are recommended. Emerging alternative Nd estimates are also considered. First, new ideas to use additional information from existing and upcoming spaceborne instruments are discussed, and second, approaches using high‐quality ground‐based observations are examined.
Plain Language Summary
Clouds have a very large influence on weather and climate. It is thus a prime task for satellite‐ and ground‐based observations to measure clouds. For satellites and many other instruments this is done by remote sensing—radiation is measured, and knowledge about clouds is inferred. Liquid water clouds consist of numerous droplets of order of 10 μm in size. A key quantity that describes clouds is the number of droplets in a given volume or the droplet number concentration. However, satellite observations of droplet number concentration are only emerging, and the quality of these observations is poorly known. This review fulfills two tasks, namely, (1) to quantify how uncertain the current way to observe droplet number concentrations from satellite is and (2) to propose ways toward better approaches. It is concluded that the current way to obtain cloud droplet number concentration works for homogeneous stratus and stratocumulus clouds, with, however, a substantial error of around 50%. For cumulus clouds the observations are substantially worse. New avenues that are proposed for a better estimate of cloud droplet concentration exploit instruments that emit light (lidar) or microwaves (radar), and measure the reflected signal, or explore the polarization of light induced by clouds.
Key Points
Satellite cloud droplet concentration uncertainties of 78% for pixel‐level retrievals and 54% for 1 by 1 degree retrievals are estimated
The effective radius retrieval is the most important aspect for improvement, and more in situ evaluation is needed
Potential improvements using passive and active satellite, and ground‐based instruments are discussed
We present a generalization of the binary-value Markovian model previously used for statistical characterization of cloud masks to a continuous-value model describing 1D fields of cloud optical ...thickness (COT). This model has simple functional expressions and is specified by four parameters: the cloud fraction, the autocorrelation (scale) length, and the two parameters of the normalized probability density function of (non-zero) COT values (this PDF is assumed to have gamma-distribution form). Cloud masks derived from this model by separation between the values above and below some threshold in COT appear to have the same statistical properties as in binary-value model described in our previous publications. We demonstrate the ability of our model to generate examples of various cloud-field types by using it to statistically imitate actual cloud observations made by the Research Scanning Polarimeter (RSP) during two field experiments.
We present, for the first time, a quantitative retrieval error-propagation study for a bistatic high spectral resolution lidar (HSRL) system intended for detailed quasi-global monitoring of aerosol ...properties from space. Our results demonstrate that supplementing a conventional monostatic HSRL with an additional receiver flown in formation at a scattering angle close to 165 degrees dramatically increases the information content of the measurements and allows for a sufficiently accurate characterization of tropospheric aerosols. We conclude that a bistatic HSRL system would far exceed the capabilities of currently flown or planned orbital instruments in monitoring global aerosol effects on the environment and on the Earth's climate. We also demonstrate how the commonly used a priori 'regularization' methodology can artificially reduce the propagated uncertainties and can thereby be misleading as to the real retrieval capabilities of a measurement system.
We present an algorithm for the retrieval of cloud droplet size distribution parameters (effective radius and variance) from the Research Scanning Polarimeter (RSP) measurements. The RSP is an ...airborne prototype for the Aerosol Polarimetery Sensor (APS), which was on-board of the NASA Glory satellite. This instrument measures both polarized and total reflectance in 9 spectral channels with central wavelengths ranging from 410 to 2260nm. The cloud droplet size retrievals use the polarized reflectance in the scattering angle range between 135° and 165°, where they exhibit the sharply defined structure known as the rain- or cloud-bow. The shape of the rainbow is determined mainly by the single scattering properties of cloud particles. This significantly simplifies both forward modeling and inversions, while also substantially reducing uncertainties caused by the aerosol loading and possible presence of undetected clouds nearby. In this study we present the accuracy evaluation of our algorithm based on the results of sensitivity tests performed using realistic simulated cloud radiation fields.
► We present an algorithm for retrieval of cloud droplet size from polarized rainbow. ► Polarized rainbow shape is determined by single scattering of cloud droplets. ► The algorithm is designed for the airborne Research Scanning Polarimeter (RSP). ► RSP measures polarized and total reflectance at 9 wavelengths between 410 to 2260nm. ► We evaluate the algorithm using tests on simulated cloud radiation fields.
In the fourth part of our “Cellular Statistical Models of Broken Cloud Fields” series we use the binary Markov processes framework for quantitative investigation of the effects of low resolution of ...idealized satellite observations on the statistics of the retrieved cloud masks. We assume that the cloud fields are Markovian and are characterized by the “actual” cloud fraction (CF) and scale length. We use two different models of observations: a simple discrete-point sampling and a more realistic “pixel” protocol. The latter is characterized by a state attribution function (SAF) which has the meaning of the probability that the pixel with a certain CF is declared cloudy in the observed cloud mask. The stochasticity of the SAF means that the cloud/clear attribution is not ideal and can be affected by external or unknown factors. We show that the observed cloud masks can be accurately described as Markov chains of pixels and use the master-matrix formalism (introduced in Part III of the series) for analytical computation of their parameters: the “observed” CF and scale length. This procedure allows us to establish a quantitative relationship (which is pixel-size dependent) between the actual and the observed cloud-field statistics. The feasibility of restoring the former from the latter is considered. The adequacy of our analytical approach to idealized observations is evaluated using numerical simulations. Comparison of the observed parameters of the simulated datasets with their theoretical expectations showed an agreement within 0.005 for the CF, while for the scale length it is within 1% in the sampling case and within 4% in the pixel case.
The cloud drop effective radius (Re) of the drop size distribution derived from passive satellite sensors is a key variable used in climate research. Validation of these satellite products has often ...taken place under stratiform cloud conditions that favor the assumption of cloud horizontal homogeneity used by the retrieval techniques. However, many studies have noted concerns with respect to significant biases in retrieved Re arising from cloud heterogeneity, for example, in cumulus cloud fields. Here, we examine data collected during the 2019 “Cloud, Aerosol and Monsoon Processes Philippines Experiment” (CAMP2Ex), which, in part, targeted the objective of providing the first detailed evaluation of Re retrieved across multiple platforms and techniques in a cumulus and congestus cloud region. Our evaluation consists of cross-comparisons of Re between the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra satellite, the Research Scanning Polarimeter (RSP) onboard the NASA P-3 aircraft, and in situ measurements from both the NASA P-3 and Learjet aircraft that are all taken in close spatiotemporal proximity to the same cloud fields. A particular advantage of our approach lies in the capability of the RSP to retrieve Re using a bi-spectral MODIS approach and a polarimetric approach, which allows for the evaluation of bi-spectral and polarimetric Re retrievals from an airborne perspective using the same samples.
Averaged over all P-3 flight segments examined here for warm clouds, the RSP polarimetric method, the in situ method, and the bias-adjusted MODIS method of Fu et al. (2019) show a comparable median (mean ± standard deviation) for the Re samples of 9.6 (10.2 ± 4.0) µm, 11.0 (13.6 ± 11.3) µm, and 10.4 (10.8 ± 3.8) µm, respectively. These values are far lower than the values of 15.1 (16.2 ± 5.5) µm and 17.2 (17.7 ± 5.7) µm from the bi-spectral retrievals of RSP and MODIS, respectively. Similar results are observed when Re is segregated by cloud-top height and in detailed case studies. The clouds sampled during CAMP2Ex consist of mostly small (mean transect length ∼ 1.4 km) and low clouds (mean cloud-top height ∼ 1 km), which had more numerous small clouds than the trade wind cumuli sampled in past field campaigns such as Rain in Shallow Cumulus over the Ocean (RICO) and the Indian Ocean Experiment (INDOEX). The overestimates of Re from the RSP bi-spectral technique compared with the polarimetric technique increased as cloud size and cloud optical depth decreased. Drizzle, cloud-top bumpiness, and solar zenith angle, however, are not closely correlated with the overestimate of bi-spectral Re. For shallow clouds that dominated the liquid cloud cover for the CAMP2Ex region and period, we show that 3-D radiative transfer and cloud heterogeneity, particularly for the optically thin and small clouds, appear to be the leading cause of the large positive biases in bi-spectral retrievals. Because this bias varies with the underlying structure of the cloud field, caution continues to be warranted in studies that use bi-spectral Re retrievals in cumulus cloud fields.
We present comparisons of cloud droplet size distributions (DSDs) retrieved from the research scanning polarimeter (RSP) data with correlative in situ measurements made during the North Atlantic ...Aerosols and Marine Ecosystems Study (NAAMES). The airborne portion of this field experiment was based out of St. John's airport, Newfoundland, Canada with the focus of this paper being on the deployment in May–June 2016. RSP was onboard the NASA C-130 aircraft together with an array of in situ and other remote sensing instrumentation. The RSP is an along-track scanner measuring the polarized and total reflectance in 9 spectral channels. Its uniquely high angular resolution allows for characterization of liquid water droplet sizes using the rainbow structure observed in the polarized reflectance over the scattering angle range from 135° to 165°. The rainbow is dominated by single scattering of light by cloud droplets, so its structure is characteristic specifically of the droplet sizes at cloud top (within unit optical depth into the cloud, equivalent to approximately 50 m). A parametric fitting algorithm applied to the polarized reflectance provides retrievals of the droplet effective radius and variance assuming a prescribed size distribution shape (gamma distribution). In addition to this, we use a non-parametric method, the Rainbow Fourier Transform (RFT), which allows us to retrieve the droplet size distribution itself. The latter is important in the case of clouds with complex microphysical structure, or multiple layers of cloud, which result in multi-modal DSDs. During NAAMES the aircraft performed a number of flight patterns specifically designed for comparisons between remote sensing retrievals and in situ measurements. These patterns consisted of two flight segments above the same straight ground track. One of these segments was flown above clouds allowing for remote sensing measurements, while the other was near the cloud top where cloud droplets were sampled. We compare the DSDs retrieved from the RSP data with in situ measurements made by the Cloud Droplet Probe (CDP). The comparisons generally show good agreement (better than 1 μm for effective radius and in most cases better than 0.02 for effective variance) with deviations explainable by the position of the aircraft within the cloud, or by the presence of additional cloud layers between the cloud being sampled by the in situ instrumentation and the altitude of the remote sensing segment. In the latter case, the multi-modal DSDs retrieved from the RSP data were consistent with the multi-layer cloud structures observed in the correlative High Spectral Resolution Lidar (HSRL) profiles. The results of these comparisons provide a rare validation of polarimetric droplet size retrieval techniques, demonstrating their accuracy and robustness and the potential of satellite data of this kind on a global scale.
•Analysis of cloud observations by the research scanning polarimeter is presented.•The RSP cloud droplet size retrievals are validated using in situ measurements.•Multimodal cloud droplet size distributions are retrieved and analyzed.•Correlative lidar profiles are used to complement RSP droplet size retrievals.
The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry ...Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.
3D effects cause substantial underestimation of cloud optical thickness (COT) in airborne and satellite retrievals based on 1D radiative transfer computations (such as in the case of widely used ...bispectral technique). For a single-layer isolated cloud we propose a simple linear correction of the retrieved COT with the renormalization factor dependent on the cloud’s aspect ratio (the ratio between vertical and horizontal dimensions of the cloud). This is an empirical assumption which we successfully test using synthetic 3D RT data. We introduce a heuristic “block model” of 3D radiative effects and show that the functional form of the renormalization factor is consistent with the process of radiation escape from cloud sides in an essentially 3D geometry. We also extend the block model to the case of single-layer broken cloud field with radiative interaction between the neighboring clouds. In this case the renormalization factor depends also on the distance between clouds.
The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances with high angular resolution. It allows for accurate characterization of ...liquid water cloud droplet sizes using the rainbow structure in the polarized reflectance. RSP's observations also provide constraints on the cumulus cloud's 2D cross section, yielding estimates of its geometric shape. In this study for the first time we evaluate the possibility to retrieve vertical profiles of microphysical characteristics along the cloud side by combining these micro- and macrophysical retrieval methods. First we constrain cloud's geometric shape, then for each point on the bright side of its surface we collect data from different scans to obtain the multi-angle polarized reflectance at that point. The rainbow structures of the reflectances from multiple points yield the corresponding droplet size distributions (DSDs), which are then combined into vertical profiles. We present the results of testing the proposed profiling algorithm on simulated data obtained using large eddy simulations and 3D radiative transfer computations. The virtual RSP measurements were used for retrieval of DSD profiles, which then were compared to the actual data from the LES-model output. A cumulus congestus cloud was selected for these tests in preparation for analysis of real measurements made during the Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex). We demonstrate that the use of the non-parametric Rainbow Fourier Transform (RFT) allows for adequate retrieval of the complex altitude-dependent bimodal structure of cloud DSDs.
•New cloud profiling algorithm for the Research Scanning Polarimeter is presented.•The proposed algorithm was tested on simulated data (3D RT on LES-generated cloud).•A cumulus congestus cloud was selected in preparation for CAMP2Ex filed campaign.•Bimodal cloud droplet size distributions were successfully retrieved and analyzed.