UBC Faculty Research and Publications

Variability of Optical Depth and Effective Radius in Marine Stratocumulus Clouds Szczodrak, Malgorzata; Austin, Philip H.; Krummel, P. B. Oct 31, 2001

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


52383-Austin_AMS_2001_JAS2912.pdf [ 439.88kB ]
JSON: 52383-1.0041783.json
JSON-LD: 52383-1.0041783-ld.json
RDF/XML (Pretty): 52383-1.0041783-rdf.xml
RDF/JSON: 52383-1.0041783-rdf.json
Turtle: 52383-1.0041783-turtle.txt
N-Triples: 52383-1.0041783-rdf-ntriples.txt
Original Record: 52383-1.0041783-source.json
Full Text

Full Text

2912 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S q 2001 American Meteorological Society Variability of Optical Depth and Effective Radius in Marine Stratocumulus Clouds M. SZCZODRAK Rosenstiel School of Marine and Atmospheric Science, University of Miami, Miami, Florida P. H. AUSTIN Atmospheric Sciences Programme, University of British Columbia, Vancouver, British Columbia, Canada P. B. KRUMMEL CSIRO Atmospheric Research, Aspendale, Victoria, Australia (Manuscript received 31 January 2000, in final form 12 March 2001) ABSTRACT Radiance measurements made by the Advanced Very High Resolution Radiometer (AVHRR) at 1-km (nadir) spatial resolution were used to retrieve cloud optical depth (t) and cloud droplet effective radius (reff) for 31 marine boundary layer clouds over the eastern Pacific Ocean and the Southern Ocean near Tasmania. In the majority of these scenes (each roughly 256 3 256 km2 in extent) t and reff are strongly correlated, with linear least squares yielding a regression curve of the form reff } t1/5. This relationship is consistent with an idealized model of a nonprecipitating layer cloud in which 1) the average cloud liquid water content increases linearly with height at some fraction of the adiabatic lapse rate in a 1 km2 vertical column, and 2) the normalized horizontal variability of the cloud liquid water path exceeds the variability of a scaled measure of the cloud droplet number concentration. In contrast, other scenes of similar horizontal extent show little or no correlation between retrieved values of t and reff. These scenes include thicker clouds in which precipitation may be occurring, as well as cloud layers with spatially distinct regions of varying reff. In situ aircraft measurements were made simultaneously with six AVHRR overpasses as part of the Southern Ocean Cloud Experiment. The clouds sampled by these flights were significantly thicker than the typically 200- m-thick eastern Pacific stratocumulus, with large vertical and horizontal variability. On five of the six flights, aircraft measurements of the cloud-top effective radius were well matched by the satellite retrievals, and in two of these layers reff } t1/5. 1. Introduction Marine stratus and stratocumulus clouds play a sig- nificant role in the planetary energy budget. Globally, boundary layer clouds act to decrease the net radiative forcing by 15 W m22 due to their large reflectivity (Hart- mann et al. 1992). This reflectivity varies with cloud parameters such as cloud fraction, column-integrated liquid water, and the mean surface area of cloud drop- lets. Knowledge of the seasonal and spatial variability of these cloud parameters is a prerequisite for under- standing feedbacks between boundary layer cloud prop- erties and natural or anthropogenic climate change. Efforts to incorporate prognostic equations for cloud liquid water content and sulphate mass in global climate models have underscored the uncertainties inherent in Corresponding author address: Dr. Philip H. Austin, Atmospheric Science Programme, University of British Columbia, 6339 Stores Road, Vancouver, BC V6T 1Z4, Canada. E-mail: paustin@eos.ubc.ca predicting the impact of aerosol concentration on cloud reflectivity (Lohmann and Feichter 1997). At the same time, a new generation of cloud-resolving models, run with horizontal and vertical resolutions of tens of meters and domain sizes of kilometers, are making detailed predictions about the distribution of fluctuations in cloud liquid water and cloud droplet size. To better con- strain predictions at both ends of the modeling spectrum, observations of cloud liquid water path, particle size, and droplet number concentration across a range of scales are needed. Several recent studies have shown that satellite-based passive remote sensing can provide information on the cloud optical depth t (related to the extinction of the direct solar beam), the surface area-weighted mean ra- dius reff (or effective radius; see appendix A for defi- nitions), the liquid water path (lwp), and a measure of the droplet number concentration in the column (see, e.g., Han et al. 1994; Nakajima and Nakajima 1995; Platnick and Valero 1995; Han et al. 1995, 1998a; Greenwald and Christopher 1999; Han et al. 1998b). 1 OCTOBER 2001 2913S Z C Z O D R A K E T A L . TABLE 1. The grid system for values of cloud base (z), cloud thickness (Dz), solar zenith angle (u0), satellite zenith angle (u), optical depth (t), and effective radius (reff) for the lookup tables used for the satellite t , reff retrievals. Quantity Grid points z (km) Dz (km) u0 (8) u (8) f (8) t reff (mm) 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0 0.1, 0.2, 0.5, 1.0, 2.0 0, 5, 10, 20, 30, 35, 40, 45, 50, 55, 60, 65, 70 0, 5, 10, 20, 30, 35, 40, 45, 50, 55, 60 0–180 (every 108) 1, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 32, 34, 50, 70 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 30 Below we use the algorithm of Nakajima and Nakajima (1995) to retrieve t and cloud-top reff from 25 relatively shallow stratocumulus cloud layers located off the Cal- ifornia coast and six deeper and more variable layers over the Southern Ocean. These scenes are represen- tative of a larger set of 404 retrievals that we have done on 63 days of satellite data. In this paper we describe the distributions of optical depth and effective radius and their correlations, and propose a simple cloud model to account for the correlations between t and reff ob- served in 19 of the 31 cloud layers. In section 2 we briefly review the Nakajima and Na- kajima (1995) algorithm and discuss the uncertainties in the retrievals of t and reff, while in section 3 we describe the datasets used in this paper. In section 4 we present the resulting retrievals for thin and thick clouds, with joint probability densities and mean layer statistics. In section 4b we use in situ aircraft measurements to validate reff retrievals in thicker clouds. Section 5 con- tains a summary and discussion of the results. 2. Retrieval algorithm and uncertainties a. Retrieval algorithm We use the algorithm of Nakajima and Nakajima (1995) to retrieve cloud optical depth and cloud droplet effective radius given radiance measurements in Ad- vanced Very High Resolution Radiometer (AVHRR) channels 1 (centered at 0.63 mm), 3 (centered at 3.74 mm), and 4 (centered at 10.8 mm). The algorithm first calculates model radiances in the same wavelength range as AVHRR channels 1 and 3 given a three layer atmosphere with specified cloud-base height, cloud thickness, visible optical depth, and droplet effective radius. The absorption coefficient for the overlying at- mosphere is specified using LOWTRAN-7, with tem- perature and vapor soundings taken from the midlatitude summer profile (Kneizys et al. 1988). The satellite-mea- sured radiances are then used to choose the best-fit set of model parameters, minimizing the difference between observed and model radiances while iteratively cor- recting for the emitted radiation in channel 3, as de- scribed in Nakajima and Nakajima (1995). The param- eter set used in the model and the view and solar zenith angles included in the calculation are listed in Table 1. b. Retrieval error estimates Uncertainty in the retrieved values of cloud optical depth and cloud droplet effective radius arises from three sources of error. 1) Approximation error: the error introduced by as- sumptions made in the forward radiative transfer model, including the assumption of a vertically uni- form profile of reff with height and of climatological vapor and temperature profiles. Nakajima et al. (1991) used radiances calculated from simulated clouds for which reff increased linearly with height to test their t, reff retrieval and found the retrieved reff to be roughly 10% less than the cloud-top reff of the vertically inhomgeneous cloud. The same retriev- al overestimated t by 1%–5%, with uncertainty in the atmospheric water vapor profile adding 1%–3% to the t and reff uncertainties. Recently Brenguier et al. (2000) performed a similar calculation with a range of cloud profiles and found their retrieved reff to decrease from between 100% to 80% of the mod- eled cloud-top value as the cloud-top reff is increased from 6 to 18 mm. 2) Independent pixel approximation (IPA) error: this error is introduced by neglecting the horizontal ra- diative transfer between pixels. For overcast scenes, Chambers et al. (1997) estimate that this error for the retrieval of t can range between 10% for 30-m Landsat pixels to less then 5% for 5.7-km pixels at a solar zenith angle of 638. As Table B1 indicates, solar zenith angles for the AVHRR scenes were gen- erally between 508 and 608, with satellite view angles typically less than 308. Given the 1-km AVHRR na- dir pixel size we anticipate that the IPA error should fall within this range of values. The IPA introduces a negative bias in retrievals of t. The effect of the IPA on reff has not been investigated for cloud re- trievals, but we anticpate that smaller channel 3 ra- diances due to inhomgeneous clouds would result in larger values of retrieved reff. The Southern Ocean clouds, which were more irregular and had lower cloud fractions, may suffer from larger IPA errors. 3) Measurements (retrieval) error: this is dominated by the uncertainty of sensor calibration and digitization, estimated for t at about 15% for our range of t and solar angle values (Pincus et al. 1995). The AVHRR 2914 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S FIG. 1. The geographical location of 25 northeastern Pacific scenes used in section 4. Markers used for each scene are explained in more detail in section 4. ▫: scenes where for which t } reff; C: scenes with bimodal distribution of Nsat (see section 4 for the definition); n: scenes with thick clouds. channel 1 radiance measurements for the five sat- ellites were calibrated following Kaufman and Hol- ben (1993). Both Han et al. (1994) and Platnick and Valero (1995) also considered the errors listed above. Platnick and Valero (1995) estimated the worse-case net uncertainty in an AVHRR retrieval of reff to be 620%, accounting for the channel 3 measurement error, an unknown in- cloud water vapor absorption and droplet size distri- bution uncertainties. Han et al. (1994) similarly esti- mated retrieval errors for their study to be between 1 and 2 mm. This is in agreement with validation results using collocated aircraft measurements. Nakajima and Nakajima (1995) and Platnick and Valero (1995) showed differences between satellite and in situ esti- mates of the effective radius differ by roughly 10%, with good qualitative agreement in regions of varying reff . In a set of comparisons with a ground-based mi- crowave radiometer and a pyranometer, satellite retriev- als of reff and liquid water path lwp by Han et al. (1995) show lwp agreement within ø20% and reff agreement to within approximately 10%. 3. Data We have selected 25 scenes of 1 km 3 1 km local area coverage AVHRR data with dimensions of ap- proximately 256 km 3 256 km that are representative of a set of roughly 400 similarly sized retrievals, drawn from 67 different orbital swaths. The 1987 data were archived as part of the First International Satellite Cloud Climatology Program (ISCCP) Regional Experiment (FIRE) (July 1987), while data from June and July 1994 and 1995 were archived and processed at the University of British Columbia. The location for these scenes is shown in Fig. 1 and Table B1; we will also discuss six Southern Hemisphere summer scenes from the Southern Ocean near Tasmania, acquired during the Southern Ocean Cloud Experiment II (SOCEX II, February 1995; see Table B2). For each scene, completely overcast and clear pixels are identified using the spatial coherence analysis tech- nique of Coakley and Bretherton (1982). We limited our retrievals to the cloudy pixels for which the standard deviation of 2 3 2 arrays of the AVHRR channel 4 brightness temperature was below 0.5. For the north- eastern Pacific scenes the fraction of cloudy pixels meet- ing this threshold exceeds 80%; for the more variable Southern Ocean scenes the fully cloudy fraction falls to between 30% and 50%. Failure to exclude pixels with larger spatial variances results in the retrieval of anom- alously large droplets. This is due to a combination of two effects: 1) lower reflectivity due to cloud inho- mogeneity and 2) an overestimate of the channel 3 ther- mal emission due to sea surface contamination of the pixel. When the erroneous thermal emission is subtract- ed from the channel 3 radiance, the resulting low re- flectivity produces a mistaken estimate of larger cloud droplets. 4. Results a. Distributions of optical depth and effective radius Figure 2 shows log–log contour plots of reff versus t for scenes 1, 4, 5, and 24 of Fig. 1. Shown in the figure is a mean regression line (solid) of the form log(r ) 5 b log(t) 1 a.eff (1) The mean regression coefficients b and a for each 256 3 256 pixel scene are found by averaging 100 individual regressions of sets of independent pixels, subsampled to remove spatial correlation. We determine the corre- lation length for the t and reff fields using the two- dimensional semivariogram (second-order structure function): 2S (Dx) 5 ^[t(x 1 Dx) 2 t(x)] &,2 (2) where Dx is the (x, y) separation vector between two pixels, and ^ & represents an average over all pixels sep- arated by Dx. The autocorrelation length scale for the scene is then taken as the distance over which S2 reaches its maxi- mum; typically this is 7–10 km for both the t and reff fields (Isaaks and Srivastava 1989; Szczodrak 1998). The random subsets are selected with pixels separated by the autocorrelation length, and a linear regression is then performed on each of these subsets of roughly 300 1 OCTOBER 2001 2915S Z C Z O D R A K E T A L . FIG. 2. Log–log contour plots of reff and t from four cloudy scenes. The contours are of the frequency density h, where hD logtD logreff denotes the number of pixels with optical depth and effective radii in the range (logt , logt9 , logt 1 D logt; logreff , log , logreff 1 D logreff). The dot–dashed, labeled lines are isolines of lwp9reff given by (3c); the dotted labeled lines are isolines of Nsat given by (3b). Also shown on each panel is the linear least squares regression line (solid), with the slope b and its uncertainty. Clouds with both high (t, reff) correlations and b ø 0.2 are denoted by a ▫ in Fig. 1. pixels yielding the individual regression coefficients that are averged. Subsampling in this way minimizes the effect of local spatial correlations on the fit, so that the resulting correlation reflects scenewide behavior. The resulting mean slope and intercept a and b with uncer- tainties are shown in Fig. 2. We also calculate the Pear- son’s linear correlation coefficient for each subsampled fit, assuming uncorrelated errors of 15% for t and 20% for reff as discussed in section 2b. The correction for the bias introduced into the regression by the log–log trans- formation is negligible for these scenes (Jansson 1985). As the solid line in each of the four panels of Fig. 2 indicates, the t, reff scatterplots in each of these four panels are well fit by a relationship of the form reff } t1/5. This is also true of 13 other scenes of the 25 shown in Fig. 1 (marked with square boxes), which all yield mean log–log regression slopes b 5 1/5 within the sam- pling uncertainty, with correlation coefficients greater than 70%. Of the 404 images analyzed for this region, 40% satisfy these two the criteria. In appendix A we review relationships between the optical depth t, the effective radius at cloud top reff, the liquid water path lwp, and the number concentration N that hold for an idealized layer cloud in which cloud liquid water content (lwc) increases linearly with height and cloud droplet number concentration N is constant with height: 1/5 22/5 1/5r 5 a N t (3a)eff 1 sat 1/2 1/2 25/2N 5 a t r ø kN/Ïb, (3b)sat 1 eff lwp 5 a tr . (3c)0 eff The definitions for the coefficients a0 and a1 and their expected range of variability are given in appendix A, along with that of k, an empirical parameter that relates reff to the volume mean radius rvol and b, which describes the rate of increase of lwc with height. The very simple cloud model described by (3) predicts that t, reff retriev- als for clouds with low variability in the quantity should exhibit the power-law relationship be-1/5 21/2a N1 sat tween t and reff observed in Fig. 2. Below we will focus on relative variability of Nsat and lwp, and use satellite- retrieved values of reff as estimates of cloud-top reff in 2916 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S FIG. 3. Histograms of t, reff, Nsat, and lwp for the four scenes of Fig. 2. (3). As mentioned in section 2b, this produces an un- derestimates of reff in clouds for which reff is increasing with height. Assuming uncorrelated errors in t and reff , the impact on the derived quantities in (3) would be a 10% underestimate lwp and a 25% overestimate in Nsat given a 10% underestimate in reff. The assumptions in (3) that number concentration is roughly constant with height in the cloud and that liquid water content increases linearly with height are consis- tent with aircraft and balloon observations of both pre- cipitating and nonprecipitating stratocumulus clouds in the north Atlantic, northeastern Pacific, and mid-Atlan- tic (e.g., Nicholls 1984; Caughey and Kitchen 1984; Austin et al. 1995; Khairoutdinov and Kogan 1999; Brenguier et al. 2000). A linear increase of liquid water content is also found in the mean lwc profiles of many numerical simulations of layer clouds, including models based on higher-order closure (e.g., Wang and Wang 1994) and two-dimensional large eddy simulations with bin-resolved microphysics (e.g., Khairoutdinov and Ko- gan 1999). There are also cases, however, in which both the assumptions of approximately constant droplet num- ber concentration and a near-linear vertical profile of liquid water content are invalid. For example, decoupled cloud layers, and clouds in which droplets at cloud top are being removed by evaporation, violate both these assumptions and break the connection between Nsat and cloud-top number concentration assumed in (3). To help characterize these cloud layers we have in- cluded isolines of Nsat 5 t1/2 and lwp 5 a0 t reff1/2 25/2a r1 eff in Fig. 2 and in the histograms for these scenes (Fig. 3), with values for a0 and a1 given in appendix A. Sev- eral studies have shown that the retrieval of lwp via (3c) from both AVHRR and Geostationary Operational En- vironmental Satellite (GOES) imager measurements of t and reff compares well with the surface microwave measurements of stratocumulus clouds of Han et al. (1995) and Greenwald et al. (1999). For example, Greenwald et al. (1999) report a root-mean-squared dif- ference of 17 g m22 between GOES imager and surface microwave data in a cloud layer with lwp values be- tween 0 and 200 g m22. In situ and satellite comparisons for Nsat are more difficult, because they require simultaneous knowledge of the cloud microphysics (N, k) and the bulk thermo- dynamics (b, a1). Han et al. (1998b) compare satellite retrievals of a related quantity (the column number con- centration Nc or product of the droplet number N and layer thickness H), for four aircraft flights in marine stratocumulus, and show 2 factor of 2 agreement be- 1 OCTOBER 2001 2917S Z C Z O D R A K E T A L . tween in situ and retrieved Nc values. In light of the uncertainties itemized in section 2b and appendix A and the assumptions underlying the relation between t, reff and the droplet concentration N in (3), we will focus here on contrasts between Nsat and lwp variability, rec- ognizing that large variations in Nsat could originate from fluctuations in cloud microphysics or departures from the assumptions underlying (3c), while small Nsat vari- ance could arise from compensating fluctuations in some or all of N, k, and b. In Fig. 3 we compare the t, reff, Nsat , and lwp distri- butions for the scenes shown in Fig. 2. The mean liquid water paths for scenes 1, 4, 5, and 24 range from 80 to 116 g m22, which are typical of thinner clouds observed by ground-based and satellite microwave radiometer measurements for this region (Han et al. 1995; Zuidema and Hartmann 1995). The Nsat distributions for all four scenes are nearly normal, with values of skewness de- fined as 3(N 2 N )sat satB 5 , (4) 3s ranging from 1023 to 1022, or 5–10 times smaller than that for the lwp distributions. The lwp and optical depth distributions for scenes 1, 4, and 5 are positively skewed, and are well fit with either gamma or lognormal distributions. In contrast, scene 24 is characterized by a unimodal Nsat distribution and bimodal distributions of t, reff, and lwp. Images of these retrieved values (not shown) show that the scene contains two distinct regions with lwp maxima of 33 and 151 g m22 but similar Nsat distributions with modes of 78 and 83 cm23, respec- tively. While the power-law relationship between t and reff shown in Fig. 2 occurs frequently, the majority of the 400 scenes for which retrievals have been performed have more complex t, reff variability. Figures 4–7 show representative examples of these more complex clouds. Figures 4a and 4b shows two scenes (9 and 10) that are similar to the bimodal scene 24 of Fig. 2d, with spatially distinct regions with differing reff distributions. To high- light this spatial separation we have shown images of the two scenes in Figs. 4c and 4d, with the pixel colors assigned according to whether the Nsat values are greater or less than 85 and 30 cm23, respectively. These two scenes differ from scene 24 in that the Nsat statistics are distinct in each subregion, as indicated by the bimodal Nsat distributions in Fig. 5c. The contrasting joint t, reff distributions are consistent with other AVHRR retrievals showing sharp contrasts between cloud properties across 100-km regions (e.g., Nakajima and Nakajima 1995). Scenes 7 and 8 (not shown) also exhibit this kind of spatially distinct variation in Nsat . Figures 6 and 7 show a third type of t, Nsat variability found in thicker clouds, defined as layers with mode optical depths larger than 20. In these four layers the mean liquid water path exceeds 100 g m22, while the mean reff for each scene is between 8.5 and 12.5 mm. Scene 14, from 16 July 1987, is also shown as a scat- terplot in Nakajima and Nakajima (1995; their Fig. 17, panel B-3). This scene has pixels with reff . 15 mm, for t , 20, but smaller reff values at larger t. As Na- kajima and Nakajima (1995) remark, the absence of pixels with both large lwp and large reff may be due to the removal of water by precipitation in thick cloud with low number concentrations, where droplets are large enough to permit collision–coalescence. In contrast, however, scene 10 (Fig. 4) shows retrieved pixels with effective radii larger than 20 mm at optical depths great- er than 20 and lwp . 270 g m22. The means of the distributions of t, reff, lwp, and Nsat for all 25 scenes are summarized in Fig. 8. We distin- guish in the figure between the 17 scenes for which the linear regression of reff , t, as described for Fig. 2, yields b 5 0.2 within the uncertainty of the fit, and other clouds (circled scene numbers in Fig. 8). The criterion that b 5 0.2 within the regression uncertainty is met by 17 of the 25 scenes in Fig. 1. Mean optical depths for these scenes range from 7.8 to 20, while mean effective radii vary from 6.6 to 13 mm, with larger mean optical depths generally corresponding with larger mean reff. In con- trast, mean values of t and reff are uncorrelated for the eight (circled) scenes which for which b ± 0.2. Figure 9 shows the means of the two derived quan- tities Nsat and lwp together with the standard deviation for each quantity. As expected, for all scenes in which regression yields reff } t1/5, the Nsat standard deviation (sN sat) normalized by the mean is more than three times the size of the normalized standard deviation of lwp. The trend for these scenes is for higher values of Nsat to occur in thinner cloud. There is little systematic geo- graphical variation in Nsat value with distance offshore; the location of scenes with Nsat . 200 cm23 ranges from 100 to more than 1000 km offshore. There is also no obvious retrieval bias with view angle or optical depth: both large and small normalized Nsat deviations occur across the entire range of t and lwp, and across a mean view angle range of 58–488 (see Table B1). b. In situ measurements In this section we present coincident aircraft and sat- ellite measurements of stratocumulus cloud layers ob- served during six flights flown to the west of Cape Grim, Tasmania, as part of SOCEX II. Boundary layer clouds during SOCEX II were typically cumulus penetrating into stratocumulus, with the solid cloud layer being ap- proximately 400–600 m thick, and overall cloud layer being approximately 1 km thick. The six analyzed SO- CEX II flights are listed in Table B2, together with the time of the matching satellite overpass. Thermodynamic profiles on most flights show discontinuities indicating decoupling between cloud and subcloud layers. All flights with the exception of flight 7 reported drizzle, with the largest rainfall rates reaching 50–60 mm day21 2918 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S FIG. 4. (a) and (b) As in Fig. 2, but for two scenes with strongly bimodal distributions of Nsat. The heavy dashed Nsat isoline in each contour plot separates large reff pixels from small reff. These pixels are then colored black (small reff) and gray (large reff) in the maps below each contour plot [(c) and (d)]. Clouds with bimodal Nsat distributions are denoted by a C in Fig. 1. and 20–30 mm day21 on flights 8 and 11, respectively. Effective radii reported here are calculated using drop- size distributions measured using the forward scattering spectrometer probe (FSSP) and 2DC probes on board the F27 aircraft as described in Boers et al. (1996). Sampling and sizing errors in the FSSP introduce un- certainties of at least 10% in droplet radius and 10% in number concentration (Baumgardner and Spowart 1990). Figure 10 shows a vertical profile of the flight path of a typical SOCEX II flight mission (flight 10); SOCEX flight missions consisted of several horizontal stacks flown at different levels within a cloud layer. Flights lasted 3–4 h, with the horizontal in-cloud legs covering 30–35 km over 7 min. The aircraft drifted with the wind, so that during a typical flight of 150 min a vertical plane 30–35 km long was repeatedly sampled. As Fig. 10 illustrates, there were typically six to eight encounters within 100–200 m of cloud top, ranging from samples of several seconds on ascent/descent to 7-min horizontal samples on level legs. The flight pattern shown in Fig. 10 allows us to obtain a picture of the vertical structure of the cloud layer over the 150-min period of a stack set. Based on this, we are able to make a rough estimate of the cloud adiabaticity b for these clouds by taking the composite vertical liquid water profile from the stacks and comparing it to an adiabatic profile taken from the mean cloud base. Figure 11 shows these estimates of b from the six SOCEX II flights plotted against total droplet number concentra- tion N (with its standard deviation), along with single sounding estimates for the FIRE region of section 4a taken from Austin et al. (1995). The error bars for the SOCEX II flights are calculated from the uncertainty in dlwc/dz given by the least squares fit to the stack profile data, neglecting the FSSP sizing errors and the uncer- tainty associated with the variable cloud base charac- teristic of the SOCEX clouds. The figure shows that the 1 OCTOBER 2001 2919S Z C Z O D R A K E T A L . FIG. 5. Histograms of t, reff, Nsat, and lwp for the two scenes of Fig. 4. FIG. 6. As in Fig. 2, but for four scenes with thicker cloud layers. Thick cloud scenes are denoted with a n in Fig. 1. 2920 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S FIG. 7. Histograms of t, reff, Nsat, and lwp for the four scenes of Fig. 6. FIG. 8. Scatterplots of the mean values of reff, t, Nsat, and lwp for all 25 northeastern Pacific scenes of Fig. 1. The eight scenes with circled numbers correspond to scenes with bimodal Nsat distributions or mean t . 20. 1 OCTOBER 2001 2921S Z C Z O D R A K E T A L . FIG. 9. As in Fig. 8, but for lwp vs Nsat. Solid lines indicate 6 one standard deviation about the mean value. FIG. 10. Aircraft vertical cloud sampling pattern (horizontal stacks) in SOCEX II flight 10. Time is in min from the start of the flight. Flight average cloud top and cloud base are indicated by the dashed lines. 2922 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S FIG. 11. Scatterplot of cloud adiabaticity b 5 (dlwc/dz)/(dlwcadiab/dz) vs cloud droplet number con- centration for the six SOCEX II flights (numbers) and for FIRE soundings (C) [using FIRE data from Austin et al. (1995)]. Vertical error bars from the least squares fit to the stack profile of lwc, horizontal error bars are 6 one standard deviation. SOCEX II b values tend to be slightly lower than those taken from FIRE, which is consistent with the obser- vations of both penetrative cumulus convection and de- coupling in the thermodynamic profiles, and drizzle re- sulting in the removal of cloud liquid water for these clouds. These values can be compared with measure- ments of 0.4 , b , 0.7 for wintertime stratocumulus in the SOCEX II region (Boers et al. 1996). Figure 12 shows the aircraft estimates of reff and t together with log–log contour plots of the coincident AVHRR retrievals. In each panel, the values of reff at or near cloud top have been denoted with circles, while aircraft measurements of reff below cloud top are shown as asterisks. We have assigned an approximate optical depth to each aircraft reff value by integrating the mean sounding produced by all the level legs in a stack. In clouds with large horizontal inhomogeneities such as these, this will be a poor estimate of t in any particular pixel. Given the irregular cloud base, penetrating cu- mulus turrets, and convective lines observed in these layers, we expect that the actual column-integrated op- tical depths might be significantly higher than those in- dicated by the time-averaged mean values presented here. Given this limitation in the aircraft t estimates, Fig. 12 does show good agreement between satellite and in situ estimates of cloud-top reff. For flights 6 and 11 the satellite retrieves large values of t with little variability in reff. For flight 6 the mean values for the in situ and satellite-estimated reff are within 0.7 mm (13.9 mm for AVHRR vs 13.2 for the aircraft). Similar agreement is found in flight 11 (14.0 vs 13.5 mm). In flights 7 and 8 the cloud-top values of reff varied between 5–10 (flight 7) and 9–12 mm (flight 8), a range that is encompassed in the satellite retrievals for these layers, although the aircraft sampled significantly thinner cloud than the sat- ellite retrievals for flight 7. For flight 9 the aircraft- measured values of reff are smaller than 10 mm, in con- trast to 12-mm radii found by the satellite in the thinest pixels subject to retrieval. In this broken cloud, as in flight 7, it is likely that some flight-leg pixels have been rejected because of failure to meet the spatial coherence criteria. Subvisible cirrus contamination is also a pos- sibility for this cloud; there are patches of cirrus away from the flight line in this image, although no high clouds are discernible in the Coakley–Bretherton dia- gram for flight 9. Overlying each of the contour plots in Fig. 12 is an Nsat isoline equal to the mean Nsat for the scene. For flights 9 and 10 this isoline is also within 1%–2% of the value of the log–log regression, which gives an av- erage slope of b 5 0.2 6 0.01 in each case. Although we expect these complicated, precipitating cloud layers to both have larger uncertainties in the retrieved values of t and reff and larger variability in b and k, we can qualitatively compare the distributions of Nsat 5 1 OCTOBER 2001 2923S Z C Z O D R A K E T A L . FIG. 12. Log–log contour plots of the retrieved t, reff frequency density h (as in Fig. 2) for the six SOCEX II flights. Open circles denote aircraft-measured reff at or near cloud top, while asterisks denote reff measurements below cloud top. Nsat isolines for the mean Nsat for each scene are also shown with the mean Nsat values indicated in each panel. t1/2 shown in Fig. 12 with the in situ measure-1/2 25/2a r1 eff ments of N and b given in Fig. 11, assuming that N 5 Nsat /k. Taking k 5 0.8, using the mean b for eachÏb flight gives satellite-retrieved values of N for flights 6, 10, and 11 of N 5 77, 176, and 76 cm23, respectively. These are within the range of N 6 sN given in Fig. 11, although the maximum in situ number concentrations observed in flights 10 and 11 (250 and 110 cm23) are lower than the maximum inferred from Nsat (350 and 181 cm23, respectively). As noted in section 4a we ex- pect satellite-retrieved reff to underestimate the cloud- top reff, which would produce overestimates of Nsat. There is little agreement between the in situ N and the Nsat estimates for the other flights, again possibly due to a sample mismatch between the Lagrangian aircraft trajectories and the location of successfully retrieved pixels. 5. Discussion Several authors [e.g., Han et al. (1994, p. 493); Na- kajima and Nakajima (1995, p. 4057); Han et al. (1998a)] have noted the tendency for satellite-retrieved reff to increase with increasing t in thinner clouds, and decrease with increasing t in optically thicker clouds. In light of Figs. 2, 4, and 6, we can be more quantitative about the variability of t and reff in thinner clouds (t , 20): we find that a power law of the form reff } t1/5 captures much of the variability between t and reff in 19 of the 31 layers we presented. Twenty-five of these scenes are drawn from a larger sample of 404 images taken from 67 satellite swaths during several summers in the northeastern Pacific. Of these 404 scenes, 40% also show a power-law relationship with an exponent equal to 1/5 within the uncertainty of the log t–log reff regression. Other scenes in this ensemble with similar values of optical depth, cloud fraction, satellite view, and solar zenith angles show larger reff variability and little correlation between t and reff. Recently, Brenguier et al. (2000) have inferred dis- tributions of cloud droplet number concentration from remotely sensed t, reff retrievals. Using an aircraft- mounted two-channel radiometer flown in the mid-At- lantic during the Second Aerosol Characterization Ex- periment, they retrieved t and reff, while simultaneously measuring cloud thickness and droplet size spectra from a second in situ aircraft. Figure 8 of their paper maps the two-channel reflectances to droplet number concen- tration assuming an adiabatic liquid water content pro- file and constant number concentration in the column beneath the radiometer; there is good agreement be- tween the number concentration inferred from this sim- 2924 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S ple model and the in situ samples, and a strong corre- lation of the two-channel radiance along isolines of con- stant number concentration. It is possible to make a similar interpretation of the satellite retrievals presented here, provided that the as- sumptions of a vertically uniform mean droplet con- centration and linearly increasing mean liquid water content hold on the 1-km horizontal scale of the satellite pixels. In that case isolines of Nsat } t1/2 map to25/2reff constant values of kN/ . We expect both precipitationÏb and entrainment to produce a negative correlation be- tween N and b (cf. Fig. 11), while Martin et al. (1994) show evidence for negative correlation between N and k, so that if the simple model holds, Nsat variability likely places only a lower bound on the variability of N in these thinner layer clouds. In the survey of northeastern Pacific scenes from which the 25 retrievals are taken we find that clouds with mean optical depths greater than 20 typically ex- hibit lower (t, reff) correlations. As SOCEX flights 9 and 10 show, however, a power law of the form reff } t1/5 can at least occasionally describe variability in thick- er clouds with significant precipitation. Simultaneous measurements of lwp, cloud-top reff, and optical depth, combined with a quantitative description of t, reff re- trieval uncertainties using three-dimensional radiative transfer models, would be very useful in determining the underlying physical constraints on the variability in lwp and cloud droplet number concentration for these layer clouds. Acknowledgments. We are grateful to Takashi Na- kajima and Teruyuki Nakajima for providing the source code for the retrieval algorithm used in this paper, along with numerous helpful discussions, and to Jorgen Jen- sen, whose comments improved this paper. The paper was improved by comments from Teruyuki Nakajima and two anonymous referees. Our research was sup- ported by grants from the Canadian Atmospheric En- vironment Service, the National Science and Engineer- ing Research Council of Canada, and the Common- wealth Scientific and Industrial Research Organisation of Australia (CSIRO). SOCEX was partly funded by the CSIRO Office for Space Science and Applications (COSSA), the Cooperative Research Centre for South- ern Hemisphere Meteorology, and CSIRO’s Multi-Di- visional Climate Change Research Program. Thanks for SOCEX work are due to Jorgen Jensen, Reinout Boers, Sunhee Lee, Dan Gogoasa, Bernard Petraitis, and other staff at CSIRO Atmospheric Research. APPENDIX A Simple Relations among t, reff, Nsat, and lwp Several authors have derived the relationship between the retrieved variables t and reff, the number concen- tration N, and the liquid water path lwp for a simple adiabatic cloud layer (e.g., Pontikis 1993; Boers and Mitchell 1994). Below we adopt the notation of Bren- guier et al. (2000), using an additional term (b) to rep- resent the departure of the liquid water content profile from adiabatic. Denoting the number concentration at height z with radii between r, r 1 dr as n(r, z)dr, the definitions for the total droplet concentration N, effective radius reff , the volume mean radius rvol , the liquid water content w, liquid water path lwp, and the optical depth t are given by ` 3n(r, z)r drE 0 r (z) 5 , (A1a)eff ` 2n(r, z)r drE 0 ` N(z) 5 n(r, z) dr, (A1b)E 0 `4 3w(z) 5 r n(r, z)pr dr, (A1c)w E3 0 1/33w(z) r (z) 5 , (A1d)vol 1 24pN(z)rw H lwp 5 w(z) dz, (A1e)E 0 H ` 2t 5 Q n(r, z)pr dr dz, (A1f)ext E E[ ] 0 0 where rw is the density of water, z is the height above cloud base, H is the cloud thickness, and Qext is the mean scattering efficiency. We will assume that the mean number concentration is roughly constant with height in the 1 km2 AVHRR pixel (N(z) 5 N), and that the mean liquid water content increases linearly with height at some fraction b of its adiabatic value w(z) 5 bC z,w (A2) where Cw, termed the moist adiabatic condensate co- efficient by Brenguier et al. (2000), is a weak function of temperature and pressure, varying between 1.8 and 2.25 3 1023 g m23 m21 in the temperature range 280– 290 K at 980 hPa. Substituting (A2) into (A1f ), inte- grating, and using the empirical relationship between and found by Martin et al. (1994) gives a slightly3 3r reff vol modified version of Eq. (12) of Brenguier et al. (2000): 3 2/3 2/3 1/3 5/3t 5 pQ A b (kN ) H , (A3)ext5 where A 5 Cw/(4/3)prw, and the coefficient k was found by Martin et al. (1994) to vary from 0.67 to 0.80 60.07 with higher values associated with maritime air masses. 1 OCTOBER 2001 2925S Z C Z O D R A K E T A L . Sensitivity tests using simple cloud models as well as comparisons with in situ aircraft measurements in- dicate that the retrieved value of reff is typically within 85%–90% of its value at cloud top (e.g., Nakajima and King 1990; Nakajima and Nakajima 1995). Making the approximation that this retrieved value, reff,sat, is equal to reff(H) and substituting k1/3rvol(H) for reff in (A3) gives 21 2 5t 5 a N r , (A4a)1 sat eff 5A a 5 , (A4b)1 3pQext with Nsat 5 kN/ .Ïb The relationship for lwp is calculated similarly, sub- stituting (A2) into (A1e), integrating, and using the re- lationship between rvol, reff, and w at cloud top to obtain lwp 5 a tr , (A5a)0 eff 10rw a 5 . (A5b)0 9Qext Values for Nsat and lwp in this paper are calculated assuming Qext 5 2 and Cw 5 Cw(temperature 5 285 K) 5 2 3 1023 g m23 m21; 285 K is a representative cloud- base temperature for both the northeastern Pacific and SOCEX clouds. APPENDIX B Satellite Scene Locations TABLE B1. Year/month/day and satellite number and orbit for the images of Fig. 1. Scene No. Date/orbit ID Lat (8N) Long (8W) Zenith (8) View (8) 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 87/07/16 N10-4294 87/07/09 N10-4194 87/07/12 N10-4237 94/07/17 N11-29943 87/06/23 N10-3951 94/07/16 N11-29929 87/06/23 N10-3951 87/07/14 N10-4265 95/06/15 N12-21220 87/07/14 N09-13324 94/07/17 N11-29943 94/07/17 N11-29943 87/07/14 N09-13324 87/07/10 N10-4208 94/07/16 N11-29929 87/07/12 N10-4237 87/07/16 N10-4294 94/07/12 N11-29802 87/07/31 N10-4507 94/06/23 N11-29802 94/07/16 N11-29830 35.03 28.66 34.89 35.38 30.61 24.40 38.40 28.19 24.86 29.23 38.27 37.93 33.84 30.60 23.63 33.53 32.11 34.89 29.29 25.54 40.32 2132.46 2124.62 2129.88 2127.93 2131.16 2121.11 2128.62 2117.65 2131.89 2126.94 2124.20 2124.93 2123.84 2121.00 2126.02 2124.71 2133.30 2129.88 2126.12 2128.83 2132.12 54.6 50.8 54.9 51.2 64.8 64.8 50.1 54.3 63.0 36.1 57.2 57.7 39.7 39.9 60.5 50.6 55.4 58.1 60.7 55.3 56.1 5.0 7.5 14.5 24.8 14.6 14.5 13.5 29.7 6.8 5.4 4.0 5.2 12.3 28.6 5.0 20.4 12.0 21.6 7.9 34.4 48.0 22 23 24 25 94/06/23 N11-29717 95/07/14 N14-20266 87/07/07 N10-4166 87/07/04 N10-4123 38.20 34.11 29.62 31.52 2132.30 2124.34 2124.81 2121.42 59.9 51.3 48.3 52.9 48.0 45.5 11.9 10.0 TABLE B2. SOCEX flight missions and coordinated satellite overpasses. Flight No. Day in 1995 Flight time (local) Satellite orbit No. Satellite time (local) 6 7 8 9 10 11 1 Feb 1 Feb 6 Feb 8 Feb 8 Feb 9 Feb 0941–1410 1537–1908 1511–1840 0945–1330 1515–1902 1003–1348 NOAA-9 52266 NOAA-14 00462 NOAA-14 00532 NOAA-9 52365 NOAA-14 00561 NOAA-9 52379 1007 1549 1455 1017 1614 1004 REFERENCES Austin, P., Y. Wang, R. Pincus, and V. Kujala, 1995: Precipitation in stratocumulus clouds: Observational and modeling results. J. Atmos. Sci., 52, 2329–2352. Baumgardner, D., and M. Spowart, 1990: Evaluation of the Forward Scattering Spectrometer Probe. Part III: Time response and laser inhomogeneity limitations. J. Atmos. Oceanic Technol., 7, 666– 672. Boers, R., and R. M. Mitchell, 1994: Absorption feedback in stra- tocumulus clouds: Influence on cloud-top albedo. Tellus, 46A, 229–241. ——, J. Jensen, P. Krummel, and H. Gerber, 1996: Microphysical 2926 VOLUME 58J O U R N A L O F T H E A T M O S P H E R I C S C I E N C E S and short-wave radiative structure of wintertime stratocumulus over the southern ocean. Quart. J. Roy. Meteor. Soc., 122, 1307– 1341. Brenguier, J.-L., H. Pawlowska, L. Schüller, R. Preusker, J. Fischer, and Y. Fouquart, 2000: Radiative properties of boundary layer clouds: Droplet effective radius versus number concentration. J. Atmos. Sci., 57, 803–821. Caughey, S. J., and C. K. Kitchen, 1984: Simultaneous measurements of the turbulent and microphysical structure of nocturnal stra- tocumulus cloud. Quart. J. Roy. Meteor. Soc., 110, 13–34. Chambers, L., B. Wielicki, and K. Evans, 1997: Accuracy of the independent pixel approximation for satellite estimates of oce- anic boundary layer cloud optical depth. J. Geophys. Res., 102, 1779–1794. Coakley, J. A., and F. P. Bretherton, 1982: Cloud cover from high- resolution scanner data: Detecting and allowing for partially filled fields of view. J. Geophys. Res., 87, 4917–4932. Greenwald, T. J., and S. A. Christopher, 1999: Daytime variation of marine stratocumulus microphysical properties as observed from geostationary satellite. Geophys. Res. Lett., 26, 1723–1726. ——, ——, J. Chou, and J. C. Liljegren, 1999: Intercomparison of cloud liquid water path derived from the GOES 9 imager and ground based microwave radiometers for continental stratocu- mulus. J. Geophys. Res., 104, 9251–9260. Han, Q., W. Rossow, and A. Lacis, 1994: Near-global survey of effective droplet radii in liquid water clouds using ISCCP data. J. Climate, 7, 465–497. ——, ——, R. Welch, A. White, and J. Chou, 1995: Validation of satellite retrievals of cloud microphysics and liquid water path using observations from FIRE. J. Atmos. Sci., 52, 4183–4195. ——, ——, J. Chou, and R. M. Welch, 1998a: Global survey of the relationships of cloud albedo and liquid water path with droplet size using ISCCP. J. Climate, 11, 1516–1528. ——, ——, ——, and ——, 1998b: Global variation of column drop- let concentration in low-level clouds. Geophys. Res. Lett., 25, 1419–1422. Hartmann, D., M. Ockert-Bell, and M. Michelsen, 1992: The effect of cloud type on earth’s energy balance: Global analysis. J. Cli- mate, 5, 1281–1304. Isaaks, E., and R. Srivastava, 1989: An Introduction to Applied Geos- tatistics. Oxford University Press, 561 pp. Jansson, M., 1985: A comparison of detransformed logarithmic re- gressions and power function regressions. Geogr. Ann., 67A, 61–70. Kaufman, Y. J., and B. N. Holben, 1993: Calibration of the AVHRR visible and near-IR bands by atmospheric scattering, ocean glint and desert. Int. J. Remote Sens., 14, 21–52. Khairoutdinov, M. F., and Y. L. Kogan, 1999: A large eddy simulation model with explicit microphysics: Validation against aircraft ob- servations of a stratocumulus-topped bounary layer. J. Atmos. Sci., 56, 2115–2131. Kneizys, F., E. Shettle, L. Abreu, J. C. Chetwynd, Jr., S. Anderson, W. Gallery, J. Selby, and S. Clough, 1988: Users guide to LOW- TRAN-7. Air Force Geophysics Laboratory Tech. Rep. AFGL- TR-88-0177, 137 pp. Lohmann, U., and J. Feichter, 1997: The atmospheric sulfur cycle in ECHAM-4 and its impact on the shortwave radiation. J. Geo- phys. Res., 102, 13 685–13 700. Martin, G. M., D. W. Johnson, and A. Spice, 1994: The measurement and parametrisation of effective radius in warm stratocumulus clouds. J. Atmos. Sci., 51, 1823–1842. Nakajima, Ta., and Te. Nakajima, 1995: Wide-area determination of cloud microphysical properties from NOAA AVHRR measure- ments for FIRE and ASTEX regions. J. Atmos. Sci., 52, 4043– 4059. Nakajima, Te., and M. King, 1990: Determination of the optical thick- ness and effective particle radius of clouds from reflected ra- diation measurements. Part I: Theory. J. Atmos. Sci., 47, 1878– 1892. ——, ——, J. Spinhirne, and L. Radke, 1991: Determination of the optical thickness and effective particle radius of clouds from reflected radiation measurements. Part II: Marine stratocumulus observations. J. Atmos. Sci., 48, 728–750. Nicholls, S., 1984: The dynamics of stratocumulus: Aircraft obser- vations and comparisons with mixed layer models. Quart. J. Roy. Meteor. Soc., 110, 783–820. Pincus, R., M. Szczodrak, P. Austin, and J. Gu, 1995: Precision and accuracy in satellite optical depth estimates. J. Climate, 8, 1453– 1462. Platnick, S., and F. Valero, 1995: A validation of satellite cloud re- trieval during ASTEX. J. Atmos. Sci., 52, 2985–3001. Pontikis, C., 1993: Parameterization of the cloud optical thickness: Influence of clear air entrainment. Geophys. Res. Lett., 20, 2655– 2658. Szczodrak, M., 1998: Variability of cloud optical depth and cloud droplet effective radius in layer clouds: Satellite based analysis. Ph.D. thesis, University of British Columbia, Vancouver, BC, Canada, 213 pp. Wang, S., and Q. Wang, 1994: Roles of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus- topped marine boundary layer. J. Atmos. Sci., 51, 1559–1576. Zuidema, P., and D. Hartmann, 1995: Satellite determination of stratus cloud microphysical properties. J. Climate, 8, 1638–1657.


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items