High-altitude Ice Research Articles

Seasonal vertical distributions of diurnal variation of ice cloud frequency by CATS measurements over a global region (51°S-51°N)

Ice clouds can strongly affect the radiation budget through the reflection (absorption) of solar shortwave (longwave) radiation in the Earth-atmosphere system. Unlike water clouds, ice clouds present a longer lifetime in the upper troposphere. Therefore, the vertical patterns of diurnal variation of ice clouds are of very importance to understand climate change in the globe. The present paper investigates the diurnal variation of ice cloud frequency (ICF) from March 2015 to October 2017 over the global scale (51 ° S and 51 ° N) by CATS (Cloud-Aerosol Transport System) LiDAR onboard the International Space Station (ISS) and CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) LiDAR measurements. Our results suggest that: (1) ICF derived from the CATS and CALIPSO LiDAR presents the similar geographical variations. (2) The vertical distribution of ICF in land area over the North Hemisphere (NH) is found larger than in the South Hemisphere (SH) likely due to the more dust aerosols as ice nuclei to influence ice cloud formation over the NH. (3) ICF at nighttime (00:00–06:00 and 18:00–23:00) exhibits higher value than daytime (06:00–18:00), especially in high altitude (more than 12 km) ice clouds in the tropics. (4) The number of ice cloud samples is higher in spring and winter compared to summer and autumn, and the highest amount of ice cloud samples in spring(winter) occurs during nighttime(daytime) for most hours. These findings provide significant insights that may help improving the ice cloud's representation in climate models.

Alluvial-lacustrine record of Mid- to Late-Holocene moisture variations trend verified by multiple proxies in the middle and lower reaches of the Hutubi River, northwest China

Alluvial-lacustrine deposits are important carriers for studying paleo-hydrological information and the evolution of paleo-environments. In this study, a total of 6 optically stimulated luminescence (OSL) samples from 2 profiles in the lower and middle reaches of the Hutubi River in the North Tianshan Mountains were dated, and a comprehensive analysis was conducted by integrating sedimentary structures, grain sizes, and magnetic susceptibility characteristics to explore the sedimentary environment patterns and regional evolution processes during the Mid to late Holocene. The results revealed the following findings: (1) The OSL signal of quartz samples indicates the predominance of fast components, allowing reliable dating of regional fluvial sediments using the appropriate condition-tested coarse-grained quartz single-aliquot regenerative dose (SAR) protocol within the 90–125 μm range. (2) Analysis of sedimentary structures, grain sizes, and magnetic susceptibility characteristics reveals lithofacies of heterotopic synchronous or homotopic asynchronous in exposed profiles, with predominantly poorly sorted fine and very fine sand. Magnetic susceptibility is primarily controlled by the content of magnetic minerals in the coarse particles, while also being influenced by the combined effects of provenance input and depositional environment, indicating a complex regional sedimentary environment influenced by variable hydrodynamic conditions driven by fluvial processes. (3) Based on OSL dating results and the analysis of various environmental proxy indicators, it is determined that the climatic environment in the Middle and Late-Holocene in the downstream areas of the northern piedmonts of the Tianshan Mountains demonstrates a fluctuating pattern of wet and dry alternations. The sedimentary environment shifted from a stable state in the Mid-Holocene to an unstable state in the Late-Holocene. The intermittent appearance of wet records may be associated with increased river activity caused presumably by the release of high-altitude ice due to rising temperatures and warming conditions.

The Mars Science Laboratory record of optical depth measurements via solar imaging

The Mars Science Laboratory Curiosity rover has monitored the Martian environment in Gale crater since landing in 2012. This study reports the record of optical depth derived from visible and near-infrared images of the Sun. Aerosol optical depth, which is mostly due to dust but also includes ice, dominates the record, with gas optical depth too small to measure. The optical depth record includes the effects of regional dust storms and one planet-encircling dust event, showing the expected peaks during southern spring and summer and relatively lower and more stable optical depth in fall and winter. The measurements show that there is a seasonally varying diurnal change in dust load, with the optical depth peaking in the morning during southern spring and summer, correlated with thermotidal pressure changes. However, there was no systematic diurnal change during autumn and winter, except after one regional storm. There were indications that the dust was relatively enhanced at high altitudes during high-optical-depth periods and that high-altitude ice was significant during winter. The observations did not provide much information about particle size or composition, but they were consistent with a smaller particle size after aphelion (in southern winter). No scattering halos were seen in associated sky images, even when there was visual evidence of ice hazes or clouds, which suggests small or amorphous ice particles. Unexpectedly, the measurement campaign revealed that the cameras collected saltating sand in their sunshades 1.97 m above the surface. As a result, the measurement strategy had to be adjusted to avoid high-elevation imaging to avoid sand covering the optics.

Ice Crystal Size Distributions in Tropical Mesoscale Convective Systems in the Vicinity of Darwin, Australia: Results from the HAIC/HIWC Campaign

Abstract Total ice water content (IWC) derived from an isokinetic evaporator probe and ice crystal particle size distributions (PSDs) measured by a two-dimensional stereo probe and precipitation imaging probe installed on an aircraft during the 2014 European High Altitude Ice Crystals–North American High IWC field campaign (HAIC/HIWC) were used to characterize regions of high IWC consisting mainly of small ice crystals (HIWC_S) with IWC ≥ 1.0 g m−3 and median mass diameter (MMD) < 0.5 mm. A novel fitting routine developed to automatically determine whether a unimodal, bimodal, or trimodal gamma distribution best fits a PSD was used to compare characteristics of HIWC_S and other PSDs (e.g., multimodality, gamma fit parameters) for HIWC_S simulations. The variation of these characteristics and bulk properties (MMD, IWC) was regressed with temperature, IWC, and vertical velocity. HIWC_S regions were most pronounced in updraft cores. The three modes of the PSD reveal different dominant processes contributing to ice growth: nucleation for maximum dimension D < 0.15 mm, diffusion for 0.15 < D < 1.0 mm, and aggregation for D > 1.0 mm. The frequency of trimodal distributions increased with temperature. The volumes of equally plausible parameters derived in the phase space of gamma fit parameters increased with temperature for unimodal distributions and, for temperatures less than −27°C, for multimodal distributions. Bimodal distributions with 0.4 mm in the larger mode were most common in updraft cores and HIWC_S regions; bimodal distributions with 0.4 mm in the smaller mode were least common in convective cores.

Comprehensive Efficiency Evaluation of Aircraft Artificial Cloud Seeding in Hunan Province, China, Based on Numerical Simulation Catalytic Method

Aircraft cloud seeding refers to the use of equipment on aircraft to release chemicals into clouds, changing their physical and chemical properties to increase rainfall or snowfall. The purpose of precipitation enhancement is to alleviate drought and water scarcity issues. Due to the complexity of the technology, the precise control of factors such as cloud characteristics and chemical release amounts is necessary. Therefore, a scientific evaluation of the potential of aircraft cloud seeding can help to improve the effectiveness of the process, and is currently a technical challenge in weather modification. This study used the mesoscale numerical model WRF coupled with a catalytic process to simulate and evaluate the seven aircraft cloud seeding operations conducted in Hunan Province in 2021. The results show that WRF can effectively evaluate the effectiveness of cloud seeding. When the water vapor conditions are suitable, the airborne dispersion of silver iodide (AgI) can significantly increase the content of large particles of high-altitude ice crystals, snow, and graupel, resulting in an increase in low-level rainwater content and, correspondingly, an increase in ground precipitation. When the water vapor conditions are insufficient, the dispersion of AgI does not trigger effective precipitation, consistent with the results of station observations and actual flight evaluations. This study provides an effective method for scientifically evaluating the potential and effectiveness of aircraft cloud seeding operations.

Dependence of Ice Crystal Size Distributions in High Ice Water Content Conditions on Environmental Conditions: Results from the HAIC-HIWC Cayenne Campaign

Abstract A new method that automatically determines the modality of an observed particle size distribution (PSD) and the representation of each mode as a gamma function was used to characterize data obtained during the High Altitude Ice Crystals and High Ice Water Content (HAIC-HIWC) project based out of Cayenne, French Guiana, in 2015. PSDs measured by a 2D stereo probe and a precipitation imaging probe for particles with maximum dimension (Dmax) > 55 μm were used to show how the gamma parameters varied with environmental conditions, including temperature (T) and convective properties such as cloud type, mesoscale convective system (MCS) age, distance away from the nearest convective peak, and underlying surface characteristics. Four kinds of modality PSDs were observed: unimodal PSDs and three types of multimodal PSDs (Bimodal1 with breakpoints 100 ± 20 μm between modes, Bimodal2 with breakpoints 1000 ± 300 μm, and Trimodal PSDs with two breakpoints). The T and ice water content (IWC) are the most important factors influencing the modality of PSDs, with the frequency of multimodal PSDs increasing with increasing T and IWC. An ellipsoid of equally plausible solutions in (No–λ–μ) phase space is defined for each mode of the observed PSDs for different environmental conditions. The percentage overlap between ellipsoids was used to quantify the differences between overlapping ellipsoids for varying conditions. The volumes of the ellipsoid decrease with increasing IWC for most cases, and (No–λ–μ) vary with environmental conditions related to distribution of IWC. HIWC regions are dominated by small irregular ice crystals and columns. The parameters (No–λ–μ) in each mode exhibit mutual dependence.

The impacts of secondary ice production on microphysics and dynamics in tropical convection

Abstract. Secondary ice production (SIP) is an important physical phenomenon that results in an increase in the ice particle concentration and can therefore have a significant impact on the evolution of clouds. In this study, idealized simulations of a mesoscale convective system (MCS) were conducted using a high-resolution (250 m horizontal grid spacing) mesoscale model and a detailed bulk microphysics scheme in order to examine the impacts of SIP on the microphysics and dynamics of a simulated tropical MCS. The simulations were compared to airborne in situ and remote sensing observations collected during the “High Altitude Ice Crystals – High Ice Water Content” (HAIC-HIWC) field campaign in 2015. It was found that the observed high ice number concentration can only be simulated by models that include SIP processes. The inclusion of SIP processes in the microphysics scheme is crucial for the production and maintenance of the high ice water content observed in tropical convection. It was shown that SIP can enhance the strength of the existing convective updrafts and result in the initiation of new updrafts above the melting layer. Agreement between the simulations and observations highlights the impacts of SIP on the maintenance of tropical MCSs in nature and the importance of including SIP parameterizations in models.

Ice particle distribution simulation in a transonic axial compressor based on Eulerian evaporation coupling model

High-altitude ice crystal icing of aero engine is a serious threat to flight safety. Previous studies on ice crystal icing focused on the influence of air flow on ice crystals, but ice crystals also affect the engine flow field, which is often ignored in the research on the influence of ice crystals on engines. Numerical simulation based on Eulerian method is adopted to realize the two-way coupling between the compressor air flow and the particles in this study. The approach is demonstrated using the NASA compressor stage 35. The changes of compressor and particle parameters with different inlet total water content, relative humidity, and median volume diameter are calculated and analyzed, and the influence of ice crystal on the compressor performance is studied. The results show that the variation of relative humidity has a great influence on the particle temperature, median volume diameter, and wet bulb temperature. Median volume diameter has a great influence on the melt ratio. The variation of total water content has little effect on particle temperature, total water content, median volume diameter, and wet bulb temperature. The particle parameters are affected by the flow field of the compressor. The parameters show that the icing is easy to occur at the leading edge of NASA stage 35 stator. By contrast, the overall compressor characteristics, after ice particle injection, the total pressure ratio, and isentropic efficiency of the compressor are increased without considering ice crystal accretion, and the chocking boundary and stall boundary are not affected.

Simulation of in-ice cosmic ray air shower induced particle cascades

We present simulations of the development of high-energy cosmic-ray air showers penetrating high-altitude ice layers that can be found at the polar regions. We use a combination of the corsika Monte Carlo code and the geant4 simulation toolkit, and focus on the particle cascade that develops in the ice to describe its most prominent features. We discuss the impact of the ice layer on the total number of particles in function of depth of the air shower, and we give a general parametrization of the charge distribution in the cascade front in function of ${X}_{\mathrm{max}}$ of the cosmic ray air shower, which can be used for analytical and semianalytical calculations of the expected Askaryan radio emission of the in-ice particle cascade. We show that the core of the cosmic ray air shower dominates during the propagation in ice, therefore creating an in-ice particle cascade strongly resembling a neutrino-induced particle cascade. Finally, we present the results of simulations of the Askaryan radio emission of the in-ice particle cascade, showing that the emission is dominated by the shower core, and discuss the feasibility of detecting the plasma created by the particle cascade in the ice using RADAR echo techniques.

A global view on stratospheric ice clouds: assessment of processes related to their occurrence based on satellite observations

Abstract. Ice clouds play an important role in regulating water vapor and influencing the radiative budget in the atmosphere. This study investigates stratospheric ice clouds (SICs) in the latitude range between ±60∘ based on the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). As polar stratospheric clouds include other particles, they are not discussed in this work. Tropopause temperature, double tropopauses, clouds in the upper troposphere and lower stratosphere (UTLS), gravity waves, and stratospheric aerosols are analyzed to investigate their relationships with the occurrence of and variability in SICs in the tropics and at midlatitudes. We found that SICs with cloud-top heights of 250 m above the first lapse rate tropopause are mainly detected in the tropics. Monthly time series of SICs from 2007 to 2019 show that high occurrence frequencies of SICs follow the Intertropical Convergence Zone (ITCZ) over time in the tropics and that SICs vary interannually at different latitudes. Results show that SICs associated with double tropopauses, which are related to poleward isentropic transport, are mostly found at midlatitudes. More than 80 % of the SICs around 30∘ N/S are associated with double tropopauses. Correlation coefficients of SICs and all the other abovementioned processes confirm that the occurrence of and variability in SICs are mainly associated with the tropopause temperature in the tropics and at midlatitudes. UTLS clouds, which are retrieved from the Atmospheric Infrared Sounder (AIRS) and used as a proxy for deep convection in the tropics and high-altitude ice cloud sources at midlatitudes, have the highest correlations with SICs in the monsoon regions and the central United States. Gravity waves are mostly related to SICs at midlatitudes, especially over Patagonia and the Drake Passage. However, the second-highest correlation coefficients show that the cold tropopause temperature, the occurrence of double tropopauses, high stratospheric aerosol loading, frequent UTLS clouds, and gravity waves are highly correlated with the SICs locally. The long-term anomaly analyses show that interannual anomalies of SICs are correlated with the tropopause temperature and stratospheric aerosols instead of the UTLS clouds and gravity waves. The overlapping and similar correlation coefficients between SICs and all processes mentioned above indicate strong associations between those processes themselves. Due to their high inherent correlations, it is challenging to disentangle and evaluate their contributions to the occurrence of SICs on a global scale. However, the correlation coefficient analyses between SICs and all abovementioned processes (tropopause temperature, double tropopauses, clouds in the upper troposphere and lower stratosphere (UTLS), gravity waves, and stratospheric aerosols) in this study help us better understand the sources of SICs on a global scale.

Microphysical processes producing high ice water contents (HIWCs) in tropical convective clouds during the HAIC-HIWC field campaign: dominant role of secondary ice production

Abstract. High ice water content (HIWC) regions in tropical deep convective clouds, composed of high concentrations of small ice crystals, were not reproduced by Weather Research and Forecasting (WRF) model simulations at 1 km horizontal grid spacing using four different bulk microphysics schemes (i.e., the WRF single‐moment 6‐class microphysics scheme (WSM6), the Morrison scheme and the Predicted Particle Properties (P3) scheme with one- and two-ice options) for conditions encountered during the High Altitude Ice Crystals (HAIC) and HIWC experiment. Instead, overestimates of radar reflectivity and underestimates of ice number concentrations were realized. To explore formation mechanisms for large numbers of small ice crystals in tropical convection, a series of quasi-idealized WRF simulations varying the model resolution, aerosol profile, and representation of secondary ice production (SIP) processes are conducted based on an observed radiosonde released at Cayenne during the HAIC-HIWC field campaign. The P3 two-ice category configuration, which has two “free” ice categories to represent all ice-phase hydrometeors, is used. Regardless of the horizontal grid spacing or aerosol profile used, without including SIP processes the model produces total ice number concentrations about 2 orders of magnitude less than observed at −10 ∘C and about an order of magnitude less than observed at −30 ∘C but slightly overestimates the total ice number concentrations at −45 ∘C. Three simulations including one of three SIP mechanisms separately (i.e., the Hallett–Mossop mechanism, fragmentation during ice–ice collisions, and shattering of freezing droplets) also do not replicate observed HIWCs, with the results of the simulation including shattering of freezing droplets most closely resembling the observations. The simulation including all three SIP processes produces HIWC regions at all temperature levels, remarkably consistent with the observations in terms of ice number concentrations and radar reflectivity, which is not replicated using the original P3 two-ice category configuration. This simulation shows that primary ice production plays a key role in generating HIWC regions at temperatures <-40 ∘C, shattering of freezing droplets dominates ice particle production in HIWC regions at temperatures between −15 and 0 ∘C during the early stage of convection, and fragmentation during ice–ice collisions dominates at temperatures between −15 and 0 ∘C during the later stage of convection and at temperatures between −40 and −20 ∘C over the whole convection period. This study confirms the dominant role of SIP processes in the formation of numerous small crystals in HIWC regions.

Экспедиция Института географии им. В.Б. Сочавы СО РАН и клуба Портулан в Восточные Саяны (горный массив Мунку-Сардык) в апреле-мае 2021 года

The results of the spring expedition of the Student and Teaching Club "Portulan" held in 2021 together with the staff of the Institute of Geography of the Siberian Branch of the Russian Academy of Sciences and the Institute of the Earth's Crust of the Siberian Branch of the Russian Acad- emy of Sciences to the area of Mount Munku-Sardyk are summed up, within the framework of which field scientific studies of nival-glacial formations, landscape-forming processes, slope forms of the relief of the territory of the basins of the White Irkut and Muguvek rivers were carried out. Traditionally, the weather, the state of the Peretolchina glacier (northern), high-altitude ice and the Active mountain stream were monitored, according to which detailed georadar profiling was carried out

Криогенные литопотоки горного массива Мунку-Сардык

The article deals with exogenous cryogenic relief-forming litopotoks: periglacial (supraglacial) zone, modern glacial zone; the thermal zone of paleoglacials of regressive stage, erosive or fluvial zone of proluvial, solifluction, high-altitude ice and permafrost-stone mountain streams.

Экспедиции клуба Портулан в район г. Мунку-Сардык в 2018 году

The results of the 2018 expeditions of the student and teacher club "Portulan" together with the staff of the Institute of Geography of the Siberian Branch of the Russian Academy of Sciences to the area of mount Munku-Sardyk are summed up, within the framework of which field scientific studies of nival-glacial formations, river and sloped ice, slope forms of relief and their processes forming along the valleys of the White Irkut, Bugovek, Muguvek rivers were carried out: regime surveys of high-altitude ice, as well as weather monitoring was traditionally carried out, study of the stone stream "Active".

Экспедиции клуба Портулан в район г. Мунку-Сардык в 2019 году

The results of the 2019 expeditions of the student and teacher club "Portulan" together with the staff of the Institute of Geography of the Siberian Branch of the Russian Academy of Sciences to the area of mount Munku-Sardyk are summed up, within the framework of which field scientific studies of nival-glacial formations, river and sloped ice, slope forms of relief and their processes forming along the valleys of the White Irkut, Bugovek, Muguvek rivers were carried out: regime surveys of high-altitude ice, as well as weather monitoring was traditionally carried out, study of the stone stream "Active".

Dependence of Ice Microphysical Properties On Environmental Parameters: Results from HAIC-HIWC Cayenne Field Campaign

AbstractHigh Ice Water Content (HIWC) regions above tropical mesoscale convective systems are investigated using data from the second collaboration of the High Altitude Ice Crystals and High Ice Water Content projects (HAIC-HIWC) based in Cayenne, French Guiana in 2015. Observations from in-situ cloud probes on the French Falcon 20 determine the microphysical and thermodynamic properties of such regions. Data from a 2-D stereo probe and precipitation imaging probe show how statistical distributions of ice crystal mass median diameter (MMD), ice water content (IWC), and total number concentration (Nt) for particles with maximum dimension (Dmax) > 55 μm vary with environmental conditions, temperature (T), and convective properties such as vertical velocity (w), MCS age, distance away from convective peak (L), and surface characteristics. IWC is significantly correlated with w, whereas MMD decreases and Nt increases with decreasing T consistent with aggregation, sedimentation and vapor deposition processes at lower altitudes. MMD typically increases with IWC when IWC < 0.5 g m-3, but decreases with IWC when IWC > 0.5 g m-3 for -15 °C ≤ T ≤ -5 °C. Trends also depend on environmental conditions, such as presence of convective updrafts that are the ice crystal source, MMD being larger in older MCSs consistent with aggregation and less injection of small crystals into anvils, and IWCs decrease with increasing L at lower T. The relationship between IWC and MMD depends on environmental conditions, with correlations decreasing with decreasing T. The strength of correlation between IWC and Nt increases as T decreases.

Evaluation of the AROME model's ability to represent ice crystal icing using in situ observations from the HAIC 2015 field campaign

AbstractSince pilots generally avoid intense convective areas, ice crystals icing (ICI) is an aeronautical weather incident that mainly occurs in the anvil of tropical deep convective clouds. Samples of favorable conditions for the occurrence of ICI and data from the High Altitude Ice Crystals (HAIC) 2015 field campaign in French Guiana are investigated and compared with simulations of the French operational mesoscale forecast system Application of Research to Operations at Mesoscales (AROME). To this end, a contextualization of convective systems into convective, stratiform, and cirriform regions is employed for both observations and AROME. General features of the microphysics of deep tropical convective systems are identified. The number concentration of crystals larger than 125 μm and total water content (TWC) are strongly correlated at each temperature level, and both decrease with increasing distance from convective cores. AROME can reproduce the general behavior of the observed microphysics, especially TWC, but seems unable to simulate extreme ICI events. Reasons are sought in the assumptions performed in the microphysical scheme ICE3, and guidelines are proposed to enhance its skills in the context of ICI. In particular, the representation of the snow particle size distribution is adjusted across observations using a generalized gamma shape. This shape is found to outperform the usual Marshall–Palmer and gamma shapes. Additionally, a temperature and snow content dependence of generalized gamma parameters is found. These changes are found to significantly improve the snow concentration diagnostic of ICE3, and these modifications open the way for improvements in the ICE3 scheme.

Microphysical processes producing high ice water contents (HIWCs) in tropical convective clouds during the HAIC-HIWC field campaign: evaluation of simulations using bulk microphysical schemes

Abstract. Regions with high ice water content (HIWC), composed of mainly small ice crystals, frequently occur over convective clouds in the tropics. Such regions can have median mass diameters (MMDs) <300 µm and equivalent radar reflectivities <20 dBZ. To explore formation mechanisms for these HIWCs, high-resolution simulations of tropical convective clouds observed on 26 May 2015 during the High Altitude Ice Crystals – High Ice Water Content (HAIC-HIWC) international field campaign based out of Cayenne, French Guiana, are conducted using the Weather Research and Forecasting (WRF) model with four different bulk microphysics schemes: the WRF single‐moment 6‐class microphysics scheme (WSM6), the Morrison scheme, and the Predicted Particle Properties (P3) scheme with one- and two-ice options. The simulations are evaluated against data from airborne radar and multiple cloud microphysics probes installed on the French Falcon 20 and Canadian National Research Council (NRC) Convair 580 sampling clouds at different heights. WRF simulations with different microphysics schemes generally reproduce the vertical profiles of temperature, dew-point temperature, and winds during this event compared with radiosonde data, and the coverage and evolution of this tropical convective system compared to satellite retrievals. All of the simulations overestimate the intensity and spatial extent of radar reflectivity by over 30 % above the melting layer compared to the airborne X-band radar reflectivity data. They also miss the peak of the observed ice number distribution function for 0.1<Dmax<1 mm. Even though the P3 scheme has a very different approach representing ice, it does not produce greatly different total condensed water content or better comparison to other observations in this tropical convective system. Mixed-phase microphysical processes at −10 ∘C are associated with the overprediction of liquid water content in the simulations with the Morrison and P3 schemes. The ice water content at −10 ∘C increases mainly due to the collection of liquid water by ice particles, which does not increase ice particle number but increases the mass/size of ice particles and contributes to greater simulated radar reflectivity.

Comparisons of Cloud In Situ Microphysical Properties of Deep Convective Clouds to Appendix D/P Using Data from the High-Altitude Ice Crystals-High Ice Water Content and High Ice Water Content-RADAR I Flight Campaigns

Comparisons of Cloud In Situ Microphysical Properties of Deep Convective Clouds to Appendix D/P Using Data from the High-Altitude Ice Crystals-High Ice Water Content and High Ice Water Content-RADAR I Flight Campaigns

A thermodynamic model for ice crystal accretion in aircraft engines: EMM-C

Ingestion of high-altitude ice particles can be hazardous to aircraft engines in flight. Ice accretion may occur in the compression system, leading to blockage of the core gas path, blade damage and/or engine flameout. Numerous engine powerloss events since 1990 have been attributed to this mechanism. An expansion in engine certification requirements to incorporate ice crystal conditions has spurred efforts to develop analytical models for phenomenon. A necessary component of a complete icing simulation is a thermodynamic accretion model. In this paper, a new model for ice crystal icing is developed through adaptation of the Extended Messinger Model (EMM) from supercooled water conditions to mixed phase conditions (ice crystals and supercooled water). The new model is termed the Extended Messinger Model -Crystals (EMM-C). Continuity and energy balances are performed using the local flow conditions and impinging mass fluxes of ice and water. For accretions grown on a warm engine surface, initially above freezing temperature, a novel three-layer (water-ice-water) accretion structure is proposed, and the underlying equations described. Finally, predictions from the EMM-C model are compared to experimental results generated in an ice crystal wind tunnel. Compared with current literature icing models, EMM-C provides a tractable framework for modelling engine icing through allowing mixed phase accretions on warm surfaces.