Cloud Radiative Transfer Research Articles

YunMa: Enabling Spectral Retrievals of Exoplanetary Clouds

In this paper, we present YunMa, an exoplanet cloud simulation and retrieval package, which enables the study of cloud microphysics and radiative properties in exoplanetary atmospheres. YunMa simulates the vertical distribution and sizes of cloud particles and their corresponding scattering signature in transit spectra. We validated YunMa against results from the literature. When coupled to the TauREx 3 platform, an open Bayesian framework for spectral retrievals, YunMa enables the retrieval of the cloud properties and parameters from transit spectra of exoplanets. The sedimentation efficiency (f sed), which controls the cloud microphysics, is set as a free parameter in retrievals. We assess the retrieval performances of YunMa through 28 instances of a K2-18 b-like atmosphere with different fractions of H2/He and N2, and assuming water clouds. Our results show a substantial improvement in retrieval performances when using YunMa instead of a simple opaque cloud model and highlight the need to include cloud radiative transfer and microphysics to interpret the next-generation data for exoplanet atmospheres. This work also inspires instrumental development for future flagships by demonstrating retrieval performances with different data quality.

Insights into 3D cloud radiative transfer effects for the Orbiting Carbon Observatory

Abstract. Clouds impose radiance perturbations upon Orbiting Carbon Observatory (OCO-2)-measured spectra. The Spherical Harmonic Discrete Ordinate radiative transfer Method (SHDOM) code is applied in both idealized bar cloud and scene-specific calculations of 1D and 3D radiances in order to understand 3D cloud effects for a wide range of gas vertical optical depths and solar- and sensor-viewing geometries for ocean and land scenes. SHDOM calculations for 36 scenes over the Amazon and the Pacific are co-analyzed with Moderate Resolution Imaging Spectroradiometer (MODIS) radiance-based cloud distance data and the OCO-2 Lite file rawXCO2 for both quality flag =0 (QF0; best quality) and quality flag =1 (QF1; poor quality) data. SHDOM calculations of the ocean and land scenes indicate that the 1D / 3D radiance intensity ratios and rawXCO2 decrease concurrently as the nearest-cloud distance decreases towards zero, especially for the ocean glint QF1 data, which provide the clearest evidence of 3D cloud effects in OCO-2 retrievals. Yearly analysis of OCO-2 O2 A-band continuum radiances indicate that 3D cloud-brightening events are predominant over cloud-shadowing events; therefore, 1D / 3D intensity ratios are predominantly less than unity. Bias corrected (bcXCO2) at cloud distances between 0 and 20 km are calculated for 20∘ latitude bands for 2015–2018. These zonal averages are used to calculate 3D-cloud-effect biases for bcXCO2 data (with a positive bias indicating that OCO-2 underestimates bcXCO2). Averages of 3D-cloud-effect biases, weighted by the number of Lite file data points in each of the nearest-cloud distance bins, in the Northern and Southern hemispheres, are 0.16 (1.31) and 0.26 (1.41) ppm (parts per million), respectively, over the ocean, and −0.13 (0.51) and −0.08 (0.47) ppm over land for QF0 (QF1) data.

Solar radiation budget in the atmosphere under broken cloudy sky. An analytical model

This study theoretically investigates the transport of broadband solar radiation in the Earth's cloudy atmosphere (extended and broken) by considering three plane parallel layers: i) thin bottom layer adjacent to the Earth's surface; ii) intermediate layer where clouds are most embedded; and iii) upper layer that is essentially cloud-free with possible occurrence of only thin clouds. These three layers are coupled by upwelling and downwelling radiation exchanges, including multiple retro-diffusion, and both direct beam and diffuse components are considered. The interrelation and relative magnitude of the solar radiation flux cascade across the atmosphere - αP, Ra, Rc, Rg and G/H – was established and explored as a tool for the present simulations. Two distinct radiation transport regimes through the inter-cloud gaps are identified, with the regime transition occurring when (AR*TAN(θ)+1)*CF = 1 where AR, θ and CF are the cloud aspect ratio, solar zenith angle, and cloud fraction, respectively. The variation of the solar radiation fluxes with CF and θ exhibit the transition between the two radiation transport regimes, as evidenced in the flux absorbed in the cloudy layer. The effect of AR upon the solar fluxes also exhibits the two radiation transport regimes, with one regime having linear dependence of the various fluxes on AR and the other one having fluxes stay at constant levels independent of AR. The shape factor SF of cloud top and bottom boundary “surfaces” is introduced to account for the self-irradiation effect for different cloud types. It is found that increased roughness of one or both boundaries increases the cloud absorptance in all cases whereas transmittance and reflectance for upwelling and downwelling can either be diminished or enhanced.With the proposed model the estimated impact of clouds on the enhancement of the shortwave radiation absorption in the atmospheric column was improved leading to results closer to the cloudy-sky absorptance observations reported in the literature. The ratio of the cloud radiative forcing at the surface to the cloud radiative forcing at the top of the atmosphere (Rnet), obtained by means of the present model, lead to Rnet = 1.3, closer to the estimate of 1.5 by Cess et al. than the values reported by other cloud radiative transfer models.

Parametrizing cloud geometry and its application in a subgrid cloud‐edge erosion scheme

AbstractTo represent the effects of unresolved cloud processes in numerical weather prediction and climate models, parametrizations of the subgrid properties of clouds are required. In this paper, we describe a method for specifying the “cloud‐edge length” within a model grid‐box, which is an important parameter for approximating the subgrid mixing of air at cloud boundaries. We begin by proposing three conceptual models that predict the cloud‐edge length using the grid‐box cloud fraction and a length‐scale to be derived empirically. The conceptual models are then evaluated using a wide range of observations and cloud‐resolving models. Based on the finding that the “effective cloud spacing” approach fits both these data best, we parametrize the effective cloud spacing as a function of pressure and model resolution. An application of this parametrization to the cloud erosion scheme in the ECMWF forecast model is then demonstrated. The effective cloud spacing approach is compared to the “effective cloud scale” approach and is shown to increase cloud fraction in stratocumulus regions, while decreasing cloud fraction in cumulus regions. These cloud changes have the overall effect of decreasing the error of the modelled top‐of‐atmosphere net short‐wave irradiance when compared to CERES observations by around 3%. Additionally, the cloud‐edge length is an important parameter for approximating subgrid radiative transfer and it is hoped that this parametrization will be useful to quantify the effect of representing 3D cloud radiative transfer in global models.

The Complex Refractive Index of Volcanic Ash Aerosol Retrieved From Spectral Mass Extinction

AbstractThe complex refractive indices of eight volcanic ash samples, chosen to have a representative range of SiO2 contents, were retrieved from simultaneous measurements of their spectral mass extinction coefficient and size distribution. The mass extinction coefficients, at 0.33–19 μm, were measured using two optical systems: a Fourier transform spectrometer in the infrared and two diffraction grating spectrometers covering visible and ultraviolet wavelengths. The particle size distribution was measured using a scanning mobility particle sizer and an optical particle counter; values for the effective radius of ash particles measured in this study varied from 0.574 to 1.16 μm. Verification retrievals on high‐purity silica aerosol demonstrated that the Rayleigh continuous distribution of ellipsoids (CDEs) scattering model significantly outperformed Mie theory in retrieving the complex refractive index, when compared to literature values. Assuming the silica particles provided a good analogue of volcanic ash, the CDE scattering model was applied to retrieve the complex refractive index of the eight ash samples. The Lorentz formulation of the complex refractive index was used within the retrievals as a convenient way to ensure consistency with the Kramers‐Kronig relation. The short‐wavelength limit of the electric susceptibility was constrained by using independently measured reference values of the complex refractive index of the ash samples at a visible wavelength. The retrieved values of the complex refractive indices of the ash samples showed considerable variation, highlighting the importance of using accurate refractive index data in ash cloud radiative transfer models.

Effects of cloud condensate vertical alignment on radiative transfer calculations in deep convective regions

Effects of cloud condensate vertical alignment on radiative transfer process were investigated using cloud resolving model explicit simulations, which provide a surrogate for subgrid cloud geometry. Diagnostic results showed that the decorrelation length Lcw varies in the vertical dimension, with larger Lcw occurring in convective clouds and smaller Lcw in cirrus clouds. A new parameterization of Lcw is proposed that takes into account such varying features and gives rise to improvements in simulations of cloud radiative forcing (CRF) and radiative heating, i.e., the peak of bias is respectively reduced by 8Wm−2 for SWCF and 2Wm−2 for LWCF in comparison with Lcw=1km.The role of Lcw in modulating CRFs is twofold. On the one hand, larger Lcw tends to increase the standard deviation of optical depth στ, as dense and tenuous parts of the clouds would be increasingly aligned in the vertical dimension, thereby broadening the probability distribution. On the other hand, larger στ causes a decrease in the solar albedo and thermal emissivity, as implied in their convex functions on τ. As a result, increasing (decreasing) Lcwleads to decreased (increased) CRFs, as revealed by comparisons among Lcw=0, Lcw=1km andLcw=∞. It also affects the vertical structure of radiative flux and thus influences the radiative heating. A better representation of στ in the vertical dimension yields an improved simulation of radiative heating. Although the importance of vertical alignment of cloud condensate is found to be less than that of cloud cover in regards to their impacts on CRFs, it still has enough of an effect on modulating the cloud radiative transfer process.

The MODIS cloud optical and microphysical products: Collection 6 updates and examples from Terra and Aqua.

The MODIS Level-2 cloud product (Earth Science Data Set names MOD06 and MYD06 for Terra and Aqua MODIS, respectively) provides pixel-level retrievals of cloud-top properties (day and night pressure, temperature, and height) and cloud optical properties (optical thickness, effective particle radius, and water path for both liquid water and ice cloud thermodynamic phases-daytime only). Collection 6 (C6) reprocessing of the product was completed in May 2014 and March 2015 for MODIS Aqua and Terra, respectively. Here we provide an overview of major C6 optical property algorithm changes relative to the previous Collection 5 (C5) product. Notable C6 optical and microphysical algorithm changes include: (i) new ice cloud optical property models and a more extensive cloud radiative transfer code lookup table (LUT) approach, (ii) improvement in the skill of the shortwave-derived cloud thermodynamic phase, (iii) separate cloud effective radius retrieval datasets for each spectral combination used in previous collections, (iv) separate retrievals for partly cloudy pixels and those associated with cloud edges, (v) failure metrics that provide diagnostic information for pixels having observations that fall outside the LUT solution space, and (vi) enhanced pixel-level retrieval uncertainty calculations. The C6 algorithm changes collectively can result in significant changes relative to C5, though the magnitude depends on the dataset and the pixel's retrieval location in the cloud parameter space. Example Level-2 granule and Level-3 gridded dataset differences between the two collections are shown. While the emphasis is on the suite of cloud optical property datasets, other MODIS cloud datasets are discussed when relevant.

A novel approach for introducing cloud spatial structure into cloud radiative transfer parameterizations

Subgrid-scale variability is one of the main reasons why parameterizations are needed in large-scale models. Although some parameterizations started to address the issue of subgrid variability by introducing a subgrid probability distribution function for relevant quantities, the spatial structure has been typically ignored and thus the subgrid-scale interactions cannot be accounted for physically. Here we present a new statistical-physics-like approach whereby the spatial autocorrelation function can be used to physically capture the net effects of subgrid cloud interaction with radiation. The new approach is able to faithfully reproduce the Monte Carlo 3D simulation results with several orders less computational cost, allowing for more realistic representation of cloud radiation interactions in large-scale models.

The Role of Clouds: An Introduction and Rapporteur Report

This paper presents an overview of discussions during the Cloud’s Role session at the Observing and Modelling Earth’s Energy FlowsWorkshop. N. Loeb and B. Soden convened this session including 10 presentations by B. Stevens, B. Wielicki, G. Stephens, A. Clement, K. Sassen, D. Hartmann, T. Andrews, A. Del Genio, H. Barker, and M. Sugi addressing critical aspects of the role of clouds in modulating Earth energy flows. Presentation topics covered a diverse range of areas from cloud microphysics and dynamics, cloud radiative transfer, and the role of clouds in large-scale atmospheric circulations patterns in both observations and atmospheric models. The presentations and discussions, summarized below, are organized around several key questions raised during the session. (1) What is the best way to evaluate clouds in climate models? (2) How well do models need to represent clouds to be acceptable for making climate predictions? (3) What are the largest uncertainties in clouds? (4) How can these uncertainties be reduced? (5) What new observations are needed to address these problems? Answers to these critical questions are the topics of ongoing research and will guide the future direction of this area of research.

Cloud optical thickness and liquid water path – does the <i>k</i> coefficient vary with droplet concentration?

Abstract. Cloud radiative transfer calculations in general circulation models involve a link between cloud microphysical and optical properties. Indeed, the liquid water content expresses as a function of the mean volume droplet radius, while the light extinction is a function of their mean surface radius. There is a small difference between these two parameters because of the droplet spectrum width. This issue has been addressed by introducing an empirical multiplying correction factor to the droplet concentration. Analysis of in situ sampled data, however, revealed that the correction factor decreases when the concentration increases, hence partially mitigating the aerosol indirect effect. Five field experiments are reanalyzed here, in which standard and upgraded versions of the droplet spectrometer were used to document shallow cumulus and stratocumulus topped boundary layers. They suggest that the standard probe noticeably underestimates the correction factor compared to the upgraded versions. The analysis is further refined to demonstrate that the value of the correction factor derived by averaging values calculated locally along the flight path overestimates the value derived from liquid water path and optical thickness of a cloudy column, and that there is no detectable relationship between the correction factor and the droplet concentration. It is also shown that the droplet concentration dilution by entrainment-mixing after CCN activation is significantly stronger in shallow cumuli than in stratocumulus layers. These various effects are finally combined to produce the today best estimate of the correction factor to use in general circulation models.

Cloud Remote Sensing Using Midwave IR CO2 and N2O Slicing Channels near 4.5 μm

Narrow channels located in the longwave IR CO2 absorption region between approximately 13.2 and 14.5 μm, the well known CO2 slicing channels, have been proven to be quite effective for the estimates of cloud heights and effective cloud amounts as well as atmospheric temperature profiles. The designs of some of the near-future multi-channel earth observing satellite sensors cannot accommodate these longwave IR channels. Based on the analysis of the multi-channel imaging data collected with the NASA Moderate Resolution Imaging SpectroRadiometer (MODIS) instrument and on theoretical cloud radiative transfer modeling, we have found that narrow channels located at the midwave IR region between approximately 4.2 and 4.55 μm, where the combined CO2 and N2O absorption effects decrease rapidly with increasing wavelength, have similar properties as the longwave IR CO2 slicing channels. The scattering of solar radiation by clouds on the long wavelength side of the 4.3 μm CO2 absorption makes only a small contribution to the upwelling radiances. In order to retain the crucial cloud and temperature sensing capabilities, future satellite sensors should consider including midwave IR CO2 and N2O slicing channels if the longwave IR channels cannot be implemented on the sensors. The hyperspectral data covering the 3.7-15.5 mm wavelength range and measured with the Infrared Atmospheric Sounding Interferometer (IASI) can be used to further assess the utility of midwave IR channels for satellite remote sensing.

The molecular interstellar medium of the Local Group dwarf NGC 6822

Do molecular clouds collapse to form stars at the same rate in all environments? In large spiral galaxies, the rate of transformation of H2 into stars (hereafter SFE) varies little. However, the SFE in distant objects (z~1) is much higher than in the large spiral disks that dominate the local universe. Some small local group galaxies share at least some of the characteristics of intermediate-redshift objects, such as size or color. Recent work has suggested that the Star Formation Efficiency (SFE, defined as the SFRate per unit H2) in local Dwarf galaxies may be as high as in the distant objects. A fundamental difficulty in these studies is the independent measure of the H2 mass in metal-deficient environments. At 490 kpc, NGC6822 is an excellent choice for this study; it has been mapped in the CO(2-1) line using the multibeam receiver HERA on the 30 meter IRAM telescope, yielding the largest sample of giant molecular clouds (GMCs) in this galaxy. Despite the much lower metallicity, we find no clear difference in the properties of the GMCs in NGC 6822 and those in the Milky Way except lower CO luminosities for a given mass. Several independent methods indicate that the total H2 mass in NGC 6822 is about 5 x 10^6 Msun in the area we mapped and less than 10^7 Msun in the whole galaxy. This corresponds to a NH2/ICO ~ 4 x 10^{21} cm^-2 /(Kkm/s) over large scales, such as would be observed in distant objects, and half that in individual GMCs. No evidence was found for H2 without CO emission. Our simulations of the radiative transfer in clouds are entirely compatible with these NH2/ICO values. The SFE implied is a factor 5 - 10 higher than what is observed in large local universe spirals.

Radiative transfer in the cloudy atmosphere

Radiative transfer in clouds is a challenging task, due to their high spatial and temporal variability which is unrivaled by any other atmospheric species. Clouds are among the main modulators of radiation along its path through the Earth’s atmosphere. The cloud feedback is the largest source of uncertainty in current climate model predictions. Cloud observation from satellites, on a global scale, with appropriate temporal and spatial sampling is therefore one of the top aims of current Earth observation missions. In this chapter three-dimensional methods for radiative transfer in cloudy atmospheres are described, which allow to study cloud-radiation interaction at the level needed to better understand the fundamental details driving climate and to better exploit remote sensing algorithms. The Monte Carlo technique is introduced which allows to handle nearly arbitrarily complex atmospheric conditions. The accuracy of the method is discussed by comparison between different models and with observations. Finally, we show some examples and discuss under which conditions three-dimensional methods are actually needed and when commonly-used one-dimensional approximations are applicable. This chapter builds upon the excellent overview of one-dimensional radiative transfer in ERCA Volume 3 [B. Pinty and M.M. Verstraete, ERCA 3, 67 (1998)].

Impact of a New Radiation Package, McRad, in the ECMWF Integrated Forecasting System

Abstract A new radiation package, “McRad,” has become operational with cycle 32R2 of the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). McRad includes an improved description of the land surface albedo from Moderate Resolution Imaging Spectroradiometer (MODIS) observations, the Monte Carlo independent column approximation treatment of the radiative transfer in clouds, and the Rapid Radiative Transfer Model shortwave scheme. The impact of McRad on year-long simulations at TL159L91 and higher-resolution 10-day forecasts is then documented. McRad is shown to benefit the representation of most parameters over both shorter and longer time scales, relative to the previous operational version of the radiative transfer schemes. At all resolutions, McRad improves the representation of the cloud–radiation interactions, particularly in the tropical regions, with improved temperature and wind objective scores through a reduction of some systematic errors in the position of tropical convection as a result of a change in the overall distribution of diabatic heating over the vertical plane, inducing a geographical redistribution of the centers of convection. Although smaller, the improvement is also seen in the rmse of geopotential in the Northern and Southern Hemispheres and over Europe. Given the importance of cloudiness in modulating the radiative fluxes, the sensitivity of the model to cloud overlap assumption (COA) is also addressed, with emphasis on the flexibility that is inherent to this new RT approach when dealing with COA. The sensitivity of the forecasts to the space interpolation that is required to efficiently address the high computational cost of the RT parameterization is also revisited. A reduction of the radiation grid for the Ensemble Prediction System is shown to be of little impact on the scores while reducing the computational cost of the radiation computations. McRad is also shown to decrease the cold bias in ocean surface temperature in climate integrations with a coupled ocean system.

Black carbon or brown carbon? The nature of light-absorbing carbonaceous aerosols

Abstract. Although the definition and measurement techniques for atmospheric "black carbon" ("BC") or "elemental carbon'' ("EC") have long been subjects of scientific controversy, the recent discovery of light-absorbing carbon that is not black ("brown carbon, Cbrown") makes it imperative to reassess and redefine the components that make up light-absorbing carbonaceous matter (LAC) in the atmosphere. Evidence for the atmospheric presence of Cbrown comes from (1) spectral aerosol light absorption measurements near specific combustion sources, (2) observations of spectral properties of water extracts of continental aerosol, (3) laboratory studies indicating the formation of light-absorbing organic matter in the atmosphere, and (4) indirectly from the chemical analogy of aerosol species to colored natural humic substances. We show that brown carbon may severely bias measurements of "BC" and "EC" over vast parts of the troposphere, especially those strongly polluted by biomass burning, where the mass concentration of Cbrown is high relative to that of soot carbon. Chemical measurements to determine "EC" are biased by the refractory nature of Cbrown as well as by complex matrix interferences. Optical measurements of "BC" suffer from a number of problems: (1) many of the presently used instruments introduce a substantial bias into the determination of aerosol light absorption, (2) there is no unique conversion factor between light absorption and "EC" or "BC" concentration in ambient aerosols, and (3) the difference in spectral properties between the different types of LAC, as well as the chemical complexity of Cbrown, lead to several conceptual as well as practical complications. We also suggest that due to the sharply increasing absorption of Cbrown towards the UV, single-wavelength light absorption measurements may not be adequate for the assessment of absorption of solar radiation in the troposphere. We discuss the possible consequences of these effects for our understanding of tropospheric processes, including their influence on UV-irradiance, atmospheric photochemistry and radiative transfer in clouds.

Radiative transfer in clouds with small‐scale inhomogeneities: Microphysical approach

Several recent publications have indicated that cloud droplets belonging to a particular size range may tend to cluster by forming small random particle groups imbedded in an otherwise uniform cloud. To analyze radiative transfer in clouds with such small‐scale inhomogeneities, we invoke the concepts of statistical electromagnetics. We show that as long as the assumptions of ergodicity and spatial uniformity hold, one can still apply the classical radiative transfer equation in which the participating extinction and phase matrices are obtained by averaging the respective single‐particle matrices over all the particles constituting the cloud. This result implies that comparisons of in situ and remote‐sensing retrievals of the cloud‐particle size distribution can be problematic and should be performed with caution.

Attenuation and impulse response for multiple scattering of light in atmospheric clouds and aerosols

Model phase functions for atmospheric clouds and aerosols typically comprise a narrow forward lobe (corona), a broad diffuse background, and a narrow backscattering peak (glory), which can reach relatively high values, especially for polyhedral scattering particles, such as hexagonal ice columns and plates. The influence of these three major components on the asymptotic and transient attenuation of the scattered light is compared for several analytic phase functions to assess the dependence of radiative transfer in clouds and aerosols on the choice of phase function. The impulse response (temporal evolution of the angular intensity distribution) is sensitive to the higher moments of the phase function and could prove to be a useful technique for inferring the optical scattering parameters of clouds and aerosols.

A New Method for Detection of Cirrus Overlapping Water Clouds and Determination of Their Optical Properties

Abstract The frequent occurrence of high cirrus overlapping low water cloud poses a major challenge in retrieving their optical properties from spaceborne sensors. This paper presents a novel retrieval method that takes full advantage of the satellite data from the Moderate Resolution Imaging Spectroradiometer (MODIS). The main objectives are identification of overlapped high cirrus and low water clouds and determination of their individual optical depths, top heights, and emissivities. The overlapped high cloud top is determined from the MODIS CO2-slicing retrieval and the underlying low cloud top is determined from the neighboring MODIS pixels that are identified as single-layer low clouds. The algorithm applies a dual-layer cloud radiative transfer model using initial cloud properties derived from the MODIS CO2-slicing channels and the visible (0.65 μm) and infrared (11 μm) window channels. An automated iterative procedure follows by adjusting the high cirrus and low water cloud optical depths until computed radiances from the dual-layer model match with observed radiances from both the visible and infrared channels. The algorithm is valid for both single-layer and dual-layer clouds with the cirrus optical depth <∼4 (emissivity <∼0.85). For more than two-layer clouds, its validity depends on the thickness of the upper-layer cloud. A preliminary validation is conducted by comparing against ground-based active remote sensing data. Pixel-by-pixel retrievals and error analyses are presented. It is demonstrated that retrievals based on a single-layer assumption can result in systematic biases in the retrieved cloud top and optical properties for overlapped clouds. Such biases can be removed or lessened considerably by applying the new algorithm.

Monte-Carlo algorithm for line-by-line calculations of thermal radiation in multiple scattering layered atmospheres

An algorithm for calculations of the longwave radiation in cloudy and aerosol slab atmospheres is described. It is based on the line-by-line and Monte-Carlo methods and is suitable for accurate treatment of both the gaseous absorption and the particulate multiple scattering in any spectral regions; other published algorithms as accurate as this can only make calculations in narrow spectral regions. It is recommended that this algorithm is well suited for radiation code validations as well as for theoretical investigations of radiative transfer in clouds and aerosols and satellite signal simulations.

Cloud modulation of zenith sky oxygen photon path lengths over Boulder, Colorado: Measurement versus model

Cloudy sky average photon path lengths were measured using the gamma band of oxygen near 628 nm on many days over Boulder, Colorado, using a zenith‐looking spectrograph with a resolution of 0.8 nm. The approach allows accurate measurement of the average photon path length. Days characterized by relatively extensive cloud cover are examined here, which exhibit very large path length variations as the cloud fields overhead evolve. The ability of a plane parallel line‐by‐line model with a single uniform cloud layer to predict the measured O2 path lengths has been tested by constraining the model to the (independently) observed surface irradiance. Overall, its performance is quite good, which demonstrates that the relationship between the diffuse transmission by clouds and the average photon path length is consistent with plane parallel radiative transfer calculations for the conditions studied. Thus large errors in the cloud radiative transfer of such models, as suggested by some recent cloud absorption measurements, are not supported by this study. Furthermore, direct observations in the blue and red spectral regions do not support even a 1% differential absorption between these spectral regions.