Differences in liquid cloud droplet effective radius and number concentration estimates between MODIS Collections 5.1 and 6 over global oceans

Differences in cloud droplet effective radius and cloud droplet number concentration (CDNC) estimates inferred from the Aqua MODIS Collections 5.1 and 6 cloud products (MYD06) are examined for warm clouds over global oceans for the year 2008. Individual pixel level retrievals for both collections are aggregated to 15 1° x 1° and compared globally and regionally for the three main spectral channel pairs used for MODIS cloud optical property retrievals. Comparisons between both collections are performed for cases in which all three effective radii retrievals are classified by the MODIS Cloud Product as valid. The contribution to the observed differences of several key MYD06 Collection 6 algorithm updates are also explored, with a focus on changes to the surface reflectance model, assumed solar irradiance, above cloud emission, cloud top pressure, and pixel registration. Global results show a neutral to positive (> 50 20 cm) change for C6-derived CDNC relative to C5.1 for the 1.6 μm and 2.1 μm channel retrievals, corresponding to a neutral to -2 μm difference in droplet effective radius. For 3.7 μm retrievals, CDNC results show a negative change in the tropics, with differences transitioning toward positive values with increasing latitude spanning -25 to +50 cm related to a +2.5 to -1 μm transition in effective radius. Cloud optical thickness differences were small relative to effective radius, and found to not significantly impact CDNC estimates. Regionally, the magnitude and behavior of the annual CDNC cycle are compared for 25 each effective radius retrieval. Results from this study indicate significant intercollection differences in aggregated values of effective radius due to changes to the pre-computed retrieval lookup tables for ocean scenes, changes to retrieved cloud top pressure, solar irradiance, or above cloud thermal emission, depending upon spectral channel. The observed differences between collections may have implications for existing MODIS derived climatologies and validation studies of effective radius and CDNC. 30


Introduction
Marine stratocumuli play a large role in the modulation of the Earth's radiative balance and hydrological cycle, owing to their persistence and large areal extent.They are maintained by moisture flux from the sea surface, longwave cooling, and dry air entrainment at the cloud top.Confined to the planetary boundary layer, most often in subtropical subsidence areas, they are susceptible to the influences of anthropogenic aerosols.Aerosol perturbation manifests changes though several mechanisms.Through the first indirect aerosol effect, anthropogenic influences potentially alter the degree of reflection of shortwave radiation by distributing a given quantity of cloud water over a larger number of droplets compared to an unperturbed cloud by supplying additional cloud condensation nuclei (CCN; Twomey, 1974).This increase in cloud albedo provides a local cooling effect as the cloud-top temperature (CTT) is near the sea surface temperature, producing a negligible longwave radiative ef-fect.With a decrease in droplet effective radius (r e ) from the addition of CCN, autoconversion efficiency can decrease, reducing precipitation efficiency (Albrecht, 1989) and potentially increasing the overall lifetime of clouds, augmenting the first indirect effect (Lohmann and Feichter, 2005).This naturally leads to a need for estimates of cloud droplet number concentration (CDNC) from observations to better understand the natural background and anthropogenic contributions to cloud droplet number concentration.
Satellites offer an opportunity to address this need through their global observational area spanning decades.The National Aeronautics and Space Administration's (NASA) Moderate Resolution Imaging Spectroradiometer (MODIS) is an instrument suited for this task.Flown aboard two Earthobserving satellites (Terra and Aqua), MODIS provides observations in the appropriate spectral bands to infer cloud optical thickness (τ ) and droplet effective radius through the bi-spectral method of Nakajima and King (1990).From the retrieved cloud optical thickness and droplet effective radius, cloud droplet number concentration can be estimated (Brenguier et al., 2000).NASA provides retrievals of the cloud optical properties through the MODIS cloud product (MOD06 and MYD06 for Terra and Aqua, respectively); note the MOD06 and MYD06 products are produced by the same algorithm (Platnick et al., 2015(Platnick et al., , 2017)).As passive imager cloud remote sensing science evolves, the MODIS cloud product undergoes periodic updates to the cloud screening and retrieval algorithms, which are collectively implemented and reprocessed as "collections".Collection 6 (C6) is the most recent release of the cloud product and includes several changes that propagate through to estimates of cloud microphysical properties relative to its predecessor, Collection 5.1 (C5.1;Platnick et al., 2015).Since there is an existing body of research based upon Collection 5.1 of effective radius (e.g., Zhang and Platnick, 2011;Painemal and Zuidema, 2011) and CDNC (e.g., Zeng et al., 2014;Ahmad et al., 2013), there is a need to understand how these retrievals differ between collections.In this study, we investigate and highlight the differences in retrieved cloud microphysical properties manifested both in effective radius and CDNC between both collections using 1 year of Aqua-MODIS observations of warm clouds over global oceans in an attempt to understand the differences in CDNC from the respective products.

Data and methods
The MODIS Collection 5.1 cloud product provides estimates of cloud optical thickness (τ ) and droplet effective radius (r e ) in addition to cloud screening, quality control flags, and ancillary data relevant to the computation of cloud droplet number concentration.For most parameters derived from visible, near, and shortwave infrared (IR) channels, the spatial resolution is 1 × 1 km at nadir.Some parameters, such as those derived from thermal infrared channels, are provided at a 5 × 5 km resolution.The C5.1 cloud product applies a Clear Sky Restoral (CSR) algorithm that identifies pixels that are expected to be only partially cloudy, and excludes them from the cloud optical properties retrieval, essentially limiting C5.1 data to cloudy scenes only.In addition, the C5.1 r e retrievals using the 1.6 and 3.7 µm channels are given as a difference relative to the 2.1 µm r e value and are therefore dependent upon the success of the 2.1 µm retrieval.
The Collection 6 cloud product offers improvements in retrievals of several relevant parameters for CDNC estimation.Changes to the optical and microphysical retrievals for C6 include improvements to the forward radiative transfer model used to create the precalculated cloud retrieval lookup tables (LUTs), re-registration of the visible and near-IR focal planes for MODIS-Aqua, and explicit reporting of all three spectral channel pair effective radius retrievals (Platnick et al., 2015); cloud-top property retrieval changes for C6 include additional cloud-top temperature and pressure retrievals at 1 × 1 km spatial resolution and improvements to the characterization of low-level clouds over water surfaces, among others (Platnick et al., 2015).With respect to the assessment of pixel quality, the CSR algorithm is still applied and its results are reported in the cloud product, but unlike C5.1 cloud optical property retrievals are attempted on pixels identified as partly cloudy and, if successful, are written to the cloud product separate from the heritage cloudy retrievals (Platnick et al., 2015).Note that pixels identified by the CSR algorithm as being not cloudy (e.g., false cloudy sunglint or thick aerosols such as smoke or dust) remain excluded from the cloud optical property retrievals.A year of MODIS-Aqua C5.1 and C6 data, corresponding to calendar year 2008, is used in this investigation.

Calculating CDNC
Relating cloud droplet number concentration to retrievals of optical thickness and droplet effective radius requires an estimate of cloud liquid water content.To attain CDNC, which is not a parameter in the MODIS cloud product, the vertical structure is assumed to follow the so-called adiabatically stratified cloud model (Brenguier et al., 2000).This model accounts for a linear increase in cloud liquid water content as a saturated air parcel rises through a cloud's vertical extent while undergoing droplet growth as the parcel cools moistadiabatically.As marine boundary layer clouds are typically shallow, the temperature is nearly constant, so it follows that the amount of condensate formed through adiabatic ascent can be assumed to be constant.The following relationships, adapted from Bennartz (2007), are presented as a summary of how CDNC relates to the adiabatic cloud model.The adiabatic liquid water content w AD at any point h above the cloud base can be expressed as where c w is the mass of condensate formed per unit volume for each meter of ascent and T is the cloud-top temperature.
The c w parameter is derived from the Clausius-Clapeyron relationship and is primarily a function of temperature and, to a lesser degree, pressure.For this study, an assumed cloudtop pressure (CTP) of 850 hPa is used, which corresponds to a geopotential height of approximately 1500 m.For marine stratocumulus regions such as the southeast Pacific, this height is near the mean top of the boundary layer (von Engeln and Teixeira, 2013) around Aqua's 13:30 local equatorial crossing time (Ho et al., 2015).
Relating the liquid water content profile to the cloud optical properties, Brenguier et al. (2000) express the adiabatic liquid water path as where the effective radius is the value at cloud top.Accounting for cloud-top entrainment of dry air, observed and modeled liquid water paths are typically near 80 % of the purely adiabatic value (Duynkerke et al., 2004;Pawlowska and Brenguier, 2003).The c w parameter in this study is therefore scaled to approximately 80 % of the purely adiabatic value, in order to avoid overestimation of the condensate observed by MODIS.
With the cloud optical parameters, Bennartz (2007) expresses CDNC as where k encapsulates the skewness and dispersion of the cloud droplet size distribution, Q is the scattering efficiency, and ρ l is the density of liquid water.Analytically, k = , where r v is the volume mean droplet radius.For marine boundary layer clouds k = 0.8 is assumed to be a representative value for this study (Brenguier et al., 2011;Martin et al., 1994), although Painemal and Zuidema (2011) suggest that k = 0.88 may be more representative for droplets near the cloud top, resulting in a narrower cloud droplet size distribution.As the size parameters for cloud droplets approach the geometric optics limit, the asymptotic value of 2 is assumed for Q.

Global-scale common pixel comparison
Common pixel scenes are those in which all three effective radius retrievals are valid for both collections with otherwise consistent selection criteria.A common pixel selection allows for an objective comparison of effective radius and CDNC estimates between both collections, free of the influence of pixel population differences due to, for example, reclassification of cloud phase due to changes to the cloud thermodynamic phase algorithm (Marchant et al., 2016), or where one or more effective radius retrievals have failed between collections (Cho et al., 2015).The global distribution of the total MODIS-Aqua common pixel count for calendar year 2008 is shown in Fig. 1.
Figure 2 shows mean annual C6 values and C6-C5.1 differences in MYD06 cloud droplet effective radius for common pixel scenes, aggregated globally to a 1 • × 1 • equal angle grid for calendar year 2008, from the 1.6, 2.1, and 3.7 µm channels.For r e,1.6 , C6 retrievals are smaller than C5.1 over most of the global oceans (differences generally less than 1 µm), with the largest differences (greater than 2 µm) occurring near the Baja California peninsula.In equatorial regions of the Indian and west Pacific basins, as well as over sea ice regions in the northern high latitudes above 60 • N, however, r e,1.6 differences are largely neutral to slightly positive.C6 retrievals of r e,2.1 are predominantly smaller than C5.1 (differences less than 1 µm), with the greatest disagreement between collections occurring in the tropics.For r e,3.7 , www.atmos-meas-tech.net/10/2105/2017/there is a latitudinal gradient in the retrieval differences, with positive inter-collection differences in the tropics (C6 > C5.1 by more than 2 µm) transitioning to slightly negative values (C6 < C5.1 by less than 2 µm) at higher latitudes.
The relationship between the inter-collection shifts in r e do not necessarily imply that the cloud liquid water path profile is steepened in C6 by the decrease in r e,1.6 and r e,2.1 and the increase of r e,3.7 at low latitudes.Contributing effective radii in this study are not required to conform to an adiabatically stratified profile and are likely subject to retrieval artifacts, complicating the determination of the liquid water content profile.This is especially the case at lower latitudes where convective clouds are more dominant and inhomogeneous relative to stratiform clouds.
Patterns similar to those of effective radius are mirrored in the respective CDNC C6 average and inter-collection difference plots of Fig. 3.As shown in Fig. 3d-e, for the N 1.6 and N 2.1 estimates, the differences are generally neutral to positive, i.e., C6 > C5.1, with values over +50 cm −3 occurring along coastal boundaries in the subtropical subsidence zones.Figure 3f shows that N 3.7 from C6 is smaller than C5.1 in the tropics and generally larger at latitudes higher that 45 • .Note that of the three spectral channels used for r e retrievals, the sensitivity of the 3.7 µm channel is weighted closest to cloud top (Platnick, 2000), and r e,3.7 is therefore considered to be the most appropriate effective radius for estimating CDNC based upon the assumption of the adiabatic model that the effective radius is the cloud-top value.Cloud optical thickness, which is another necessary retrieved parameter for estimation of CDNC, is neutral to slightly higher in aggregated retrievals for Collection 6 (not shown).For regions in which stratocumuli are often present, generally in the subtropical subsidence zones, C6 cloud optical thickness is within 0.25 of C5.1.This translates to a relative overestimation by C6 of roughly 2 %.Considering the relatively weak sensitivity of CDNC to changes in cloud optical thickness relative to droplet effective radius, this results in an increase of C6 CDNC by less than 1 %.

Algorithmic and calibration factors in inter-collection differences
Due to the strong sensitivity to effective radius (N ∝ r −5/2 e ) in Eq. (3), larger effective radii translate into lower relative CDNC (and vice versa), thus the r e differences shown in Fig. 2 largely explain the CDNC differences observed in Fig. 3.As previously noted, the use of common pixels in the analysis in Sect.3.1 removes the influence of pixel population differences between the two collections.The large effective radius differences between C5.1 and C6 observed in Fig. 2 are therefore the result of either changes to the MYD06 optical property retrieval algorithm itself or changes to upstream or ancillary products such as the Level 1b calibrated radiometric data or the MYD06 cloud-top property retrieval algorithm.
For C6, changes to the liquid phase cloud optical property retrieval algorithm were, for the case of r e,1.6 and r e,2.1 , limited primarily to the precomputed retrieval lookup tables.The change having the most significant impact on the r e retrievals is the use of new LUTs for retrievals over ocean, which incorporate an ocean surface bidirectional reflectance model directly into the LUT that uses the Cox-Munk wind speed and direction-dependent wave-slope distribution (Cox and Munk, 1954b, a).Previously in C5.1 the same LUT was used over land and ocean scenes, with the ocean surface assumed to be a Lambertian reflector with 5 % albedo (Platnick et al., 2017).Changes to input datasets, such as Aqua Level 1b re-registration of the visible and near-IR focal plane and updates to the cloud-top property retrievals, are also found to significantly impact the r e retrievals.Figures 5, 6, and 7 show a single granule case study illustrating the impacts of these changes.An Aqua-MODIS granule of open and closed cell marine stratocumuli acquired off the coast of Baja California on 9 June 2014 (21:55 UTC) was selected.The true color images, 1 km CTP, τ , and r e,2.1 , are shown in Fig. 5ad.A region of sunglint is visible slightly west of nadir below 15 • N. Also, a glory is present east of nadir above 23 • N, which is evident in the difference plots in Fig. 6.The C5.1 to C6 r e,1.6 , r e,2.1 , and r e,3.7 changes, defined here as the difference r e,λ (C6) -r e,λ (C5.1), are shown in Fig. 6a-c.Note that purple colors arise in regions where adjacent pixels have positive (red) and negative (blue) differences in effective radii.Subsequent rows in Fig. 6 show the individual impacts of the C6 changes listed above.
The impacts of the new LUT for ocean scenes on r e,1.6 , r e,2.1 , and r e,3.7 are shown in Fig. 6d-f.This change primarily impacts r e retrievals for optically thinner clouds (τ roughly < 2 to 3) and those over sunglint, lowering r e by up to 1 µm.The strength of this effect is dependent on the orthogonality of the solution space, and thus the impacts are largest for r e,1.6 and smallest for r e,3.7 .These r e decreases due to LUT changes can be offset or enhanced, depending on the location of the pixel relative to cloud illumination or shadowing, by the effects of the Aqua Level 1b re-registration of the visible and near-IR focal plane, shown in Fig. 6g-i that are again strongest on r e,1.6 due to enhanced LUT non-orthogonality (r e,3.7 is largely unaffected).The largest contribution to the r e,2.1 changes observed in Fig. 2 result from the updates to the cloud-top pressure retrievals, which are discussed in de- tail in Baum et al. (2012).The impacts on each channel are shown in Fig. 6j-l.For the liquid phase clouds in this granule, C6 CTP is largely higher than in C5.1 with increases of 100 hPa or more (not shown).Higher CTP, i.e., lower cloudtop height, results in larger above-cloud atmospheric gaseous absorption corrections, brighter top of cloud reflectance, and thus smaller r e from all three spectral channels, by up to 1 µm or more.The global yearly mean CTP change for 2008 is shown in Fig. 4a, again illustrating that C6 CTP is largely higher than in C5.1, particularly over the marine stratocumulus regions where CTP has been increased to better place the height of the capping inversions (Baum et al., 2012).In contrast, in regions where the liquid cloud occurrence is low (e.g., the tropics; see Fig. 1) C6 CTP slightly decreases.
While the above C6 changes largely explain the intercollection differences for r e1.6 and r e2.1 , algorithm updates for the r e3.7 retrieval were more extensive, involving fixes to known shortcomings in the C5.1 algorithm, and their net contribution to the dissimilar difference patterns shown in Fig. 2 are unclear.Retrieved r e3.7 decreases resulting from the CTP changes are enhanced by a correction to the assumed band-averaged solar irradiance (F 0 ) at 3.7 µm.For C6, F 0 is smaller than for C5.1 and earlier versions (10.9 W m −2 µm −1 for C6 versus 11.7 W m −2 µm −1 for C5.1), resulting in a decrease of r e,3.7 by 1 µm or more, as shown in Fig. 7a.However, above-cloud atmospheric emission, previously ignored in C5.1, is now accounted for in C6 and yields a sizably larger r e,3.7 , on the order of 1 µm or more, as shown in Fig. 7b.The cumulative effects of the C6 changes discussed here appear to explain the observed granule-level inter-collection r e,3.7 differences (Fig. 6c), which for the case study granule are strongly positive (C6 larger) or negative (C6 smaller) and are seemingly dependent on view angle.Nevertheless, it is difficult to determine how the granule-level r e,3.7 differences Figure 6.C6-C5.1 r e retrieval differences for the granule shown in Fig. 5 for r e1.6 , r e,2.1 , and r e,3.7 effective radii (a-c) in total and contributions from (d-f) lookup table differences, (g-i) pixel reaggregation, and (j-l) cloud-top pressure.
translate to the global aggregated differences in Fig. 2. In addition, given the extensive C6 changes and the known shortcomings of the C5.1 r e,3.7 retrievals, it is of limited benefit to further assess their inter-collection differences.
A final consideration regarding the inter-collection CDNC differences shown in Fig. 3 is the fact that CDNC estimates from the retrieved cloud optical properties also depend on the retrieved cloud-top temperature through the c w parameter.As discussed above and in Sect.2, a number of algorithm improvements were implemented in the C6 MYD06 cloud-top property product, including cloud-top retrievals at 1 km resolution in addition to the 5 km cloud-top retrievals that have heritage to C5.1 and earlier collections (Baum et al., 2012).For the C6 CDNC estimates in this study, the 1 km CTT is used.To assess the impact of differences in CTT on CDNC, the CTT values from C5.1 at 5 km and C6 at 1 km were aggregated to 1 • × 1 • for the common pixel population, shown in Fig. 4b.The differences in CTT reveal a relatively lower temperature (3 K) from C6 in subsidence regions at low latitudes, neutral to slightly positive differences at higher latitudes, and largely positive differences in the intertropical convergence zone.Although the CTP increases in subsidence regions between collections, CTP placement is independent of CTT changes between collections, which creates a seemingly contradictory situation when viewed together.To determine how the disagreement in CTT translates to CDNC apart from the cloud-top retrieval impacts to r e , the c w was also aggregated to 1 • × 1 • .The differences in c w,6 and c w,5.1 translate to a 0.5 % relative decrease in globally averaged N for C6, provided all other parameters are held constant.The largest observed difference is in the subtropical subsidence regions, where the decrease is near 4 %.Therefore, c w has only a marginal impact on Nfor any of the three absorption channels, a result that provides further evidence in favor of attributing inter-collection N differences largely to changes in the effective radius retrievals.

Regional common pixel comparison
Four maritime regions were selected for additional examination primarily for their proximity to sources of anthropogenic aerosols and generally high number of liquid cloud observations within their domains: southern Africa (SAF), South America (SAM), North America (NAM), and Southeast Asia (SEA).The boundaries of these regions are shown as white outlines in Fig. 1.
The annual CDNC cycles from C5.1 and C6 are given in Fig. 8a for the southern Africa region.Similar to the global results, the 1.6 and 2.1 µm CDNC from C6 are systematically higher (> 10 cm −3 ), while the C6 3.7 µm values vary from being lower (< 10 cm −3 ) to near agreement with C5.1.During the southern African dry biomass burning season between June and October (Roberts et al., 2009), the 1.6 µm C6 CDNC presents the greatest differences between collections, exceeding 30 cm −3 with a maximum CDNC of over 180 cm −3 .Interestingly, there is no discernable difference at 3.7 µm for either collection over the biomass burning season; however, it is quite muted relative to N 1.6 ; peak values during the season are roughly in agreement with Austral summer values between 90 and 100 cm −3 .For the region, absorbing aerosols from biomass burning often overlie the cloud deck, decreasing the reflected radiation observed by MODIS.With increasing aerosol optical depth, this can result in decreased retrieved cloud optical thickness and a general decrease in retrieved effective radius (Haywood et al., 2004).The sensitivity to overlying aerosols on the retrieved quantities depends on the degree of orthogonality of the retrieval solution space.Since CDNC estimates are significantly more sensitive to r e changes than optical thickness, the decrease in r e may result in an overestimation of CDNC if overlying absorbing aerosols are present.Results from the SAM region in Fig. 8b show a stratification and general response between collections, which is similar to SAF.For this region, a generally consistent aerosol emission rate gives rise to little variability (Huneeus et al., 2006) in CDNC, with annual cycle ampli- tudes limited to a range of approximately 20 cm −3 across all channels for both collections.For NAM, the common pixel comparison also yields similar results to the previous regions as shown in Fig. 8c.For this region, both the 1.6 and 2.1 µm channels also exhibit a 20-30 cm −3 difference for Collection 6 retrievals while 3.7 µm is in near agreement with Collection 5.1.All three channels have a uniform response to increases in CDNC coincident with a relatively early wildfire season that occurred during the spring and summer of 2008 in California (Brioude et al., 2009).
In contrast to the other regions in the study, the SEA domain yields a different stratification of CDNC shown in Fig. 8d, where N 3.7 > N 1.6 > N 2.1 for both collections with general agreement between N 3.7 retrievals.Disagreements between collections are less than 15 cm −3 for each channel.The different stratification may be a result of differing cloud dynamics in this region due to a weaker inversion strength, resulting in greater entrainment compared to the other regions which lie in subsidence zones.There is a strong seasonal response to aerosols visible in each channel coincident with the advection of aerosols by the Asian winter monsoon (Bennartz et al., 2011).
The adiabatic model assumes the retrieved effective radius is the near-cloud-top value.For an ideal adiabatically stratified cloud r e,3.7 > r e,2.1 > r e,1.6 , implying that N 1.6 > N 2.1 > N 3.7 , as each effective radius retrieval has a varying vertical penetration depth.None of the focus regions in this study exhibit this behavior.In addition to entrainment processes, retrieved droplet effective radius and therefore estimated CDNC can be influenced by drizzle processes (Suzuki et al., 2010) or 3-D radiative transfer effects and cloud field inhomogeneity at the sub-pixel level that violate the plane-parallel cloud assumption used in the optical properties retrieval (Zhang and Platnick, 2011).Since values in this study are aggregated to 1 × 1 • and averaged monthly with no screening criteria with respect to vertical stratification, it is likely these processes are contributing to the nonadiabatic stratification observed.

Conclusions
It has previously been shown, under the condition of adiabatic clouds, that estimates of the CDNC of warm marine liquid-phase clouds can be derived from passive satellite remote sensing observations.To this end, observations from MODIS on the Terra and Aqua satellites have seen wide use.The operational MODIS cloud product (MYD06 for Aqua) provides the cloud property retrievals necessary for computing CDNC estimates, namely cloud-top temperature, cloud optical thickness, and effective droplet radius (r e ).Collection 6, the most recent release of MYD06, includes numerous updates relative to its predecessor, Collection 5.1, that can propagate through to estimates of CDNC.Using 1 year ( 2008) of global MODIS-Aqua observations, inter-collection differences are investigated for CDNC derived from three independent spectral r e retrievals, namely from the 1.6, 2.1, and 3.7 µm channels.For the pixel population having successful r e retrievals from all three spectral channels in both collections, C6 r e,1.6 and r e,2.1 retrievals are smaller than those of C5.1, with the greatest differences generally found for r e,1.6 .These inter-collection r e differences result in relatively larger estimates of C6 CDNC for both channels.Nevertheless, CDNC from these two spectral channels offer similar annual cycles for both C5.1 and C6, regardless of the geographic region.N 3.7 differences in stratocumulus regions are generally smaller and, unlike N 1.6 and N 2.1 , the signs of the differences are rarely uniform throughout the annual cycle.Moreover, the global inter-collection differences in r e,3.7 and N 3.7 exhibit quite different behavior than those derived from the 1.6 and 2.1 µm channels.These differences are inherently attributable to the more wholesale C6 changes that addressed known shortcomings in the C5.1 r e,3.7 retrieval.But because the individual effects of the C6 changes can be quite large, are often of opposite sign, and have angular and potentially other unknown dependencies, it remains unclear how granule-level r e,3.7 differences translate to the global aggregated differences shown here.Additional research is necessary in order to quantify the contribution of these dependencies to the observed inter-collection differences of r e,3.7 and N 3.7 .Furthermore, it is recommended that quantitative use of the C5.1 r e,3.7 retrieval be avoided given its known shortcomings.
Among the updates for C6 is the independent reporting of r e for each wavelength, rather than as differences with respect to r e,2.1 , as well as the inclusion of retrievals of the so-called partly cloudy pixels previously discarded in C5.1.These changes offer additional options in the analysis of CDNC from MODIS.Several permutations in data screening scenarios beyond those used in this study are now available, which were not possible in C5.1, and will likely be explored in future studies.Caution will be warranted when interpreting the results of these and other future studies of r e and CDNC intercomparisons which may be subject to different screening criteria available in C6 which allows for the inclusion of more potentially problematic retrievals.The differences observed here are only for overcast scenes for pixels in which all three effective radius retrievals are valid for liquid phase clouds.Different screening criteria may exacerbate the effects of aggregated retrieval artifacts, which have not been considered in this study.

Figure 1 .
Figure 1.Annual total of Aqua-MODIS liquid phase cloud observations where all three effective radii retrievals are valid per 1 × 1 • grid box.Focus regions are overlain.

Figure 4 .
Figure 4. Cloud-top pressure (a) and temperature (b) differences between aggregated 1 km resolution C6 and 5 km resolution C5.1 for the common pixel population.

Figure 7 .
Figure 7. C6-C5.1 r e,3.7 retrieval differences due to (a) the change of the assumed solar irradiance (F 0 ) at 3.7 µm and (b) correctly accounting for above-cloud atmospheric emission previously ignored in C5.1.

Figure 8 .
Figure 8. Annual CDNC cycle plots for the (a) southern African, (b) South American, (c) North American, and (d) Southeast Asian regions.