1. Introduction
The question of how atmospheric turbulence impacts cloud microphysical properties and cloud evolution has benefited from renewed interest in the past several years with many recent advances summarized by Shaw (2003) and Vaillancourt and Yau (2000). Long-lived vortex tubes on spatial scales close to the Kolmogorov microscale (which typically is ∼1 mm or so) have been hypothesized to cause preferential concentration of droplets and supersaturation intermittency (e.g., Shaw et al. 1998; Shaw 2000), thereby promoting the growths of “favored” droplets and the onset of active coalescence. The ability of turbulence to enhance droplet collision rates has also been demonstrated to be effective at relatively low Reynolds numbers (Vohl et al. 1999; Shaw 2003; Pinsky and Khain 2004), thus raising the question of the efficiency of these processes in real clouds. Recent studies have also begun to look at the effect of the cloud microphysical structure on the statistics of the turbulent flow field. The modification of the turbulent energy spectrum by solid particles with diameters smaller than the Kolmogorov microscale was investigated with direct numerical simulation (DNS) methods by Ferrante and Elgobashi (2003) and Elgobashi and Truesdell (1993). They found that solid particles increase the fluid turbulence energy at high wavenumbers. The influence of water droplets in a multiphase flow proves to be a more challenging problem, in large part because the evaporation or condensation of vapor from or onto droplets strongly affects the local buoyancy. However, because of the small volume and mass fraction the impact of cloud droplets on the flow is supposed to be negligible (Vaillancourt and Yau 2000).
One of the first airborne measurements of vertical wind velocities in clouds was presented in Telford and Warner (1962) indicating a significant increase of turbulence within clouds. The first power spectral densities S( f ) of the liquid water content (LWC) were shown in Ackerman (1967) with a resolution of better than 100 m. MacPherson and Isaac (1977) presented the first airborne measurements of S( f ) of all three wind velocity components in cumulus clouds and estimated a mean energy dissipation rate of
Today, single droplet measurements on the millimeter scale using the Forward Scattering Spectrometer Probe (FSSP), Fast-FSSP (Brenguier 1993; Brenguier et al. 1998), and the modified Fast-FSSP (M-Fast FSSP; Schmidt et al. 2004) are available. Airborne data of the LWC with the Particle Volume Monitor (PVM-100A; Gerber et al. 2001) and of static air temperature in clouds using the ultrafast thermometers (UFT; Haman et al. 2001) are also on-hand and deliver data with centimeter resolution. Vaillancourt and Yau (2000) pointed out that turbulence measurements with much higher resolution (than several meters) would be needed to investigate the intermittent structure of clouds. Furthermore, all the discussions concerning droplet–turbulence interaction are based on the assumption that turbulence in clouds can be described with the same methods as turbulence in a droplet-free flow. However, the validity of this hypothesis has only been established for spatial scales on the order of several meters since in situ measurements of the turbulent velocity field in clouds with higher resolution were simply not available so far (see discussion in Vaillancourt and Yau 2000).
To overcome the limited spatial resolution (due to typically high TAS of research aircraft) slower-flying helicopters (e.g., Muschinski et al. 2001) have been used yielding spatial resolutions of the measurements below one meter. Tethered balloons have also been used to further increase the spatial resolution. Such balloon measurements were reported by Muschinski et al. (2004), who estimated local ε values and temperature structure parameters. These measurements were based on fine-wire sensors, which are difficult to use under cloudy conditions. Other balloon-borne measurements were reported by Gerber (1986), Blanc et al. (1989), and Plant et al. (1998). These measurements did not include high-resolution turbulence sensors. Tethered balloon–borne turbulence measurements in small cumulus clouds were reported by Kitchen and Caughey (1981), where power spectral densities of the vertical wind at frequencies up to 5 Hz were reported.
To investigate the spatial structure of turbulence in clouds on smaller scales, Siebert et al. (2003) developed the tethered balloon–borne payload Airborne Cloud Turbulence Observation System (ACTOS). The spatial resolution that can be achieved with the measurements of ACTOS (wind vector components, temperature, and LWC) is in the order of 10 cm, which is the highest possible resolution so far.
The main objective of this work is to test the hypothesis that basic ideas of classical turbulence theory (Kolmogorov 1941, 1962; Oboukhov 1962; Batchelor 1953) in boundary layer clouds are valid even to the decimeter scale. For this purpose, two datasets with measurements in different types of boundary layer clouds are analyzed. For temperature measurements in clouds this hypothesis has already been discussed (e.g., Haman et al. 2001), therefore, here we focus on high-resolution wind velocity measurements.
After a brief overview of the instrumental setup (section 2), the two datasets used in this study are introduced (section 3). In section 4 the structure of the turbulent flow, down to centimeter scales, is explored in more detail, with a particular emphasis on local (8 m) estimates of ε and the associated pdfs of the wind.
2. Instrumentation
High-frequency measurements of the three-dimensional wind vector in boundary layer clouds were made with ACTOS carried by the tethered balloon Mobiles Autonomes Positionierungs System (MAPS-Y; Siebert et al. 2003). ACTOS is an instrumental payload that, in addition to measuring the three-dimensional wind vector, also measures static temperature, and humidity each with a sampling frequency of at least 100 Hz. A PVM-100A (Gerber 1991; Gerber et al. 1994) measures the LWC; cloud droplet size distributions are derived from measurements with the M-Fast-FSSP. A photo of the current version of ACTOS is shown in Fig. 1. ACTOS consists of four shielded 19-in racks that include power supply, real-time data acquisition, and further sensor electronics. The turbulence sensors are fixed to a 1.5-m-long outrigger made of lightweight carbon fiber. Two tail units keep ACTOS aligned with the mean wind. A navigation unit including a differential global positioning system (GPS), an inertial navigation system (INS), and a high-resolution barometer allows reconstruction of the ACTOS attitude, position, and motion. The total mass of ACTOS is 120 kg.
This work is mainly based on data from the ultrasonic anemometer (hereafter called sonic) and the navigation unit. The LWC and the static air temperature measured with a shielded fine-wire resistance thermometer (UFT-B; Haman et al. 2000) are used to give an overview of the general conditions only.
a. Sonic measurement
The sonic used in this study is a three-dimensional Solent-Research Horizontally Symmetric (HS), manufactured by Gill Instruments, Lymington, United Kingdom. The maximum sampling frequency is 100 Hz; that is, the Nyquist frequency fNy is 50 Hz. From laboratory investigations (Siebert and Muschinski 2001) the spectral noise floor S(n)ui for ui (i = 1, 2, and 3 for the three wind velocity components) is estimated to about 10−5 m2 s−2 Hz−1, the standard deviation due to uncorrelated noise is σ(n)ui =
b. Corrections for payload attitude and motion
The measurements of the wind vector with the sonic refer to a payload-fixed coordinate system and, therefore, have to be corrected for attitude and motion of ACTOS. The correction has been performed following Lenschow (1986) and Edson et al. (1998).
In the standard configuration of ACTOS, the roll and pitch angles are directly derived from the INS data. The GPS system (TANS Vector, Trimble Ltd., Sunnyvale, California) is a multiantenna array that measures all three attitude angles using differential carrier phase measurements between the four antennas. Only the heading is used to support the INS in the low frequency range. The resolution of roll and pitch from the INS is about 0.0006° and 0.2° for the heading (including stabilization with GPS heading). The measure of the translational velocity vector of ACTOS is also a combination of INS and GPS data.
The three attitude angles can also be derived directly from the GPS, though with lower resolution (0.3° for roll and pitch and 0.2° for the heading angles). In this case, the wind measurements can only be corrected for attitude since the accuracy and data availability of the translational velocity components from the GPS is incomplete. Therefore, calculated wind velocity data show a peak around the pendulum frequency at 0.1 Hz, most significant for the lateral wind component. Only high-pass-filtered data with a cutoff frequency >0.1 Hz can then be used for further analysis.
c. ACTOS and MAPS-Y
ACTOS is lifted with the tethered balloon MAPS-Y up to a height of around 1.5 km depending on wind conditions. The maximum wind speed for the system is 15 m s−1. MAPS-Y has a volume of 400 m3 and can be filled either with hydrogen or helium. The maximum climbing rate is around 1.5 m s−1 for moderate wind speeds and up to 3 m s−1 for wind speeds below 5 m s−1. The payload was fixed 25 m below MAPS-Y to overcome the influence of flow distortions induced by the balloon itself (cf. Wills and Cole 1986). More details and a photo of MAPS-Y can be found in Siebert et al. (2003).
3. Data overview
In this work data from two field campaigns are considered. The first experiment was the Baltex Bridge Campaign (BBC2; Crewell et al. 2004), which was conducted near Cabauw, Netherlands, in May 2003. The second was called Influence of Clouds on Spectral Actinic Flux in the Lower Troposphere (INSPECTRO2) and took place in southern Germany in May 2004. From both campaigns, a time series recorded during constant level flights is analyzed. Measurements of the attitude-corrected longitudinal wind velocity U, vertical velocity W, static air temperature Ta, and LWC are presented. Actual values are indicated by capitals, whereas the fluctuating part is described by lower case letters: X(t) =
a. BBC2 campaign
A 20-min-long record taken on 21 May 2003 between 0947 and 1006 UTC was analyzed in detail. More than 50% of the flight time was in clouds. The flight height was 760 ± 20 m. Shallow cumulus clouds (Cu med) were observed by a horizontal-looking video camera onboard ACTOS. The cloud base was around 500 m as estimated from a ground-based ceilometer (CT 75, Vaisala, Finland). Figure 2 shows two vertical profiles of static air temperature Ta, relative humidity RH, horizontal wind speed U, and wind direction D as derived from radiosonde data (soundings from 0830 and 1121 UTC). The temperature at cloud base was 8°C as derived from Fig. 2. A small temperature inversion with a significant decrease of relative humidity was observed in the first profile around 1200-m altitude in association with cloud top. However, the picture given by the two soundings is not completely consistent. In particular, the relative humidity profiles differ significantly, which is most likely caused by the local cloud situation for the individual sounding. In the measurement height of ACTOS (marked with a horizontal line in Fig. 2), no significant temperature inversion or change in the wind direction happened in between the soundings. It was not possible to measure meaningful vertical profiles in clouds with ACTOS because of the limited horizontal extension of the clouds.
Figure 3 shows the time series of U, W, LWC, and Ta measured with ACTOS. During this flight, the INS was not running because of a technical problem; therefore, the wind velocity components were corrected for roll and pitch angle derived from GPS. In Fig. 3 both wind components are low-pass filtered (second-order Butterworth, fcut = 0.05 Hz) to reduce the oscillation due to the pendulum motion. The mean horizontal wind was
b. INSPECTRO2 campaign
From the INSPECTRO2 campaign a time series from the 16 May 2004 is analyzed (cf. Fig. 4). This record is 33 min long and was also taken at a nearly constant altitude of 1540 ± 40 m. The sky was covered with a layer of scattered shallow fair weather cumuli (Cu humilis), which run flat and disappeared shortly after these measurements. Using the onboard camera, the position of ACTOS was estimated to be in the lower third of the clouds.
During this period the navigation unit worked properly and, therefore, the wind vector was completely corrected for attitude and translational motion as described in section 2b. The mean horizontal wind was 7.4 m s−1 (σU = 1.4 m s−1). The vertical wind velocity showed strong fluctuations with positive and negative values outside the cloud, whereas inside the cloud W was almost always positive. A few small cloud holes corresponded well with downdrafts. After the cloud section W fluctuated around zero with an amplitude of about 1 m s−1. The temperature was near zero with lowest values of around −1°C.
The structure of the LWC field was very inhomogeneous even inside the cloud. A 100-s-long subsequence of the LWC time series is enlarged at the top panel of Fig. 4. Part A of this subsequence shows a more homogeneous LWC field, whereas part B exhibits more fluctuations of the LWC. These clouds are much less dynamic compared with the BBC2 clouds and are in the stage of decaying.
4. Data analysis
a. Spectral analysis and mean energy dissipation rates 







Figure 5 shows S( f ) for all parameters presented in Figs. 3 and 4 (LWC, W, U, Ta). All S( f ) of the BBC2 data are divided by a factor of 10 to clearly distinguish from INSPECTRO2 data. For frequencies higher 0.1 Hz all spectra closely follow a −5/3 slope. The spectra of u and w for the BBC2 data show a slight peak around 0.1 Hz due to the pendulum motion that could not completely be removed. The S( f ) of both velocity components of the INSPECTRO2 data show a slight flattening for frequencies higher 30 Hz. This behavior was observed with the same type of sonic during a previous ground-based experiment (Siebert and Muschinski 2001).




b. Local isotropy

c. Local energy dissipation rates
The small-scale structure of the turbulence is investigated by means of so-called local energy dissipation rates ετ where the index τ indicates the averaging time over which the parameter is estimated. In this context, based on the work by Kolmogorov (1962) and Oboukhov (1962), ετ is used as a random variable compared to
- (i) the inertial dissipation method based on Eq. (3), which was applied to the complete record to estimate
can be used with shorter averaging times τ to estimate ετ; - (ii) the direct method based on the definition of ε (local isotropy and the validity of Taylor's hypothesis are assumed):Here, 〈x〉τ indicates an average of x over the time τ. The direct method requires measurements with a spatial resolution in the order of the Kolmogorov length scale of ∼1 mm (Oncley et al. 1996; Champagne et al. 1977), which is not possible with sonics.
- (iii) The relation between the variance and S( f ) [cf. Eq. (4)] can be solved for the case that the complete frequency range is within the inertial subrange:in the frequency interval [ f0; fNy] and f0 ≪ fNy:with f0 = 1/T0 = 1/(2T); and, finally
- (iv) ετ can be derived from the second-order structure parameter (C2u)τ = 2ε2/3τ of u:and, therefore,with the second-order structure function Dτu(t′) = 〈(u(t + t′) − u(t))2〉τ and the time lag t′.
In Fig. 8 the time series of the local ετ is presented for BBC2 data. Additionally, a running average over 10 points is included. The values of local ετ cover a range of nearly four orders of magnitude.
For the BBC2 data, minimum values of ετ down to 10−5 m2 s−3 were found in the cloud-free areas, absolute minima are often found just outside the clouds (cf. marked regions i and ii in Fig. 8). The mean value is

For the INSPECTRO2 data (Fig. 9) the correlation between local ετ and LWC was not as obvious. The smallest values were found inside homogeneous parts of clouds, such as part A. Higher values were found in regions of high fluctuations of LWC, especially at cloud edges (see, e.g., dashed box), where lateral entrainment of dry air created stronger fluctuations due to evaporative cooling.
d. PDFs of wind velocity increments and ετ
To analyze the statistics of wind velocity fluctuations in an intermittent flow it is common to calculate the PDFs of the wind velocity increments Δui instead of ui. The increment of a variable x(t) is defined as Δx ≡ x(t + τ) − x(t). Here, the time lag τ = 10 ms is used, which is the sampling time of our measurements. In Fig. 10, the PDFs of Δu and Δw for the BBC2 data are presented. To compare these PDFs with model distributions, a Gaussian fit, which is often used as a first-order assumption, is included. All data show a smooth and symmetric PDF, only some scatter can be observed for high absolute values of Δu and Δw. The difference between the Gaussian assumption and the measured PDF is obvious. The shape of the measured PDF is typical for intermittent time series where high values are more frequent than predicted by a Gaussian PDF (Böttcher et al. 2003; Frisch 1995). The probability of the highest values of the measured velocity increments is at least one order of magnitude higher than predicted by a Gaussian distribution that emphasizes the intermittent character of the wind velocity field.
A widely used model assumption for the PDFs of ετ in high Reynolds number flows is a lognormal distribution (e.g., Shaw and Oncley 2001). A discussion of different models for PDFs of ετ can be found in the review article by Sreenivasan and Antonia (1997). Figure 11 shows the PDFs of the natural logarithm of ετ for the BBC2 (top) and INSPECTRO2 (bottom) datasets. Both PDFs are quantitatively well described by a Gaussian distribution; attempts to more quantitatively describe the distributions were hindered by the limited sample.
In Fig. 12 the PDFs of ετ conditionally sampled on LWC are presented for BBC2 (top) and INSPECTRO2 (bottom) data. In this presentation the differences between the two types of clouds become more evident. For the BBC2 data the PDFs of the cloud and the cloud-free data are clearly shifted to each other. The mean ετ inside of clouds is 1 × 10−2 m2 s−3 in contrast to 2 × 10−3 m2 s−3 in the cloud-free regions. A significant increase of the PDF(ετ) for cloud-free data coincides with the mean of PDF(ετ) for cloud data. This indicates that there are regions associated with updrafts and increased turbulence but not necessarily resulting in clouds. For the INSEPCTRO2 data, this difference is much less obvious (5 × 10−3 m2 s−3 in clouds and 8 × 10−3 m2 s−3 for the cloud-free regions).
5. Discussion
The observed wide range of ετ has consequences for the interaction between turbulence and microphysical processes: “The utilization of an averaged value of the dissipation rate can lead to an underestimation of the effects of turbulence on collisions and on precipitation formation. It is quite possible that regions of effective collisions should be well correlated with zones of enhanced turbulence” (Pinsky and Khain 2004).
The key question is whether there are turbulence-induced local accelerations of fluid parcels resulting in different velocities of droplets as a function of their inertia. Pinsky and Khain (2004) pointed out that “the effect of Lagrangian acceleration is mainly caused by the increase in the collision efficiency that is highly sensitive even to weak variations of interdrop relative velocity.”


For this discussion, several points have to be considered. First, the observed clouds are only fair weather clouds and much higher mean and local values of ε can be expected for more convective cumulus clouds. Second, the local ετ are averaged over a period τ, which is at least two orders of magnitude higher than the Kolmogorov time scale τK = (ν/
The observation of high Taylor Reynolds number resulting in high intermittency of ετ shows also the general limitation for direct numerical simulations (DNS) for investigations of turbulence–particle interactions. The largest DNS of turbulent flows reached Reλ ∼103 (Kaneda et al. 2003) whereas DNS studies of particle-laden turbulence are limited to Reλ ∼ 102 (e.g., Ferrante and Elgobashi 2003; Squires and Eaton 1991) which is about two orders of magnitude smaller than observed in our data. Therefore, results concerning particle–turbulence interaction obtained from DNS have to be interpreted with caution.
6. Summary
From earlier airborne studies of turbulence in clouds only data with a resolution of several meters are available. Therefore, only mean energy dissipation rates averaged over a complete flight lag were reported. In this work, turbulence data with a resolution on decimeter scales are presented for two measurement campaigns (BBC2 and INSPECTRO2). The two datasets represents two types of small cumulus fields at different stages of development. The aim of this work was to test the hypotheses whether turbulence in clouds can be described with the same methods as used for intermittent (but cloud-free) flows.
Power spectra show an inertial subrange behavior down to length scales of about 0.2 m. For length scales between 0.8 and 80 m the ratios of S( f ) of longitudinal and transversal wind velocity components are close to the theoretical value of 4/3 indicating local isotropy. Deviations for smaller scales are due to general limitations of the measurements.
For the BBC2 data, high values of ετ were found in clouds in association with high turbulent activity due to vertical updrafts. For the observations during INSPECTRO2, highest ετ values were found in cloud gaps and mostly small values were observed in homogeneous parts inside the clouds. For this case, the vertical wind speed showed no significant correlation with the LWC indicating an “old” advected cloud. The higher ετ values between the clouds or in regions with cloud holes may be a signature of lateral mixing at cloud edges.
These differences between the two datasets become most evident in PDFs derived from data conditionally sampled on the LWC. For the BBC2 data, the average of ετ measured in clouds is about 5 times higher than for cloud-free regions, whereas there is no significant difference for the INSPECTRO2 data. However,
Probability density functions of wind velocity increments are derived. The PDFs show significant deviations from a Gaussian distribution consistent with strong intermittency. The shape of the ετ PDFs could be parameterized with a lognormal distribution that is in good agreement with theory (Kolmogorov 1962) and earlier experiments made in the cloud-free boundary layer (Shaw and Oncley 2001). From these data we conclude that small-scale turbulence in clouds at least to decimeter scales can be well described with the same methods as used for cloud-free conditions.
We thank G. Sanftleben and his crew from the German Bundeswehr Technical Centre for Ships and Naval Weapons for the balloon operation. D. Schell from the enviscope GmbH (Frankfurt/Main) is acknowledged for his technical support during both experiments. Thanks are also due to R. Shaw, O. Hellmuth, and the two anonymous reviewers for their helpful comments on an earlier version of the manuscript. Finally, we thank KNMI, Netherlands, for hosting the BBC2 campaign and providing the radiosonde data.
REFERENCES
Ackerman, B., 1967: The nature of the meteorological fluctuations in clouds. J. Appl. Meteor, 6 , 61–71.
Batchelor, G. K., 1953: The Theory of Homogenous Turbulence. Cambridge University Press, 197 pp.
Biltoft, C. A., 2001: Some thoughts on local isotropy and the 4/3 lateral to longitudinal velocity spectrum ratio. Bound.-Layer Meteor, 100 , 393–404.
Blanc, T. V., , W. J. Plant, , and W. C. Keller, 1989: The Naval Research Laboratory's air–sea interaction blimp experiment. Bull. Amer. Meteor. Soc, 70 , 354–364.
Böttcher, F., , C. Renner, , H-P. Waldl, , and J. Peinke, 2003: On the statistics of wind gusts. Bound.-Layer Meteor, 108 , 163–173.
Brenguier, J-L., 1993: Observations of cloud microstructure at the centimeter scale. J. Appl. Meteor, 32 , 783–793.
Brenguier, J-L., , T. Bourrianne, , A. A. Coelho, , J. Isbert, , R. Peytavi, , D. Trevarin, , and P. Weschler, 1998: Improvements of droplet size distribution measurements with the fast-FSSP (forward scattering spectrometer probe). J. Atmos. Oceanic Technol, 15 , 1077–1090.
Champagne, F. H., , C. A. Friehe, , J. C. LaRue, , and J. Wyngaard, 1977: Flux measurements, flux estimation techniques, and fine-scale turbulence measurements in the unstable surface layer over land. J. Atmos. Sci, 34 , 515–530.
Crewell, S., and Coauthors, 2004: The BALTEX Bridge Campaign: An integrated approach for a better understanding of clouds. Bull. Amer. Meteor. Soc, 85 , 1565–1584.
Edson, J. B., , A. A. Hinton, , K. E. Prada, , J. E. Hare, , and C. W. Fairall, 1998: Direct covariance flux estimates from mobile platforms at sea. J. Atmos. Oceanic Technol, 15 , 547–562.
Elgobashi, S., , and C. Truesdell, 1993: On the two-way interaction between homogeneous turbulence and dispersed solid particles. I: Turbulence modification. Phys. Fluids, 5 , 1790–1801.
Ferrante, A., , and S. Elgobashi, 2003: On the physical mechanisms of two-way coupling in particle laden isotropic turbulence. Phys. Fluids, 15 , 315–329.
Frehlich, R., , Y. Meillier, , M. L. Jensen, , and B. Balsley, 2004: A statistical description of small-scale turbulence in the low-level nocturnal jet. J. Atmos. Sci, 61 , 1079–1085.
Frisch, U., 1995: Turbulence—The Legacy of A. N. Kolmogorov. Cambridge University Press, 296 pp.
Gerber, H., 1986: Tethered balloon measurements at San Nicolas Island (Oct. 1984): Instrumentation, data summary, preliminary data interpretation. NRL Tech. Rep. 8972, Naval Research Laboratory, Atmospheric Physics Branch, Space Science Division, 29 pp.
Gerber, H., 1991: Direct measurement of suspended particulate volume concentration and far-infrared extinction coefficient with a laser-diffraction instrument. Appl. Opt, 30 , 4824–4831.
Gerber, H., , B. G. Arends, , and A. S. Ackerman, 1994: New microphysics sensor for aircraft use. Atmos. Res, 31 , 235–252.
Gerber, H., , J. B. Jensen, , A. B. Davis, , A. Mashak, , and W. J. Wiscombe, 2001: Spectral density of cloud liquid water content at high frequencies. J. Atmos. Sci, 58 , 497–503.
Haman, K., , S. P. Malinowski, , A. Makulski, , B. D. Struś, , R. Busen, , A. Stefko, , and H. Siebert, 2000: A family of ultrafast aircraft thermometers for warm and supercooled clouds and various types of aircraft. Preprints, 13th Int. Conf. on Clouds and Precipitation, Reno, NV, International Commission on Clouds and Precipitation, 224–227.
Haman, K. E., , S. P. Malinowski, , B. D. Struś, , R. Busen, , and A. Stefko, 2001: Two new types of ultrafast aircraft thermometer. J. Atmos. Oceanic Technol, 18 , 117–134.
Hill, R. J., 2002: Scaling of acceleration in locally isotropic turbulence. J. Fluid Mech, 452 , 361–370.
Kaimal, J. C., , and J. J. Finnigan, 1994: Atmospheric Boundary Layer Flows. Oxford University Press, 289 pp.
Kaimal, J. C., , J. C. Wyngaard, , and D. A. Haugen, 1968: Deriving power spectra from a three-component sonic anemometer. J. Appl. Meteor, 7 , 827–837.
Kaimal, J. C., , J. C. Wyngaard, , Y. Izumi, , and O. R. Cote, 1972: Spectral characteristics of surface-layer turbulence. Quart. J. Roy. Meteor. Soc, 98 , 563–589.
Kaneda, Y., , T. Ishihara, , M. Yokokawa, , K. Itakura, , and A. Uno, 2003: Energy dissipation rate and energy spectrum in high resolution direct numerical simulations of turbulence in a periodic box. Phys. Fluids, 15 .L21–L24, doi:10.1063/1.1539855.
Kitchen, M., , and S. J. Caughey, 1981: Tethered-balloon observations of the structure of small cumulus clouds. Quart. J. Roy. Meteor. Soc, 107 , 853–874.
Kolmogorov, A. N., 1941: The local structure of turbulence in incompressible viscous fluid for very large Reynolds numbers. Dokl. Akad. Nauk SSSR, 30 , 301–304.
Kolmogorov, A. N., 1962: A refinement of previous hypotheses concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number. J. Fluid Mech, 13 , 82–85.
Lenschow, D. H., 1986: Aircraft measurements in the boundary layer. Probing the Atmospheric Boundary Layer, D. H. Lenschow, Ed., Amer. Meteor. Soc., 239 pp.
Lenschow, D. H., , M. Zhou, , X. Zeng, , L. Chen, , and X. Xu, 2000: Measurements of fine-scale structures at the top of marine stratocumulus. Bound.-Layer Meteor, 97 , 331–357.
MacPherson, J. I., , and G. A. Isaac, 1977: Turbulent characteristics of some Canadian cumulus clouds. J. Appl. Meteor, 16 , 81–90.
Muschinski, A., , R. G. Frehlich, , M. L. Jensen, , R. Hugo, , A. M. Hoff, , F. Eaton, , and B. B. Balsley, 2001: Fine-scale measurements of turbulence in the lower troposphere: An intercomparison between a kite- and balloon-borne and a helicopterborne measurement system. Bound.-Layer Meteor, 98 , 219–250.
Muschinski, A., , R. G. Frehlich, , and B. B. Balsley, 2004: Small-scale and large-scale intermittency in the nocturnal boundary layer and the residual layer. J. Fluid Mech, 515 , 319–351.
Nicholls, S., 1984: The dynamics of stratocumulus: Aircraft observations and comparisons with a mixed layer model. Quart. J. Roy. Meteor. Soc, 110 , 783–820.
Nicholls, S., , and J. Leighton, 1986: An observational study of the structure of stratiform cloud sheets: Part I. Structure. Quart. J. Roy. Meteor. Soc, 112 , 431–460.
Oboukhov, A. M., 1962: Some specific features of atmospheric turbulence. J. Fluid Mech, 13 , 77–81.
Oncley, S., , C. A. Friehe, , J. C. Larue, , J. A. Businger, , E. C. Itsweire, , and S. S. Chang, 1996: Surface-layer fluxes, profiles, and turbulence measurements over uniform terrain under near-neutral conditions. J. Atmos. Sci, 53 , 1029–1043.
Pinsky, M. B., , and A. P. Khain, 2004: Collision of small drops in a turbulent flow. Part II: Effects of flow accelerations. J. Atmos. Sci, 61 , 1926–1939.
Plant, W. J., , W. C. Keller, , V. Hesany, , K. Hayes, , K. W. Hoppel, , and T. V. Blanc, 1998: Measurements of the marine boundary layer from an airship. J. Atmos. Oceanic Technol, 15 , 1433–1458.
Pope, S. B., 2000: Turbulent Flows. Cambridge University Press, 771 pp.
Saddoughi, S. G., , and S. V. Veeravalli, 1994: Local isotropy in turbulent boundary layers at high Reynolds numbers. J. Fluid Mech, 268 , 333–372.
Schmidt, S., , K. Lehmann, , and M. Wendisch, 2004: Minimizing instrumental broadening of the drop size distribution with the M-Fast-FSSP. J. Atmos. Oceanic Technol, 21 , 1855–1867.
Shaw, R. A., 2000: Supersaturation intermittency in turbulent clouds. J. Atmos. Sci, 57 , 3452–3456.
Shaw, R. A., 2003: Particle-turbulence interactions in atmospheric clouds. Annu. Rev. Fluid Mech, 35 , 183–227.
Shaw, R. A., , and S. P. Oncley, 2001: Acceleration intermittency and enhanced collision kernels in turbulent clouds. Atmos. Res, 59–60 , 77–87.
Shaw, R. A., , W. C. Reade, , L. R. Collins, , and J. Verlinde, 1998: Preferential concentration of cloud droplets by turbulence: Effects on the early evolution of cumulus cloud droplet spectra. J. Atmos. Sci, 55 , 1965–1976.
Siebert, H., , and U. Teichmann, 2000: The behaviour of an ultrasonic under cloudy conditions. Bound.-Layer Meteor, 94 , 165–169.
Siebert, H., , and A. Muschinski, 2001: Relevance of a tuning-fork effect for temperature measurements with the Gill Solent HS ultrasonic anemometer-thermometer. J. Atmos. Oceanic Technol, 18 , 1367–1376.
Siebert, H., , M. Wendisch, , T. Conrath, , U. Teichmann, , and J. Heintzenberg, 2003: A new tethered balloon-borne payload for fine-scale observations in the cloudy boundary layer. Bound.-Layer Meteor, 106 , 461–482.
Smith, S. A., , and P. P. Jonas, 1995: Observations of the turbulent fluxes in fields of cumulus clouds. Quart. J. Roy. Meteor. Soc, 121 , 1185–1208.
Squires, K. D., , and J. K. Eaton, 1991: Preferential concentration of particles by turbulence. Phys. Fluids, 3 , 1169–1178.
Sreenivasan, K. R., , and R. A. Antonia, 1997: The phenomenology of small-scale turbulence. Annu. Rev. Fluid Mech, 29 , 435–472.
Telford, J. W., , and J. Warner, 1962: On the measurement from aircraft of buoyancy and vertical air velocity in cloud. J. Atmos. Sci, 19 , 415–423.
Vaillancourt, P. A., , and M. K. Yau, 2000: Review of particle–turbulence interactions and consequences for cloud physics. Bull. Amer. Meteor. Soc, 81 , 285–298.
Vohl, O., , S. K. Mitra, , S. Wurzler, , and H. R. Pruppacher, 1999: A wind tunnel study of the effect of turbulence on the growth of cloud drops by collision and coalescence. J. Atmos. Sci, 56 , 4088–4099.
Willis, G. E., , and J. W. Deardorff, 1976: On the use of Taylor's translation hypothesis for diffusion in the mixed layer. Quart. J. Roy. Meteor. Soc, 102 , 817–822.
Wills, J. A. B., , and L. R. Cole, 1986: Model studies of wind flow distortion on an AI600 airship. Tech. Rep. 45046, British Maritime Technology, Middlesex, United Kingdom, 15 pp.

The turbulence payload ACTOS with sonic, UFT-B, PVM-100A, and M-Fast-FSSP. Also shown are the inlets for humidity and aerosol particle measurements.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

Two vertical profiles of (left) air temperature (Ta) and RH and (right) horizontal wind speed (U) and wind direction (D). The profiles were derived from radiosondes that were launched close to the balloon site. The first radiosonde was launched at 0830 UTC, the second at 1121 UTC. A horizontal line marks the height of the balloon between 0947 and 1006 UTC.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

Time series of LWC, vertical wind velocity W, horizontal wind velocity U, and temperature Ta as measured with ACTOS at a height of around 760 m AGL on 21 May 2003 during the BBC2 campaign.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

Time series of LWC, vertical wind velocity W, horizontal wind velocity U, and temperature Ta as measured with ACTOS at a height of around 1540 m AGL on 16 May 2004 during the INSPECTRO campaign.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

Power spectral densities S( f ) of the same data as presented in Figs. 3 and 4. All spectra are in units of their variance per frequency; spectra of BBC data are divided by a factor of 10 for better resolution. For the top panel the frequencies are converted into wavelength assuming a constant horizontal wind speed of 8 m s−1.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

Spectral ratios Sυ/Su and Sw/Su for (top) BBC2 data and (bottom) INSPECTRO2 data. A 4/3 ratio is plotted for reference, indicating the value predicted by classical theory for isotropic turbulence.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

(top) Correlation of εSu vs εσ3 and εC2u, respectively. Each point represents an ετ averaged over a 1-s period (τ = 1 s). (bottom) A short time series of ετ calculated with the three methods described in the text.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

(top) Time series of local energy dissipation rate ετ and (bottom) LWC of BBC2 data. The integration time τ for ετ is 1 s; a running average over 10 points is included.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

(top) Time series of local energy dissipation rate ετ and (bottom) LWC of INSPECTRO2data. The integration time τ for ετ is 1 s; a running average over 10 points is included.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

Semilogarithmic plot of the PDFs of the increments Δu and Δw for the BBC2 data. A Gauss fit is included for reference.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

PDF of natural logarithm of local energy dissipation rates ετ. A Gauss fit is included for reference.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1

PDF of natural logarithm of energy dissipation rates of the (top) BBC2 data and (bottom) INSPECTRO2 data inside of clouds and outside of clouds. The energy dissipation rates are conditionally sampled on the LWC.
Citation: Journal of the Atmospheric Sciences 63, 5; 10.1175/JAS3687.1