Impacts of the Three-dimensional Radiative Effects on Cloud Droplet Number Concentration Retrieval and Aerosol Cloud Interaction Analysis
Abstract. Cloud droplet number concentration (ππ) in warm liquid clouds play a crucial role in understanding cloud microphysical processes and the influence of aerosol–cloud interactions (ACI) on Earth’s climate. ππ from satellite-retrieved cloud properties such as the cloud optical thickness (τ) and cloud droplet effective radius (ππ) can be biased due to the three-dimensional (3D) radiative transfer (RT) effects. Using Large-Eddy Simulation (LES) cloud fields and RT simulations, this study investigates how biases in cloud property retrievals caused by 3D-RT effects impact the derived ππ and subsequent ACI analyses. Our sensitivity studies confirm that the bi-spectral retrievals using the 3.7 ππ channel—whose ππ retrieval is closest to cloud top— shows better agreement with ππ from our LES models, compared to results based on the 1.6 and 2.1 ππ retrievals. At native LES resolution, ππ across all absorbing channels is strongly impacted by the 3D-effects, with the magnitude depending on the solar zenith angles (SZAs); on average, for high/low sun conditions ππ under 3D-RT underestimates/overestimates its 1D-RT counterpart, which indicates dominant darkening/brightening effects. At coarser satellite-like resolutions, average statistics between 1D and 3D retrievals agree better, indicating compensation between 3D and plane-parallel effects. Furthermore, the impact of 3D-effects on ACI analyses produced similar results across all spectral band pairings, with minimal disagreement between 1D and 3D at coarse spatial resolution. Together, these results indicate that 3D retrieval artifacts in bi-spectral ππ retrievals do not seem to drive uncertainties associated with radiative impact applications, resulting in reliable ACI and flux-related analyses.
Competing interests: One of the (co) authors is a member of the editorial board of ACP journal.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Review of egusphere-2025-4169
General Comments
This study examines uncertainties in the most widely used remote sensing retrievals of droplet number concentrations and their implications of albedo susceptibility calculations, which are a core part of the analysis of aerosol cloud interactions and highly relevant to our understanding of climate. The authors make a valuable novel contribution in this study by using Large Eddy Simulations and 3D radiative transfer (3DRT) models and so can provide end-to-end analysis of errors. They also claim to have analysed albedo susceptibility but there are significant flaws in this analysis for which I recommend major revisions. I believe that once these are addressed that this article will make a valuable contribution to ACP.
Specific Comments
βACI analysisβ in the title & abstract etc. is far too broad for what is covered in this paper. This should be narrowed to βalbedo susceptibility estimatesβ.
Here, I will be blunt, if only to impress upon the editor the importance of the issues I raise. This study neglects the anisotropy of radiation while making claims about a radiative flux (and its susceptibility) in Section 2.4. There is also a lack of a reference calculation of albedo susceptibility (discussed further below). These omissions render the results in this section largely meaningless.
A large part of the 3DRT effect on TOA irradiances is changes in the anisotropy of radiation with optical depth. For example, for near overhead sun, nadir reflectance saturates more quickly for finite clouds (e.g., a cube) than the oblique views. In other words, the angular pattern of the susceptibility changes from 1DRT to 3DRT. The effects of 3DRT on albedo susceptibility cannot be estimated while neglecting the anisotropy of radiation. The burden of proof is on the authors to demonstrate that the anisotropy susceptibility is unimportant for their analysis.
This study presents several different estimates of albedo susceptibility and compares them in Section 2.4. For example, a 1DRT+regression estimate applied to measurements generated using 1DRT and measurements generated using 3DRT (at 100 m or 800 m resolution) are compared. None of these four estimates are the truth and so it is not possible to tell whether agreement between estimates indicates good or poor performance. Without a clear reference, this makes it easy for a reader to come to an incorrect conclusion about what these results imply about the accuracy of using 1DRT to estimate albedo susceptibility.
The derivative of TOA albedo (averaged over some reference area) with respect to optical depth in each column has an exact expression in terms of linearized 3DRT (Doicu and Efremenko, 2019) that should be used as the reference for evaluating all approximations of albedo susceptibility.
These elements need to be present or other supporting evidence included to justify the claims presented here. Either the claims about the efficacy of estimating albedo susceptibility using 1DRT need to be removed or appropriate evidence needs to be added.
I note that the calculation of reference-quality albedo susceptibility will require a greater application of technical effort and computational power than the entire contents of the manuscript, which is why I suggest that it be removed.
Other specific flux issues:
Line 211-212: CERES is an estimate of F3D, it is not measuring a fundamentally different quantity. Yes, there will be errors. By construction an empirical ADM method with a one-dimensional scene ID will produce a high-variance estimate of the sceneβs contribution to the flux. Note that the factor mapping from flux to radiance is the anisotropy factor, not the ADM. The ADM is a model for the anisotropy factor (or correction).
Line 213: This does not make sense to me. CERES flux is computed over on the order of 20 km FOVs, CERES does not measure at the scales of pixels in the LES simulations. Moreover, each CERES flux estimates a contribution to the flux at TOA from a localized region (the ~20 km FOV), the radiation from that region reaches TOA spread over a ~hundred km region. See, for example, the difference between the MISR restrictive and expansive albedos. F3D is the flux at TOA, and so of course it is spatially smoothed as the upwelling flux at each horizontal location is sourced from a huge area on the surface, due to integration of upwelling radiance over angle.
Line 219-223: It sounds like the CERES-type flux referred to in this paper is a non-standard quantity that warrants a clear precise definition.
Line 223-224: βThe ADM is the same for all cloudy pixels.β This is incorrect.
The ratio between flux and radiance (the anisotropy factor) varies everywhere in space and angle and is defined everywhere within the 3DRT domain. It is not constant over a scene or small LES domain. CERES approximates the anisotropy factor as constant over a large (order ~20 km) region. The definition for the CERES flux relationship (i.e., between a localized region in a FOV and a flux spread over the TOA) can be applied at higher spatial resolution. This is exemplified in the MISR restrictive albedo and local albedo products, the latter of which is a TOA flux product at 2.2 km resolution (though not broadband). The difficulties here are in the geometric definition of the reference levels for registering TOA flux contributions to, and the ability to produce a model for the anisotropy. This becomes increasingly difficult at high resolution (e.g., ~100 m) because there is no one-size-fits-all definition of a reference level for which we can construct useful low variance models for anisotropy. These issues are irrelevant when the cloud is specified and the albedo susceptibility can be calculated exactly using appropriate numerical modelling.
Line 733-734: Its not flux, itβs a broadband radiance. Given fixed profiles of absorbing gases (or their absence), it would be weird if broadband radiance and bi-spectral radiances werenβt highly correlated.
Line 737-740: These statements are not justified based on the evidence presented either in the paper or with appropriate references.
Line 829-831: This is unjustified, again, as 3D effects on albedo susceptibility have not actually been calculated.
Reference (Zuidema et al., 2008) and discuss.
My challenge to the authors is to better explain their choice of definition for droplet concentration retrieved from bispectral retrievals, as it differs from how such retrievals are used in practice. So, for general readers to get useful information from the study, I think it would be beneficial to directly engage with this issue and address it clearly in the text.
Eq. 1: Why should anyone care about this definition of vertically-averaged Nd? As noted above this Equation in the text, the reference Nd for in situ studies is the average of Nd across the vertical extent of the cloud (or the Nd at cloud base). The vertical averaging of Nd is also the reference used for field validation of remote sensing retrievals of droplet number concentration (Gryspeerdt et al., 2022; Painemal and Zuidema, 2011). Most of the theory used in ACI analysis uses a simplified adiabatic framework in which Nd is constant with height (not accounting for expansion with height i.e., geometrically thin layers). The reasons for these definitions are simple; to isolate the effects of droplet concentrations from other effects (e.g., dynamical and thermodynamic), even if through an incomplete characterization of the 3D Nd field using assumptions that introduce errors. The isolation of droplet concentration is critical to ACI analysis. By using an Nd weighted by extinction coefficient as a reference, this utility is lost, so by choosing this reference there is no traceability to existing ACI analysis.
So, why should I worry whether retrievals agree with the extinction weighted Nd_LES, or not?
If there is a need to change the definition/interpretation of remotely-sensed droplet concentration from quasi-adiabatic to retrieving Eq. 1, then this is an important point that needs to be clearly communicated to the communities that use these retrievals. If this is the case, then I suggest making that argument be the focus of the paper.
For example, the argument might go: βDue to our inability to specify appropriate microphysical profile information, we are unable to usefully bound uncertainty in droplet concentration and therefore must restrict the definition to an extinction-weighted quantity and limit the theoretical insight that is available.β
If the authors are not taking a strong position on how droplet concentration βretrievalsβ are interpreted, it would be useful for them to cleave to the existing definitions and thereby provide useful error characteristics for people using bispectral retrievals to infer droplet concentration. As is, the presented results have limited utility for refining the interpretation of existing studies. The same limitation holds for Kokkola et al. (2025). Gonzalez et al. (2025) take a different approach and should be discussed.
References: (Gonzalez et al., 2025; Kokkola et al., 2025) (see bottom for full citation).
Line 159: Cumulus rising into stratocumulus is rather complex microphysically. Should this cloud field be considered representative? How can the veracity of the cloud microphysical properties in the simulations be evaluated? Some further supporting evidence or commentary on this aspect should be included.
Line 288: How is the CBH and CTH defined? i.e., how do you separate between and cloudy and clear volumes in the model? I ask because the order of magnitude difference between the extinction-weighted and pure droplet concentration averages in the appendix suggest a very strict criterion for clear sky. Commonly, these definitions are designed to pursue coherence between models and in-situ measurements and are based on the sensitivity of the in-situ measurements; e.g., 0.01 g/kg or 10#/cm^3 etc. as in other studies that examine microphysical effects on susceptibility equations (Feingold et al., 2022) (This reference should be discussed in the text). I suggest that this point be examined in more detail as the order-of-magnitude differences seem larger than the errors in Eq. 3 estimated from other LES or from in-situ measurements (Feingold et al., 2022; Grosvenor et al., 2018).
Line 350-352: Nd_LES is not the reference point for adiabatic behaviour so disagreement between Nd_1D and Nd_LES are not necessarily evidence of deviations from adiabatic conditions. If the cloud is adiabatic then Nd_1D will agree with both the vertical average of the droplet concentration (no extinction weighting) and Nd_LES. However, agreement between Nd_1D and Nd_LES does not necessarily mean the cloud is adiabatic. This is just one of the complications of the chosen definition.
Several figures (e.g., Fig. 4&5) show comparisons of retrievals but not the reference (regardless of its definition). Please include the reference/truth in the figures so we can understand which errors compensate.
There are exact calculations for the partial derivatives Β (whatever the definition of Nd). Regression across columns and a cloud field does not provide an exact reference because it can include unwanted covariates, e.g., a correlation between mean droplet concentration and droplet profile shape, which will bias the power-law relationship. Please compare the estimates to the exact calculation for each column.
Β
Technical Comments.
Line 53: Perhaps remote sensing retrievals rather than instruments.
Line 59: Twomey & Seton 1980 is the earlier paper where bispectral sensitivity is introduced.
Line 70: Should βMostβ operational retrievals be βAllβ?
Line 81: I suggest adding Zinner et al. 2010; Cornet & Davies
Line 90: VΓ‘rnai and Davies (1999) did not analyse retrieval accuracy.
Line 86-93: At which wavelength are changes in cloud reflectance being used to define βbrighteningβ and βdarkeningβ. Changes in the relative amount of brightening and darkening across visible and absorbing wavelengths affects retrieval errors.
Line 93-99: Iβm not sure what the message is. With the βAlthough XXX, previous studies have shown YYYβ construction, is it saying that out of the several factors (variability in cloud top height, solar and view geometry, and photons leaking from optically thick to optically thin cloudy regions) that only solar-viewing geometry is significant? I suggest changing βphotons leaking β¦β to βlight scattering from optically thick to optically thin regionsβ.
Line 111: Gryspeerdt et al. 2019 and associated sentence does not seem to mix with the rest of the paragraph, which are about retrieval errors, not ACI in general.
Line 167: What is the range of the droplet size bins?
Line 196: Does this mean there is integration of multiple monochromatic reflectances over the band? Averaging of optical properties over the band? Or just a single monochromatic calculation at each nominal wavelength.
Line 191-205: What is the domain top of the RT calculations? Is it the LES top of 3 km? or based on a particular atmospheric profile?
Line 315-316: No need to cite Loveridge & Di Girolamo (2024) here as that study just employs simplified versions of microphysics based on other studies.
Line 531-532: Arenβt PPHA biases caused by sub-pixel heterogeneity?
Line 533-534: It shouldnβt be regardless. There should be logic motivating the analysis that shouldnβt be waved aside.
Line 547: And plane-parallel optical depth bias contribution?
Line 839-840: Sinclair et al. 2019 does not use active instruments.
Fig. 6&7: The dashed and solid lines for the PDF and the mean do not correspond to each other in the legend, which is confusing.
References
Doicu, A. and Efremenko, D. S.: Linearizations of the Spherical Harmonic Discrete Ordinate Method (SHDOM), Atmosphere, 10, 292, https://doi.org/10.3390/atmos10060292, 2019.
Feingold, G., Goren, T., and Yamaguchi, T.: Quantifying albedo susceptibility biases in shallow clouds, Atmospheric Chemistry and Physics, 22, 3303β3319, https://doi.org/10.5194/acp-22-3303-2022, 2022.
Gonzalez, J., Dipu, S., Jimenez, G., Camps-Valls, G., and Quaas, J.: Machine Learning-Based Retrieval of Cloud Droplet Number Concentration and Liquid Water Path From Satellite Spectral Data, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 18, 21910β21922, https://doi.org/10.1109/JSTARS.2025.3601981, 2025.
Gryspeerdt, E., McCoy, D. T., Crosbie, E., Moore, R. H., Nott, G. J., Painemal, D., Small-Griswold, J., Sorooshian, A., and Ziemba, L.: The impact of sampling strategy on the cloud droplet number concentration estimated from satellite data, Atmospheric Measurement Techniques, 15, 3875β3892, https://doi.org/10.5194/amt-15-3875-2022, 2022.
Painemal, D. and Zuidema, P.: Assessment of MODIS cloud effective radius and optical thickness retrievals over the Southeast Pacific with VOCALS-REx in situ measurements, Journal of Geophysical Research: Atmospheres, 116, https://doi.org/10.1029/2011JD016155, 2011.
Grosvenor, D. P., Sourdeval, O., Zuidema, P., Ackerman, A., Alexandrov, M. D., Bennartz, R., Boers, R., Cairns, B., Chiu, J. C., Christensen, M., Deneke, H., Diamond, M., Feingold, G., Fridlind, A., HΓΌnerbein, A., Knist, C., Kollias, P., Marshak, A., McCoy, D. T., Merk, D., Painemal, D., Rausch, J., Rosenfeld, D., Russchenberg, H., Seifert, P., Sinclair, K., Stier, P., van Diedenhoven, B., Wendisch, M., Werner, F., Wood, R., Zhang, Z., and Quaas, J.: Remote Sensing of Droplet Number Concentration in Warm Clouds: A Review of the Current State of Knowledge and Perspectives, Reviews of Geophysics, 56, 409β453, https://doi.org/10.1029/2017RG000593, 2018.
Kokkola, H., Tonttila, J., CalderΓ³n, S. M., Romakkaniemi, S., Lipponen, A., PerΓ€korpi, A., Mielonen, T., Gryspeerdt, E., Virtanen, T. H., Kolmonen, P., and Arola, A.: Model analysis of biases in the satellite-diagnosed aerosol effect on the cloud liquid water path, Atmospheric Chemistry and Physics, 25, 1533β1543, https://doi.org/10.5194/acp-25-1533-2025, 2025.
Zuidema, P., Xue, H., and Feingold, G.: Shortwave Radiative Impacts from Aerosol Effects on Marine Shallow Cumuli, Journal of the Atmospheric Sciences, 65, 1979β1990, https://doi.org/10.1175/2007JAS2447.1, 2008.