Supercooled Ice Research Articles

Engulfment Avalanches and Thermal Hysteresis for Antifreeze Proteins on Supercooled Ice.

Antifreeze proteins (AFPs) bind to the ice-water surface and prevent ice growth at temperatures below 0 °C through a Gibbs-Thomson effect. Each adsorbed AFP creates a metastable depression on the surface that locally resists ice growth, until ice engulfs the AFP. We recently predicted the susceptibility to engulfment as a function of AFP size, distance between AFPs, and supercooling [ J. Chem. Phys. 2023, 158, 094501]. For an ensemble of AFPs adsorbed on the ice surface, the most isolated AFPs are the most susceptible, and when an isolated AFP gets engulfed, its former neighbors become more isolated and more susceptible to engulfment. Thus, an initial engulfment event can trigger an avalanche of subsequent engulfment events, leading to a sudden surge of unrestrained ice growth. This work develops a model to predict the supercooling at which the first engulfment event will occur for an ensemble of randomly distributed AFP pinning sites on an ice surface. Specifically, we formulate an inhomogeneous survival probability that accounts for the AFP coverage, the distribution of AFP neighbor distances, the resulting ensemble of engulfment rates, the ice surface area, and the cooling rate. We use the model to predict thermal hysteresis trends and compare with experimental data.

Supercooled Jahn-Teller ice

When the spins on the frustrated pyrochlore lattice obey the celebrated 2-$in$-2-$out$ ice rule, they stay in a correlated disordered phase and break the third law of thermodynamics. Similarly, if the atomic ions on the pyrochlore lattice move in and outward of the tetrahedra, they may obey a constraint resembling the ice rule. We discover that a model for pyrochlore molybdates $A_2$Mo$_2$O$_7$ ($A=$Y, Dy, Tb) exhibits a "supercooled ice" state of the displacement degrees of freedom of Mo$^{4+}$ ions, when we take account of the Jahn-Teller (JT) effect. The JT effect occurs when the lattice distortions reduce the symmetry of the local crystal field, resulting in the orbital-energy-splitting that causes the local energy gain. Unlike the standard JT effect that leads to periodic long range ordering, the displacements of Mo$^{4+}$ ions are disordered following the ice-like rule. We microscopically derive a model that describes this situation by having the 2nd and 3rd neighbor interactions between in-out lattice displacements comparably as strong as the nearest neighbor interactions of standard ice. There, the well-known nearly flat energy landscape of the ice state is altered to a metastable highly quasi-degenerate ice-like liquid state coexisting with a crystalline-like ground state. Our Monte Carlo simulations show that this liquid remains remarkably stable down to low temperatures by avoiding the putative first order transition. The relaxation in the supercooled JT ice state exhibits glassy dynamics with a plateau structure. They fit the feature of a "good glassformer" very often found in molecular liquids, but that has never been observed in material solids. The high glass-forming ability of the interacting lattice degrees of freedom will play a key role in the spin-glass transition of the material.

Detection of supercooled liquid water containing clouds with ceilometers: development and evaluation of deterministic and data-driven retrievals

Abstract. Cloud and aerosol lidars measuring backscatter and depolarization ratio are the most suitable lidars to detect cloud phase (liquid, ice, or mixed phase). However, such instruments are not widely deployed as part of operational networks. In this study, we propose a new algorithm to detect supercooled liquid water containing clouds (SLCC) based on ceilometers measuring only co-polarization backscatter. We utilize observations collected at Davis, Antarctica, where low-level, mixed-phase clouds, including supercooled liquid water (SLW) droplets and ice crystals, remain poorly understood due to the paucity of ground-based observations. A 3-month set of observations were collected during the austral summer of November 2018 to February 2019, with a variety of instruments including a depolarization lidar and a W-band cloud radar which were used to build a two-dimensional cloud phase mask distinguishing SLW and mixed-phase clouds. This cloud phase mask is used as the reference to develop a new algorithm based on the observations of a single polarization ceilometer operating in the vicinity for the same period. Deterministic and data-driven retrieval approaches were evaluated: an extreme gradient boosting (XGBoost) framework ingesting backscatter average characteristics was the most effective method at reproducing the classification obtained with the combined radar–lidar approach with an accuracy as high as 0.91. This study provides a new SLCC retrieval approach based on ceilometer data and highlights the considerable benefits of these instruments to provide intelligence on cloud phase in polar regions that usually suffer from a paucity of observations. Finally, the two algorithms were applied to a full year of ceilometer observations to retrieve cloud phase and frequency of occurrences of SLCC: SLCC was present 29 ± 6 % of the time for T19 and 24 ± 5 % of the time for G22-Davis over that annual cycle.

Response of Mixed-Phase Cloud Microphysical Properties to Cloud-Seeding Near Cloud Top Over Hebei, China

On November 29, 2019, an aircraft observation during the period of cloud-seeding was carried out for a mixed-phase cloud over Xingtai, Hebei Province, China. This study investigates the response of mixed-phase cloud microphysical properties to cloud-seeding near cloud top. Before cloud seeding, the cloud droplet concentration from fast cloud droplet probe (NC_FCDP) presented a multi-peak vertical distribution structure, with a maximum concentration of 192 cm−3 at a height of 3,322 m; the maximum concentration of ice crystals from cloud imaging probe (NC_CIP) was 10 L−1, which appeared at 4,500 m in the upper part of cloud; and the peak value of liquid water content (LWC) in the cloud also appeared at 4,500 m, with a value of 0.15 g/m3. The coexistence of supercooled liquid water and ice crystals implies that they are particularly suitable for cloud seeding at the height of 4,550 m. About 7–8 min later after cloud seeding at this height, the average NC_FCDP decreased from 160.3 to 129 cm−3, and the average NC_CIP increased from 7.1 to 10 L−1. Moreover, after cloud seeding, high NC_CIP as well as larger and more ice crystals appeared almost in the same areas within the cloud, and LWC presented an obvious decreasing trend. In contrast, the concentration of small cloud droplets and LWC decreased obviously after seeding. The findings suggest that the cloud microphysical properties showed obvious responses to the artificial introduction of silver iodide, which is important for human weather modification.

Satellite altimetry detection of ice-shelf-influenced fast ice

Abstract. The outflow of supercooled Ice Shelf Water from the conjoined Ross and McMurdo ice shelf cavity augments fast ice thickness and forms a thick sub-ice platelet layer in McMurdo Sound. Here, we investigate whether the CryoSat-2 satellite radar altimeter can consistently detect the higher freeboard caused by the thicker fast ice combined with the buoyant forcing of a sub-ice platelet layer beneath. Freeboards obtained from CryoSat-2 were compared with 4 years of drill-hole-measured sea ice freeboard, snow depth, and sea ice and sub-ice platelet layer thicknesses in McMurdo Sound in November 2011, 2013, 2017 and 2018. The spatial distribution of higher CryoSat-2 freeboard concurred with the distributions of thicker ice-shelf-influenced fast ice and the sub-ice platelet layer. The mean CryoSat-2 freeboard was 0.07–0.09 m higher over the main path of supercooled Ice Shelf Water outflow, in the centre of the sound, relative to the west and east. In this central region, the mean CryoSat-2-derived ice thickness was 35 % larger than the mean drill-hole-measured fast ice thickness. We attribute this overestimate in satellite-altimeter-obtained ice thickness to the additional buoyant forcing of the sub-ice platelet layer which had a mean thickness of 3.90 m in the centre. We demonstrate the capability of CryoSat-2 to detect higher Ice Shelf Water-influenced fast ice freeboard in McMurdo Sound. Further development of this method could provide a tool to identify regions of ice-shelf-influenced fast ice elsewhere on the Antarctic coastline with adequate information on the snow layer.

Designing Self-Sustainable Icephobic Layer by Introducing a Lubricating Un-Freezable Water Hydrogel from Sodium Polyacrylate on the Polyolefin Surface

A convenient, environment-friendly, and cost-effective method to keep anti-icing for a long time was highly desirable. Slippery lubricant layers were regarded to be effective and promising for anti-icing on different surfaces, but the drought-out of lubricants and the possible detriments to the environment were inevitable. By combining super-high molecular weight sodium polyacrylate (H-PAAS) with polyolefin through a one-pot method, a self-sustainable lubricating layer with extremely low ice adhesion of un-freezable water hydrogel was achieved at subzero conditions. The lubricant hydrogel layer could auto-spread and cover the surface of polyolefin after encountering supercooled water, frost, or ice. Due to the reduction of storage modulus in the interface, the ice adhesion of the specimen surfaces was far below 20 kPa, varying from 5.13 kPa to 18.95 kPa. Furthermore, the surfaces could preserve the fairly low adhesion after icing/de-icing cycles for over 15 times and thus exhibited sustainable durability. More importantly, this method could be introducing to various polymers and is of great promise for practical applications.

Modeling the vertical structure of the ice shelf–ocean boundary current under supercooled condition with suspended frazil ice processes: A case study underneath the Amery Ice Shelf, East Antarctica

In contrast with the severe thinning of ice shelves along the coast of West Antarctica, large ice shelves (specifically, the Filchner–Ronne and Amery Ice Shelves) with deep grounding lines gained mass during the period 1994–2012. This positive mass budget is potentially associated with the marine ice production, which originates from the supercooled Ice Shelf Water plume carrying suspended frazil ice along the ice shelf base. In addition, the outflow of this supercooled plume from beneath the ice shelf arguably exerts a significant impact on the properties of Antarctic Bottom Water, as well as its production. However, knowledge of this buoyant and supercooled shear flow is still limited, let alone its structure that is generally assumed to be vertically uniform. In this study we extended the vertical one-dimensional model of ice shelf–ocean boundary current from Jenkins (2016) by incorporating a frazil ice module and a fairly sophisticated turbulence closure (i.e., k-ε model) with the effects of density stratification. On the basis of this extended model, the study reproduced the measured thermohaline properties of a perennially-prominent supercooled ice shelf–ocean​ boundary current underneath the Amery Ice Shelf in East Antarctica, and conducted extensive sensitivity runs to a wide range of factors, including advection of scalar quantities, far-field geostrophic currents, basal slope, and the distribution of frazil ice crystal size. Based on the simulation results, the following conclusions can be drawn: Firstly, it can be difficult to reasonably reproduce the vertical structure of the ice shelf–ocean boundary current using a constant eddy viscosity/diffusivity near the ice shelf base. Secondly, although there are no direct observations of the size of frazil ice crystals beneath the ice shelves, the size of the finest ice crystals that play an important role in controlling the ice shelf–ocean boundary current is strongly suggested. Lastly, but most importantly, the ice shelf–ocean boundary layer response to the vertical gradient of frazil ice concentration will significantly reduce the level of turbulence. Therefore, this study highlights the importance of the strong interaction between frazil ice formation and the hydrodynamics and thermodynamics of ice shelf–ocean boundary layer. This interaction must not only be included, but also be resolved at high resolutions in three-dimensional coupled ice shelf–ocean models applied to cold ice cavities, which will have a potential impact on the overall ice shelf mass balance and the Antarctic Bottom Water production.

Platelet Ice Under Arctic Pack Ice in Winter

AbstractThe formation of platelet ice is well known to occur under Antarctic sea ice, where subice platelet layers form from supercooled ice shelf water. In the Arctic, however, platelet ice formation has not been extensively observed, and its formation and morphology currently remain enigmatic. Here, we present the first comprehensive, long‐term in situ observations of a decimeter thick subice platelet layer under free‐drifting pack ice of the Central Arctic in winter. Observations carried out with a remotely operated underwater vehicle (ROV) during the midwinter leg of the MOSAiC drift expedition provide clear evidence of the growth of platelet ice layers from supercooled water present in the ocean mixed layer. This platelet formation takes place under all ice types present during the surveys. Oceanographic data from autonomous observing platforms lead us to the conclusion that platelet ice formation is a widespread but yet overlooked feature of Arctic winter sea ice growth.

Variability in the Distribution of Fast Ice and the Sub‐ice Platelet Layer Near McMurdo Ice Shelf

AbstractVariability in the volume of supercooled Ice Shelf Water outflow in McMurdo Sound is reflected in the thickness and distribution of fast ice and the sub‐ice platelet layer beneath. Ground‐based electromagnetic induction and drill hole surveys of the distribution and thickness of ice shelf‐influenced fast ice and the sub‐ice platelet layer in McMurdo Sound were carried out in late spring of 2011, 2013, 2016, and 2017. In 2011 and 2017, thicker sub‐ice platelet layers of up to 7.5 and 6 m were observed, respectively. Fast ice formation throughout the winters of 2011 and 2017 was influenced by a higher occurrence of strong southerly wind events and resultant activity of the Ross Sea Polynya. In contrast, lower wind conditions in 2016 led to largely undisturbed sea ice growth and anomalously extensive fast ice coverage. A thinner sub‐ice platelet layer of up to 4 m was observed in 2016. In 2011 and 2017, substantial and variable sub‐ice platelet layers were detected in a region of exchange of water masses between the Ross Sea and the McMurdo‐Ross ice shelf cavity, which were not observed in 2013 and 2016. We hypothesize that a higher frequency of strong southerly wind events, resultant polynya activity, and High Salinity Shelf Water production over winter accelerates circulation and increases melting in the proximal shallow McMurdo Ice Shelf and the deeper Ross Ice Shelf regions of the conjoined cavity. The outflow of supercooled Ice Shelf Water and sub‐ice platelet layer formation in McMurdo Sound are consequently promoted.

Oceanographic observations in supercooled water: Protocols for mitigation of measurement errors in profiling and moored sampling

Supercooled water is liquid water that is colder than its pressure- and salinity-dependent freezing temperature. As supercooling is a metastable state, it is sensitive to small environmental changes that may force switches between solid and liquid phases, including disruption from instruments introduced to the environment. Here we present a thoroughly-tested set of observational protocols that have been developed over 15 years working in supercooled Ice Shelf Water in McMurdo Sound, Antarctica. We identify the four issues most commonly encountered in observations of supercooled ocean water as (i) the requirement for an extended process of thermal equilibration, (ii) contamination by ice within conductivity cells, (iii) anchor ice accumulation on suspended equipment, and (iv) ice coating and resultant sensing of artificially-produced latent heat. We present protocols for identifying and managing these potential modes of corruption, which include determining and allowing appropriate time for thermal equilibration, identifying and eliminating data that have been corrupted by ice capture and growth, and practical steps that can be taken to minimise accumulation of anchor ice on suspended equipment. We stress the need for planned data redundancy to enable identification and elimination of data contaminated by ice capture and phase change. We also include a complete procedure for CTD profiling which includes a 10-min soak (for an instrument initially warm relative to the ocean) and exercises data redundancy in the form of multiple up- and down-casts. Finally, we discuss how these protocols may be extended to other platforms since there is increasing scope for ice-related contamination of data as access to ice-covered oceans, and the use of autonomous vehicles in the polar oceans, increases.

Supercooling extends preservation time of human livers.

The inability to preserve vascular organs beyond several hours contributes to the scarcity of organs for transplantation1,2. Standard hypothermic preservation at +4°C1,3 limits liver preservation to less than 12 hours. Our group previously showed that supercooled ice free storage at −6°C can extend viable preservation of rat livers4,5 However, scaling supercooling preservation to human organs is intrinsically limited because of volume-dependent stochastic ice formation. Here, we describe an improved supercooling protocol that averts freezing of human livers by minimizing favourable sites of ice nucleation and homogeneous preconditioning with protective agents during machine perfusion. We show that human livers can be stored at −4°C with supercooling followed by subnormothermic machine perfusion, effectively extending the ex vivo life of the organ by 27 hours. We show that viability of livers before and after supercooling is unchanged, and that after supercooling livers can withstand the stress of simulated transplantation by ex vivo normothermic reperfusion with blood.

Interaction of supercooled droplets and nonspherical ice crystals with a solid body in a mixed cloud

Mathematical models of nonspherical particles’ motion in nonuniform flow were proposed for continuous flow approach as well as for single particles’ motion simulation. Lagrangian approach was used to simulate aerosol flow with nonspherical particles around transversal cylinder. Peculiarities of supercooled water crystallization were investigated via cryostat, infrared spectroscopy, and original devices for supercooled droplets’ impact simulation. The results of the experimental study of the supercooled water crystallization are presented. The relationship between the parameters of intermolecular interaction with the parameters of macroscopic phenomena accompanying the crystallization of supercooled water was determined. Mathematical models of the supercooled liquid crystallization and the expression for the minimum rate for the beginning of the supercooled droplets crystallization on impact with the surface of the aircraft were developed. The numerical and experimental evaluations of physical parameters characterizing the supercooled liquid crystallization in the icing problem are performed. Two-phase boundary layer formation near nanostructutred surface was investigated via molecular dynamics method.

Multiphase Reactive Transport and Platelet Ice Accretion in the Sea Ice of McMurdo Sound, Antarctica

AbstractSea ice seasonally to interannually forms a thermal, chemical, and physical boundary between the atmosphere and hydrosphere over tens of millions of square kilometers of ocean. Its presence affects both local and global climate and ocean dynamics, ice shelf processes, and biological communities. Accurate incorporation of sea ice growth and decay, and its associated thermal and physiochemical processes, is underrepresented in large‐scale models due to the complex physics that dictate oceanic ice formation and evolution. Two phenomena complicate sea ice simulation, particularly in the Antarctic: the multiphase physics of reactive transport brought about by the inhomogeneous solidification of seawater, and the buoyancy driven accretion of platelet ice formed by supercooled ice shelf water onto the basal surface of the overlying ice. Here a one‐dimensional finite difference model capable of simulating both processes is developed and tested against ice core data. Temperature, salinity, liquid fraction, fluid velocity, total salt content, and ice structure are computed during model runs. The model results agree well with empirical observations and simulations highlight the effect platelet ice accretion has on overall ice thickness and characteristics. Results from sensitivity studies emphasize the need to further constrain sea ice microstructure and the associated physics, particularly permeability‐porosity relationships, if a complete model of sea ice evolution is to be obtained. Additionally, implications for terrestrial ice shelves and icy moons in the solar system are discussed.

Airfoil performance degradation by coupling effects of supercooled raindrop icing and heavy rainfall

Icing and rainfall are the two critical meteorological factors that threaten the aircraft flight safety. A number of previous studies have shown their individual influences on aircraft aerodynamics; however, to date no studies on their coupling effects exist. In this paper, a numerical study is conducted to focus on the aerodynamic performance of a NACA 23012 airfoil exposed to heavy rain in the presence of an ice accretion by supercooled raindrop on the leading edge. An Eulerian–Lagrangian two-phase flow approach developed for rain simulation in our previous work is adopted here with some improvement in the turbulence model. A series of new phenomena about the aerodynamic performance and wake characteristics under the coupling effects of ice accretion and rainfall are found and discussed. These results do not seem to have been published previously and should be of significance to the aircraft industry for improved aircraft design and pilot training.

Icing test capabilities at DGA Aero-engine Testing

Since the 1960s, DGA Aero-engine Testing has built facilities to test and certify parts and equipment of aircraft and rotorcraft in icing conditions. For aerofoil or blade testing, a 2D test set-up has been developed. For engine tests, a free jet configuration is used. Until recently, the standard conditions for icing were associated with the presence of clouds composed of supercooled small droplets. Nevertheless, recent accidents and incidents have highlighted the danger of other icing conditions such as supercooled large drops (SLD), ice crystals and mixed phase which is composed of ice crystals and small droplets. So, current and future challenges include the ability to simulate these conditions. To generate droplets larger than 100 micrometres, a new atomiser has been developed and calibration methods have been modified in consequence. The next planned step is the generation of ice crystals and mixed phase.

New particle-dependent parameterizations of heterogeneous freezing processes: sensitivity studies of convective clouds with an air parcel model

Abstract. Based on the outcome of laboratory results, new particle-dependent parameterizations of heterogeneous freezing were derived and used to improve and extend a two-dimensional spectral microphysics scheme. They include (1) a particle-type-dependent parameterization of immersion freezing using the numbers of active sites per mass, (2) a particle-type and size-resolved parameterization of contact freezing, and (3) a particle-type-dependent description of deposition freezing. The modified microphysical scheme was embedded in an adiabatic air parcel model with entrainment. Sensitivity studies were performed to simulate convective situations and to investigate the impact of ice nuclei concentrations and types on ice formation. As a central diagnostic parameter, the ice water fraction (IWF) was selected, which is the relation of the ice water content to the total amount of water in the condensed form. The following parameters were varied: initial aerosol particle number size distributions, types of ice nucleating particles, final temperature, and the fractions of potential ice nucleating particles. Single and coupled freezing processes were investigated. The results show that immersion freezing seems to be the most efficient process. Contact freezing is constrained by the collision kernel between supercooled drops and potential ice nucleating particles. The importance of deposition freezing lies in secondary ice formation; i.e., small ice particles produced by deposition nucleation trigger the freezing of supercooled drops by collisions. Thus, a broader ice particle spectrum is generated than that by immersion and contact freezing. During coupled immersion–contact and contact–deposition freezing no competition was observed, and both processes contribute to cloud ice formation but do not impede each other. As already suggested in the literature, mineral dust particles seem to be the most important ice nucleating particles. Biological particles are probably not involved in significant ice formation. The sensitive parameters affecting cloud properties are temperature, aerosol particle composition and concentration, and particle size distribution.

Observed platelet ice distributions in Antarctic sea ice: An index for ocean‐ice shelf heat flux

AbstractAntarctic sea ice that has been affected by supercooled Ice Shelf Water (ISW) has a unique crystallographic structure and is called platelet ice. In this paper we synthesize platelet ice observations to construct a continent‐wide map of the winter presence of ISW at the ocean surface. The observations demonstrate that, in some regions of coastal Antarctica, supercooled ISW drives a negative oceanic heat flux of −30 Wm−2 that persists for several months during winter, significantly affecting sea ice thickness. In other regions, particularly where the thinning of ice shelves is believed to be greatest, platelet ice is not observed. Our new data set includes the longest ice‐ocean record for Antarctica, which dates back to 1902 near the McMurdo Ice Shelf. These historical data indicate that, over the past 100 years, any change in the volume of very cold surface outflow from this ice shelf is less than the uncertainties in the measurements.

The sub-ice platelet layer and its influence on freeboard to thickness conversion of Antarctic sea ice

Abstract. This is an investigation to quantify the influence of the sub-ice platelet layer on satellite measurements of total freeboard and their conversion to thickness of Antarctic sea ice. The sub-ice platelet layer forms as a result of the seaward advection of supercooled ice shelf water from beneath ice shelves. This ice shelf water provides an oceanic heat sink promoting the formation of platelet crystals which accumulate at the sea ice–ocean interface. The build-up of this porous layer increases sea ice freeboard, and if not accounted for, leads to overestimates of sea ice thickness from surface elevation measurements. In order to quantify this buoyant effect, the solid fraction of the sub-ice platelet layer must be estimated. An extensive in situ data set measured in 2011 in McMurdo Sound in the southwestern Ross Sea is used to achieve this. We use drill-hole measurements and the hydrostatic equilibrium assumption to estimate a mean value for the solid fraction of this sub-ice platelet layer of 0.16. This is highly dependent upon the uncertainty in sea ice density. We test this value with independent Global Navigation Satellite System (GNSS) surface elevation data to estimate sea ice thickness. We find that sea ice thickness can be overestimated by up to 19%, with a mean deviation of 12% as a result of the influence of the sub-ice platelet layer. It is concluded that within 100 km of an ice shelf this influence might need to be considered when undertaking sea ice thickness investigations using remote sensing surface elevation measurements.

Evolution of a supercooled Ice Shelf Water plume with an actively growing subice platelet matrix

AbstractWe use new observations in Western McMurdo Sound, combined with longitudinal hydrographic transects of the sound, to identify a northward‐flowing Ice Shelf Water (ISW) plume exiting the cavity of the McMurdo‐Ross Ice Shelf. We estimate the plume's net northward transport at 0.4 ± 0.1 Sv, carving out a corridor approximately 35 km wide aligned with the Victoria Land Coast. Basal topography of the McMurdo Ice Shelf is such that the plume is delivered to the surface without mixing with overlying warmer water, and is therefore able to remain below the surface freezing temperature at the point of observation beneath first‐year ice. Thus, the upper ocean was supercooled, by up to 50 mK at the surface, due to pressure relief from recent rapid ascent of the steep basal slope. The 70 m thick supercooled layer supports the growth and maintenance of a thick, semirigid, and porous matrix of platelet ice, which is trapped by buoyancy at the ice‐ocean interface. Continued growth of individual platelets in supercooled water creates significant brine rejection at the top of the water column which resulted in convection over the upper 200 m thick, homogeneous layer. By examining the diffusive nature of the intermediate water between layers of ISW and High Salinity Shelf Water, we conclude that the ISW plume must have originated beneath the Ross Ice Shelf and demonstrate that it is likely to expand eastward across McMurdo Sound with the progression of winter.

Development and Theoretical Analysis of an Aircraft Supercooled Icing Model

In the last two decades, the phenomena of rapid ice growth and severe runback were found in several serious aircraft-icing crashes, which were believed by many researchers to be due to the supercooled-large-droplet environment. However, the mechanism of the phenomena was not well acknowledged from the conclusion of the icing researches. The dendritic ice growth is thought to greatly influence the ice-growth velocity, and is rarely mentioned in aircraft-icing researches. In this paper, an aircraft supercooled icing model for ice-shape simulation is developed, based on the properties of dendritic ice growth in supercooled water. And the potential consequence of dendritic ice growth to ice-accretion rate on aircraft is discussed theoretically. The supercooled icing model is compared with the Messinger model through the analysis of icing film and the ice-shape simulation. It concludes that the supercooled icing model can well predict the ice-accretion rate, the ice shape, and the surface temperature. Compared with the Messinger model, the supercooled icing model has obvious advantages in predicting rapid ice growth and severe runback phenomena. Further investigations based on the supercooled-water physical features on aircraft will shed more light on the icing mechanism.