Next Article in Journal
Predicting Carbohydrate Concentrations in Avocado and Macadamia Leaves Using Hyperspectral Imaging with Partial Least Squares Regressions and Artificial Neural Networks
Previous Article in Journal
Classification of Small Targets on Sea Surface Based on Improved Residual Fusion Network and Complex Time–Frequency Spectra
Previous Article in Special Issue
Long-Term Evaluation of Aerosol Optical Properties in the Levantine Region: A Comparative Analysis of AERONET and Aqua/MODIS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Considering the Effects of Horizontal Heterogeneities in Satellite-Based Large-Scale Statistics of Cloud Optical Properties

by
Tamás Várnai
1,2,* and
Alexander Marshak
2
1
Goddard Earth Sciences Technology and Research (GESTAR) II, University of Maryland Baltimore County, Baltimore, MD 21250, USA
2
NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(18), 3388; https://doi.org/10.3390/rs16183388
Submission received: 11 August 2024 / Revised: 7 September 2024 / Accepted: 10 September 2024 / Published: 12 September 2024

Abstract

:
This paper explores a new approach to improving satellite measurements of cloud optical thickness and droplet size by considering the radiative impacts of horizontal heterogeneity in boundary-layer cumulus clouds. In contrast to the usual bottom-up approach that retrieves cloud properties for individual pixels and subsequently compiles large-scale statistics, the proposed top-down approach first determines the effect of 3D heterogeneity on large-scale cloud statistics and then distributes the overall effects to individual pixels. The potential of this approach is explored by applying a regression-based scheme to a simulated dataset containing over 3000 scenes generated through large eddy simulations. The results show that the new approach can greatly reduce the errors in widely used bispectral retrievals that assume horizontal homogeneity. Errors in large-scale mean values and cloud variability are typically reduced by factors of two to four for 1 km resolution retrievals—and the reductions remain significant even for a 4 km resolution. The calculations also reveal that over vegetation heterogeneity-caused droplet size retrieval biases are often opposite to the biases found over oceans. Ultimately, the proposed approach shows potential for improving the accuracy of both old and new satellite datasets.

1. Introduction

Clouds constitute a critically important element of the climate system. Their complexity, however, poses formidable challenges, and clouds are widely recognized to be a key source of uncertainty in our understanding and predictions of climate [1,2,3]. As the Intergovernmental Panel on Climate Change 6th assessment report released in 2021 states, “…clouds remain the largest contribution to overall uncertainty in climate feedbacks” [4].
Due to their extensive spatial coverage, satellite observations of clouds can help to reduce current uncertainties, as they provide a rich archive for characterizing the Earth’s cloud cover [5], for evaluating the realism of clouds simulated by climate models [6], and for improving our understanding and the simulation capabilities of cloud-related processes such as aerosol–cloud interactions [7].
The critical cloud information provided by satellites includes the widely used measurements of cloud optical depth and cloud droplet effective radius. These two parameters can help, among others, to better understand cloud radiative effects [8] and to estimate both the water content of clouds and the number of cloud droplets [9,10]—which, in turn, enables further analyses, for example, of aerosol–cloud interactions [11].
Most satellite measurements of cloud optical thickness and droplet size use a bispectral technique developed in the 1970s and 1980s [12,13]. This technique combines observations of reflected sunlight at a visible or near-infrared wavelength where clouds do not absorb sunlight, and at a near-infrared wavelength where clouds absorb light at a rate that increases with droplet size. The technique, however, is based on one-dimensional (1D) radiative transfer theory and relies on the assumption of homogeneous, plane-parallel clouds. This means that, in interpreting the radiance data from a particular pixel, the method considers neither horizontal photon transport (that is, radiative interactions between nearby areas that have different cloud properties), nor the effects of unresolved small-scale variability within the pixel. Since the first studies in the 1970s [14,15], as described below, numerous theoretical and observational studies have indicated that the 1D assumption can cause substantial errors (including biases) in satellite retrievals of both optical thickness and effective radius.
Cloud optical thickness retrievals were shown to underestimate the mean optical thickness for high sun (more than ≈30° above the horizon) and overestimate it for low sun (less than ≈30° above the horizon) in theoretical studies. For high sun, thick clouds scatter light into thin areas where they can easily pass through the cloud layer [16], thus reducing the area-averaged reflectance and, hence, the retrieved optical depth. For low sun, interception by cloud sides prevents some incoming sunlight from simply passing through cloud-free gaps, thereby increasing reflectance and the retrieved optical depth [14,17,18]. Statistical analyses of satellite observations have shown that this introduces a strong solar elevation dependence into the mean optical thicknesses retrieved by the Advanced Very-High-Resolution Radiometer (AVHRR) and MODerate-resolution Imaging Spectroradiometer (MODIS) [19,20,21]. Naturally, such biases can also greatly affect the perceived latitude (or even the local time) dependence of cloud statistics (e.g., Level 3 products) obtained from polar-orbiting or geostationary satellites [22]. Finally, we note that three-dimensional (3D) radiative processes not considered in 1D retrievals distort not only the mean, but also the perceived spatial distribution of cloud optical thickness. Photon diffusion can make clouds appear too smooth [23,24,25], whereas 3D shadowing and side illumination effects can make clouds appear too rough [26,27], asymmetric [28,29,30], and anisotropic [27].
Cloud droplet effective radius retrievals were also shown to be affected by horizontal cloud heterogeneity. Theoretical studies explored these effects in detail [31,32], and satellite data analyses found systematic effects, for example, by linking cloud heterogeneity to differences between droplet sizes estimated using different MODIS wavelengths [33]. An analysis of global multiangle satellite data also indicated large biases in particle size retrievals for heterogeneous clouds [34]. Finally, we note that cloud heterogeneity distorts our view of not only mean, but also the spatial distribution of droplet radius values. Satellite-retrieved effective radii were found to be systematically larger when the cloud top was tilted away from the sun and not toward it [29,30].
Recognizing the limitations of current bispectral retrievals, researchers have explored various alternatives [35], most of which can also consider information on spatial and angular variations and even polarization. One approach seeks to iteratively adjust the estimated cloud structure in order to reduce the discrepancy between the observed cloud reflectances and the cloud reflectances simulated by a 3D radiative transfer model [36], often as part of cloud tomography [37,38,39,40,41,42]. Another approach estimates cloud properties using a neural network, machine learning, or artificial intelligence [43,44,45,46,47]. Such new approaches have yielded promising results and, indeed, revealed potential paths for great improvements in our cloud retrieval capabilities. At present, however, none of these methods have reached the mature phase of being used operationally or as a tool in studies that do not focus on retrieval methodologies, but just need accurate cloud properties. Additionally, some of these techniques require multiangle or polarization observations, which may limit their applicability to advanced satellite sensors with such capabilities.
The present study explores a new approach to improving satellite measurements of cloud optical depth and droplet size for heterogeneous clouds. The primary goal is to help improve current estimates of large-scale cloud statistics, such as the mean and standard deviation of the cloud optical thickness and effective radius over areas that are several dozen kilometers along each side. This could help to remove systematic biases from various cloud studies and could improve the accuracy of some Level 3 cloud products such as the MODIS and VIIRS (Visible Infrared Imaging Radiometer Suite) Level 3 Cloud Properties Products provided at a 1° by 1° latitude–longitude resolution [48].
Given this primary goal, the proposed approach reverses the flow of past efforts that have focused on retrievals for individual pixels, with any large-scale statistics being compiled subsequently. The proposed approach first determines the effect of 3D heterogeneity on 1D estimates of large-scale statistics of cloud parameters and then distributes the overall effects to individual pixels. Compared to the traditional “bottom-up” approach, the proposed top-down approach offers advantages such as higher computational efficiency and a better-posed inversion problem. It promises improved cloud statistics without requiring us to untangle the complex web of 3D radiative interactions that, for decades, have stymied the retrievals of heterogeneous cloud properties for individual pixels. (The terms “bottom-up” and “top-down” seek to capture the following difference: in the “bottom-up” approach, using a physics-based or empirical method, we obtain results for each pixel and then build statistics using the individual pixel values as our foundation, whereas in the “top-down” approach, we start by obtaining overall top-level statistics and then estimate the pixel values underlying these estimated statistics.)
We note that the proposed approach is inspired, in part, by past efforts to parameterize the impact of 3D radiative processes (on radiative fluxes or radiances), not for specific locations, but for heterogeneous cloud fields as a whole [49,50,51,52,53]. The closest predecessor approaches, which also considered biases in cloud retrievals, were described in [54,55].
The outline of this paper is as follows. First, Section 2 describes the dataset used in this study. Next, Section 3 describes the proposed retrieval technique, and Section 4 presents the results. Finally, Section 5 provides a brief summary and offers some concluding remarks.

2. Dataset

The proposed technique seeks to adjust (and improve) current 1D retrieval results using statistical relationships between the optimal adjustments and satellite-observable quantities. To establish and test these statistical relationships, we need three types of data, as follows: (i) heterogeneous cloud fields specified in detail, (ii) radiation fields for these clouds, and (iii) look-up tables for the 1D retrievals whose results we seek to adjust. These three types of data are described in the following subsections.

2.1. Cloud Fields

This study used the 3D liquid water content fields of liquid-phase trade cumulus clouds generated using 52 Large Eddy Simulations (LESs). The LES simulations were performed by the LES model developed by Bjorn Stevens in 2005 [56] and upgraded a few years later [57]. One half (26) of the 52 simulations were performed by Frank Evans and were described in detail in [58]. The other half of the simulations were repeat runs performed by Guoyong Wen using the same setup and inputs, but with different random perturbations to the initial atmospheric fields, which yielded completely different cloud arrangements during the LES simulation period.
Each simulation followed cloud evolution in a 20 by 20 km region through an 8 h long period. The horizontal resolution was 62.5 m, and the vertical resolution was 40 m or finer. From each 8 h simulation, we used cloud structures at 17 instances separated by half-hour intervals, during which, the arrangement of cloud elements could change drastically. At the selected instances, the cloud cover ranged from 0 to 60%, the scene-average optical thickness ranged from 0 to 15, and the cloud top altitude ranged up to 4 km. To illustrate the variety of clouds in our database, Figure 1 shows a few sample scenes.
We note that the LES-generated scenes (e.g., Figure 1) fell into the cloud organization pattern termed “sugar” [59], which occurs frequently over both land and ocean [60].
Because LES simulations are computationally expensive, we enhanced the variety of cloud structures in our database and considered four versions of each scene, with each version being rotated by 0°, 90°, 180°, or 270°. This rotation helped because, from the perspective of the incoming sunlight, the cloud arrangements (and radiative heterogeneity effects) became very different after each rotation—for example, a point in shadow may have become sunlit after a rotation. Naturally, this rotation would not make a difference for exact overhead sun, but for the solar zenith angles (≥15°) considered in this study, the rotation made a significant difference: The results shown throughout the paper rarely feature closely grouped data points that likely come from differently rotated versions of the same LES scene.
Ultimately, after all completely cloud-free instances were eliminated, the dataset contained 3176 LES scenes for use in this study.

2.2. Radiation Fields

The next step in creating our dataset was to simulate the radiation fields a satellite would observe over the LES-generated cloud fields. Specifically, the simulations provided 1 km resolution radiation fields at three wavelengths often used in cloud retrievals over land, as follows: 0.65 µm, 2.13 µm, and 11 µm. For example, these wavelengths and this resolution are fairly similar to the parameters of the MODIS and VIIRS observations used in cloud optical retrievals.
Since this study seeks only to explore a possible new approach and gauge its potential usefulness, the simulations did not need the most realistic, detailed, and all-encompassing set of scene and observation parameters. As a result, the simulations provided radiation fields only for a 0° viewing zenith angle (Θ) and three solar zenith angles (Θ0), 15°, 45°, and 75°. Also, the simulations ignored factors such as the curvature of the Earth, gaseous absorption, and atmospheric aerosols, and used simple approximations for other environmental effects. Rayleigh scattering was considered using the formulas in [61], and a Lambertian underlying surface was assumed with albedos of 0.1, 0.3, and 0.0 at 0.65, 2.13, and 11 µm wavelengths, respectively.
As the LES simulations provided only liquid water content values and gave no information on droplet size distribution, we assumed a modified gamma drop size distribution with an effective radius (re) of 8 µm. We note that using a constant drop size distribution made the analysis much simpler, as we did not need to deal with ambiguities in the “correct” result that would occur if re varied within a pixel. (For example, vertical variations in drop size cause even perfectly accurate methods to yield different results if they weigh differently the drops occurring at various altitudes [62,63].)
The wavelength-dependent cloud droplet extinction efficiencies, single-scattering albedos, and scattering phase functions were calculated using Mie theory [64]. The 3D distribution of 0.65 µm volume extinction coefficients (β0.65) was calculated using the formula,
β 0.65 x , y , z , t = 3 2 L W C x , y , z , t ρ r e
where LWC is the liquid water content provided by the LES simulation (as a function of the coordinates x, y, and z and time t) and ρ is the density of water. The extinction coefficients for other wavelengths, β λ , were then calculated using the equation,
β λ = β 0.65 Q λ Q 0.65
where λ identifies the wavelength (2.13 µm or 11 µm) and Q is the extinction efficiency calculated using Mie theory [64]. The cloud droplet single-scattering albedos and scattering phase functions were calculated using Mie theory at all three wavelengths.
Once all the scene parameters were set, 3D radiative transfer calculations provided the radiation fields that a hypothetical satellite would observe at the three wavelengths. The radiative simulations were performed using an updated version of our Monte Carlo model that was validated through the Intercomparison of 3D Radiation Codes (I3RC) project [65,66] and used in several earlier studies [67,68]. The radiative processes were simulated using the 62.5 m resolution LES clouds, but the simulated fields of nadir radiances were degraded to a 1 km horizontal resolution, which is representative of satellite sensors such as MODIS or VIIRS. The 0.65 and 2.13 µm radiances were converted to reflectances, also called bidirectional reflectance factors (BRFs), defined as B R F = π I / μ 0 F 0 , where I is the nadir radiance, µ0 is the cosine of the solar zenith angle, and F0 is the solar irradiance. The 11 µm radiances were converted first to equivalent brightness temperatures using Planck’s law, and then to apparent effective altitudes, assuming a 20 °C surface temperature and a 6.5 °K/km vertical temperature gradient.
The calculations simulated 400,000 photon paths for each 1 km size pixel (using the method of local estimates [69]), which ensured a sufficient computational accuracy for our purposes. The mean of the largest single-pixel uncertainty values encountered within each scene was less than 0.013 and 0.03 for the 0.65 and 2.13 µm reflectances, respectively, and about 15 m for the effective altitudes.
Finally, we note that our computational resources allowed radiation simulations for all 3176 scenes only for the 45° solar zenith angle; for the 15° and 75° solar zenith angles, we performed simulations for a randomly selected set of 1200 scenes (≈37% of the 3176 scenes).

2.3. Retrieval Tables

After simulating what a satellite would see over the LES scenes, our next task was to examine how the widely used bispectral retrieval method (based on 1D radiative transfer) would interpret the simulated observations. In practice, this meant creating look-up tables that provided cloud optical thickness (τ) and droplet effective radius (re) values as a function of 0.65 µm and 2.13 µm reflectances. These tables—created using the Discrete Ordinates Radiative Transfer (DISORT) model [70]—are illustrated in Figure 2.
As discussed in the next two paragraphs, the figure reveals that droplet size retrieval biases caused by subpixel variability were often opposite over vegetated areas than in the oceanic regions that were the focus of earlier studies on cloud heterogeneity effects [31,32,71,72].
Although the figure considers vegetated regions, the points marked A, B, and C illustrate the nonlinearity-caused biases that were found earlier in oceanic regions (often simulated assuming a dark, non-reflecting ocean surface). For this, let us consider a retrieval performed for a pixel that is clear in one half and is covered by a thick cloud in the other half (with corresponding reflectances identified by points A and B, respectively). A retrieval for the entire pixel will then be based on the average 0.65 µm and 2.13 µm reflectances of points A and B, which is identified as point C. The figure shows that, while the dark blue color at point B indicates that the true re of the cloud is 4 µm, the lighter color at point C shows that the re retrieved for the pixel will be closer to 8 µm. This illustrates that, in the presence of subpixel variability, the curved (nonlinear) pattern of colors identifying the retrieved re values will cause overestimations of the effective radius values. While 3D radiative effects such as side illumination can reverse these errors for individual pixels, overall biases always meant re overestimations over oceans [72].
Over vegetation, however, the situation is different. While even large, absorbing cloud droplets increase the 2.13 µm reflectance over dark oceans, they reduce the 2.13 µm reflectance over vegetation, which is highly reflective at this wavelength. Naturally, the reduction is stronger for thicker clouds, but, as Figure 2 shows, it levels off for very thick clouds. This makes the nonlinearity effect shift the other way. If a sample pixel contains an equal mix of a clear region and a cloud that contains large droplets with re = 30 µm (reflectances identified by points A and D, respectively), the pixel’s average 0.65 µm and 2.13 µm reflectances (at point E) will yield an underestimated effective radius value of around 19 µm. This implies that, over vegetated land, cloud heterogeneity can cause systematic re underestimations for clouds containing large droplets. Moreover, Figure 2f shows that, for very oblique sun, underestimation is expected, even for clouds containing small droplets. As cloud reflection is concentrated into forward scattering directions, the nadir 2.13 µm reflectance of the bright vegetation is reduced by the presence of even mildly absorbing small cloud droplets.
While bispectral retrievals based on 1D radiative transfer can provide reasonable (even if biased) results for relatively flat clouds such as stratocumulus, Figure 3 shows that, for cumulus clouds, this is often not the case. Three-dimensional cloud heterogeneity not only causes retrieval errors, but altogether prevents retrievals for a large portion of pixels. This is because retrievals are usually constrained to the 4–30 µm range regarding droplet effective radius (marked by yellow lines in Figure 3), and 3D radiative effects often move these (simulated) observations outside the range of reflectances corresponding to 4 and 30 µm droplet sizes.
It is important to point out that many of the outlying pixels are cloud-free and have too-high or too-low reflectances only because of the shadowing or side-scattering enhancement effect of nearby clouds. Nonetheless, the yellow text in the bottom-right corner of each figure panel reveals that a large fraction of cloudy pixels also do not allow 1D retrievals.
Finally, we note that the relatively mild curvature of the 8 µm lines in Figure 3a,b shows that the retrieval biases due to subpixel variability (discussed for Figure 2) are fairly small for the 8 µm droplet size used in this study, while Figure 3c shows that this is not the case for very oblique sun (Θ0 = 75°). This means that, for not-too-oblique sun (Θ0 = 15°, 45°), most errors in the re values provided by the 1D bispectral retrievals come from 3D radiative effects, and not from subpixel heterogeneity.

3. Proposed Technique

As discussed in the introduction, the primary goal of this study is to explore the potential help a new technique may offer in improving the current estimates of large-scale cloud statistics—namely, the mean and standard deviation of the cloud optical thickness and effective radius over areas dozens of kilometers on each side. This could help to remove systematic biases from cloud studies and improve the accuracy of Level 3 cloud products based on bispectral retrievals (e.g., MODIS or VIIRS products).
Specifically, we explore reversing the usual data processing flow, which estimates the cloud properties for each pixel and then compiles large-area statistics based on the individual pixel results. Instead, the proposed approach first determines the effect of 3D heterogeneity on 1D estimates of large-area statistics and then distributes this effect among individual pixels. The main advantage of this top-down approach is a better-posed inversion problem. It promises improved cloud statistics without requiring us to untangle the complex web of 3D radiative interactions for individual pixels.
For an initial exploration of the top-down approach, this study used a simple regression-based method to estimate the 3D effects based on readily observable scene characteristics. We considered 30 such potentially relevant scene characteristics: (i) cloud fraction; (ii) the mean and standard deviation of the cloud optical thickness and the ratio of these two; (iii) the χ-parameter defined as χ = e l n τ ¯ / τ ¯ , with overbars indicating averaging over all cloudy pixels [73,74]; (iv) the mean differences in along-sun or cross-sun directions over distances of one, two, or three pixels in cloudy pixels with a 0.65 µm reflectance, optical thickness, or effective cloud top height (based on brightness temperature); (v) the mean relative difference in the 0.65 µm reflectance or optical thickness of cloudy pixels neighboring each other in the along-sun direction; (vi) the standard deviation of the difference in the 0.65 µm reflectance or optical thickness of cloudy pixels neighboring each other in the along-sun direction; (vii) the ratio of the mean differences between the 0.65 µm reflectances of pixels separated by one pixel and three pixels in the cross-sun direction; and (viii) the difference between the mean reflectance, optical thickness, or effective radius over slopes facing toward and away from the sun (as identified by brightness temperature gradients) [28]. Admittedly, this was an arbitrary set chosen based on intuition to capture and take advantage of key aspects of cloud variability such as scale invariance [75]; future studies may very well identify more helpful characteristics.
We note that, in calculating these parameters, we considered pixels to be cloudy if either their true or 1D-retrieved τ-value exceeded 0.4, which is near the cloud detection sensitivity in MODIS cloud products [76]. We also note that, in the future processing of actual satellite images, one may skip an area and jump to the next one if the area is cloud-free or otherwise unsuitable for the proposed method, for example, if cloud height retrievals indicate high (ice) clouds.
Next, we tested which subset of characteristics allowed us to best estimate 3D effects using the algorithm described below (and outlined in Figure 4). These tests compared the results from using all possible combinations of up to 4 characteristics (out of the 30), and then checked whether adding any 5th characteristic to the best 4 would help. We found that, in all cases, four characteristics were sufficient, and including additional ones did not significantly improve the results. However, the exact set of the four most helpful characteristics varied with solar elevation and was different for estimations of cloud optical thickness and particle size, and even for estimations of mean or standard deviation values. Oftentimes, though, the accuracy differences between the most helpful set and a few other sets of four characteristics were minimal. As an example, Table 1 lists the characteristics that were most helpful for estimating the influence of 3D heterogeneity on the scene-averaged values of 1D-retrieved cloud optical thickness ( τ ¯ ) values. The table reveals that, as the sun becomes more oblique and photons travel larger horizontal distances, it becomes increasingly helpful to consider pixel-to-pixel differences over larger distances in the along-sun direction. For very oblique sun, even the asymmetry (A) between the mean τ-values retrieved over slopes facing toward and away from the sun becomes quite helpful.
Step 1 of the proposed top-down technique estimates—given a set of Nc number of characteristics out the 30 described above (Nc ≤ 5)—the sought-after large-scale statistics such as the average cloud optical thickness of (20 km)2 scenes. This involves two sub-steps.
First, a preliminary estimate ( τ ¯ p ) is obtained through linear regression using,
τ ¯ p = τ ¯ 1 D + τ ¯ p = τ ¯ 1 D + a + i = 1 N c b i x i
where τ ¯ 1 D is the scene’s average τ provided by 1D retrievals, τ ¯ p is our preliminary adjustment for 3D effects, x i is the ith observed scene characteristic (one of the 30 mentioned above), and a and b i are empirically determined regression coefficients.
Subsequently, a final estimate ( τ ¯ f ) is obtained by another regression that can capture some nonlinearities by using a second-order polynomial in the regression,
τ ¯ f = τ ¯ 1 D + τ ¯ f = τ ¯ 1 D +   c + d τ ¯ p + e τ ¯ p 2
where τ ¯ f is our final adjustment for 3D effects, and c, d , and e are empirically determined regression coefficients. (All regression coefficients are determined by minimizing the root-mean-square (RMS) difference between the true LES τ ¯ values and the estimates from Equation (3) or (4)).
Step 2 of the proposed top-down technique takes the area-averaged 3D effect ( τ ¯ f ) and distributes it among individual pixels. In practice, this is performed in two sub-steps.
The first sub-step obtains pixel-by-pixel estimates of each pixel’s τ (and Δτ), similarly to Equations (3) and (4), but without the area averaging indicated by overbars, and using pixel-specific local (and not scene-averaged) reflectance, τ, or altitude gradients and variabilities as input characteristics.
The second sub-step normalizes these pixel-to-pixel estimates so their average over the (20 km)2 area matches the result we obtained in Step 1 using Equation (4). We found that the optimal method for this normalization (i.e., rescaling or shifting) depends on solar elevation; the best normalized τ-values ( τ f * ) for individual pixels come from,
τ f * = τ 1 D + τ f * = τ 1 D + τ f τ ¯ f τ f ¯    for   Θ 0 45 °
and from,
τ f * = τ 1 D + τ f * = τ 1 D + τ f + τ ¯ f τ f ¯    for   Θ 0 = 75 °
In these equations, τ f is the pre-normalization estimate of the individual pixel 3D effect, τ f ¯ is the area average of these pre-normalization individual pixel estimates (i.e., τ f values), and τ ¯ f is the area-averaged 3D effect provided by Equation (4).
We note that we prefer using τ ¯ 1 D as a starting point in the estimation process (instead of estimating τ ¯ p directly, without involving τ ¯ 1 D ) because the calculation of τ ¯ 1 D using an existing data-processing algorithm can take care of complicating issues such as Earth curvature, gaseous absorption, or instrument spectral response function, and not dealing with these issues can make the 3D adjustment much simpler. However, relying on the result of 1D retrievals raises the question of how to deal with 3D effects moving pixels outside the range of reflectances that allow valid 1D retrievals (i.e., outside the 1D retrieval space bounded by the r e = 4 µm and r e = 30 µm lines in Figure 3). In such cases, we force 1D retrievals to yield r e = 4 µm or r e = 30 µm and the τ-value corresponding to the “observed” 0.65 µm reflectance along the r e = 4 µm or r e = 30 µm lines in Figure 3. If the “observed” 0.65 µm reflectance lies outside the range of 1D reflectances for 0 < τ < 150, τ 1 D is set to 0 (cloud free) or 150.
Finally, we also note that our technique works the same way for estimating re instead of τ, and for estimating standard deviation instead of mean value—but for brevity, we do not include the equations for those estimations here.

4. Results

For an initial exploration of the technique described above, we divided our LES dataset (described in Section 2) into two parts: 70% of the scenes were used to determine the regression coefficients in Equations (3) and (4), and the remaining 30% of scenes were used to test the performance of the method. Accordingly, each dot in Figure 5, Figure 6, Figure 7 and Figure 8 will represent one of the scenes in the 30% used for testing; 396 testing scenes were used for Θ 0 = 15° and 75°, and 1049 testing scenes were used for Θ 0   = 45°.
Figure 5 shows that the proposed method was effective in removing systematic biases from the scene-averaged cloud optical thickness values estimated via traditional 1D retrievals.
Figure 6 shows that the proposed method also had considerable skill in estimating the error of 1D retrievals—that is, the adjustment needed to correct the scene-averaged cloud optical thickness values. Panels (a) and (b) show that, for moderate sun angles, the method captured and correctly accounted for over 80% of the scene-to-scene variability in 1D retrieval errors (R2 > 0.8) and reduced 4-fold the mean average errors (MAEs) to a level below 0.7. Panel c shows that, at very oblique sun, the error reduction was smaller, but the method yielded even more accurate estimates (smaller MAE) than for less oblique sun.
We note that, while the proposed method fully eliminates systematic biases for this dataset, it may not perform quite so well for other datasets that include statistically different cloud structures. Even so, the method can be expected to reduce the overall biases to somewhere between a perfect reduction to zero and a reduction by the same factor that MAE is reduced by in this figure. (This assumes that the distributions of cloud variability in a future dataset and the current dataset will not differ from each other by more than the distributions of cloud variability in individual scenes and in the dataset as a whole tend to differ from each other—or that, if this assumption does not hold, the regression coefficients will be recalculated using an improved dataset.)
To put the results in Figure 6 in perspective, it is helpful to test how the accuracy of the proposed top-down approach compares to that of the traditional bottom-up approach, which first estimates the 1D retrieval errors for each 1 km size pixel and then averages the results over large areas. For these comparisons, we used the same regression-based technique in both approaches: For the top-down approach, we used Equations (3) and (4), and for the bottom-up approach, we the same equations, but without the overbars and using local (as opposed to scene-wide) gradients and variabilities as inputs.
Table 2 shows that, given similar inputs and (regression-based) methodology, the proposed top-down approach outperforms the traditional bottom-up approach. Its higher R2 values show that it can explain a larger fraction of variability in the 1D retrieval errors that we tried to correct, and its lower MAE values show that it yields more accurate results.
As mentioned earlier, most calculations in this paper considered only pixels where cloud detection could be expected if actual (and not simulated) satellite observations were used (i.e., τ 1 D > 0.4). However, Table 3 shows that the proposed method had comparable skill in obtaining the average optical thickness of all clouds, even including the pixels where instruments such as MODIS would not detect a cloud. This suggests that the approach could also help researchers to consider the statistical impacts of thin, small, or shadowed clouds that tend to be missed by cloud detection algorithms.
After considering the scene-averaged τ-values, we can also check whether the proposed approach could also help in improving characterizing the τ variability within each scene. Specifically, Figure 7 demonstrates that the proposed approach can significantly improve our estimates not only for mean cloud properties, but for cloud variability as well. Namely, the proposed approach reduces by a factor from two to four the errors that 3D effects cause in the estimated standard deviation of the cloud optical thicknesses encountered within a scene. We note that this standard deviation is part of some Level 3 satellite data products such as the MODIS or VIIRS cloud products [77,78].
Having explored the potential for improving optical thickness statistics, Figure 8 examines the potential for improving the retrieved mean values of droplet effective radii, r e . (Standard deviation of r e values within each scene is not examined, as our dataset contains constant r e = 8 µm throughout all scenes, implying σ r e = 0 , which would provide a lopsided view revealing only increases but no decreases in σ r e ). The figure shows that the proposed approach can significantly reduce the current retrieval biases not only for optical thickness, but for cloud droplet size as well. It captures 65–90% of the scene-by-scene variability in the 3D effects on the 1D r e retrievals and reduces the mean average error to less than 0.6 µm.
The results, so far, point toward a potentially useful skill in improving the retrievals of cloud properties based on 1 km resolution observations representative of data from polar-orbiting instruments such as MODIS or VIIRS. Thus, it can be of interest to test whether the proposed approach could also help in interpreting coarser-resolution observations, such as those of the Advanced Baseline Imager (ABI) on GOES 16 and 17, the Advanced Himawari Imager (AHI) on Himawari-8, or the Flexible Combined Imager (FCI) on Meteosat-12. To check this, we degraded the 1 km resolution simulated reflectance fields to 2 km and 4 km resolutions and re-tested the performance of the proposed approach for the coarse-resolution data.
Table 4, Table 5 and Table 6 suggest that, while the new method worked best for 1 km resolution satellite data (which can provide detailed statistical information of cloud variability), the proposed method can bring about significant improvements even for satellite data with a coarser (2 km or even 4 km) resolution. All biases were eliminated by design, R2 values dropped only modestly ( r ¯ e at high sun was the one exception), and all MAE values remained significantly lower for the proposed method than for the 1D retrievals.
We note that one factor that helped the MAE values to remain small or even drop for coarser resolutions in Table 4 and Table 5 is that, during spatial averaging, error increases due to losses of information were counteracted by drops in the true τ ¯ and στ values of cloudy pixels. The true mean τ values decrease for coarser resolutions in Table 4 because spatial averaging increases cloud cover. In calculating the average τ of cloudy pixels at a coarser resolution, we divide the unchanged amount of scatterers by a larger cloudy area. (The cloudy area increases because the entire area of a coarse-resolution pixel is considered as cloudy, even if clouds cover only a part of the pixel.) The true στ values decrease for coarser resolutions in Table 5 both because the τ values decrease in general (as mentioned above) and because averaging washes out small-scale variability.
Although the main purpose of the proposed approach is to improve large-scale cloud statistics, the final test of this study explores the second step of the proposed approach: distributing the estimated area-averaged 3D effects among individual pixels. Figure 9 illustrates that, while the new method offers substantial improvements for estimating large-area cloud statistics (as shown in preceding figures and tables), it provides more modest improvements for estimating the cloud properties at individual pixels. This is precisely the reason that motivated the development of the new top-down approach. Correctly estimating radiative heterogeneity effects is a more difficult endeavor for individual pixels than for large areas. We emphasize, however, that further improvements are not at all out of the question—and that even as is, the method removes most systematic biases from the individual pixel estimates and (as indicated by R2 ≈ 0.5) correctly accounts for about half of the variability in the errors of the 1D τ retrievals.

5. Summary and Conclusions

This paper presented a study that was motivated by two key points, as follows: (i) satellite measurements of cloud optical properties are critically important for reducing key uncertainties in climate prediction and (ii) current methods for such satellite measurements use the one-dimensional (1D) approximation of horizontal homogeneity, which causes errors for both individual pixels and large-scale cloud statistics. The specific goal was to explore a new approach to improve the satellite measurements of cloud optical depth and droplet size by considering the radiative impacts of horizontal heterogeneity. The main focus was on improving large area statistics, such as Level 3 satellite products, although the study also briefly explored the issue of improving the accuracy for individual pixels in satellite images.
To improve large-area cloud statistics, this study explored a new approach to the remote sensing of heterogeneous clouds. While past efforts have used a bottom-up approach that starts with retrievals for individual pixels and then subsequently compiles large-scale statistics, this study reversed the order. We sought to first determine the effect of 3D heterogeneity on 1D estimates of large-scale cloud statistics and then distribute the overall effects to individual pixels. The main advantage of this top-down approach is a better-posed inversion problem. We can seek accurate statistics without untangling the complex web of 3D radiative interactions that has, for decades, stymied the retrievals of cloud properties for individual pixels in heterogeneous scenes.
Although this study only considered adjustments to bispectral retrievals based on a single instrument (not unlike MODIS or VIIRS or geostationary satellite imagers), the top-down approach (in which the key inputs are statistical parameters) may also be easier to apply in situations where data are combined from sensors with different view directions, spatial resolutions, or observation times.
To explore this top-down approach, we implemented it in a simple algorithm and tested its performance using a dataset containing over 3000 scenes of boundary-layer cumulus clouds generated through large eddy simulations (LESs). The dataset included both 3D cloud structures and the corresponding radiation fields simulated by a 3D Monte Carlo radiative transfer model. These radiation fields consisted of 0.65 µm, 2.13 µm, and 11 µm radiances that a 1 km resolution satellite looking straight down would observe.
Using this dataset, we tested the accuracy of our simple regression-based algorithm that uses statistical properties such as mean horizontal gradients to adjust the scene-averaged cloud optical thickness and effective radius values provided by 1D retrievals. The main results from these tests were as follows.
  • The results showed that the algorithm greatly reduced the heterogeneity-caused biases in 1 km resolution satellite retrievals of cloud optical thickness and effective radius. For individual (20 km)2 cloud fields, the algorithm reduced the typical mean errors of 1D retrievals to below 0.7 in optical thickness and 0.6 µm in effective radius. In most cases, this meant a reduction by a factor of four or more. Overall biases and errors for larger individual areas are expected to drop even further and be close to zero.
  • The results indicated that, in addition to removing biases from scene-averaged cloud parameters, the new algorithm also improved our estimates of the cloud variability within each cloud field. Specifically, it significantly reduced (by a factor ranging from two to four) the errors that 3D radiative effects cause in current estimates of the standard deviations of cloud optical thicknesses. We note that this standard deviation is also part of some (e.g., MODIS) Level 3 satellite data products.
  • The data revealed that, while the algorithm worked best for 1 km resolution satellite data, it brought about almost as large improvements for 2 km resolution satellite images. Even for 4 km resolution images, the algorithm reduced current errors by factors ranging from about two to six. This implies that the algorithm can help to improve even long-term time series that include data from past or future coarse-resolution (e.g., geostationary) satellites.
  • The results indicated that the proposed top-down approach provided more accurate cloud statistics than comparable bottom-up estimations of large-scale statistics based on individual pixel retrievals.
  • The results also yielded scientific insights into the impact of cloud heterogeneities in satellite remote sensing. Most importantly, they revealed that 1D bispectral satellite retrievals of cloud droplet size often have biases with opposite signs over vegetated surfaces and over oceanic areas. While nonlinearities tend to result in overestimations of cloud droplet sizes over oceans, they cause underestimations over vegetation if the cloud droplet effective radius is large (≥15 µm) or the sun is low above the horizon.
The current study sought only to provide initial insights about the potential of a proposed approach, and it is very important to recognize its limitations. These limitations need to be overcome if this approach is pursued in future studies or perhaps even operational applications. First, future studies should consider a wider variety of scene and observation parameters, for example, by including more view directions and surface albedos, 3D variations in droplet effective radius, ice crystals in addition to liquid water droplets, and scenes from several LES models and atmospheric conditions—including larger scenes that, for simulations over water, should include cloud organization patterns other than the “sugar” pattern used here. Due to the large computational needs of 3D radiative transfer simulations, such expansions in scene and observation parameters will likely need to rely heavily on interpolation (e.g., in view direction or surface albedo), and perhaps take advantage of emerging emulators for quick estimations of the radiation fields in heterogeneous cloud scenes [79]. Real-world applications would also need to develop criteria for identifying the scenes where the proposed approach can be applied. Establishing such criteria may involve a satellite-level verification of any improvement in retrieval capabilities, which is not yet available for the proposed method and some other promising techniques such as neural networks. Finally, future studies may seek improvements to the regression-based estimations of scene statistical parameters, perhaps by using more complex nonlinear regressions based on physical insights or additional input scene characteristics, and perhaps by using machine learning techniques other than the straightforward regressions used here.
Overall, our main conclusion from this study is that the proposed top-down approach does indeed show potential to improve the accuracy of both old and new satellite datasets.

Author Contributions

Conceptualization, T.V.; software, T.V.; investigation, T.V. and A.M.; resources, T.V. and A.M.; writing—original draft preparation, T.V.; writing—review and editing, A.M.; funding acquisition, T.V and A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Aeronautics and Space Administration under grant 80NSSC21K0187.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Acknowledgments

We are grateful to Guoyong Wen and K. Frank Evans for allowing us to use the cloud fields they generated through Large Eddy Simulations.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Bony, S.; Colman, R.; Kattsov, V.M.; Allan, R.P.; Bretherton, C.S.; Dufresne, J.-L.; Hall, A.; Hallegatte, S.; Holland, M.M.; Ingram, W.; et al. How Well Do We Understand and Evaluate Climate Change Feedback Processes? J. Clim. 2006, 19, 3445–3482. [Google Scholar] [CrossRef]
  2. Slingo, J.; Palmer, T. Uncertainty in weather and climate prediction. Phil. Trans. R. Soc. 2011, 369, 4751–4767. [Google Scholar] [CrossRef]
  3. Zelinka, M.D.; Randall, D.A.; Webb, M.J.; Klein, S.A. Clearing clouds of uncertainty. Nature Clim. Chang. 2017, 7, 674–678. [Google Scholar] [CrossRef]
  4. Arias, P.A.; Bellouin, N.; Coppola, E.; Jones, R.G.; Krinner, G.; Marotzke, J.; Naik, V.; Palmer, M.D.; Plattner, G.-K.; Rogelj, J.; et al. 2021: Technical Summary. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S.L., Pe, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M.I., et al., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2021; pp. 33–144. [Google Scholar] [CrossRef]
  5. King, M.D.; Platnick, S.E.; Menzel, W.P.; Ackerman, S.A.; Hubanks, P.A. Spatial and temporal distribution of clouds observed by MODIS onboard the Terra and Aqua satellites. IEEE Trans. Geosci. Remote Sens. 2013, 51, 3826–3852. [Google Scholar] [CrossRef]
  6. Pincus, R.; Hubanks, P.A.; Platnick, S.; Meyer, K.; Holz, R.E.; Botambekov, D.; Wall, C.J. Updated observations of clouds by MODIS for global model assessment. Earth Syst. Sci. Data 2023, 15, 2483–2497. [Google Scholar] [CrossRef]
  7. Quaas, J. Approaches to Observe Anthropogenic Aerosol-Cloud Interactions. Curr. Clim. Chang. Rep. 2015, 1, 297–304. [Google Scholar] [CrossRef]
  8. Stephens, G.L.; O’Brien, D.; Webster, P.J.; Pilewski, P.; Kato, S.; Li, J. The albedo of Earth. Rev. Geophys. 2015, 53, 141–163. [Google Scholar] [CrossRef]
  9. Stephens, G.L. Remote Sensing of the Lower Atmosphere; Oxford University Press: Oxford, UK, 1994; p. 523. [Google Scholar]
  10. Grosvenor, D.P.; Sourdeval, O.; Zuidema, P.; Ackerman, A.; Alexandrov, M.D.; Bennartz, R.; Boers, R.; Cairns, B.; Chiu, J.C.; Christensen, M.; et al. Remote sensing of droplet number concentration in warm clouds: A review of the current state of knowledge and perspectives. Rev. Geophys. 2018, 56, 409–453. [Google Scholar] [CrossRef]
  11. McCoy, D.T.; Bender, F.A.; Mohrmann, J.K.C.; Hartmann, D.L.; Wood, R.; Grosvenor, D.P. The global aerosol-cloud first indirect effect estimated using MODIS, MERRA, and AeroCom. J. Geophys. Res. Atmos. 2017, 122, 1779–1796. [Google Scholar] [CrossRef]
  12. Twomey, S.; Seton, K.J. Inferences of gross microphysical properties of clouds from spectral reflectance measurements. J. Atmos. Sci. 1980, 37, 1065–1069. [Google Scholar] [CrossRef]
  13. Nakajima, T.; King, M.D. Determination of the optical depth and effective particle radius of clouds from reflected solar radiation measurements. Part I: Theory. J. Atmos. Sci. 1990, 47, 1878–1893. [Google Scholar] [CrossRef]
  14. Wendling, P. Albedo and reflected radiance of horizontally inhomogeneous clouds. J. Atmos. Sci. 1977, 34, 642–650. [Google Scholar] [CrossRef]
  15. Reynolds, D.W.; McKee, T.B.; Danielson, K.S. Effects of cloud size and cloud particles on satellite-observed reflected brightness. J. Atmos. Sci. 1978, 35, 160–164. [Google Scholar] [CrossRef]
  16. Cannon, J.C. Line transfer in two dimensions. Astrophys. J. 1970, 161, 255–264. [Google Scholar] [CrossRef]
  17. Várnai, T.; Marshak, A. A method for analyzing how various parts of clouds influence each other’s brightness. J. Geophys. Res. 2003, 108, 4706. [Google Scholar] [CrossRef]
  18. Fauchez, T.; Platnick, S.; Sourdeval, O.; Wang, C.; Meyer, K.; Cornet, C.; Szczap, F. Cirrus horizontal heterogeneity and 3D radiative effects on cloud optical property retrievals from MODIS near to thermal infrared channels as a function of spatial resolution. J. Geophys. Res. 2018, 123, 11141–11153. [Google Scholar] [CrossRef]
  19. Loeb, N.G.; Várnai, T.; Winker, D.M. Influence of subpixel-scale cloud-top structure on reflectances from overcast stratiform cloud layers. J. Atmos. Sci. 1998, 55, 2960–2973. [Google Scholar] [CrossRef]
  20. Loeb, N.G.; Coakley, J.A. Inference of marine stratus cloud optical depths from satellite measurements: Does 1D theory apply? J. Clim. 1998, 11, 215–233. [Google Scholar] [CrossRef]
  21. Seethala, C.; Horváth, A. Global assessment of AMSR-E and MODIS cloud liquid water path retrievals in warm oceanic clouds. J. Geophys. Res. 2010, 115, D13202. [Google Scholar] [CrossRef]
  22. Grosvenor, D.P.; Wood, R. The effect of solar zenith angle on MODIS cloud optical and microphysical retrievals within marine liquid water clouds. Atmos. Chem. Phys. 2014, 14, 7291–7321. [Google Scholar] [CrossRef]
  23. Cahalan, R.F.; Snider, J.B. Marine stratocumulus structure. Remote Sens. Environ. 1989, 28, 95–107. [Google Scholar] [CrossRef]
  24. Marshak, A.; Davis, A.; Wiscombe, W.; Cahalan, R.F. Radiative smoothing in fractal clouds. J. Geophys. Res. 1995, 100, 26247–26261. [Google Scholar] [CrossRef]
  25. Davis, A.; Marshak, A.; Cahalan, R.F.; Wiscombe, W.J. The Landsat scale break in stratocumulus as a three-dimensional radiative transfer effect: Implications for cloud remote sensing. J. Atmos. Sci. 1997, 54, 241–260. [Google Scholar] [CrossRef]
  26. Oreopoulos, L.; Marshak, A.; Cahalan, R.; Wen, G. Cloud three-dimensional effects evidenced in Landsat spatial power spectra and autocorrelation function. J. Geophys. Res. 2000, 105, 14777–14788. [Google Scholar] [CrossRef]
  27. Várnai, T. Influence of three-dimensional radiative effects on the spatial distribution of shortwave cloud reflection. J. Atmos. Sci. 2000, 57, 216–229. [Google Scholar] [CrossRef]
  28. Várnai, T.; Marshak, A. Observations of three-dimensional radiative effects that influence MODIS cloud optical depth retrievals. J. Atmos. Sci. 2002, 59, 1607–1618. [Google Scholar] [CrossRef]
  29. Várnai, T.; Marshak, A. Observations of three-dimensional radiative effects that influence satellite retrievals of cloud properties. Quart. J. Hung. Meteor. Serv. 2002, 106, 265–278. [Google Scholar]
  30. Vant-Hull, B.; Marshak, A.; Remer, L.; Li, Z. The effects of scattering angle and cumulus cloud geometry on satellite retrievals of cloud drop effective radius. IEEE Geosci. Remote Sens. Lett. 2007, 45, 1039–1045. [Google Scholar] [CrossRef]
  31. Marshak, A.; Platnick, S.; Várnai, T.; Wen, G.; Cahalan, R.F. Impact of 3D radiative effects on satellite retrievals of cloud droplet sizes. J. Geophys. Res. 2006, 111, DO9207. [Google Scholar] [CrossRef]
  32. Zhang, Z.; Werner, F.; Cho, H.-M.; Wind, G.; Platnick, S.; Ackerman, A.S.; Di Girolamo, L.; Marshak, A.; Meyer, K. A framework based on 2-D Taylor expansion for quantifying the impacts of subpixel reflectance variance and covariance on cloud optical thickness and effective radius retrievals based on the bispectral method. J. Geophys. Res. 2016, 121, 7007–7025. [Google Scholar] [CrossRef]
  33. Zhang, Z.; Platnick, S. An assessment of differences between cloud effective particle radius retrievals for marine water clouds from three MODIS spectral bands. J. Geophys. Res. 2011, 116, D20215. [Google Scholar] [CrossRef]
  34. Liang, L.; Di Girolamo, L.; Sun, W. Bias in MODIS cloud drop effective radius for oceanic water clouds as deduced from optical thickness variability across scattering angle. J. Geophys. Res. 2015, 120, 7661–7681. [Google Scholar] [CrossRef]
  35. Marshak, A.; Knyazikhin, Y.; Várnai, T. A new spectrally-invariant approach to the remote sensing of inhomogeneous clouds. Front. Remote Sens. 2024, 5, 1392596. [Google Scholar] [CrossRef]
  36. Martin, W.; Cairns, B.; Bal, G. Adjoint methods for adjusting three-dimensional atmosphere and surface properties to fit multi-angle/multi-pixel polarimetric measurements. J. Quant. Spectrosc. Rad. Transf. 2014, 144, 68–85. [Google Scholar] [CrossRef]
  37. Levis, A.; Schechner, Y.Y.; Davis, A.B.; Loveridge, J. Multi-View Polarimetric Scattering Cloud Tomography and Retrieval of Droplet Size. Remote Sens. 2020, 12, 2831. [Google Scholar] [CrossRef]
  38. Forster, L.; Davis, A.B.; Diner, D.J.; Mayer, B. Toward Cloud Tomography from Space Using MISR and MODIS: Locating the “Veiled Core” in Opaque Convective Clouds. J. Atmos. Sci. 2021, 78, 155–166. [Google Scholar] [CrossRef]
  39. Doicu, A.; Efremenko, D.; Trautmann, T. Cloud tomographic retrieval algorithms. I: Surrogate minimization method. J. Quant. Spectrosc. Radiat. Transf. 2022, 277, 107954. [Google Scholar] [CrossRef]
  40. Doicu, A.; Efremenko, D.; Trautmann, T. Cloud tomographic retrieval algorithms. II: Adjoint method. J. Quant. Spectrosc. Rad. Transf. 2022, 285, 108177. [Google Scholar] [CrossRef]
  41. Loveridge, J.; Levis, A.; Di Girolamo, L.; Holodovsky, V.; Forster, L.; Davis, A.B.; Schechner, Y.Y. Retrieving 3D distributions of atmospheric particles using Atmospheric Tomography with 3D Radiative Transfer—Part 1: Model description and Jacobian calculation. Atmos. Meas. Tech. 2023, 16, 1803–1847. [Google Scholar] [CrossRef]
  42. Loveridge, J.; Levis, A.; Di Girolamo, L.; Holodovsky, V.; Forster, L.; Davis, A.B.; Schechner, Y.Y. Retrieving 3D distributions of atmospheric particles using Atmospheric Tomography with 3D Radiative Transfer—Part 2: Local optimization. Atmos. Meas. Tech. 2023, 16, 3931–3957. [Google Scholar] [CrossRef]
  43. Faure, T.; Isaka, H.; Guillemet, B. Neural network retrieval of cloud parameters from high-resolution multispectral radiometric data—A feasibility study. Remote Sens. Environ. 2002, 80, 285–296. [Google Scholar] [CrossRef]
  44. Faure, T.; Isaka, H.; Guillemet, B. Neural network retrieval of cloud parameters of inhomogeneous and fractional clouds: Feasibility study. Remote Sens. Environ. 2001, 77, 123–138. [Google Scholar] [CrossRef]
  45. Cornet, C.; Isaka, H.; Guillemet, B.; Szczap, F. Neural network retrieval of cloud parameters of inhomogeneous clouds from multispectral and multiscale radiance data: Feasibility study. J. Geophys. Res. 2004, 109, D12203. [Google Scholar] [CrossRef]
  46. Evans, K.F.; Marshak, A.; Várnai, T. The potential for improved cloud optical depth retrievals from the multiple directions of MISR. J. Atmos. Sci. 2008, 65, 3179–3196. [Google Scholar] [CrossRef]
  47. Okamura, R.; Iwabuchi, H.; Schmidt, K.S. Feasibility study of multi-pixel retrieval of optical thickness and droplet effective radius of inhomogeneous clouds using deep learning. Atmos. Meas. Tech. 2017, 10, 4747–4759. [Google Scholar] [CrossRef]
  48. Platnick, S.; Platnick, S.; Meyer, K.G.; Hubanks, P.; Holz, R.; Ackerman, S.A.; Heidinger, A.K. VIIRS Atmosphere L3 Cloud Properties Product. Version-1.1; NASA Level-1 and Atmosphere Archive & Distribution System (LAADS) Distributed Active Archive Center (DAAC); Goddard Space Flight Center: Greenbelt, MD, USA, 2019. [Google Scholar] [CrossRef]
  49. Schmetz, J. On the parameterization of the radiative properties of broken clouds. Tellus A Dyn. Meteorol. Oceanogr. 1984, 36A, 417–432. [Google Scholar] [CrossRef]
  50. Kobayashi, T. Parameterization of Reflectivity for Broken Cloud Fields. J. Atmos. Sci. 1988, 45, 3034–3045. [Google Scholar] [CrossRef]
  51. Welch, R.M.; Wielicki, B.A. Reflected Fluxes for Broken Clouds over a Lambertian Surface. J. Atmos. Sci. 1989, 46, 1384–1395. [Google Scholar] [CrossRef]
  52. Bréon, F. Reflectance of Broken Cloud Fields: Simulation and Parameterization. J. Atmos. Sci. 1992, 49, 1221–1232. [Google Scholar] [CrossRef]
  53. Barker, H.W.; Liu, D. Inferring Optical Depth of Broken Clouds from Landsat Data. J. Clim. 1995, 8, 2620–2630. [Google Scholar] [CrossRef]
  54. Iwabuchi, H. Retrieval of cloud optical thickness and effective radius using multispectral remote sensing and accounting for 3D effects. In Light Scattering Reviews 2; Kokhanovsky, A.A., Ed.; Springer Praxis Books; Springer: Berlin/Heidelberg, Germany, 2007; pp. 97–124. [Google Scholar] [CrossRef]
  55. Iwabuchi, H.; Hayasaka, T. Effects of cloud horizontal inhomogeneity on the optical thickness retrieved from moderate-resolution satellite data. J. Atmos. Sci. 2002, 59, 2227–2242. [Google Scholar] [CrossRef]
  56. Stevens, B.; Moeng, C.-H.; Ackerman, A.S.; Bretherton, C.S.; Chlond, A.; de Roode, S.; Edwards, J.; Golaz, J.-C.; Jiang, H.; Khairoutdinov, M.; et al. Evaluation of Large-Eddy Simulations via Observations of Nocturnal Marine Stratocumulus. Mon. Weather Rev. 2005, 133, 1443–1462. [Google Scholar] [CrossRef]
  57. Savic-Jovcic, V.; Stevens, B. The Structure and Mesoscale Organization of Precipitating Stratocumulus. J. Atmos. Sci. 2008, 65, 1587–1605. [Google Scholar] [CrossRef]
  58. Wen, G.; Marshak, A.; Várnai, T.; Levy, R. Testing the two-layer model for correcting near-cloud reflectance enhancement using LES/SHDOM-simulated radiances. J. Geophys. Res. Atmos. 2016, 121, 9661–9674. [Google Scholar] [CrossRef] [PubMed]
  59. Bony, S.; Schulz, H.; Vial, J.; Stevens, B. Sugar, gravel, fish and flowers: Dependence of mesoscale patterns of trade-wind clouds on environmental conditions. Geophys. Res. Lett. 2020, 47, e2019GL085988. [Google Scholar] [CrossRef]
  60. Rasp, S.; Schulz, H.; Bony, S.; Stevens, B. Combining Crowdsourcing and Deep Learning to Explore the Mesoscale Organization of Shallow Convection. Bull. Amer. Meteor. Soc. 2020, 101, E1980–E1995. [Google Scholar] [CrossRef]
  61. Iqbal, M. An Introduction to Solar Radiation; Academic Press: Cambridge, MA, USA, 1983; p. 390. [Google Scholar]
  62. Platnick, S. Vertical photon transport in cloud remote sensing problems. J. Geophys. Res. 2000, 105, 22919–22935. [Google Scholar] [CrossRef]
  63. Bréon, F.M.; Doutriaux-Boucher, M. A comparison of cloud droplet radii measured from space. IEEE Trans. Geosci. Remote Sens. 2005, 43, 1796–1805. [Google Scholar] [CrossRef]
  64. Wiscombe, W.J. Improved Mie scattering algorithms. Appl. Opt. 1980, 19, 1505–1509. [Google Scholar] [CrossRef]
  65. Cahalan, R.F.; Oreopoulos, L.; Marshak, A.; Evans, K.F.; Davis, A.B.; Pincus, R.; Yetzer, K.H.; Mayer, B.; Davies, R.; Ackerman, T.P.; et al. The International Intercomparison of 3D Radiation Codes (I3RC): Bringing together the most advanced radiative transfer tools for cloudy atmospheres. Bull. Amer. Meteor. Soc. 2005, 86, 1275–1293. [Google Scholar] [CrossRef]
  66. Oreopoulos, L.; Marshak, A.; Cahalan, R.F.; Várnai, T.; Davis, A.B.; Macke, A. New Directions in 3D Atmospheric Radiative Transfer. EOS 2006, 87, 52. [Google Scholar] [CrossRef]
  67. Várnai, T. Multiyear statistics of 2D shortwave radiative effects at three ARM sites. J. Atmos. Sci. 2010, 67, 3757–3762. [Google Scholar] [CrossRef]
  68. Várnai, T.; Marshak, A.; Yang, W. Multi-satellite aerosol observations in the vicinity of clouds. Atmos. Chem. Phys. 2013, 13, 3899–3908. [Google Scholar] [CrossRef]
  69. Marchuk, G.I.; Mikhailov, G.A.; Nazaraliev, M.A.; Darbinjan, R.A.; Kargin, B.A.; Elepov, B.S. The Monte Carlo Methods in Atmospheric Optics; Springer: New York, NY, USA, 1980; p. 208. [Google Scholar]
  70. Stamnes, K.; Tsay, S.-C.; Wiscombe, W.; Jayaweera, K. Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media. Appl. Opt. 1988, 27, 2502–2509. [Google Scholar] [CrossRef] [PubMed]
  71. Zhang, Z.; Ackerman, A.S.; Feingold, G.; Platnick, S.; Pincus, R.; Xue, H. Effects of cloud horizontal inhomogeneity and drizzle on remote sensing of cloud droplet effective radius: Case studies based on large-eddy simulations. J. Geophys. Res. 2012, 117, D19208. [Google Scholar] [CrossRef]
  72. Fu, D.; Di Girolamo, L.; Liang, L.; Zhao, G. Regional biases in MODIS marine liquid water cloud drop effective radius deduced through fusion with MISR. J. Geophys. Res. 2019, 124, 13182–13196. [Google Scholar] [CrossRef]
  73. Cahalan, R.F.; Ridgway, W.; Wiscombe, W.J.; Bell, T.L.; Snider, J.B. The albedo of fractal stratocumulus clouds. J. Atmos. Sci. 1994, 51, 2434–2455. [Google Scholar] [CrossRef]
  74. Oreopoulos, L.; Cahalan, R.F. Cloud Inhomogeneity from MODIS. J. Clim. 2005, 18, 5110–5124. [Google Scholar] [CrossRef]
  75. Davis, A.; Marshak, A.; Gerber, H.; Wiscombe, W. Horizontal Structure of Marine Boundary-Layer Clouds from Cm– to Km–Scales. J. Geophys. Res. 1999, 104, 6123–6144. [Google Scholar] [CrossRef]
  76. Ackerman, S.A.; Holz, R.E.; Frey, R.; Eloranta, E.W.; Maddux, B.C.; McGill, M. Cloud Detection with MODIS. Part II: Validation. J. Atmos. Ocean. Technol. 2008, 25, 1073–1086. [Google Scholar] [CrossRef]
  77. Platnick, S.; King, M.D.; Ackerman, S.A.; Menzel, W.P.; Baum, B.A.; Riédi, J.C.; Frey, R.A. The MODIS cloud products: Algorithms and examples from Terra. IEEE Trans. Geosci. Remote Sens. 2003, 41, 459–473. [Google Scholar] [CrossRef]
  78. Platnick, S.; Meyer, K.; Wind, G.; Holz, R.E.; Amarasinghe, N.; Hubanks, P.A.; Marchant, B.; Dutcher, S.; Veglio, P. The NASA MODIS-VIIRS Continuity Cloud Optical Properties Products. Remote Sens. 2021, 13, 2. [Google Scholar] [CrossRef]
  79. Chiu, J.Y.C.; Yang, K. Three-dimensional radiative transfer emulation using machine learning techniques. In Proceedings of the 104th Annual Meeting of the American Meteorological Society, Baltimore, MD, USA, 28 January–1 February 2024. [Google Scholar]
Figure 1. Side and top views of a few scenes that illustrate the variety of clouds in the used LES dataset. The plots show horizontally or vertically integrated liquid water content values. The images were provided by K. Frank Evans.
Figure 1. Side and top views of a few scenes that illustrate the variety of clouds in the used LES dataset. The plots show horizontally or vertically integrated liquid water content values. The images were provided by K. Frank Evans.
Remotesensing 16 03388 g001
Figure 2. A depiction of the look-up tables used for simulated 1D retrievals of (ac) cloud optical thickness τ and (df) effective radius re. The solar zenith angle Θ0 in the three rows is 15° (top), 45° (middle), and 75° (bottom). τ-values can range up to 150; the color scale in panels (ac) saturates at τ = 30 only to highlight variations at low values. Marks A–E are explained in the text.
Figure 2. A depiction of the look-up tables used for simulated 1D retrievals of (ac) cloud optical thickness τ and (df) effective radius re. The solar zenith angle Θ0 in the three rows is 15° (top), 45° (middle), and 75° (bottom). τ-values can range up to 150; the color scale in panels (ac) saturates at τ = 30 only to highlight variations at low values. Marks A–E are explained in the text.
Remotesensing 16 03388 g002
Figure 3. Two-dimensional histograms of 0.65 and 2.13 µm reflectances simulated for all 1 km size pixels. The color scale saturates at 30 only to make the variations at lower pixels numbers easier to see; the actual histogram values reach well above 30. The solar zenith angle Θ0 in the three figure panels is (a) 15°, (b) 45°, and (c) 75°. Yellow lines mark the 1D look-up table reflectances for 4, 8, and 30 µm droplet sizes, with τ ranging from 0 to 150 along each line. The plots also specify what fraction of cloudy pixels (with τ > 0.4) lie outside the 1D retrieval space bounded by the 4 µm and 30 µm lines. The size of the short plum-colored lines near the origins illustrates the average (over all scenes) of the scene-maximum Monte Carlo simulation uncertainty. Panel b is based on 1,270,400 pixels in 3176 scenes; panels a and c are based on 480,000 pixels in 1200 scenes.
Figure 3. Two-dimensional histograms of 0.65 and 2.13 µm reflectances simulated for all 1 km size pixels. The color scale saturates at 30 only to make the variations at lower pixels numbers easier to see; the actual histogram values reach well above 30. The solar zenith angle Θ0 in the three figure panels is (a) 15°, (b) 45°, and (c) 75°. Yellow lines mark the 1D look-up table reflectances for 4, 8, and 30 µm droplet sizes, with τ ranging from 0 to 150 along each line. The plots also specify what fraction of cloudy pixels (with τ > 0.4) lie outside the 1D retrieval space bounded by the 4 µm and 30 µm lines. The size of the short plum-colored lines near the origins illustrates the average (over all scenes) of the scene-maximum Monte Carlo simulation uncertainty. Panel b is based on 1,270,400 pixels in 3176 scenes; panels a and c are based on 480,000 pixels in 1200 scenes.
Remotesensing 16 03388 g003
Figure 4. Flowchart of the proposed algorithm.
Figure 4. Flowchart of the proposed algorithm.
Remotesensing 16 03388 g004
Figure 5. Comparison of actual scene-averaged cloud optical thickness values ( τ ¯ ) with estimates obtained through 1D retrievals (red) or the proposed method for considering 3D radiative effects (blue). The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Figure 5. Comparison of actual scene-averaged cloud optical thickness values ( τ ¯ ) with estimates obtained through 1D retrievals (red) or the proposed method for considering 3D radiative effects (blue). The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Remotesensing 16 03388 g005
Figure 6. Comparisons of the estimated and actual errors of 1D retrievals in obtaining the average cloud optical thickness of each scene in our dataset. τ ¯ t r u e is the true scene-averaged τ value, while τ ¯ 1 D is the estimate based on 1D retrievals. Each dot represents a scene, the position along the vertical axis represents τ ¯ f estimated in Equation (4), and the position along the horizontal axis represents the actual error of 1D retrievals that we seek to determine. If the method worked perfectly, all dots would be along the 1to-1 line. The specified statistics include the correlation coefficient (R), the coefficient of determination (R2), the mean absolute error (MAE), and the bias (mean error, with negative values meaning underestimation). No subscript means the values are for the proposed method, while the subscript “1D” identifies the MAE or bias of assuming that 1D retrievals provide accurate results. The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Figure 6. Comparisons of the estimated and actual errors of 1D retrievals in obtaining the average cloud optical thickness of each scene in our dataset. τ ¯ t r u e is the true scene-averaged τ value, while τ ¯ 1 D is the estimate based on 1D retrievals. Each dot represents a scene, the position along the vertical axis represents τ ¯ f estimated in Equation (4), and the position along the horizontal axis represents the actual error of 1D retrievals that we seek to determine. If the method worked perfectly, all dots would be along the 1to-1 line. The specified statistics include the correlation coefficient (R), the coefficient of determination (R2), the mean absolute error (MAE), and the bias (mean error, with negative values meaning underestimation). No subscript means the values are for the proposed method, while the subscript “1D” identifies the MAE or bias of assuming that 1D retrievals provide accurate results. The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Remotesensing 16 03388 g006
Figure 7. Same as Figure 6, but for the σ τ standard deviation (and not for the mean) of cloud optical thickness values within each scene. The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Figure 7. Same as Figure 6, but for the σ τ standard deviation (and not for the mean) of cloud optical thickness values within each scene. The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Remotesensing 16 03388 g007
Figure 8. Same as Figure 6, but for the scene average droplet effective radius r ¯ e (instead of the scene average cloud optical thickness τ ¯ ). The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Figure 8. Same as Figure 6, but for the scene average droplet effective radius r ¯ e (instead of the scene average cloud optical thickness τ ¯ ). The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Remotesensing 16 03388 g008
Figure 9. Same as Figure 6, but showing results for individual (1 km)2 size pixels instead of average values over (20 km)2 size scenes. Because of the abundance of individual pixel values, not all data points look identical; instead, colors indicate the number of pixels at each coordinate. The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Figure 9. Same as Figure 6, but showing results for individual (1 km)2 size pixels instead of average values over (20 km)2 size scenes. Because of the abundance of individual pixel values, not all data points look identical; instead, colors indicate the number of pixels at each coordinate. The three figure panels are for solar zenith angles (Θ0) of (a) 15°, (b) 45°, and (c) 75°.
Remotesensing 16 03388 g009
Table 1. List of scene characteristics most helpful in estimating 3D effects on τ ¯ . [Notation: Θ0: solar zenith angle, σ: standard deviation, Δ: difference, z: cloud top height, τ: cloud optical thickness, overbar: averaging over cloudy pixels in area, as: along-sun, cs: cross-sun, numeric subscripts: distance in units of pixels, CF: cloud fraction, and A: asymmetry between cloud slopes facing toward and away from the sun].
Table 1. List of scene characteristics most helpful in estimating 3D effects on τ ¯ . [Notation: Θ0: solar zenith angle, σ: standard deviation, Δ: difference, z: cloud top height, τ: cloud optical thickness, overbar: averaging over cloudy pixels in area, as: along-sun, cs: cross-sun, numeric subscripts: distance in units of pixels, CF: cloud fraction, and A: asymmetry between cloud slopes facing toward and away from the sun].
Θ0Most Helpful Parameters
iiiiiiiv
15° σ τ z ¯ a s , 1 τ ¯ c s , 1 CF
45° σ τ z ¯ a s , 1 τ ¯ c s , 1 z ¯ a s , 2
75° σ τ z ¯ a s , 1 τ ¯ a s , 3 Aτ
Table 2. Comparison of two approaches for estimating the error of scene-averaged cloud optical thickness (τ) values provided by 1D retrievals: (1) the traditional bottom-up approach of averaging 1 km resolution pixel-by-pixel error estimations over each scene and (2) the proposed top-down approach of estimating scene-average errors based on scene variability statistics. [Notation: Θ0: solar zenith angle, R2: coefficient of determination (square of the correlation between true and estimated errors of 1D retrievals), and MAE: mean absolute error].
Table 2. Comparison of two approaches for estimating the error of scene-averaged cloud optical thickness (τ) values provided by 1D retrievals: (1) the traditional bottom-up approach of averaging 1 km resolution pixel-by-pixel error estimations over each scene and (2) the proposed top-down approach of estimating scene-average errors based on scene variability statistics. [Notation: Θ0: solar zenith angle, R2: coefficient of determination (square of the correlation between true and estimated errors of 1D retrievals), and MAE: mean absolute error].
Θ0Bottom-Up ApproachTop-Down Approach
R2MAER2MAE
15°0.530.990.810.65
45°0.580.850.850.53
75°0.100.630.430.43
Table 3. Accuracy in estimating the impact of 3D effects on the average optical thickness τ ¯ only for cloudy pixels with τ 1 D > 0.4 (where clouds are likely to be detected by instruments such as MODIS), or for entire scenes (including even the thin-cloud or cloud-free pixels where instruments such as MODIS are unlikely to detect any cloud presence, thus capturing the effect of 3D heterogeneities even over undetected clouds). Higher R2 values and smaller MAE/MAE1D ratios indicate better performance. [Notation: Θ0: solar zenith angle, R2: coefficient of determination (square of the correlation between the true and estimated errors of 1D retrievals), and MAE: mean absolute error. No subscript means the values are for the proposed method, while the subscript “1D” identifies the MAE of assuming that 1D retrievals provide accurate results].
Table 3. Accuracy in estimating the impact of 3D effects on the average optical thickness τ ¯ only for cloudy pixels with τ 1 D > 0.4 (where clouds are likely to be detected by instruments such as MODIS), or for entire scenes (including even the thin-cloud or cloud-free pixels where instruments such as MODIS are unlikely to detect any cloud presence, thus capturing the effect of 3D heterogeneities even over undetected clouds). Higher R2 values and smaller MAE/MAE1D ratios indicate better performance. [Notation: Θ0: solar zenith angle, R2: coefficient of determination (square of the correlation between the true and estimated errors of 1D retrievals), and MAE: mean absolute error. No subscript means the values are for the proposed method, while the subscript “1D” identifies the MAE of assuming that 1D retrievals provide accurate results].
Θ0 τ 1 D > 0.4 OnlyEntire Scene
R2MAE/MAE1DR2MAE/MAE1D
15°0.810.250.880.33
45°0.850.240.830.41
75°0.430.640.690.20
Table 4. Accuracy statistics for estimating the errors of the τ   ¯ scene-averaged cloud optical thickness values provided by 1D retrievals. [Notation: Θ0: solar zenith angle, R2: coefficient of determination (square of the correlation between true and estimated errors), and MAE: mean absolute error. No subscript means the values are for the proposed method, while the subscript “1D” identifies the MAE or bias of assuming that 1D retrievals provide accurate results].
Table 4. Accuracy statistics for estimating the errors of the τ   ¯ scene-averaged cloud optical thickness values provided by 1D retrievals. [Notation: Θ0: solar zenith angle, R2: coefficient of determination (square of the correlation between true and estimated errors), and MAE: mean absolute error. No subscript means the values are for the proposed method, while the subscript “1D” identifies the MAE or bias of assuming that 1D retrievals provide accurate results].
Resolution
[km]
Θ0 τ ¯
R2MAEBiasMAE1DBias1DTrue Mean τ ¯
115°0.810.650.02.61−2.614.65
45°0.850.530.02.17−2.16
75°0.430.430.00.670.38
215°0.740.530.01.70−1.703.23
45°0.810.570.01.76−1.75
75°0.380.430.00.59−0.42
415°0.710.530.01.46−1.462.74
45°0.770.530.01.38−1.36
75°0.630.510.00.92−0.85
Table 5. Same as Table 4, but for estimating the error in the standard deviation (and not the mean) of cloud optical thickness values within each scene (στ).
Table 5. Same as Table 4, but for estimating the error in the standard deviation (and not the mean) of cloud optical thickness values within each scene (στ).
Resolution
[km]
Θ0στ
R2MAEBiasMAE1DBias1DTrue Mean στ
115°0.841.180.04.70−4.706.24
45°0.811.140.04.56−4.56
75°0.660.920.02.061.53
215°0.740.710.02.17−2.173.39
45°0.790.820.02.64−2.64
75°0.740.590.01.11−0.17
415°0.750.510.01.27−1.271.87
45°0.700.530.01.32−1.32
75°0.700.450.00.78−0.61
Table 6. Same as Table 4, but for the scene average droplet effective radius r ¯ e instead of the scene average cloud optical thickness τ ¯ .
Table 6. Same as Table 4, but for the scene average droplet effective radius r ¯ e instead of the scene average cloud optical thickness τ ¯ .
Resolution
[km]
Θ0 r ¯ e
R2MAEBiasMAE1DBias1DTrue Mean r ¯ e
115°0.700.580.03.523.528.00
45°0.650.560.00.91−0.08
75°0.900.290.02.52−2.46
215°0.470.620.02.782.778.00
45°0.531.150.04.724.69
75°0.860.400.02.32−2.14
415°0.260.550.02.342.348.00
45°0.512.100.012.312.3
75°0.830.500.02.56−2.23
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Várnai, T.; Marshak, A. Considering the Effects of Horizontal Heterogeneities in Satellite-Based Large-Scale Statistics of Cloud Optical Properties. Remote Sens. 2024, 16, 3388. https://doi.org/10.3390/rs16183388

AMA Style

Várnai T, Marshak A. Considering the Effects of Horizontal Heterogeneities in Satellite-Based Large-Scale Statistics of Cloud Optical Properties. Remote Sensing. 2024; 16(18):3388. https://doi.org/10.3390/rs16183388

Chicago/Turabian Style

Várnai, Tamás, and Alexander Marshak. 2024. "Considering the Effects of Horizontal Heterogeneities in Satellite-Based Large-Scale Statistics of Cloud Optical Properties" Remote Sensing 16, no. 18: 3388. https://doi.org/10.3390/rs16183388

APA Style

Várnai, T., & Marshak, A. (2024). Considering the Effects of Horizontal Heterogeneities in Satellite-Based Large-Scale Statistics of Cloud Optical Properties. Remote Sensing, 16(18), 3388. https://doi.org/10.3390/rs16183388

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop