Ice formation and development in aged, wintertime cumulus over the UK: observations and modelling

In situ high resolution aircraft measurements of cloud microphysical properties were made in coordination with ground based remote sensing observations of a line of small cumulus clouds, using Radar and Lidar, as part of the Aerosol Properties, PRocesses And InfluenceS on the Earth's climate (APPRAISE) project. A narrow but extensive line (~100 km long) of shallow convective clouds over the southern UK was studied. Cloud top temperatures were observed to be higher than −8 °C, but the clouds were seen to consist of supercooled droplets and varying concentrations of ice particles. No ice particles were observed to be falling into the cloud tops from above. Current parameterisations of ice nuclei (IN) numbers predict too few particles will be active as ice nuclei to account for ice particle concentrations at the observed, near cloud top, temperatures (−7.5 °C). The role of mineral dust particles, consistent with concentrations observed near the surface, acting as high temperature IN is considered important in this case. It was found that very high concentrations of ice particles (up to 100 L−1) could be produced by secondary ice particle production providing the observed small amount of primary ice (about 0.01 L−1) was present to initiate it. This emphasises the need to understand primary ice formation in slightly supercooled clouds. It is shown using simple calculations that the Hallett-Mossop process (HM) is the likely source of the secondary ice. Model simulations of the case study were performed with the Aerosol Cloud and Precipitation Interactions Model (ACPIM). These parcel model investigations confirmed the HM process to be a very important mechanism for producing the observed high ice concentrations. A key step in generating the high concentrations was the process of collision and coalescence of rain drops, which once formed fell rapidly through the cloud, collecting ice particles which caused them to freeze and form instant large riming particles. The broadening of the droplet size-distribution by collision-coalescence was, therefore, a vital step in this process as this was required to generate the large number of ice crystals observed in the time available. Simulations were also performed with the WRF (Weather, Research and Forecasting) model. The results showed that while HM does act to increase the mass and number concentration of ice particles in these model simulations it was not found to be critical for the formation of precipitation. However, the WRF simulations produced a cloud top that was too cold and this, combined with the assumption of continual replenishing of ice nuclei removed by ice crystal formation, resulted in too many ice crystals forming by primary nucleation compared to the observations and parcel modelling.

Observations of supercooling and frazil ice formation in the Laptev Sea coastal polynya

This paper examines a hydrographic response to the wind‐driven coastal polynya activity over the southeastern Laptev Sea shelf for April–May 2008, using a combination of Environmental Satellite (Envisat) advanced synthetic aperture radar (ASAR) and TerraSAR‐X satellite imagery, aerial photography, meteorological data, and SBE‐37 salinity‐temperature‐depth and acoustic Doppler current profiler land‐fast ice edgemoored instruments. When ASAR observed the strongest end‐of‐April polynya event with frazil ice formation, the moored instruments showed maximal acoustical scattering within the surface mixed layer, and the seawater temperatures were either at or 0.02°C below freezing. We also find evidence of the persistent horizontal temperature and salinity gradients across the fast ice edge to have the signature of geostrophic flow adjustment as predicted by polynya models.

Study on transient ice formation of laminar flow inside externally supercooled rectangular duct

The transient process of solidification of laminar liquid flow (water) submitted to super-cooling was investigated both theoretically and experimentally. In this study an alternative analytical formulation and numerical approach were adopted resulting in the unsteady model with temperature dependent thermophysical properties in the solid region. The proposed model is based upon the fundamental equations of energy balance in the solid and liquid regions as well as across the solidification front. The basic equations and the associated boundary and initial conditions were made dimensionless by using the Landau transformation to immobilize the moving front and render the problem to a fixed plane type problem. A laminar velocity profile is admitted in the liquid domain and the resulting equations were discretized using the finite difference approach. The numerical predictions obtained were compared with the available results based on other models and concepts such as Neumann analytical model, the apparent thermal capacity model due to Bonacina and the conventional fixed grid energy model due to Goodrich. To obtain further comparisons and more validation of the model and the numerical solution, an experimental rig was constructed and instrumented permitting very well controlled experimental measurements. The numerical predictions were compared with the experimental results and the agreement was found satisfactory.

Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite

Most ice in nature forms because of impurities which boost the exceedingly low nucleation rate of pure supercooled water. However, the microscopic details of ice nucleation on these substances remain largely unknown. Here, we have unraveled the molecular mechanism and the kinetics of ice formation on kaolinite, a clay mineral playing a key role in climate science. We find that the formation of ice at strong supercooling in the presence of this clay is about 20 orders of magnitude faster than homogeneous freezing. The critical nucleus is substantially smaller than that found for homogeneous nucleation and, in contrast to the predictions of classical nucleation theory (CNT), it has a strong two-dimensional character. Nonetheless, we show that CNT describes correctly the formation of ice at this complex interface. Kaolinite also promotes the exclusive nucleation of hexagonal ice, as opposed to homogeneous freezing where a mixture of cubic and hexagonal polytypes is observed.

The effects of rotation and ice shelf topography on frazil-laden ice shelf water plumes

A model of the dynamics and thermodynamics of a plume of meltwater at the base of an ice shelf is presented. Such ice shelf water plumes may become supercooled and deposit marine ice if they rise (because of the pressure decrease in the in situ freezing temperature), so the model incorporates both melting and freezing at the ice shelf base and a multiple-size-class model of frazil ice dynamics and deposition. The plume is considered in two horizontal dimensions, so the influence of Coriolis forces is incorporated for the first time. It is found that rotation is extremely influential, with simulated plumes flowing in near-geostrophy because of the low friction at a smooth ice shelf base. As a result, an ice shelf water plume will only rise and become supercooled (and thus deposit marine ice) if it is constrained to flow upslope by topography. This result agrees with the observed distribution of marine ice under Filchner–Ronne Ice Shelf, Antarctica. In addition, it is found that the model only produces reasonable marine ice formation rates when an accurate ice shelf draft is used, implying that the characteristics of real ice shelf water plumes can only be captured using models with both rotation and a realistic topography.

The influence of ocean flow on newly forming sea ice

The heat and mass balance of the Arctic Ocean is very sensitive to the growth and decay of sea ice and the interaction between the heat and salt fields in the oceanic boundary layer. The hydraulic roughness of sea ice controls the detailed nature of turbulent fluxes in the boundary layer and hence is an important ingredient in model parameterizations. We describe a novel mechanism for the generation of corrugations of the sea ice–ocean interface, present a mathematical analysis elucidating the mechanism, and present numerical calculations for geophysically relevant conditions. The mechanism relies on brine flows developing in the sea ice due to Bernoulli suction by flow of ocean past the interface. For oceanic shears at the ice interface of 0.2 s−1, we expect the corrugations to form with a wavelength dependent upon the permeability structure of the sea ice which is described herein. The mechanism should be particularly important during sea ice formation in wind-maintained coastal polynyas and in leads. This paper applies our earlier analyses of the fundamental instability to field conditions and extends it to take account of the anisotropic and heterogeneous permeability of sea ice.

Validity of the ice shelf water plume concept beneath Filchner-Ronne Ice Shelf

In winter, brine rejection from sea ice formation and export in the Weddell Sea, offshore of Filchner-Ronne Ice Shelf (FRIS), leads to the formation of High Salinity Shelf Water (HSSW). This dense water mass enters the cavity beneath FRIS by sinking southward down the sloping continental shelf towards the grounding line. Melting occurs when the HSSW encounters the ice shelf, and the meltwater released cools and freshens the HSSW to form a water mass known as Ice Shelf Water (ISW). If this ISW rises, the ‘ice pump’ is initiated (Lewis and Perkin, 1986), whereby the ascending ISW becomes supercooled and deposits marine ice at shallower locations due to the pressure increase in the in-situ freezing temperature. Sandh¨ager et al. (2004) were able to infer the thickness patterns of marine ice deposits at the base of FRIS (figure 1), so the primary aim of this work is to try to understand the ocean flows that determine these patterns. The plume model we use to investigate ISW flow is described fully by Holland and Feltham (accepted) so only a relatively brief outline is presented here. The plume is simulated by combining a parameterisation of ice shelf basal interaction and a multiplesize- class frazil dynamics model with an unsteady, depth-averaged reduced-gravity plume model. In the model an active region of ISW evolves above and within an expanse of stagnant ambient fluid, which is considered to be ice-free and has fixed profiles of temperature and salinity. The two main assumptions of the model are that there is a well-mixed layer underneath the ice shelf and that the ambient fluid outside the plume is stagnant with fixed properties. The topography of the ice shelf that the plume flows beneath is set to the FRIS ice shelf draft calculated by Sandh¨ager et al. (2004) masked with the grounding line from the Antarctic Digital Database (ADD Consortium, 2002). To initiate the plumes, we assume that the intrusion of dense HSSW initially causes melting at the points on the grounding line where the glaciological tributaries feeding FRIS go afloat.

Development of frazil modelling in ice shelf water plumes

Abstract Preliminary results are presented from a modelling study directed at the spatial variation of frazil ice formation and its effects on flow underneath large ice shelves. The chosen plume and frazil models are briefly introduced, and results from two simplified cases are outlined. It is found that growth and melting dominate the frazil model in the short term. Secondary nucleation converts larger crystals into several nuclei due to crystal collisions (microattrition) and fluid shear and therefore governs the ice crystal dynamics after the initial supercooling has been quenched. Frazil formation is found to have a significant depth-dependence in an idealised study of an Ice Shelf Water plume. Finally, plans for more extensive and realistic studies are discussed.

Sea ice production and water mass modification in the eastern Laptev Sea

A simple polynya flux model driven by standard atmospheric forcing is used to investigate the ice formation that took place during an exceptionally strong and consistent western New Siberian (WNS) polynya event in 2004 in the Laptev Sea. Whether formation rates are high enough to erode the stratification of the water column beneath is examined by adding the brine released during the 2004 polynya event to the average winter density stratification of the water body, preconditioned by summers with a cyclonic atmospheric forcing (comparatively weakly stratified water column). Beforehand, the model performance is tested through a simulation of a well‐documented event in April 2008. Neglecting the replenishment of water masses by advection into the polynya area, we find the probability for the occurrence of density‐driven convection down to the bottom to be low. Our findings can be explained by the distinct vertical density gradient that characterizes the area of the WNS polynya and the apparent lack of extreme events in the eastern Laptev Sea. The simple approach is expected to be sufficiently rigorous, since the simulated event is exceptionally strong and consistent, the ice production and salt rejection rates are likely to be overestimated, and the amount of salt rejected is distrusted over a comparatively weakly stratified water column. We conclude that the observed erosion of the halocline and formation of vertically mixed water layers during a WNS polynya event is therefore predominantly related to wind‐ and tidally driven turbulent mixing processes.

Daily thin-ice thickness maps from modis thermal infrared imagery

Accurate knowledge of ice-production rates within the marginal ice zones of the Arctic Ocean requires monitoring of the thin-ice distribution within polynyas. The thickness of the ice layer controls the heat loss and hence the new-ice formation. An established thinice algorithm using high-resolution MODIS data allows deriving the ice-thickness distribution within polynyas. The average uncertainty is ±4.7 cm for ice thicknesses below 0.2 m. In this study, the ice-thickness distributions within the Laptev Sea polynya for the two winter seasons 2007/08 and 2008/09 are calculated. Then, a new method is applied to determine a daily MODIS thin-ice product.

Spatio-temporal variability of polynya dynamics and ice production in the Laptev Sea between the winters of 1979/80 and 2007/08

Polynyas in the Laptev Sea are examined with respect to recurrence and interannual wintertime ice production.We use a polynya classification method based on passive microwave satellite data to derive daily polynya area from long-term sea-ice concentrations. This provides insight into the spatial and temporal variability of open-water and thin-ice regions on the Laptev Sea Shelf. Using thermal infrared satellite data to derive an empirical thin-ice distribution within the thickness range from 0 to 20 cm, we calculate daily average surface heat loss and the resulting wintertime ice formation within the Laptev Sea polynyas between 1979 and 2008 using reanalysis data supplied by the National Centers for Environmental Prediction, USA, as atmospheric forcing. Results indicate that previous studies significantly overestimate the contribution of polynyas to the ice production in the Laptev Sea. Average wintertime ice production in polynyas amounts to approximately 55 km39 27% and is mostly determined by the polynya area, wind speed and associated large-scale circulation patterns. No trend in ice production could be detected in the period from 1979/80 to 2007/08.

Improvement and sensitivity analysis of thermal thin-ice thickness retrievals

Considering the sea ice decline in the Arctic during the last decades, polynyas are of high research interest since these features are core areas of new ice formation. The determination of ice formation requires accurate retrieval of polynya area and thin-ice thickness (TIT) distribution within the polynya.We use an established energy balance model to derive TITs with MODIS ice surface temperatures (Ts) and NCEP/DOE Reanalysis II in the Laptev Sea for two winter seasons. Improvements of the algorithm mainly concern the implementation of an iterative approach to calculate the atmospheric flux components taking the atmospheric stratification into account. Furthermore, a sensitivity study is performed to analyze the errors of the ice thickness. The results are the following: 1) 2-m air temperatures (Ta) and Ts have the highest impact on the retrieved ice thickness; 2) an overestimation of Ta yields smaller ice thickness errors as an underestimation of Ta; 3) NCEP Ta shows often a warm bias; and 4) the mean absolute error for ice thicknesses up to 20 cm is ±4.7 cm. Based on these results, we conclude that, despite the shortcomings of the NCEP data (coarse spatial resolution and no polynyas), this data set is appropriate in combination with MODIS Ts for the retrieval of TITs up to 20 cm in the Laptev Sea region. The TIT algorithm can be applied to other polynya regions and to past and future time periods. Our TIT product is a valuable data set for verification of other model and remote sensing ice thickness data.

Relative roles of biogenic emissions and Saharan dust as ice nuclei in the Amazon basin

Some aerosol particles, known as ice nuclei, can initiate ice formation in clouds, thereby influencing precipitation, cloud dynamics and the amount of incoming and outgoing solar radiation. In the absence of biomass burning, aerosol mass concentrations in the Amazon basin are low(1). Tropical forests emit primary biological particles directly into the atmosphere; secondary organic aerosols form from the emission and oxidation of biogenic gases(2). In addition, particles derived from biomass burning in central Africa, marine aerosols, and windblown dust from North Africa(3-5) often reach the central part of the Amazon basin during the wet season. The contribution of these aerosol sources to ice nucleation in the region is uncertain. Here we present observations of the concentration and elemental composition of ice nuclei in the Amazon basin during the wet season. Using transmission electron microscopy combined with energy-dispersive X-ray spectroscopy, we show that ice nuclei are primarily composed of carbonaceous material and dust. We show that biological particles dominate the carbonaceous fraction, whereas import of Saharan dust explains the intermittent appearance of dust-containing nuclei. We conclude that ice-nucleus concentration and abundance can be explained almost entirely by local emissions of biological particles supplemented by import of Saharan dust. Using a simple model, we tentatively suggest that the contribution of local biological particles to ice nucleation is increased at higher atmospheric temperatures, whereas the contribution of dust particles is increased at lower temperatures.

Changes in the intermediate water mass formation rates in the global ocean for the Last Glacial Maximum, mid-Holocene and pre-industrial climates

The paleoclimate version of the National Center for Atmospheric Research Community Climate System Model version 3 (NCAR-CCSM3) is used to analyze changes in the water formation rates in the Atlantic, Pacific, and Indian Oceans for the Last Glacial Maximum (LGM), mid-Holocene (MH) and pre-industrial (PI) control climate. During the MH, CCSM3 exhibits a north-south asymmetric response of intermediate water subduction changes in the Atlantic Ocean, with a reduction of 2 Sv in the North Atlantic and an increase of 2 Sv in the South Atlantic relative to PI. During the LGM, there is increased formation of intermediate water and a more stagnant deep ocean in the North Pacific. The production of North Atlantic Deep Water (NADW) is significantly weakened. The NADW is replaced in large extent by enhanced Antarctic Intermediate Water (AAIW), Glacial North Atlantic Intermediate Water (GNAIW), and also by an intensified of Antarctic Bottom Water (AABW), with the latter being a response to the enhanced salinity and ice formation around Antarctica. Most of the LGM intermediate/mode water is formed at 27.4 < sigma(theta) < 29.0 kg/m(3), while for the MH and PI most of the subduction transport occurs at 26.5 < sigma(theta) < 27.4 kg/m(3). The simulated LGM Southern Hemisphere winds are more intense by 0.2-0.4 dyne/cm(2). Consequently, increased Ekman transport drives the production of intermediate water (low salinity) at a larger rate and at higher densities when compared to the other climatic periods.