The structure of turbulence and mixed-phase cloud microphysics in a highly supercooled altocumulus cloud

Abstract

Observations of vertically resolved turbulence and cloud microphysics in a mixed-phase altocumulus cloud are presented using in situ measurements from an instrumented aircraft. The turbulence spectrum is observed to have an increasingly negative skewness with distance below cloud top, suggesting that long-wave radiative cooling from the liquid cloud layer is an important source of turbulence kinetic energy. Turbulence measurements are presented from both the liquid cloud layer and ice virga below. Vertical profiles of both bulk and microphysical liquid and ice cloud properties indicate that ice is produced within the liquid layer cloud at a temperature of $-\mathrm{30}$ ${}^{\circ }$C. These high-resolution in situ measurements support previous remotely sensed observations from both ground-based and space-borne instruments and could be used to evaluate numerical model simulations of altocumulus clouds at spatial scales from eddy-resolving models to global numerical weather prediction models and climate simulations.

1 Introduction

Mixed-phase layer clouds are common in the Earth's atmosphere , from the tropics where detrainment from convection forms long-lived altocumulus layers to the mid-latitudes where humidity is brought to the mid-troposphere by cyclonic activity . Upwards air motion associated with gravity waves may also generate altocumulus cells.

observed that mid-latitude altocumulus layer clouds were of mixed-phase composition on more than two-thirds of occasions and that mixed-phase conditions were observed within a few tens of metres of observable cloud top and extended down through the cloud. Peak liquid water content (LWC) was found at cloud top, ice water content (IWC) reached a maximum in the lower half of the cloud system and the similarity to Arctic boundary layer mixed-phase stratocumulus was noted. More than half of the observed clouds are thinner than 500 m, with a mean LWC of 0.14 g m${}^{-\mathrm{3}}$ and with thinner clouds being correlated with lower temperatures. found, for altocumulus in the mid-latitudes, that cloud systems can consist of single and multiple layers. The maintenance of altocumulus clouds is the result of a complex network of processes relating supercooled water to ice through long-wave radiative cooling (LWRC), turbulence, underlying aerosol properties and entrainment, similar to the network in Arctic clouds described by .

The glaciation of a liquid cloud has significant consequences for fractional cloud coverage and albedo. A liquid or mixed-phase altocumulus cloud may have large areal coverage and significant optical depth, although the amount of condensed water may be relatively low, with 50 % of clouds having liquid water path (LWP) $\le$ 100 g m${}^{-\mathrm{2}}$ , as observed by in situ instrumented aircraft. Radiative transfer calculations performed by suggest that the radiative impact of the liquid layer is extremely significant. Once glaciated, the coverage can be greatly reduced and the optical depth of the ice cloud much lower and therefore understanding the processes involved in the production of ice particles is crucial for being able to quantify the radiative balance of the global climate system . showed that maintenance of the supercooled liquid layer in a large eddy simulation of mixed-phase altocumulus was dependant on good representation of the distribution of vertical velocity fluctuations, as derived from ground-based radar and lidar. Previous in situ observations show the range of turbulent fluctuations in altocumulus in the UK to be typically $\pm \mathrm{1}$ m s${}^{-\mathrm{1}}$, with a root-mean-square value of 0.5 m s${}^{-\mathrm{1}}$ in the middle of the (liquid) cloud . Similar results were found by . presented ground-based remote-sensing observations of the life cycle of a tropical altocumulus cloud system above Cabo Verde ($T=-\mathrm{34}$ ${}^{\circ }$C), which began as a liquid layer cloud, and prior to the development of ice had vertical velocity fluctuations, $w^{\prime }\pm \mathrm{1.0}$ m s${}^{-\mathrm{1}}$, with a standard deviation of 0.44 m s${}^{-\mathrm{1}}$. A warmer mid-latitude cloud with $T=-\mathrm{6}$ ${}^{\circ }$C observed from the ground by had a similar range $-\mathrm{1.5}$ m s${}^{-\mathrm{1}}<w^{\prime }<\mathrm{1.0}$ m s${}^{-\mathrm{1}}$.

found that the supercooled water in mid-level clouds is drastically underestimated in GCMs (general circulation models) and NWP (numerical weather prediction) models. The GCM simulations are found to have too little cloud in the mid-levels, resulting in a warm bias in sea surface temperatures, one of the largest of which is found in the Southern Ocean . Observations made by found that warm conveyor belt type clouds, characterised by warm advection and moderate to strong vertical ascent, i.e. altocumulus, were partly responsible for the bias. However, optically thin mid-level clouds were shown to be of low global significance by . Systematic GCM deficiencies are often attributed to low vertical resolution in the mid-levels, poor mixed-phase microphysics, and a lack of subgrid-scale processes such as cellular convection . The lifetime and albedo of the clouds are found to be extremely sensitive to the properties of aerosols and ice nucleation processes in climate models .

The problem of resolution is compounded at higher altitudes, where model levels are more widely spaced, as thinner clouds are correlated with lower temperatures , something which also makes in situ measurements of mid-latitude altocumulus clouds difficult. Many of the existing studies took place before the impact of shattering of ice on cloud physics instrumentation was fully appreciated and so there is some doubt cast on the quality of the microphysical measurements.

The aim of this paper is to report on improved cloud and aerosol microphysical observations from a detailed case study and place them in the context of new highly detailed turbulence measurements in order to better understand the processes that maintain mixed-phase altocumulus clouds in the mid-latitudes. The results should also be relevant to mixed-phase layer clouds in other geographical locations but similar temperature regimes. The following section details the instrumentation and methods of data processing. A case study of altocumulus is presented in Sect. 3. Results are presented in Sect. 4, including the thermodynamics, turbulence structure and cloud microphysics, followed by discussion and conclusions.

2 The case study

A mid-latitude cyclone was centred off the southern tip of Greenland on 2 February 2012. The warm front extended eastwards across Iceland towards Scandinavia, whilst the cold front extended roughly north–south over the Atlantic Ocean a few hundred kilometres west of Ireland. The broad warm sector covered the north-west of the UK, with high pressure to the south and east. Figure  shows a 10.8 $\mathrm{\mu }$m infrared satellite image of the North Atlantic region from the AVHRR (Advanced Very High Resolution Radiometer) for the north-eastern Atlantic, including the western UK and Ireland, Iceland, and part of Greenland. Extensive layer cloud was observed within the warm sector. Close to the low-pressure centre, west of 10${}^{\circ }$ W, the cloud top temperature was colder than $-\mathrm{50}$ ${}^{\circ }$C, and to the south-east the cloud top temperature was between $-\mathrm{18}$ and $-\mathrm{29}$ ${}^{\circ }$C, as estimated using Met Office products derived from MSG (Meteosat Second Generation) and NWP output. Discussion regarding this technique can be found in .

Figure 1

Infrared satellite imagery (10.8 $\mathrm{\mu }$m) from the AVHRR at 14:40 UTC, 2 February 2012. Layer cloud associated with the warm sector is visible to the west of the UK. The red box indicates the location of the airborne sampling.

[Figure omitted. See PDF]

In situ data were collected from the FAAM BAe146 Atmospheric Research Aircraft from single-layer mixed-phase altocumulus clouds within the warm sector of a mid-latitude cyclone, sampled on 2 February 2012 (FAAM Flight B674). Data collection began in the afternoon (16:00 UTC) with the end of the measurement period being flown in twilight conditions (19:00 UTC). The choice of flight track was restricted by air traffic considerations and so it was impossible to advect with the cloud in a Lagrangian sampling strategy. Fortunately the north–south flight track was closely aligned to the direction of the mean wind, which was predominantly from the south and ranged in strength from 6 m s${}^{-\mathrm{1}}$ at the southern end of the flight track to 8 m s${}^{-\mathrm{1}}$ in the north. There was some degree of wind shear above the cloud-containing layer, with mean wind direction being close to northwesterly. This shear may have resulted in production of turbulence in the layer below.

An optically opaque liquid cloud with cellular structure and an areal coverage of between 6 and 7 Oktas was observed upon arrival in the area, along with ice virga which extended below the liquid cloud base, as shown in the photograph in Fig. . Visual observation from the flight deck and real-time inspection of in situ data revealed that the cloud-containing layer was capped by a weak temperature inversion ($\approx \mathrm{1}$ K), below which sat the liquid cloud layer. Liquid cloud top sloped from 5800 m ($T=-\mathrm{31}$ ${}^{\circ }$C) at the southern extent of the sampling region to 5400 m ($T=-\mathrm{27}$ ${}^{\circ }$C) at the northern end. Constant-altitude flight segments interspersed with slant profiles were flown between 4500 and 7500 m. The slant profiles indicated a pseudo-adiabatic liquid water structure with larger liquid water content towards cloud tops, similar to observations by .

Figure 2

Photograph of the altocumulus layer from underneath the liquid layer, within the precipitating ice virga layer, taken on arrival in the operating area at 16:17 UTC, 2 February 2012. The broken cellular structure of the liquid altocumulus layer cloud is visible above the dashed red line. Ice virga can be seen below the liquid cloud, advected with the mean wind from right to left (south to north), between the dashed red and yellow lines. Boundary layer stratocumulus cloud is visible below the dashed yellow line. Photo: S. Abel.

[Figure omitted. See PDF]

3 Instrumentation and methods

The BAe146 carries a scientific payload capable of measuring meteorological and thermodynamic conditions and bulk and microphysical cloud properties, described in part by and .

Data were selected from times when the aircraft was flying a constant heading using the limits for rate of change of heading and roll in Table . Slant profiles were flown at a vertical rate of change of 5 m s${}^{-\mathrm{1}}$, whereas a straight and level run (SLR) had a rate of change of altitude less than 0.5 m s${}^{-\mathrm{1}}$. During sampling the aircraft typically had a nominal indicated airspeed (IAS) of 210 kn, approximately 140 m s${}^{-\mathrm{1}}$ in the mid-troposphere.

Table 1

Parameters and limits used in determining when aircraft is in a straight and level run (SLR) and a slant profile (Pro).

Segment type	Parameter	Units	Limit
SLR/Pro	$\frac{\mathrm{d}}{\mathrm{d}t}$ heading	${}^{\circ }$ s${}^{-\mathrm{1}}$	$\pm \mathrm{1.0}$
SLR/Pro	Roll	${}^{\circ }$	$\pm \mathrm{3.0}$
SLR	$\frac{\mathrm{d}}{\mathrm{d}t}$ altitude	m s${}^{-\mathrm{1}}$	$\le \pm \mathrm{0.5}$
Pro	$\frac{\mathrm{d}}{\mathrm{d}t}$ altitude	m s${}^{-\mathrm{1}}$	$\pm \mathrm{5.0}$

Turbulent wind components were sampled in three dimensions at 32 Hz using a five-port turbulence probe located on the nose of the aircraft (e.g. ). No problems due to icing of the turbulence probe pressure ports were observed during data collection.

Temperature measurements were recorded at 32 Hz by a non-de-iced Goodrich Type 102 platinum resistance thermometer and reported at 1 Hz. There was no evidence of contamination on this sensor housing due to the presence of condensed liquid water when compared against the deiced sensor. Humidity data were provided at 0.4 Hz by a near-infrared tunable diode laser (TDL) absolute humidity spectrometer, a WVSS2 (Water Vapour Sensing System Mk. 2) fitted to a flush-mounted inlet , allowing computation of dew point temperature and water vapour mixing ratio, $q_{\mathrm{v}}$.

Bulk LWC measurements were made using a Nevzorov hot-wire probe . Data from the Nevzorov total water sensor (LWC $+$ IWC) were also interrogated to identify regions of cloud-free air that would permit use of out-of-cloud aerosol particle observations; for details, see Appendix . A Lyman-alpha total water content ($q_{\mathrm{t}}$) instrument was employed to estimate total water mixing ratio, $q_{\mathrm{t}}=q_{\mathrm{v}}$ $+$ LWC $+$ IWC.

Liquid cloud particle size and number concentration were measured using a CDP (cloud droplet probe) with anti-shatter tips fitted . The performance of the probe was monitored using glass spheres of known diameter following the method presented in . The CDP was calibrated prior to the field campaign, and additional glass bead checks were routinely performed throughout in order to ensure that there was no drift in performance. The integrated LWC from CDP compared well with that measured by the Nevzorov hot-wire probe (Pearson correlation coefficient $=$ 0.98), with slight over-reading at larger sizes, where the calibration is more uncertain. Ice particle number concentrations were measured by two optical array probes (OAPs), i.e. the Cloud Imaging Probes (CIPs), with CIP15 for diameters 30 $\mathrm{\mu }$m $\le D\le \mathrm{960}$ $\mathrm{\mu }$m in 15 $\mathrm{\mu }$m increments and CIP100 (100 $\mathrm{\mu }$m $\le D\le \mathrm{6.2}$ mm) in 100 $\mathrm{\mu }$m increments . Data were processed using the SODA2 package (System for OAP Data Analysis version 2; ) from NCAR (National Centre for Atmospheric Research). Maximum observed ice crystal diameters were often less than 1 mm in this study and so the impact of shattered ice fragments on microphysical measurements is expected to be low, and careful inspection of particle imagery from the CIP probes supports this. Additionally, an algorithm to remove shattering artefacts based on the inter-arrival times of particles was applied. Imagery of cloud particles was provided by a Cloud Particle Imager (CPI) probe . Cloud fraction was computed for both liquid (CDP) and ice particle (CIP15) observations by examination of the ratio of 1 s periods with and without cloud particles in a given altitude range. Full details are given in Appendix .

Aerosol particle number concentration and size were sampled in the nominal size range 0.1 $\mathrm{\mu }$m $\le D\le \mathrm{3.0}$ $\mathrm{\mu }$m with a Passive Cavity Aerosol Spectrometer Probe (PCASP) , which was calibrated at FAAM using the method of . Full details of the data processing are given in Appendix .

In the absence of an ice nucleus counter (INC) an assessment of the ice-nucleating particle (INP) concentration was made by measuring the number concentration of aerosols with diameter, $D\ge \mathrm{0.5}$ $\mathrm{\mu }$m and applying empirically based temperature-dependent parameterisations from for general aerosol particles and for forest emissions of INP.

When flying above cloud the cloud top height (CTH) was measured using data from a downward facing Leosphere ALS450 backscatter lidar, using an integration time of 2 s and vertical resolution of 1.5 m . During a preliminary measurement leg, flying above cloud, it was observed that the cloud top height sloped along the flight track. This was coupled with the fact that air traffic restrictions prevented Lagrangian cloud sampling in the horizontal. It was therefore decided to develop an air-relative coordinate system, anchored to cloud top height (or inversion altitude) in the vertical and neglecting the horizontal dimension. Cloud top was constrained at an isentropic surface by a thermodynamic inversion which was crossed nine times between 17:10 and 17:35 UTC and a further six times later in the measurement period, between 18:25 and 18:40 UTC, in the southern half of the region. Many of these profiles contained cloud. The cloud layer below the inversion was identified from measurements of total water mixing ratio. The cloud layer was measured at $q_{\mathrm{t}}=\mathrm{0.725}$ g m${}^{-\mathrm{3}}$, as compared to the free troposphere above, where $q_{\mathrm{t}}=\mathrm{0.45}$ g m${}^{-\mathrm{3}}$. Within the cloud layer wind speed was greater in the north, and a linear fit was produced to give wind speed as a function of latitude. This allowed for an estimate of the air-relative horizontal coordinate and hence two new vertical coordinates: the vertical position of the aircraft relative to the spatially and temporally varying (i) CTH, $\mathrm{\Delta }{z}_{\left(\mathrm{CTH}\right)}$, and (ii) the inversion (top) altitude, $\mathrm{\Delta }{z}_{(\mathrm{inv}.)}$. Estimated time series of both inversion and cloud top altitude using these spot measurements were generated by interpolating in space and time throughout the measurement period.

To enable use of wind data from slant profiles, the three-dimensional wind components were filtered using a four-pole Butterworth high-pass filter, following the studies of and others (e.g. ), who have used this method for analysis of stable boundary layers. All data from both slant profiles and geometrically level flight segments were included outside of turns. Investigation of the length scales of turbulence in the system were investigated by applying a range of Butterworth filters (see Appendix ) with filter lengths between 1.5 and 16 km (Table ). Turbulence kinetic energy (TKE) was computed from the residual high-frequency fluctuations in three dimensions, for each length of filter, using 1$$\mathrm{TKE}={\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}\sqrt{{u^{\prime }}^{\mathrm{2}}+{v^{\prime }}^{\mathrm{2}}+{w^{\prime }}^{\mathrm{2}}}.$$

Application of Taylor's frozen-turbulence hypothesis permits the compositing of data from throughout the measurement period on to a “virtual meteorological mast” floating in the free troposphere.

4 Results

Figure  shows the flight profile on a latitude–altitude cross section. The potential temperature ($\mathit{\theta }$) measurements, although not uniformly distributed throughout the section, were interpolated in latitude–altitude space to give an overview of the thermodynamic environment in which the cloud was found. The liquid water content from the CDP and ice number concentration ($N_{\mathrm{I}}$) from the CIP15 are plotted and show the spatial cellular structure of the liquid water and vertical structure of the cloud system.

Figure 3

Aircraft altitude as a function of latitude from flight B674 (2 February 2012; black line), LWC from CDP (blue circles) and $N_{\mathrm{I}}$ from CIP15 (red triangles). Potential temperature ($\mathit{\theta }$) observations are interpolated in latitude–altitude space (greyscale contours).

[Figure omitted. See PDF]

The CTH derived from profile observations as a function of time and latitude was extrapolated back in time to allow comparison with the lidar-derived CTH observations from the earlier above-cloud run. Residual differences between the derived time series of CTH and lidar CTH were typically less than 50 m and in all cases less than $\approx \mathrm{100}$ m (see Table ). De-trending the lidar CTH data at two scales, 30 and 3 km, shows that the range of CTH was of a similar order to the uncertainty. The standard deviation of CTH at the 30 km scale was of the order of 25 m and at a 3 km scale was close to 12 m, with the two scales perhaps corresponding to the altocumulus cell scale (3 km) and gravity wave structures (30 km).

Table 2

Cloud top height derived from CDP LWC, inversion altitude derived from temperature profiles, lidar-derived CTH (and with 30 and 3 km trends removed), and difference (residual) between the measured inversion altitude and the derived inversion altitude, shown in terms of mean, standard deviation and range.

Parameter	Mean	SD	Range
	$[$m$]$	$[$m$]$	$[$m$]$
Cloud top	5656	150	432
Inv. top	5666	146	439
Lidar CTH	5553	58	267
Lidar CTH (30 km)	$<\mathrm{1}$	25	137
Lidar CTH (3 km)	$<\mathrm{1}$	12	91
Inv. top diff.	30	25	112

With distance from inversion altitude as the vertical coordinate it can be seen that there was an inversion in potential temperature of the order of 3 ${}^{\circ }$C that extended over less than 100 m (Fig. ). Below the inversion was a seemingly well-mixed layer in terms of potential temperature, for at least 600 m below the inversion. The air was much drier above the inversion where relative humidity w.r.t. liquid fell from saturation at the top of the layer to $\approx \mathrm{50}$ % above the inversion (Fig. c). Ice supersaturation extended down through the top 500 m of the cloud layer. The liquid cloud at the top of the layer had a depth of just over 200 m with a maximum coverage of 0.8 just below the inversion (Fig. d). Ice cloud fraction was close to zero at the top of the liquid cloud layer, increasing with depth below this to a similar value close to 0.8 some 300 m below cloud top.

Figure 4

Vertical profiles of statistics w.r.t. inversion altitude, $\mathrm{\Delta }{z}_{(\mathrm{inv}.)}$, for (a) dew point (black), as well as mean and standard deviation of temperature (red), and (b) potential temperature and (c) relative humidity w.r.t. ice (blue) and liquid (black). The boxes indicate the interquartile range (two-dimensional), the whiskers indicate the 5th and 95th percentiles, and the filled triangles indicate the 1st and 99th percentiles. (d) Areal cloud fraction calculated from number concentration of liquid (CDP) (solid black) and ice particles (CIP15) (solid blue) over 40 m altitude bins, with dashed lines showing the volume mean equivalent fraction over a 600 m deep layer (black is liquid and blue is ice) and calculated mixed-phase cloud fraction in this layer assuming maximum overlap (dotted–dashed black line).

[Figure omitted. See PDF]

The current generation of operational NWP models have a typical vertical resolution of the order of 600 m in the mid-troposphere . If the cloud fraction (cf) observations were to be volume-averaged over the full depth of the cloud layer, cf${}^{\mathrm{grid}}$ (similar to a typical NWP grid box depth), then the volume cloud fractions of the liquid and ice parts of the cloud are cf${}_{\mathrm{liquid}}^{\mathrm{grid}}\le \mathrm{0.25}$ and cf${}_{\mathrm{ice}}^{\mathrm{grid}}\le \mathrm{0.60}$. The resulting mixed-phase cloud fraction is therefore cf${}_{\mathrm{mixed}}^{\mathrm{grid}}\le \mathrm{0.15}$. The maximum value assumes that both liquid and ice are evenly distributed throughout the volume, something the observations do not support: liquid is concentrated in the top third of the cloud layer, and ice is concentrated in the lower two-thirds. This greatly reduces the volume in which liquid and ice are expected to co-exist, which has important consequences for mixed-phase microphysical processes in NWP simulations. Microphysical aspects of the cloud are discussed in detail in Sect. .

In the following section the analysis begins to focus on the vertical structure of the observed vertical velocity fluctuations in the cloud layer system described above.

Figure 5

Two-dimensional histograms for 9 km filtered data against distance from inversion altitude for (a) vertical velocity fluctuations, $w^{\prime }$ and (b) TKE, with a PDF (probability distribution function) of 32 Hz data (greyscale) and percentiles (yellow) also shown. (a) Mean value of $w^{\prime }$ (green), standard deviation of $w^{\prime }$ (magenta) and skewness of $w^{\prime }$ (cyan). (b) Mean value of TKE from the three component winds (solid green) and only the vertical component (see text for details) (dashed green).

[Figure omitted. See PDF]

4.1 Vertical profiles of TKE and vertical velocity distribution

Figure  shows profiles of vertical wind fluctuations, $w^{\prime }$ (Fig. 5a) and TKE (Fig. 5b) constructed from raw data with a high-pass Butterworth filter that is 9 km in length (see Appendix ). The percentiles are shown in 20 vertical levels, with each representing 25 km of measurements. Vertical velocity fluctuations (Fig. a) ranged between $-\mathrm{3}$ and $+\mathrm{2}$ m s${}^{-\mathrm{1}}$, and the variability, characterised by the standard deviation of $w^{\prime }$, increased from the inversion top to a maximum of 0.5 m s${}^{-\mathrm{1}}$ at $-\mathrm{150}$ m and remained reasonably constant to $-\mathrm{400}$ m before diminishing by $-\mathrm{600}$ m. The skewness of the distribution of vertical velocity fluctuations is increasingly negative with distance below cloud top, as the spectra of $w^{\prime }$ becomes increasingly dominated by downdraughts. The TKE profile (Fig. b) shows median values close to 0.5 m${}^{\mathrm{2}}$ s${}^{-\mathrm{2}}$ in a layer from $-\mathrm{150}$ to $-\mathrm{400}$ m, below which they reduce. Horizontal high-pass filtered wind components, $u^{\prime }$ and $v^{\prime }$, have higher-magnitude residual mean values close to the inversion compared to $w^{\prime }$, a result of slight wind shear across the inversion and contamination by the filter. This will result in spuriously large values for TKE in this region. By substituting the vertical velocity fluctuations, $w^{\prime }$, in place of $u^{\prime }$ and $v^{\prime }$ in the calculation of TKE (Eq. ) it can be seen that the turbulence kinetic energy tends towards zero and increases steadily below cloud top towards the local maxima at $-\mathrm{150}$ m. Below this altitude the estimate of TKE using only $w^{\prime }$ is biased high as compared to the full three-dimensional estimate.

Vertical velocity fluctuation probability distribution functions (PDFs) as a function of distance from inversion altitude are presented in Fig.  for five vertical ranges. The highest level, within cloud tops, centred on $\mathrm{\Delta }{z}_{(\mathrm{inv}.)}=-\mathrm{50}$ m, had a narrow $w^{\prime }$ distribution with low skewness ($-\mathrm{0.23}$), and had a range that was predominantly $\pm \mathrm{1.0}$ m s${}^{-\mathrm{1}}$. Within the bulk of the liquid cloud, centred on $-\mathrm{150}$ m, the range was larger with peak updraughts of 1.5 m s${}^{-\mathrm{1}}$ offset by stronger maximum downdraughts of $-\mathrm{2.0}$ m s${}^{-\mathrm{1}}$, resulting in low skewness of $-\mathrm{0.20}$. Peak updraughts were weaker in the next vertical level down, centred on $-\mathrm{250}$ m and roughly corresponding to the cloud base layer and start of sub-cloud virga, where negative skewness developed to $-\mathrm{0.62}$. Deep within the virga layer at $-\mathrm{400}$ m there will still be moderate updraughts generated, but the magnitude of the skewness was largest at $-\mathrm{0.79}$, relaxing to $-\mathrm{0.44}$ at $-\mathrm{550}$ m. The PDF at the lowest level showed that most fluctuations were of low magnitude, similar to cloud top, but occasional strong downdraughts up to $-\mathrm{1.5}$ m s${}^{-\mathrm{1}}$ were present and still influenced by LWRC from cloud top.

Figure 6

Vertical velocity fluctuations PDFs from five altitudes below inversion altitude with levels given in the legend. The first value gives the centre altitude of the level, with the second number giving skewness of the vertical velocity distribution. Standard deviations of depth in each level were between 15 and 50 m.

[Figure omitted. See PDF]

The impact of filter length on the measurements is explored to provide information about the length scales that are operating within the cloud system (Fig. ). Close to cloud top at $-\mathrm{50}$ m, both TKE (standard deviation) and skewness are low and vary little as a function of length scale. Within the liquid cloud layer at $-\mathrm{150}$ m, the TKE is observed to peak and is dominated by longer length scale circulations. There is low-magnitude negative skewness to the distribution of vertical velocity fluctuation, consistent for all filter lengths: here the induced updraughts are almost sufficient to offset the negative skewness introduced by the downdraughts. Just below the liquid cloud is where a separation of skewness as a function of length scale becomes apparent, whilst TKE is generally similar to the level above. By 250 m below cloud top the magnitude of the skewness has doubled in general, with a slightly stronger increase for longer filter lengths. Below $-\mathrm{400}$ m the downdraughts begin to dominate the spectra, with the magnitude of skewness with a 9 km being 50 % larger than for a 1.5 km filter, although the overall magnitude of turbulence intensity tends to reduce. The lower-magnitude skewness for shorter length scales at this altitude implies that resultant updraughts occur at smaller scales than the downdraughts driving the circulation. Updraughts are occurring on a scale of the order of a few hundred metres, compared to a couple of kilometres for downdraughts. The strongest difference between long and short filter lengths is apparent at $-\mathrm{500}$ m, although the overall magnitudes of skewness are lower, implying that the few downdraughts that penetrate this low do not have as large an impact in driving updraughts, possibly due to reduced convergence between ever less frequent cold pools to this depth. Turbulence intensity is a factor of 4 lower than the peak.

Figure 7

(a) Standard deviation of $w^{\prime }$ as a function of filter length, normalised to the layer maximum value for five altitudes below the inversion altitude. Length scales longer than the chosen 9 km filter length are shown with dotted line to indicate how much of the variance is at scale that pass through this filter. (b) Skewness of distribution of $w^{\prime }$ as a function of filter length, normalised to layer absolute maximum value, the same as in (a). (c) Vertical profile of standard deviation of $w^{\prime }$ relative to the inversion altitude for a range of filter lengths from 1.5 km (blue) to 14 km (green) and the limit of 16 km (dotted–dashed black); 9 km is highlighted with dotted red. (d) Vertical profile of skewness of the distribution of $w^{\prime }$ relative to inversion altitude, as (c).

[Figure omitted. See PDF]

4.2 Cloud and aerosol microphysical observations

The updraughts are responsible for the maintenance of the supercooled liquid cloud layer. The liquid water content profile as measured by the Nevzorov LWC sensor is shown in Fig.  with data plotted with respect to $\mathrm{\Delta }{z}_{\left(\mathrm{CTH}\right)}$. The data show a sub-adiabatic profile of liquid water content for the majority of observations but with the extreme largest values being close to adiabatic.

Figure 8

LWC profile percentiles relative to CTH for the mixed-phase cloud layer with data from Nevzorov LWC sensor. Theoretical adiabatic ascents are shown, starting from three altitudes (see legend for details), which correspond to potential cloud bases (see text for details). Statistics are the same as in Fig. .

[Figure omitted. See PDF]

Theoretical, undiluted adiabatic LWC profiles were calculated by assuming an ascent of a saturated air parcel from three initial altitudes corresponding to potential cloud bases. The first, from the minimum liquid cloud base at $-\mathrm{223}$ m, shows that peak observed cloud top LWC values compare well with this theoretical estimate. An ascent from $-$167 m peaks close to the 75th percentile of cloud top LWC. A third ascent from $-\mathrm{97}$ m has a peak LWC close to the 50th percentile at cloud top. Whilst entrainment of dry air from aloft at the inversion may be non-zero, these calculations demonstrate that the non-uniform cloud base may have contributed to the observed in-cloud variability in LWC at a given level. This suggests a range of turbulent eddies and updraught depths contributed to the overall spectrum of in-cloud liquid water contents.

Now the full mixed-phase cloud system is considered. Statistics were calculated for vertically resolved liquid and ice particle in-cloud number concentrations (Fig. a), effective radius (Fig. b), and condensed water content (Fig. c). Vertical layers of ice supersaturation frequency are shown to aid interpretation.

Figure 9

Vertical profiles relative to cloud top of (a) cloud particle number concentrations of liquid, $N_{\mathrm{d}}$, (CDP, cm${}^{-\mathrm{3}}$; blue) and ice, $N_{\mathrm{i}}$, (CIP15, L${}^{-\mathrm{1}}$; orange) and (CIP100, L${}^{-\mathrm{1}}$; light brown), (b) effective radius for liquid cloud particles, $R_{\mathrm{eff}}$(liquid), (CDP; blue) and ice particles $R_{\mathrm{eff}}$(ice), (CIP100; light brown). (c) Condensed water content (CWC) with LWC (CDP; blue) and IWC (CIP15; orange). The altitude of observed frequency (areal fraction) of ice supersaturation is shown for 100 % (dashed black line), 50 % (dashed–dotted black line) and 25 % (dotted black line).

[Figure omitted. See PDF]

Cloud drop numbers were constant in-cloud at $\approx \mathrm{30}$ cm${}^{-\mathrm{3}}$, implying nucleation at cloud base (Fig. a). Cloud top liquid particle effective radius was typically observed in the range 11 to 15 $\mathrm{\mu }$m (interquartile range), with smaller particles close to cloud base (Fig. b). As previously shown (Fig. ), the largest LWC values were at cloud top, peaking close to 0.1 g m${}^{-\mathrm{3}}$. Liquid cloud particle size distributions were as expected, with cloud particles growing in size towards cloud tops (Fig. ). These composite distributions, segregated by distance from cloud top (Table ) show that clouds bases typically have high number concentrations of small particles below the mode of $\approx \mathrm{11}$ $\mathrm{\mu }$m and very few particles larger than 20 $\mathrm{\mu }$m. Within the mid-cloud layer there are greater concentrations of larger particles and a similar mode. At cloud tops the largest particles are about 30 $\mathrm{\mu }$m and concentrations of particles smaller than 10 $\mathrm{\mu }$m are reduced compared to the layers below, whilst the mode remains constant.

Figure 10

(a) Mean liquid cloud particle size distributions for three vertical levels as Table , with cloud top level (black), mid-cloud (blue) and cloud bases (red). Error bars represent standard deviations. (b) Ice particle size distribution from CIP15 (dashed line) and CIP100 (solid line) with colours the same as in (a) and vertical levels the same as Table . Example particle images are taken from CPI probe in ice virga and from within the mixed-phase cloud regions, and the scale bar is shown for 100 $\mathrm{\mu }$m length, along with the time range of the imagery.

[Figure omitted. See PDF]

Table 3

Vertical level (w.r.t. CTH) ranges for liquid particle size distribution (PSD) from CDP and ice PSD from CIP probes and layer temperature ranges, with the distance spent sampling at that level along with time spent sampling mixed-phase (MP) clouds, number of cloud drops larger than given diameter and size-resolved ice production rates. ${}^{*}$ Lower level is deeper for CIP data.

				Sample size				$N_{\mathrm{D}}>D$ $[$m${}^{-\mathrm{3}}]$				Ice production rate $[$m${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}]$
	Altitude $[$m$]$			$[$km$]$	$[$km$]$	$[$s$]$		where $D$ $>$ $[$$\mathrm{\mu }$m$]$				for $D$ range $[$$\mathrm{\mu }$m$]$
Name	Base	Top	$T$ $[$K$]$	CDP	CIP	MP		18.8	19.8	21.9		60–105	120–225	240–345
Upper	$-\mathrm{50}$	0	$\mathrm{243.8}\pm \mathrm{0.7}$	44	29	31		84 456	10 944	1073		$\mathrm{1.60}\pm \mathrm{0.48}$	$\mathrm{0.59}\pm \mathrm{0.11}$	$\mathrm{0.32}\pm \mathrm{0.12}$
Intermediate	$-\mathrm{150}$	$-\mathrm{50}$	$\mathrm{243.8}\pm \mathrm{0.6}$	114	106	358		7411	687	152		$\mathrm{0.84}\pm \mathrm{0.14}$	$\mathrm{0.62}\pm \mathrm{0.06}$	$\mathrm{0.43}\pm \mathrm{0.04}$
Lower	$-\mathrm{250}$	$-\mathrm{150}$	$\mathrm{244.2}\pm \mathrm{0.5}$	58	205	244		896	149	74		$\mathrm{0.55}\pm \mathrm{0.11}$	$\mathrm{0.63}\pm \mathrm{0.06}$	$\mathrm{0.72}\pm \mathrm{0.08}$
	($-{\mathrm{450}}^{*}$)

Median ice particle number concentrations measured by CIP15 were of the order of 0.5 to 5 L${}^{-\mathrm{1}}$, peaking just below the liquid cloud base and remaining approximately constant below this whilst reducing with altitude towards the top of the liquid cloud (Fig. a). Data from CIP100 showed a similar trend, with a peak concentration just below cloud base and the concentrations falling off more rapidly towards cloud top than data from the CIP15, where particles tend to be smaller. Ice particle effective radius (Fig. b; shown for CIP100, which has a larger sample volume) increased with depth from liquid cloud top to a maximum value close to 300 m below cloud top. Peak IWC was slightly greater than peak LWC (Fig. c), whilst median values were lower in magnitude. Ice particles were not found at the very top of the liquid layer and their concentration peaked in magnitude below the liquid cloud bases. Ice supersaturation frequency (Fig. ) fell below 100 % at a level of 350 m below cloud top. IWC reduced below 600 m, where ice supersaturation frequency fell below 25 %. Effective diameter decreased below 500 m, where the ice supersaturation frequency was less than 50 %. Ice properties remained consistent, even when ice supersaturation frequency reduced, suggesting that the ice particles initially fell within branches of descending air with larger values of ice supersaturation.

Observations of ice number, effective radius and IWC are consistent with ice nucleation occurring within the liquid cloud. Following this, ice particles continue to grow and sediment through the ice supersaturated layers beneath. Mean ice particle size distributions are plotted in Fig. b according to the vertical levels in Table  to show the properties of the ice virga. Number concentrations of the smallest particles are much larger for the very top level in cloud, perhaps identifying the ice particle generation zone. Below this cloud top layer, the cumulative particle volume distribution has a constant slope up to a size of 300 $\mathrm{\mu }$m, implying that the same process is responsible for the ice production, i.e. diffusional growth. Particle imagery from CPI (Fig. b) showed the dominant ice particle type to be complex polycrystalline structures (also observed by ), which likely grew as single crystals, followed by aggregation to larger sizes. The particles share many features of “assemblages of large plates” that are present at temperatures close to $-\mathrm{30}$ ${}^{\circ }$C above water saturation (see Fig. 5 in ). Some evidence of riming was present within the liquid layer.

The ice production rate for these altocumulus clouds was calculated in a manner similar to , with the addition that the rate was computed between adjacent pairs of size channels (which have a width of 15 $\mathrm{\mu }$m) for CIP15 data up to 400 $\mathrm{\mu }$m (full details in Appendix ). Data were selected from mixed-phase regions ($N_{\mathrm{d}}>\mathrm{5}$ cm${}^{-\mathrm{3}}$ and $N_{\mathrm{i}}>\mathrm{0.1}$ L${}^{-\mathrm{1}}$ for particles with diameters larger than 100 $\mathrm{\mu }$m) for the initial cloud sampling period between 17:00 and 17:30 UTC and in the three vertical levels defined in Table . The results are shown in Fig. . Peak rates are notably greater for the level nearest to the cloud top and for sizes smaller than 140 $\mathrm{\mu }$m. At this cloud top level there are very little data for sizes larger than 230 $\mathrm{\mu }$m. For the two lower levels within the liquid cloud there is general similarity between the computed ice production rates. A summary of these results can be found in Table , which presents weighted mean ice production rates for groupings of adjacent size channels on CIP15, in the size ranges 60–105, 120–225 and 240–345 $\mathrm{\mu }$m. For the smallest particles ($D<\mathrm{105}$ $\mathrm{\mu }$m), the rate of ice production near cloud top is double that at the level below and a factor of 3 larger than at cloud base, which supports the suggestion that ice nucleation occurs at cloud top.

Figure 11

Ice production rate for three levels (Table ) in altocumulus cloud with $-\mathrm{50}$ m $<$ $z$ $<$ 0 m (black), $-\mathrm{150}$ m $<$ $z$ $<$ $-\mathrm{50}$ m (blue) and $-\mathrm{250}$ m $<$ $z$ $<$ $-\mathrm{150}$ m (red), computed between adjacent channels of CIP15 data between 45 and 400 $\mathrm{\mu }$m. Weighted mean values and errors are shown for the smallest particle size range (as Table ) at all three levels with squares and errors bars.

[Figure omitted. See PDF]

As mentioned in Sect. , ice-nucleating particles were not observed directly. Total aerosol number concentration (accumulation mode) and concentrations of the larger particles are shown in Fig. a for two layers, one above cloud and one within the cloud layer, when out of cloud and precipitation. The observed number concentrations of aerosol particles larger than 0.5 $\mathrm{\mu }$m from PCASP (Fig. a) were used as inputs into the parameterisations of for general aerosol and for forest sources of aerosol (Fig. b). Back trajectories (not shown) indicate a source region for the air mass over the boreal forests of North America with large-scale ascent in the mid-latitude storm tracks providing a potential uplift mechanism. A nucleation temperature of $-\mathrm{30}$ ${}^{\circ }$C was applied, assuming nucleation occurred at cloud top. Number concentrations of aerosol particles larger than 0.5 $\mathrm{\mu }$m are of the order of 0.07 cm${}^{-\mathrm{3}}$, both above and within the cloud layer. Resulting INP concentrations are $N_{\mathrm{INP}}=\mathrm{1.0}$ L${}^{-\mathrm{1}}$ and 0.60 L${}^{-\mathrm{1}}$ . Observed ice number concentrations are plotted in Fig. c from CIP15 and CIP100. Observed ice concentrations compare closely with the INC-corrected INP concentrations, particularly for the CIP15, which is able to observe smaller ice particles, and for the forest emissions parameterisation of . Whilst the observations show that there are larger concentrations of accumulation mode aerosols within the cloud layer than above it, the same cannot be said for number concentrations of aerosols larger than 0.5 $\mathrm{\mu }$m, limiting the ability of the data in this study to distinguish the source of the INP particles. Peak accumulation mode aerosol particle number concentrations are below 10 cm${}^{-\mathrm{3}}$, thus some particles smaller than 0.1 $\mathrm{\mu }$m must have acted as CCN (cloud condensation nuclei) to produce a cloud with droplet concentrations of 30 cm${}^{-\mathrm{3}}$.

Figure 12

(a) Percentiles of aerosol particle number concentrations larger than 0.1 $\mathrm{\mu }$m (brown) and 0.5 $\mathrm{\mu }$m (teal) from PCASP for the cloud layer (left) and above 6500 m (right). (b) INP concentrations (L${}^{-\mathrm{1}}$) using (red) and (green). (c) Ice particle number concentrations $N_{\mathrm{i}}$ from CIP15 (light blue) and CIP100 (dark blue).

[Figure omitted. See PDF]

5 Discussion

The measurements presented show a highly supercooled altocumulus layer cloud, with precipitating ice virga. Along with previously published results from the literature, the life cycle of altocumulus clouds will now be considered. Supersaturation and hence liquid cloud formation in the mid-levels of the troposphere may be achieved through large-scale ascent with upwards air motion accentuated through wind shear or gravity wave activity or convective detrainment at a stable interface (Fig. a), where the distribution of relative humidity with respect to liquid permits. showed that for tropical altocumulus clouds the liquid phase is always present before the ice phase, and found that for mid-latitude clouds observed from the ground that liquid cloud layers are situated above ice clouds in the majority of cases warmer than $-\mathrm{27}$ ${}^{\circ }$C.

Figure 13

Schematic of processes that control altocumulus clouds as a function of normalised cloud layer depth below cloud top altitude, $Z$. The lower boundary is typically diffuse, unlike for stratocumulus, which has the surface as a lower boundary (a) cloud formation through large-scale ascent, possibly augmented by wind shear or gravity wave activity acting at a stable interface in potential temperature ($\mathit{\theta }$: black trace), with a suitable relative humidity distribution (solid blue: RH w.r.t. liquid; dashed blue: RH w.r.t. ice; arrows indicate distribution up to 100 %) just below the temperature inversion. (b) Clouds (blue parcels) cool through long-wave radiative cooling from cloud top, a process which imparts turbulence to the layer below with narrow strong downdraughts (blue arrows) and broad weaker updraughts (red arrows). (b-i) to (b-v) See text for details of the resulting turbulence structure. (c) Turbulence acts to maintain the liquid water content profile (blue line) in a pseudo-adiabatic form (dashed black line). (d) At temperatures lower than 0 ${}^{\circ }$C, ice nucleation may occur, located within the liquid cloud, resulting in ice precipitating from the liquid cloud layer. (e) The turbulence kinetic energy peak occurs below the liquid cloud layer, with the skewness of the vertical velocity distribution occurring lower still.

[Figure omitted. See PDF]

Once the formation of liquid has occurred, cooling from cloud top through emission of long-wave radiation to space will generate negatively buoyant air parcels that will tend to descend through the cloud layer. Descending air parcels occur in narrow downdraughts and through mass continuity force upward motion that, dependent on humidity, may continue to support the production of supersaturation and hence liquid water. A more detailed discussion is presented below. Figure  shows that the observed liquid water content followed a pseudo-adiabatic profile; see Fig. c.

At some later stage, where the liquid clouds reside at temperatures colder than 0 ${}^{\circ }$C, the production of ice may occur (Fig. d). The computation of ice production rate suggests a tendency towards ice production close to cloud top (Fig. ), as depicted in Fig. d-i, but the increase in both ice number concentration and ice mass with distance below cloud top (Fig. a and c) indicates that the production of ice occurs at all levels within the liquid cloud. Number concentrations of ice particles remain constant in the ice-supersaturated layer below the liquid cloud (Fig. a), indicating that ice nucleation is not active in that region (Fig.  d-ii). The negatively skewed vertical velocity distribution (Fig. ) is similar to that found in nocturnal stratocumulus . (Fig. b) found a similar profile of skewness for LWRC-driven nocturnal stratocumulus clouds using ground-based measurements, in contrast with the profiles of skewness that were obtained when surface-heating-driven cumulus clouds were overhead. The turbulence kinetic energy spectra are also similar to those found in stratocumulus (Fig. e). The eventual dissipation of the altocumulus clouds may occur through erosion of humidity within the cloud layer through precipitation or an increase in subsidence, which would act to warm the layer . Likewise, turbulent mixing and radiation can also be important.

The turbulence observations presented in Sect.  will now be placed into context by comparing them with the remote-sensing observations and derived conceptual model of of the circulation structures within altocumulus clouds. Other observations from ground-based studies conducted in both the tropics and mid-latitudes, as well as previously published in situ observations, will be drawn upon where appropriate. At the top of the cloud, where $\mathrm{\Delta }{z}_{(\mathrm{inv}.)}=-\mathrm{50}$ m (Fig. ), these in situ observations show a narrow distribution of $w^{\prime }$ with a low skewness value, features which are likely to be characteristic of shallow cloud top eddies, similar to those resolved by the remote-sensing observations in Fig. 7 of : at this level, close to the inversion, the broader and deeper circulations do not have the opportunity to develop, as depicted in Fig. b-i.

The PDF of $w^{\prime }$ at $\mathrm{\Delta }{z}_{(\mathrm{inv}.)}=-\mathrm{150}$ m has positive and negative fluctuations of almost equal magnitude and hence low skewness (Fig. ), which therefore corresponds to the “cloud layer” where Rayleigh–Bénard type cells were observed (Fig. b-ii). also found a symmetrical distribution of vertical velocities on a horizontal scale comparable to or just larger than the liquid cloud depth. The previous in situ observations of turbulence in altocumulus clouds were made within the liquid layer and therefore all correspond to these liquid layer observations and show similar magnitudes. and both presented remote-sensing results for the liquid portion of altocumulus clouds with similar magnitudes and distributions.

Altocumulus clouds often occur in multiple layers . These could form when the resulting updraughts produce liquid supersaturation at a level below the cloud layer driving the turbulence (Fig. b-iii). It is possible that long-lived altocumulus clouds such as those observed by and were maintained through a particular set of conditions, including vertical stability and the distribution of humidity, which permitted the formation of secondary “daughter cells” at the same altitude as the driving clouds.

Within the ice virga, between $-\mathrm{250}$ and $-\mathrm{400}$ m below cloud top (Fig. ), the observed turbulence properties differed from those found in the liquid cloud layer with a shift to increasingly negative skewness; these results are consistent with the “Subcloud Layer 1” of , corresponding to Fig. b-iv, with downwards penetrating mammatus type features. This connection between the liquid cloud layer and below was observed by to be on a longer length scale than the in-cloud circulations of the order of 1 km. presented turbulence observation from within the ice virga layer (using radar sensitive to the larger ice particles) for late-stage development clouds, which showed a shift to negative skewness, $-\mathrm{1.75}$ m s${}^{-\mathrm{1}}<w^{\prime }<\mathrm{0.75}$ m s${}^{-\mathrm{1}}$, and an increase in standard deviation of w${}^{\prime }$, comparable to the direct turbulence observations presented here. Unlike remote-sensing-derived observations of turbulence, these in situ measurements are not reliant on, or biased by, observations of, or assumptions about the ice particle size distribution. At the lowest level of the in situ observations ($\mathrm{\Delta }{z}_{(\mathrm{inv}.)}=-\mathrm{500}$ m, Fig. ), the magnitude of both skewness and TKE is lower, and the width of the PDF of $w^{\prime }$ is as narrow as at cloud top, sharing characteristics with the “Subcloud Layer 2” of (Fig. b-v).

There is observational evidence in the literature of a link between turbulence and the production of ice in supercooled layer clouds. For example, found that ice production in two altocumulus clouds, both at similar temperatures to each other and to the clouds in this study, differed along with the magnitude of turbulence and liquid water content. One of the clouds had peak LWC $=$ 0.05 g m${}^{-\mathrm{3}}$ with vertical velocity values of $\pm \mathrm{0.75}$ m s${}^{-\mathrm{1}}$, and produced virga with IWC $=$ 0.08 g m${}^{-\mathrm{3}}$. The second cloud had peak LWC $=$ 0.02 g m${}^{-\mathrm{3}}$, and negligible ice concentrations, whilst the vertical velocity fluctuations were $<\mathrm{0.25}$ m s${}^{-\mathrm{1}}$. Crucially, that cloud with ice had liquid particles of 15 $\mathrm{\mu }$m at cloud top, whereas the low LWC ice-free cloud had higher drop concentrations and a mean liquid particle diameter of less than 6 $\mathrm{\mu }$m.

saw ice enhancement in Arctic stratiform clouds close to cloud tops when liquid cloud particles were larger than $\approx \mathrm{20}$ $\mathrm{\mu }$m. Table  presents integrated number concentrations of liquid cloud drops, $N_{\mathrm{D}}$, larger than three particular diameters, close to 20 $\mathrm{\mu }$m for the three vertical cloud levels. For a particular minimum size the concentrations are greater towards the top of the cloud along with the observed ice production rate. There is a steep gradient in concentrations as a function of size at each level in this size range, in the case of the top cloud level, up to 3 orders of magnitude across the three size bins. If indeed the production of ice is dependent on the presence of “large” liquid drops then it could be particularly sensitive to the size distribution of liquid particles. Taking the cloud top maximum ice production rate of 1.6 m${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$ and numbers of large drops given in Table , the production of ice could be sustained over 14 h, assuming that the droplets larger than 18.8 $\mathrm{\mu }$m form ice. The steep slope of the PSD in these clouds results in a strong sensitivity to this calculation of ice production timescale. Close to 2 h of ice production could be maintained for ice-forming droplets larger than 19.8 $\mathrm{\mu }$m but only 11 min if ice formation required liquid particles larger than 21.9 $\mathrm{\mu }$m.

later showed that the CCN budget was capable of modulating the ice phase but rejected the notion that the inhibition of freezing of small drops due to dissolved solute (noted by ) was important in the clouds, as the cloud particles typically exceeded 10 $\mathrm{\mu }$m, thus negating the effect. postulated that the ice may have been formed during partial evaporation of cloud droplets during mixing of dry air from aloft. Later work defined such a process as contact nucleation inside out (CNIO), similar to the INP production through evaporation mechanism that demonstrated could explain ice production in Arctic mixed-phase clouds in eddy-permitting cloud simulations. Increased turbulence resulting from additional cooling in higher-LWC clouds may promote saturation and evaporation cycling of cloud and aerosol particles, which would lead to an increase in the rate of CNIO events for a given temperature range. This mechanism may be able to explain many of the supercooled layer cloud observations, including those presented here.

Many of the in situ and remote-sensing observations discussed above show similar properties to the turbulence structure, seemingly independent of temperature and geographical location. Given that there is a great variety of altocumulus clouds , this is somewhat surprising and therefore additional observations and studies using eddy permitting simulations should be used to determine the turbulence budget of altocumulus clouds as a function of parameters including temperature, LWP and CCN. Such studies could then inform representation or parameterisation of weakly forced turbulent layer clouds in NWP and climate models.

New instrumentation including holographic imaging probes may prove useful in trying to identify the conditions and locations in which the ice first forms in liquid layer clouds. Observations are required that can positively identify small ice ($\ge \mathrm{10}$ $\mathrm{\mu }$m) in the presence of large liquid drops ($\approx \mathrm{25}$ $\mathrm{\mu }$m) or that could identify small ice at cloud edges in regions of evaporation and entrainment. Supercooled stratocumulus clouds may provide a more accessible natural laboratory.

6 Conclusions

Highly supercooled single-layer mixed-phase mid-level altocumulus clouds with a cloud top temperature of $-\mathrm{30}$ ${}^{\circ }$C with precipitating ice virga were measured by an instrumented aircraft. There have been few detailed observational studies of this kind mainly because the clouds are usually transient in nature. The in situ observations presented in this paper show the turbulence structure and microphysical properties of this cloud system at high vertical resolution.

Cloud top was found to slope significantly and thus cloud top and inversion height were estimated as a function of time and space. This allowed the data to be studied purely as a function of the vertical dimension, as though collected by a virtual meteorological mast floating in the free troposphere as the clouds advected past. Comparison of thermodynamic and cloud microphysics data from both geometrically level flight segments and slant profiles was therefore possible, increasing the quantity of data available. These observations are a snapshot in time of a dynamic environment and there is an implicit assumption that the turbulence is frozen in an extension of Taylor's hypothesis for these airborne measurements . Application of the same technique to existing layer cloud flight data may be useful. The method was accurate to a scale of approximately 30 m in the vertical, as confirmed by above-cloud lidar observations.

The key features observed are turbulence sustained by long-wave radiative cooling, resulting in a vertical velocity distribution where negative skewness increases with distance below cloud top and a PDF in the range $\pm \mathrm{2}$ m s${}^{-\mathrm{1}}$ but typically $\pm \mathrm{1}$ m s${}^{-\mathrm{1}}$. Direct measurement of the turbulence in the ice virga layer below the supercooled liquid cloud layer are presented for the first time. These in situ observations support the observation-based conceptual model of turbulence in altocumulus shown by , as well as other remotely sensed ground-based observations of turbulence and the cloud life cycle observations of . The measurements are also similar to previous in situ observations made in the liquid layers of altocumulus clouds .

The adiabatic-type liquid water content profiles were shown to be consistent with updraughts from variable cloud bases, consistent with turbulence driven from cloud top. The deepest liquid cells were just over 200 m in depth. Cloud drop number concentration was constant in height at 30 cm${}^{-\mathrm{3}}$, with the modal and maximum size increasing with altitude, suggesting little entrainment of dry air from aloft. Production of ice within the liquid cloud layer was shown to occur at a rate of $\mathrm{0.84}\pm \mathrm{0.14}$ m${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$ in the middle of the liquid cloud, with evidence that freezing events may be enhanced close to cloud top by a factor of 2 or more to $\mathrm{1.60}\pm \mathrm{0.48}$ m${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$. Ice particle concentrations appear to be well represented by INP concentrations derived from in situ observations of aerosol particle size distributions and parameterisations of and .

Appendix A Processing of data from microphysics probes

Optical array probe data from the CIP15 and CIP100 were processed using SODA 2 software . Specific SODA2 settings used were (i) to not reconstruct particles that are at the edge of the diode array and (ii) not apply corrections for out-of-focus liquid drops. OAP probes occasionally suffer from “stuck bits” where one pixel remains constantly “on”. Inspection of the imagery did not indicate that this problem occurred for any of the pixels and so the option to correct for stuck bits was turned off. Counts per bin were corrected for size-dependent depth of field and ice particle number concentration calculated using measured true airspeed (TAS) and

A1$$N_{\mathrm{i}}=\left({\displaystyle \frac{\sum _{\mathrm{bin}=\mathrm{1}}^{n_{\mathrm{chn}}}\mathrm{count}\left(\mathrm{bin}\right)}{\mathrm{DoF}\left(\mathrm{Fn}\right(\mathrm{bin}\left)\right)}}\right)/(\mathrm{SA}\times \mathrm{TAS}).$$ The errors in sample volume can be up to 100 % for the smallest size bins. To reduce counting uncertainty, individual bins were combined at larger sizes: CIP15 at 100 to 420 $\mathrm{\mu }$m (2 bins), 435 to 600 $\mathrm{\mu }$m (4 bins), and $>\mathrm{615}$ $\mathrm{\mu }$m (8 bins) and CIP100 at 700 to 2800 $\mathrm{\mu }$m (2 bins), 2900 to 4000 $\mathrm{\mu }$m (4 bins), and $>\mathrm{4100}$ $\mathrm{\mu }$m (8 bins). CIP15 data and the CIP100 data agree well where sizes overlap. IWC was calculated from the integrated size distribution using the mass-dimensional relation from .

Following , the data from bins adjacent to gain-stage crossovers in the PCASP were summed into a single wider bin, namely bins 4 and 5 and bins 15 and 16. The lowest size channel was rejected as the lower bound is unknown. A composite error was calculated by combining fractional errors from each contributor: bin size, counts and flow rate. It is apparent that, even for low bin counts, the major contribution to the uncertainty in PCASP was from the sizing.

PCASP aerosol data are only valid when out of cloud due to contamination by a break-up of cloud and precipitation particles. Cloud-free regions were determined using the standard deviation of raw power on the Nevzorov total water content (TWC) probe (TWC $=$ LWC $+$ IWC), exploiting the observation that variability in-cloud is very different from that out-of-cloud. A 1 Hz time series of the standard deviation of electrical power was computed from the 32 Hz data record, and a threshold of 2.0 mW ($\approx \mathrm{1}\times {\mathrm{10}}^{-\mathrm{4}}$ g m${}^{-\mathrm{3}}$) was placed on this parameter to partition the data into cloud-contaminated and clear-sky time periods with a 2 s “safety window” to account for in-cloud variability, potential timing offsets between individual data logging system clocks and cloud edges where the condensed water content may fall below the sensitivity of the probe. Typical out-of-cloud variability was between 0.5 and 3 mW, and when cloud microphysics probes reported cloud particles the Nevzorov TWC reported orders of magnitude greater variability of up to 1 W.

Appendix B Butterworth filters

Turbulent fluctuations within the planetary boundary layer are typically calculated using Reynolds decomposition from vertical stacks of level flight legs by removing a linear trend from the data record. Cloud boundaries of mid-tropospheric layer clouds tend to follow isentropic surfaces, removed from the surface of the Earth. This leads to difficulty when attempting to use this filtering method in altocumulus cloud layers, as geometrically level flight legs are effectively slant profiles through thermodynamic space. Here the high-frequency fluctuations in horizontal and vertical wind parameters were extracted from the data record by filtering using high-pass four-pole Butterworth filters and taking data from both level flight and slant profiles.

A range of filter lengths were applied to investigate scales of motion within the cloud system (Table ). It was assumed that TAS $=$ 140 m s${}^{-\mathrm{1}}$ throughout the measurement period as the true value varied by less than 5 %.

Table B1

Butterworth filter lengths and frequencies used to investigate turbulence.

Filter length $($km$)$	Frequency $($Hz$)$
1.5	0.0933
2.0	0.0700
2.5	0.0560
3.0	0.0467
3.5	0.0400
4.0	0.0350
4.5	0.3111
5.0	0.0280
6.0	0.0234
7.0	0.2000
8.0	0.0175
9.0	0.0156
12.0	0.0117
14.0	0.0100
16.0	0.0088

Assuming that the synoptic and turbulence scales are spectrally distinct, a perfect high-pass filter would result in the mean of the parameter in question being equal to zero at all locations. Throughout the depth of the layer the residual mean of the high-pass filtered vertical velocity fluctuations is $\le \mathrm{0.01}$ m s${}^{-\mathrm{1}}$ for filter lengths below 5 km and only as large as $\le \mathrm{0.02}$ m s${}^{-\mathrm{1}}$ for a 12 km filter length. The greatest magnitude for any given filter occurred close to the inversion altitude as a result of smoothing of the vertical wind shear across the inversion, thus contaminating the residual turbulence fluctuations.

Loss of variance (TKE) occurs for filter lengths of up to 9 km: the resolved TKE falls below 90 % of the maximum value for filter lengths shorter than 9 km, below 75 % at 6 km and below 50 % at 2.5 km. The reduction in overall TKE indicates that scales of this order were contributing to the turbulent fluxes in the vicinity. As a compromise, the 9 km filter was taken to be sufficient to separate the turbulence scales of motion from synoptic and other larger scales of motion, and is set as the default filter.

Appendix C Calculation of ice production rate

Calculations of ice production rate from observations were first performed by (hereafter HHC87), originally to calculate the rime-splintering secondary ice production (SIP) rate in cumulus clouds within the Hallett–Mossop (HM) temperature range , and have been since repeated by . Those same calculations were applied here to measure primary ice production rate (PIP) in these mixed-phase altocumulus layer clouds.

HHC87 calculated the ice production rate, $P$, from the difference in the cumulative size distribution of measured ice particles, $C$, between two size thresholds, $L_{\mathrm{1}}$ and $L_{\mathrm{2}}$, C1$$P=\left[C\left(L_{\mathrm{2}}\right)-C\left(L_{\mathrm{1}}\right)\right]/t_{\mathrm{21}},$$ where the growth time is given by C2$$t_{\mathrm{21}}=\left(L_{\mathrm{2}}-L_{\mathrm{1}}\right)/G\left(T\right),$$ and $G\left(T\right)$ is the average ice particle growth rate. Using the growth curve from for the water-saturated “Region A” (temperature-dependent growth rate for plate-like particles) and the measured layer temperatures (Table ), the growth rates were found to be of the order of $G\left(R\right)=\mathrm{0.38}\pm \mathrm{0.02}$ $\mathrm{\mu }$m s${}^{-\mathrm{1}}$. It was assumed that the growing ice particles remained in a water-saturated environment and that in the early stages growth through riming and aggregation was negligible and diffusional growth dominated. Ice particles were assumed to have negligible size upon production as the maximum liquid cloud particles had diameters smaller than $\le \mathrm{30}$ $\mathrm{\mu }$m.

Only ice particles larger than 137 $\mathrm{\mu }$m were considered in HHC87, but with the current generation of shadow imaging probes and the maximum cloud droplet size in these altocumulus clouds of only 30 $\mathrm{\mu }$m, this restriction can be reduced to smaller sizes. The ice production rate was calculated for particles with diameters larger than 30 $\mathrm{\mu }$m and smaller than 400 $\mathrm{\mu }$m.

In order to validate this extension to the technique, the same method was applied to calculations of ice production in two cumulus clouds at $T=-\mathrm{5}$ and $T=-\mathrm{8}$ ${}^{\circ }$C previously sampled by on 24 November 2013. Both clouds were in the SIP temperature range. Liquid particles as large as 100 $\mathrm{\mu }$m were observed and thus the calculations considered only ice particles with sizes larger than 100 $\mathrm{\mu }$m and smaller than 200 $\mathrm{\mu }$m. Weighted mean values are $P^{[-\mathrm{5}{}^{\circ }\mathrm{C}]}=\mathrm{31.44}\pm \mathrm{4.05}$ and $P^{[-\mathrm{8}{}^{\circ }\mathrm{C}]}=\mathrm{11.86}\pm \mathrm{2.38}$ m${}^{-\mathrm{3}}$ s${}^{-\mathrm{1}}$, more than an order of magnitude larger than those found in these altocumulus clouds and comparable to previous measurements made by HHC87 and (Fig. ).

Figure C1

Ice production rate calculated between adjacent OAP bins smaller than 400 $\mathrm{\mu }$m for the altocumulus clouds in this study at $-\mathrm{30}$ ${}^{\circ }$C (black) and also for passes through two cumulus clouds in the HM temperature range from flight B816 at $-\mathrm{5}$ ${}^{\circ }$C (blue) $-\mathrm{8}$ ${}^{\circ }$C (red). Previous measurements in HM zone are shown from HHC87 (green) and (grey).

[Figure omitted. See PDF]

For the altocumulus clouds the computed rate is a lower bound because aggregation would act to reduce the observed rate, as would growth through riming, which was observed on some CPI particle imagery. present results from previous observational and modelling studies that indicate a range of minimum riming dimension thresholds, $D_{\mathrm{thresh}}$, for various ice particle habits, between 35 $\mathrm{\mu }$m for hexagonal columns and 200 $\mathrm{\mu }$m for broad-branched plates.

Code and data availability

Data from the FAAM research flights are available from the Centre for Environmental Data Analysis, https://catalogue.ceda.ac.uk/uuid/829a7d5425424584aa8455ca512ad0da (last access: 17 February 2020) . Code for processing of core data from FAAM research flights is available from https://github.com/ncasuk/decades-pp (last access: 17 February 2020) . Other bespoke processing code is available on request.

Author contributions

PAB designed and conducted the research flight, analysed the data, and wrote the paper. AB and PRAB were involved in the development of the analysis methods. AB, PRAB and SJA reviewed the written document.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

Airborne data were obtained using the BAe146-301 Atmospheric Research Aircraft operated by Directflight Ltd and managed by the FAAM, which is a joint entity of NERC and the Met Office. The staff of the Met Office, FAAM, Directflight Ltd, Avalon Engineering and BAE Systems are thanked for their work that made the research flights a success. Franco Marenco (Met Office) is thanked for his assistance in providing the cloud top height product from the lidar data. James Dorsey (University of Manchester) is thanked for providing the data and imagery from the CPI probe. We thank two anonymous reviewers for their comments which helped us to improve the manuscript.

Review statement

This paper was edited by Timothy Garrett and reviewed by two anonymous referees.