Validating ECMWF forecasts for the occurrence of ice supersaturation using visual observations of persistent contrails and radiosonde measurements over England

One of the largest uncertainties in quantifying the impact of aviation on climate concerns the formation and spreading of persistent contrails. The inclusion of a cloud scheme that allows for ice supersaturation into the integrated forecast system (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF) can be a useful tool to help reduce these uncertainties. This study evaluates the quality of the ECMWF forecasts with respect to ice super saturation in the upper troposphere by comparing them to visual observations of persistent contrails and radiosonde measurements of ice supersaturation over England. The performance of 1- to 3-day forecasts is compared including also the vertical accuracy of the supersaturation forecasts. It is found that the operational forecasts from the ECMWF are able to predict cold ice supersaturated regions very well. For the best cases Peirce skill scores of 0.7 are obtained, with hit rates at times exceeding 80% and false-alarm rates below 20%. Results are very similar for comparisons with visual observations and radiosonde measurements, the latter providing the better statistical significance.

Contemporary (2001) and ‘Little Ice Age’ glacier extents in the Buordakh Massif, Cherskiy Range, north east Siberia

The Buordakh Massif of the Cherskiy Range of sub-arctic north east Siberia, Russia has a cold continental climate and supports over 80 glaciers. Despite previous research in the region, a georeferenced map of the glaciers has only recently been completed and an enhanced version of it is reproduced in colour here. The mountains of this region reach heights in excess of 3,000 m and the glaciers on their slopes range in size from 0.1 to 10.4 km2. The mapping has been compiled through the interpretation of Landsat 7 ETM+ satellite imagery from August 2001 which has been augmented by data from a field campaign undertaken at the same time. The glaciers of the region are of the cold, ‘firn-less’ continental type and their mass balance relies heavily on the formation of superimposed ice. Moraines which lie in front of the glaciers by up to a few kilometres are believed to date from the Little Ice Age (ca. 1550-1850 AD). Over half of the glaciers mapped have shown marked retreat from these moraines.

Doppler lidar measurements of oriented planar ice crystals falling from supercooled and glaciated layer clouds

The properties of planar ice crystals settling horizontally have been investigated using a vertically pointing Doppler lidar. Strong specular reflections were observed from their oriented basal facets, identified by comparison with a second lidar pointing 4° from zenith. Analysis of 17 months of continuous high-resolution observations reveals that these pristine crystals are frequently observed in ice falling from mid-level mixed-phase layer clouds (85% of the time for layers at −15 °C). Detailed analysis of a case study indicates that the crystals are nucleated and grow rapidly within the supercooled layer, then fall out, forming well-defined layers of specular reflection. From the lidar alone the fraction of oriented crystals cannot be quantified, but polarimetric radar measurements confirmed that a substantial fraction of the crystal population was well oriented. As the crystals fall into subsaturated air, specular reflection is observed to switch off as the crystal faces become rounded and lose their faceted structure. Specular reflection in ice falling from supercooled layers colder than −22 °C was also observed, but this was much less pronounced than at warmer temperatures: we suggest that in cold clouds it is the small droplets in the distribution that freeze into plates and produce specular reflection, whilst larger droplets freeze into complex polycrystals. The lidar Doppler measurements show that typical fall speeds for the oriented crystals are ≈ 0.3 m s−1, with a weak temperature correlation; the corresponding Reynolds number is Re ∼ 10, in agreement with light-pillar measurements. Coincident Doppler radar observations show no correlation between the specular enhancement and the eddy dissipation rate, indicating that turbulence does not control crystal orientation in these clouds. Copyright © 2010 Royal Meteorological Society

Rate-dependent fracture toughness of polycrystalline ice

A series of three-point bend tests using single edge notched testpieces of pure polycrystalline ice have been performed at three different temperatures (–20°C, –30°C and –40°C). The displacement rate was varied from 1 mm/min to 100 mm/min, producing the crack tip strain rates from about 10–3 to 10–1 s–1. The results show that (a) the fracture toughness of pure polycrystalline ice given by the critical stress intensity factor (K IC) is much lower than that measured from the J—integral under identical conditions; (b) from the determination of K IC, the fracture toughness of pure polycrystalline ice decreases with increasing strain rate and there is good power law relationship between them; (c) from the measurement of the J—integral, a different tendency was appeared: when the crack tip strain rate exceeds a critical value of 6 × 10–3 s–1, the fracture toughness is almost constant but when the crack tip strain rate is less than this value, the fracture toughness increases with decreasing crack tip strain rate. Re-examination of the mechanisms of rate-dependent fracture toughness of pure polycrystalline ice shows that the effect of strain rate is related not only to the blunting of crack tips due to plasticity, creep and stress relaxation but also to the nucleation and growth of microcracks in the specimen.

Texture development in ice cream - effects of adding cocoa solids

Dynamic rheological techniques can aid the understanding of the factors contributing to ice cream structure, though the data obtained differs from that deduced from destructive techniques. Studies have shown that ice cream systems are both strain- and frequency-dependent. Chocolate ice cream is normally more viscous than the equivalent vanilla ice cream during mix preparation and has more body on freezing. Ice creams were prepared with and without cocoa solids and frequency sweeps were made from 0.1 to 100 Hz at 0.1% strain. With rapidly frozen ice creams, both G' and G" increased in the presence of cocoa solids. Comparison of mixes made with and without low-fat cocoa powder or non-gelatinizing starch demonstrated a similar relationship, with higher apparent viscosities in those mixes containing either cocoa powder or the starch. The results were consistent with the cocoa particles adding to the effect of the fat globules in increasing viscosity.

Viability of Salmonella enterica subsp. enterica during the preparation and cold storage of Egyptian soft cheeses and ice-cream

Changes occurring in the viability of Salmonella enterica subsp. enterica during the preparation and cold storage of Domiati cheese, Kariesh cheese and ice-cream were examined. A significant decrease in numbers was observed after whey drainage during the manufacture of Domiati cheese, but Salmonella remained viable for 13 weeks in cheeses prepared from milks with between 60 and 100 g/L NaCl; the viability declined in Domiati cheese made from highly salted milk during the later stages of storage. The method of coagulation used in the preparation of Kariesh cheese affected the survival time of the pathogen, and it varied from 2 to 3 weeks in cheeses made with a slow-acid coagulation method to 4-5 weeks for an acid-rennet coagulation method. This difference was attributed to the higher salt-in-moisture levels and lower pH values of Kariesh cheese prepared by the slow-acid coagulation method. A slight decrease in the numbers of Salmonella resulted from ageing ice-cream mix for 24 h at 0degreesC, but a greater reduction was evident after one day of frozen storage at -20degreesC. The pathogen survived further frozen storage for four months without any substantial change in numbers.

Combined CloudSat-CALIPSO-MODIS retrievals of the properties of ice clouds

In this paper, data from spaceborne radar, lidar and infrared radiometers on the “A-Train” of satellites are combined in a variational algorithm to retrieve ice cloud properties. The method allows a seamless retrieval between regions where both radar and lidar are sensitive to the regions where one detects the cloud. We first implement a cloud phase identification method, including identification of supercooled water layers using the lidar signal and temperature to discriminate ice from liquid. We also include rigorous calculation of errors assigned in the variational scheme. We estimate the impact of the microphysical assumptions on the algorithm when radiances are not assimilated by evaluating the impact of the change in the area-diameter and the density-diameter relationships in the retrieval of cloud properties. We show that changes to these assumptions affect the radar-only and lidar-only retrieval more than the radar-lidar retrieval, although the lidar-only extinction retrieval is only weakly affected. We also show that making use of the molecular lidar signal beyond the cloud as a constraint on optical depth, when ice clouds are sufficiently thin to allow the lidar signal to penetrate them entirely, improves the retrieved extinction. When infrared radiances are available, they provide an extra constraint and allow the extinction-to-backscatter ratio to vary linearly with height instead of being constant, which improves the vertical distribution of retrieved cloud properties.

Ice nucleation on a model hexagonal surface

The adsorption of water on a model hexagonal surface has been studied using accurate intermolecular potentials. The structure and binding energies of single molecules, clusters, and adlayers are obtained. The limiting case of weak, nondirectional surface-water interactions presented here is compared with other cases involving water-water and water-surface interactions of a similar magnitude (partial templating) and dominating water-surface interactions (perfect templating) from the literature. None of these models is conducive to the nucleation of ice, each for different reasons.Wecommenton the requirements for a good ice-nucleating surface.

Using Doppler radar with a simple explicit microphysics model to diagnose problems with ice sublimation depth-scales in forecast models

Several previous studies have attempted to assess the sublimation depth-scales of ice particles from clouds into clear air. Upon examining the sublimation depth-scales in the Met Office Unified Model (MetUM), it was found that the MetUM has evaporation depth-scales 2–3 times larger than radar observations. Similar results can be seen in the European Centre for Medium-Range Weather Forecasts (ECMWF), Regional Atmospheric Climate Model (RACMO) and Météo-France models. In this study, we use radar simulation (converting model variables into radar observations) and one-dimensional explicit microphysics numerical modelling to test and diagnose the cause of the deep sublimation depth-scales in the forecast model. The MetUM data and parametrization scheme are used to predict terminal velocity, which can be compared with the observed Doppler velocity. This can then be used to test the hypothesis as to why the sublimation depth-scale is too large within the MetUM. Turbulence could lead to dry air entrainment and higher evaporation rates; particle density may be wrong, particle capacitance may be too high and lead to incorrect evaporation rates or the humidity within the sublimating layer may be incorrectly represented. We show that the most likely cause of deep sublimation zones is an incorrect representation of model humidity in the layer. This is tested further by using a one-dimensional explicit microphysics model, which tests the sensitivity of ice sublimation to key atmospheric variables and is capable of including sonde and radar measurements to simulate real cases. Results suggest that the MetUM grid resolution at ice cloud altitudes is not sufficient enough to maintain the sharp drop in humidity that is observed in the sublimation zone.

Advancements in the estimation of ice particle fall speeds using laboratory and field measurements

Accurate estimates for the fall speed of natural hydrometeors are vital if their evolution in clouds is to be understood quantitatively. In this study, laboratory measurements of the terminal velocity vt for a variety of ice particle models settling in viscous fluids, along with wind-tunnel and field measurements of ice particles settling in air, have been analyzed and compared to common methods of computing vt from the literature. It is observed that while these methods work well for a number of particle types, they fail for particles with open geometries, specifically those particles for which the area ratio Ar is small (Ar is defined as the area of the particle projected normal to the flow divided by the area of a circumscribing disc). In particular, the fall speeds of stellar and dendritic crystals, needles, open bullet rosettes, and low-density aggregates are all overestimated. These particle types are important in many cloud types: aggregates in particular often dominate snow precipitation at the ground and vertically pointing Doppler radar measurements. Based on the laboratory data, a simple modification to previous computational methods is proposed, based on the area ratio. This new method collapses the available drag data onto an approximately universal curve, and the resulting errors in the computed fall speeds relative to the tank data are less than 25% in all cases. Comparison with the (much more scattered) measurements of ice particles falling in air show strong support for this new method, with the area ratio bias apparently eliminated.

Sea-ice anomalies in the Sea of Okhotsk and the relationship with storm tracks in the Northern Hemisphere during winter

Previous studies have shown that sea-ice in the Sea of Okhotsk can be affected by local storms; in turn, the resultant sea-ice changes can affect the downstream development of storm tracks in the Pacific and possibly dampen a pre-existing North Atlantic Oscillation (NAO) signal in late winter. In this paper, a storm tracking algorithm was applied to the six hourly horizontal winds from the National Centers for Environmental Prediction (NCEP) reanalysis data from 1978(9) to 2007 and output from the atmospheric general circulation model (AGCM) ECHAM5 forced by sea-ice anomalies in the Sea of Okhotsk. The life cycle response of storms to sea-ice anomalies is investigated using various aspects of storm activity—cyclone genesis, lysis, intensity and track density. Results show that, for enhanced positive sea-ice concentrations in the Sea of Okhotsk, there is a decrease in secondary cyclogenesis, a westward shift in cyclolysis and changes in the subtropical jet are seen in the North Pacific. In the Atlantic, a pattern resembling the negative phase of the NAO is observed. This pattern is confirmed by the AGCM ECHAM5 experiments driven with above normal sea-ice anomalies in the Sea of Okhotsk

Observations of ice multiplication in a weakly convective cell embedded in supercooled mid-level stratus

Simultaneous observations of cloud microphysical properties were obtained by in-situ aircraft measurements and ground based Radar/Lidar. Widespread mid-level stratus cloud was present below a temperature inversion (~5 °C magnitude) at 3.6 km altitude. Localised convection (peak updraft 1.5 m s−1) was observed 20 km west of the Radar station. This was associated with convergence at 2.5 km altitude. The convection was unable to penetrate the inversion capping the mid-level stratus. The mid-level stratus cloud was vertically thin (~400 m), horizontally extensive (covering 100 s of km) and persisted for more than 24 h. The cloud consisted of supercooled water droplets and small concentrations of large (~1 mm) stellar/plate like ice which slowly precipitated out. This ice was nucleated at temperatures greater than −12.2 °C and less than −10.0 °C, (cloud top and cloud base temperatures, respectively). No ice seeding from above the cloud layer was observed. This ice was formed by primary nucleation, either through the entrainment of efficient ice nuclei from above/below cloud, or by the slow stochastic activation of immersion freezing ice nuclei contained within the supercooled drops. Above cloud top significant concentrations of sub-micron aerosol were observed and consisted of a mixture of sulphate and carbonaceous material, a potential source of ice nuclei. Particle number concentrations (in the size range 0.1Ice crystal concentrations in the cloud were constant at around 0.2 L−1. It is estimated that entrainment of aerosol particles into cloud cannot replenish the loss of ice nuclei from the cloud layer via precipitation. Precipitation from the mid-level stratus evaporated before reaching the surface, whereas rates of up to 1 mm h−1 were observed below the convective feature. There is strong evidence for the Hallett-Mossop (HM) process of secondary ice particle production leading to the formation of the precipitation observed. This includes (1) Ice concentrations in the convective feature were more than an order of magnitude greater than the concentration of primary ice in the overlaying stratus, (2) Large concentrations of small pristine columns were observed at the ~−5 °C level together with liquid water droplets and a few rimed ice particles, (3) Columns were larger and increasingly rimed at colder temperatures. Calculated ice splinter production rates are consistent with observed concentrations if the condition that only droplets greater than 24 μm are capable of generating secondary ice splinters is relaxed. This case demonstrates the importance of understanding the formation of ice at slightly supercooled temperatures, as it can lead to secondary ice production and the formation of precipitation in clouds which may not otherwise be considered as significant precipitation sources.