Cloud-clear Air Research Articles

Microphysical timescales and local supersaturation balance at a warm cloud top boundary

Recent results have shown that there is an acceleration in the spread of the size distribution of droplet populations in the region bordering the cloud and undersaturated ambient. We have analyzed the supersaturation balance in this region, which is typically a highly intermittent shearless turbulent mixing layer, under a condition where there is no mean updraft. We have investigated the evolution of the cloud–clear air interface and of the droplets therein via direct numerical simulations. We have compared horizontal averages of the phase relaxation, evaporation, reaction, and condensation times within the cloud–clear air interface for the size distributions of the initial monodispersed and polydisperse droplets. For the monodisperse population, a clustering of the values of the reaction, phase, and evaporation times, that is around 20–30 s, is observed in the central area of the mixing layer, just before the location where the maximum value of the supersaturation turbulent flux occurs. This clustering of values is similar for the polydisperse population but also includes the condensation time. The mismatch between the time derivative of the supersaturation and the condensation term in the interfacial mixing layer is correlated with the planar covariance of the horizontal longitudinal velocity derivatives of the carrier air flow and the supersaturation field, thus suggesting that a quasi-linear relationship may exist between these quantities.

A linear stability analysis of two-layer moist convection with a saturation interface

The linear convective instability of a mixture of dry air, water vapour and liquid water, with a stable unsaturated layer residing on an unstable saturated layer, is studied. It may serve as a prototype model for understanding the instability that causes mixing at the top of stratocumulus cloud or fog. Such a cloud-clear air interface is modelled as an infinitely thin saturation interface where radiative and evaporative cooling take place. The interface position is determined by the Clausius–Clapeyron equation, and can undulate with the evolution of moisture and temperature. In the small-amplitude regime two physical mechanisms are revealed. First, the interface undulation leads to the undulation of the cooling source, which destabilizes the system by superposing a vertical dipole heating anomaly on the convective cell. Second, the evolution of the moisture field induces non-uniform evaporation at the interface, which stabilizes the system by introducing a stronger evaporative cooling in the ascending region andvice versain the descending region. These two mechanisms are competing, and their relative contribution to the instability is quantified by theoretically estimating their relative contribution to buoyancy flux tendency. When there is only evaporative cooling, the two mechanisms break even, and the marginal stability curve remains the same as the classic two-layer Rayleigh–Bénard convection with a fixed cooling source.

Droplet dynamics and fine-scale structure in a shearless turbulent mixing layer with phase changes

Three-dimensional direct numerical simulations of a shearless mixing layer in a small fraction of the cloud–clear air interface are performed to study the response of an ensemble of cloud water droplets to the turbulent entrainment of clear air into a cloud filament. The main goal of this work is to understand how mixing of cloudy and clear air evolves as turbulence and thermodynamics interact through phase changes, and how the cloud droplets respond. In the main simulation case, mixing proceeds between a higher level of turbulence in the cloudy filament and a lower level of turbulence in the clear air environment – the typical shearless mixing layer set-up. Fluid turbulence is driven solely by buoyancy, which incorporates feedbacks from the temperature, the vapour content and the liquid water content fields. Two different variations on the core set of shearless mixing layer simulations are discussed, a simulation in a larger domain and a simulation with the same turbulence level inside the filament and its environment. Overall, it is found that, as evaporation occurs for the droplets that enter subsaturated clear air regions, buoyancy comes to dominate the subsequent evolution of the mixing layer. The buoyancy feedback leads initially to downdraughts at the cloudy–clear air interface and to updraughts in the bulk regions. The strength of the turbulence after initial transients depends on the domain size, showing that the range of scales is an important parameter in the shearless mixing layer set-up. In contrast, the level of turbulence in the clear air is found to have little effect on the evolution of the mixing process. The distributions of cloud water droplet size, supersaturation at the droplet positions and vertical velocity are more sensitive to domain size than to the details of the turbulence profile, suggesting that the evolution of cloud microphysics is more sensitive to large-scale as opposed to small-scale properties of the flow.

Energy and water vapor transport across a simplified cloud-clear air interface

We consider a simplified physics of the could interface where condensation, evaporation and radiation are neglected and momentum, thermal energy and water vapor transport is represented in terms of the Boussinesq model coupled to a passive scalar transport equation for the vapor. The interface is modeled as a layer separating two isotropic turbulent regions with different kinetic energy and vapor concentration. In particular, we focus on the small scale part of the inertial range of the atmospheric boundary layer as well as on the dissipative range of scales which are important to the micro-physics of warm clouds. We have numerically investigated stably stratified interfaces by locally perturbing at an initial instant the standard temperature lapse rate at the cloud interface and then observing the temporal evolution of the system. When the buoyancy term becomes of the same order of the inertial one, we observe a spatial redistribution of the kinetic energy which produce a concomitant pit of kinetic energy within the mixing layer. In this situation, the mixing layer contains two interfacial regions with opposite kinetic energy gradient, which in turn produces two intermittent sublayers in the velocity fluctuations field. This changes the structure of the field with respect to the corresponding non-stratified shearless mixing: the communication between the two turbulent region is weak, and the growth of the mixing layer stops. These results are discussed with respect to Large Eddy Simulations data for the Planetary Boundary Layers.

Modeling of subgrid-scale cloud-clear air turbulent mixing in Large Eddy Simulation of cloud fields

This paper presents computational approach to locally predict the homogeneity of subgrid-scale turbulent mixing between a cloud and its environment in large-eddy simulation of warm (ice-free) shallow convective clouds applying a double-moment bulk microphysics scheme. The term homogenity of mixing refers to the change of the mean droplet size associated with evaporation of cloud water due to entrainment. The two contrasting limits are the homogeneous mixing, where the radius of all droplets is reduced and the concentration does not change during microscale homogenization, and the extremely inhomogeneous mixing, where the microscale homogenization leads to complete evaporation of some droplets and does not affect the rest. The novel approach is applied to simulations of shallow convective cloud field. The results show that locally the homogeneity of mixing can vary significantly because of the spatial variability of the intensity of turbulence and the mean droplet radius. On average, however, the mixing becomes more homogeneous with height because of higher turbulence intensities and larger droplet sizes aloft.

Numerical Simulation of Cloud–Clear Air Interfacial Mixing: Homogeneous versus Inhomogeneous Mixing

Abstract This note presents an analysis of several dozens of direct numerical simulations of the cloud–clear air mixing in a setup of decaying moist turbulence with bin microphysics. The goal is to assess the instantaneous relationship between the homogeneity of mixing and the ratio of the time scales of droplet evaporation and turbulent homogenization. Such a relationship is important for developing improved microphysical parameterizations for large-eddy simulation of clouds. The analysis suggests a robust relationship for the range of time scale ratios between 0.5 and 10. Outside this range, the scatter of numerical data is significant, with smaller and larger time scale ratios corresponding to mixing scenarios that approach the extremely inhomogeneous and homogeneous limits, respectively. This is consistent with the heuristic argument relating the homogeneity of mixing to the time scale ratio.

Laboratory and modeling studies of cloud–clear air interfacial mixing: anisotropy of small-scale turbulence due to evaporative cooling

Small-scale mixing between cloudy air and unsaturated clear air is investigated in numerical simulations and in a laboratory cloud chamber. Despite substantial differences in physical conditions and some differences in resolved scales of motion, results of both studies indicate that small-scale turbulence generated through cloud–clear air interfacial mixing is highly anisotropic. For velocity fluctuations, numerical simulations and cloud chamber observations demonstrate that the vertical velocity variance is up to a factor of two larger than the horizontal velocity variance. The Taylor microscales calculated separately for the horizontal and vertical directions also indicate anisotropy of turbulent eddies. This anisotropy is attributed to production of turbulent kinetic energy (TKE) by buoyancy forces due to evaporative cooling of cloud droplets at the cloud–clear air interface. Numerical simulations quantify the effects of buoyancy oscillations relative to the values expected from adiabatic and isobaric mixing, standardly assumed in cloud physics. The buoyancy oscillations result from microscale transport of liquid water due to the gravitational sedimentation of cloud droplets. In the particular modeling setup considered here, these oscillations contribute to about a fifth of the total TKE production.

Numerical Simulation of Cloud–Clear Air Interfacial Mixing: Effects on Cloud Microphysics

This paper extends the previously published numerical study of Andrejczuk et al. on microscale cloud–clear air mixing. Herein, the primary interest is on microphysical transformations. First, a convergence study is performed—with well-resolved direct numerical simulation of the interfacial mixing in the limit—to optimize the design of a large series of simulations with varying physical parameters. The principal result is that all conclusions drawn from earlier low-resolution (Δx = 10−2 m) simulations are corroborated by the high-resolution (Δx = 0.25 × 10−2 m) calculations, including the development of turbulent kinetic energy (TKE) and the evolution of microphysical properties. This justifies the use of low resolution in a large set of sensitivity simulations, where microphysical transformations are investigated in response to variations of the initial volume fraction of cloudy air, TKE input, liquid water mixing ratio in cloudy filaments, relative humidity (RH) of clear air, and size of cloud droplets. The simulations demonstrate that regardless of the initial conditions the evolutions of the number of cloud droplets and the mean volume radius follow a universal path dictated by the TKE input, RH of clear air filaments, and the mean size of cloud droplets. The resulting evolution path only weakly depends on the progress of the homogenization. This is an important conclusion because it implies that a relatively simple rule can be developed for representing the droplet-spectrum evolution in cloud models that apply parameterized microphysics. For the low-TKE input, when most of the TKE is generated by droplet evaporation during mixing and homogenization, an inhomogeneous scenario is observed with approximately equal changes in the dimensionless droplet number and mean volume radius cubed. Consistent with elementary scale analysis, higher-TKE inputs, higher RH of cloud-free filaments, and larger cloud droplets enhance the homogeneity of mixing. These results are discussed in the context of observations of entrainment and mixing in natural clouds.

Mixing of cloud and clear air in centimeter scales observed in laboratory by means of Particle Image Velocimetry

Cloudy air, containing small (∅ 14 μm) water droplets, undergoes mixing with the unsaturated environment inside the cloud chamber in the process resembling smallest scales of entrainment and mixing in real clouds. Particle Image Velocimetry (PIV) applied to images from the chamber interior is used to investigate dynamics of the process in scales from 1.2 mm to few centimeters. A special algorithm, allowing for investigation of droplets motion, is developed, tested and adapted to the experimental data. Two velocity components retrieved in the vertical cross-section through the chamber interior indicate anisotropy of small-scale turbulent motions, with the preferred vertical direction. This result confirms earlier numerical studies, indicating that evaporation of cloud droplets at the cloud–clear air interface may substantially influence the small-scale turbulence in clouds.

Numerical Simulation of Cloud–Clear Air Interfacial Mixing

Abstract This paper discusses results from a series of direct numerical simulations of the microscale cloud–clear air mixing, set forth in the idealized scenario of decaying moist turbulence. In the moist case, kinetic energy of microscale motions comes not only from the classical downscale energy cascade, but it can also be generated internally due to the evaporation of cloud droplets. Three sets of numerical simulations are performed for three intensities of initial large-scale eddies. In each set, a control dry simulation is performed, as well as two moist simulations applying either bulk or detailed representation of cloud microphysics. Model results suggest that, as far as the evolutions of enstrophy and turbulent kinetic energy are concerned, the most significant impact of moist processes occurs at the low intensity of initial large-scale eddies (the input turbulent kinetic energy of 2 × 10−4 m2 s−2 resulting in the maximum eddy dissipation rate of 5 × 10−4 m2 s−3). In such a case, mixing and homoge...

Influence of Solar Wind on the Global Electric Circuit, and Inferred Effects on Cloud Microphysics, Temperature, and Dynamics in the Troposphere

There are at least three independent ways in which the solar wind modulates the flow of current density (Jz) in the global electric circuit. These are (A) changes in the galactic cosmic ray energy spectrum, (B) changes in the precipitation of relativistic electrons from the magnetosphere, and (C) changes in the ionospheric potential distribution in the polar caps due to magnetosphere-ionosphere coupling. The current density J z flows between the ionosphere and the surface, and as it passes through conductivity gradients it generates space charge concentrations dependent on J z and the conductivity gradient. The gradients are large at the surfaces of clouds and space charge concentrations of order 1000 to 10,000 elementary charges per cm3 can be generated at cloud tops. The charge transfers to droplets, many of which are evaporating at the cloud-clear air interface. The charge remains on the residual evaporation nuclei with a lifetime against leakage of order 1000 sec, and for a longer period the nuclei also retain coatings of sulfate and organic compounds adsorbed by the droplet while in the cloud. The charged evaporation nuclei become well mixed with more droplets in many types of clouds with penetrative mixing. The processes of entrainment and evaporation are also efficient for these clouds. The collection of such nuclei by nearby droplets is greatly increased by the electrical attraction between the charge on the particle and the image charge that it creates on the droplet. This process is called electroscavenging. Because the charge on the evaporation nuclei is derived from the original space charge, it depends on J z, giving a rate of electroscavenging responsive to the solar wind inputs. There may be a number of ways in which the electroscavenging has consequences for weather and climate. One possibility is enhanced production of ice. The charged evaporation nuclei have been found to be good ice forming nuclei because of their coatings, and so in supercooled clouds droplet freezing can occur by contact ice nucleation, as the evaporation nuclei are electroscavenged. Although quantitative models for the all the cloud microphysical processes that may be involved have not yet been produced, we show that for many clouds, especially those with broad droplet size distributions, relatively high droplet concentrations, and cloud top temperatures just below freezing, this process is likely to dominate over other primary ice nucleation processes. In these cases there are likely to be effects on cloud albedo and rates of sedimentation of ice, and these will depend on J z. For an increase in ice production in thin clouds such as altocumulus or stratocumulus the main effect is a decrease in albedo to incoming solar radiation, and in opacity to outgoing longwave radiation. At low latitudes the surface and troposphere heat, and at high latitudes in winter they cool. The change in meridional temperature gradient affects the rate of cyclogenesis, and the amplitude of planetary waves. For storm clouds, as in winter cyclones, the effect of increased ice formation is mainly to increase the rate of glaciation of lower level clouds by the seeder-feeder process. The increase in precipitation efficiency increases the rate of transfer of latent heat between the air mass and the surface. In most cyclones this is likely to result in intensification, producing changes in the vorticity area index as observed. Cyclone intensification also increases the amplitude of planetary waves, and shifts storm tracks, as observed. In this paper we first describe the production of space charge and the way in which it may influence the rate of ice nucleation. Then we review theory and observations of the solar wind modulation of J z, and the correlated changes in atmospheric temperature and dynamics in the troposphere. The correlations are present for each input, (A, B, and C), and the detailed patterns of responses provide support for the inferred electrical effects on the physics of clouds, affecting precipitation, temperature and dynamics.

Properties of the turbulent cloud-clear air interface observed in the laboratory experiment

Two-dimensional horizontal and vertical cuts through patterns created in the process of turbulent mixing of a cloud with unsaturated clear air are imaged in the laboratory experiment. The extracted cloud-clear air interface is then investigated in terms of a modified box-counting analysis. Results indicate that the observed interface scales in the range 2×10−2 to 5×10−1, and is not self-similar in scales below 2×10−2 m. The detailed analysis shows that the scaling exponents in the horizontal and vertical are not equal. An interface scales rather like a product of independent sets than according to the additive codimension rule, which is the case for isoconcentration surfaces of passive scalars in typical turbulent flows. Based on the above results alternative estimates of the area of the cloud-clear interface are proposed.

Laboratory investigation of the droplet concentration at the cloud-clear air interface

Preliminary results of laboratory observations of the scattered light on cloud droplets at an interface between the cloud and clear air are presented. Results suggest a weak dependence 'of the absolute value of the brightness gradient (related to the gradient of droplet concentration and liquid water content) at the interface on the interface orientation and turbulence level. There is an evidence of anisotropy (vertical direction preferred) of filaments created in the process of turbulent mixing.

Laboratory Observations of Cloud–Clear Air Mixing at Small Scales

Abstract Cloud–clear air mixing at scales from 1 mm to 1 m is observed in a laboratory chamber. Cross sections through the volume in which the mixing takes place are obtained by illuminating a planar sheet of cloud with laser light (λ = 0.488 μm, 1.2-mm thickness); the light is scattered by cloud droplets and photographed. Images indicate that complicated filament-like structures are created during mixing. Due to the properties of Mie scattering, this technique is in principle more sensitive to the larger cloud drops, and volumes with the small droplets may be underrepresented in the images. After digitization of the images, an interface between cloudy and clear-air filaments is investigated. Preliminary results indicate that at the scale of 2 cm the nature of the interface changes: at larger scales it exhibits self-similar properties, whereas at smaller scales it has a simple geometrical structure.

Fractal Analyses of High-Resolution Cloud Droplet Measurements

Abstract Fractal analyses of individual cloud droplet distributions using aircraft measurements along one-dimensional horizontal cross sections through clouds are performed. Box counting and cluster analyses are used to determine spatial scales of inhomogeneity of cloud droplet spacing. These analyses reveal that droplet spatial distributions do not exhibit a fractal behavior. A high variability in local droplet concentration in cloud volumes undergoing mixing was found. In these regions, thin filaments of cloudy air with droplet concentration close to those observed in cloud cores were found. Results suggest that these filaments may be anisotropic. Additional box counting analyses performed for various classes of cloud droplet diameters indicate that large and small droplets are similarly distributed, except for the larger characteristic spacing of large droplets. A cloud-clear air interface defined by a certain threshold of total droplet count (TDC) was investigated. There are indications that this inte...

On the Surface of Clouds

Abstract The cloud–environment mixing process is considered in terms of fractal properties of the cloud-clear air interface. The fractal dimension of the cloud surface is estimated from high-resolution airborne data. The value obtained is D = 2.55 in a range of scales from at least 10 m to over 1000 m with the possibility of even greater extension. This differs significantly from values obtained in shear-generated, well-developed, and homogeneous turbulence. The distribution of filament sizes of cloudy and clear air and estimates of the cloud surface and characteristic time of mixing process are given.

Saturation Point Analysis of Moist Convective Overturning

A unified approach to the thermodynamics of cloudy air, cloud-clear air mixing processes, atmospheric thermodynamic equilibrium structure and instability is formulated, using a new concept: the Saturation Point. This permits the representation of mixing processes and virtual potential temperature isopleths for clear and cloudy air on a thermodynamic diagram (a tephigram is used here), and their comparison with the atmospheric stratification. Illustrative examples will be given for evaporative mixing instability and convective equilibrium structure for stratocumulus, cumulus and cumulonimbus convection and convection in the incipient severe storm atmosphere.