Ice Number Research Articles

Revisiting Contrail Ice Formation: Impact of Primary Soot Particle Sizes and Contribution of Volatile Particles.

Aircraft contrails, formed largely on soot particles in current flights, are important for aviation's non-CO2 climate impact. Here we show that the activation of nonvolatile soot particles during contrail formation is likely determined by the sizes of primary soot particles rather than the effective sizes of soot aggregates as assumed in previous studies, which can explain less-than-unity fractions of soot particles forming contrail ice particles as recently observed during ECLIF (Emission and CLimate Impact of alternative Fuels) campaigns. The smaller soot primary sizes compared to aggregate sizes delay the onset of contrail ice formation, increase the maximum plume supersaturation reached in the contrail plume, and thus increase the probability of small volatile particles contributing to the total contrail ice particle number. This study suggests that the range of conditions for volatile plume particles to contribute significantly to the contrail ice number budget is wider than previously thought. As the aviation industry is moving toward sustainable aviation fuel and/or lean-burning engine technology, which is expected to reduce not only the emission index of nonvolatile soot particles but also the sizes of primary soot particles, this study highlights the need to better understand how the combined changes may affect contrail formation, contribution of volatile particles, and climate impacts.

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Pathways of ice multiplication in nimbostratus clouds during the Indian summer monsoon

The present study illustrates the microphysical parameters of nimbostratus clouds during the Indian summer monsoon using airborne and radar observations conducted as part of the Cloud Aerosol Interaction and Precipitation Enhancement Experiment (CAIPEEX), and simulations with the Weather Research and Forecasting (WRF) model. The study also analyzes the impact of secondary ice production (emphasizing the microphysical pathways of Hallett-Mossop (HM) as represented in the models) on cloud processes using model sensitivity simulations. The effect of possible pathways of riming during the HM process is the focus of the sensitivity simulations. Aircraft observations showed higher ice particle concentrations in the Hallett–Mossop zone (−3 °C to - 8 °C) along with the existence of smaller and larger cloud droplets, rimed needles, columns, and snow particles. Observations strongly suggested the active presence of the HM process in this cloud. The observed mean values of microphysical parameters including liquid water content, ice number concentrations, and maximum vertical velocity, agreed well with model simulations. The number concentration and water content of ice crystals and snow decreased in the HM zone during the HM inactive case. A reduction in rainfall is also observed in the absence of HM. HM increases convection as it causes a sudden increase in the number concentration of small ice particles, which causes rainwater to freeze quickly and water vapour to be consumed by diffusional growth. Latent heating is produced by both effects, which energizes the convection.

Differences in microphysical properties of cirrus at high and mid-latitudes

Abstract. Despite their proven importance for the atmospheric radiative energy budget, the effect of cirrus on climate and the magnitude of their modification by human activity is not well quantified. Besides anthropogenic pollution sources on the ground, aviation has a large local effect on cirrus microphysical and radiative properties via the formation of contrails and their transition to contrail cirrus. To investigate the anthropogenic influence on natural cirrus, we compare the microphysical properties of cirrus measured at mid-latitude (ML) regions (<60∘ N) that are often affected by aviation and pollution with cirrus measured in the same season in comparatively pristine high latitudes (HLs; ≥60∘ N). The number concentration, effective diameter, and ice water content of the observed cirrus are derived from in situ measurements covering ice crystal sizes between 2 and 6400 µm collected during the CIRRUS-HL campaign (Cirrus in High Latitudes) in June and July 2021. We analyse the dependence of cirrus microphysical properties on altitude and latitude and demonstrate that the median ice number concentration is an order of magnitude larger in the measured mid-latitude cirrus, with 0.0086 cm−3, compared to the high-latitude cirrus, with 0.001 cm−3. Ice crystals in mid-latitude cirrus are on average smaller than in high-latitude cirrus, with a median effective diameter of 165 µm compared to 210 µm, and the median ice water content in mid-latitude cirrus is higher (0.0033 g m−3) than in high-latitude cirrus (0.0019 g m−3). In order to investigate the cirrus properties in relation to the region of formation, we combine the airborne observations with 10 d backward trajectories to identify the location of cirrus formation and the cirrus type, i.e. in situ or liquid origin cirrus, depending on whether there is only ice or also liquid water present in the cirrus history, respectively. The cirrus formed and measured at mid-latitudes (M–M) have a particularly high ice number concentration and low effective diameter. This is very likely a signature of contrails and contrail cirrus, which is often observed in the in situ origin cirrus type. In contrast, the largest effective diameter and lowest number concentration were found in the cirrus formed and measured at high latitudes (H–H) along with the highest relative humidity over ice (RHi). On average, in-cloud RHi was above saturation in all cirrus. While most of the H–H cirrus were of an in situ origin, the cirrus formed at mid-latitudes and measured at high latitudes (M–H) were mainly of liquid origin. A pristine Arctic background atmosphere with relatively low ice nuclei availability and the extended growth of few nucleated ice crystals may explain the observed RHi and size distributions. The M–H cirrus are a mixture of the properties of M–M and H–H cirrus (preserving some of the initial properties acquired at mid-latitudes and transforming under Arctic atmospheric conditions). Our analyses indicate that part of the cirrus found at high latitudes is actually formed at mid-latitudes and therefore affected by mid-latitude air masses, which have a greater anthropogenic influence.

Arctic Cloud‐Base Ice Precipitation Properties Retrieved Using Bayesian Inference

AbstractCloud‐climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long‐lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This Arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at Arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud‐base properties a more appealing target for inference and evaluation of model simulations. Here I describe an inverse model for the estimation of cloud base ice precipitation properties over Utqiagvik, North Slope of Alaska, using the synthesis of ground‐based radar and lidar measurements. By leveraging a Markov Chain Monte Carlo algorithm as the core of the inverse model, a wide range of particle size distributions are sampled, and different combinations of ice habit models are examined, both of which are typically fixed in other retrieval methods. Results show intriguing links between different cloud base thermodynamic and ice precipitation properties. Apparent ice number concentration enhancements at temperatures of −5 and −15°C suggest possible secondary ice production (SIP). The analysis alludes to an overestimation of SIP occurrence and intensity, especially in studies relying only on radar or lidar measurements. Finally, reflectivity‐dependent ice precipitation rate and ice water content parameterizations are presented.

Investigating the vertical extent and short-wave radiative effects of the ice phase in Arctic summertime low-level clouds

Abstract. Low-level (cloud tops below 2 km) mixed-phase clouds are important in amplifying warming in the Arctic region through positive feedback in cloud fraction, water content and phase. In order to understand the cloud feedbacks in the Arctic region, good knowledge of the vertical distribution of the cloud water content, particle size and phase is required. Here we investigate the vertical extent of the cloud-phase and ice-phase optical properties in six case studies measured in the European Arctic during the ACLOUD campaign. Late spring- and summertime stratiform clouds were sampled in situ over pack ice, marginal sea ice zone and open-ocean surface, with cloud top temperatures varying between −15 and −1.5 ∘C. The results show that, although the liquid phase dominates the upper parts of the clouds, the ice phase was frequently observed in the lower parts down to cloud top temperatures as warm as −3.8 ∘C. In the studied vertical cloud profiles, the maximum of average liquid phase microphysical properties, droplet number concentration, effective radius and liquid water content, varied between 23 and 152 cm−3, 19 and 26 µm, 0.09 and 0.63 g m−3, respectively. The maximum of average ice-phase microphysical properties varied between 0.1 and 57 L−1 for the ice number concentration, 40 and 70 µm for the effective radius, and 0.005 and 0.08 g m−3 for the ice water content. The elevated ice crystal number concentrations and ice water paths observed for clouds, with cloud top temperatures between −3.8 and −8.7 ∘C can be likely attributed to secondary ice production through rime splintering. Low asymmetry parameters between 0.69 and 0.76 were measured for the mixed-phase ice crystals with a mean value of 0.72. The effect of the ice-phase optical properties on the radiative transfer calculations was investigated for the four cloud cases potentially affected by secondary ice production. Generally the choice of ice-phase optical properties only has a minor effect on the cloud transmissivity and albedo, except in a case where the ice phase dominated the upper cloud layer extinction. In this case, cloud albedo at solar wavelengths was increased by 10 % when the ice phase was given its measured optical properties instead of treating it as liquid phase. The presented results highlight the importance of accurate vertical information on cloud phase for radiative transfer and provide a suitable data set for testing microphysical parameterizations in models.

Orographic Cirrus and Its Radiative Forcing in NCAR CAM6

AbstractCirrus clouds play an important role in the Earth's radiative energy budget, thereby affecting the climate state and climate change. Orographic gravity wave (OGW)‐induced sub‐grid scale vertical velocity (i.e., cooling rate) is not resolved by large‐scale models and its impact on ice formation in cirrus clouds is not well quantified. In this study, one sub‐grid scale OGW scheme (e.g., McFarlane) is used in the Community Atmosphere Model version 6 (CAM6) to generate vertical velocity variance (σw) for cirrus formation. Results from the default model and simulations with the OGW‐induced σw are evaluated against the DOE ARM Small Particles in Cirrus (SPARTICUS) campaign observations. The OGW based on the McFarlane scheme increases the sub‐grid scale σw over mountains compared to the default model and improves the model agreement with the SPARTICUS observations. Larger σw due to OGWs can trigger more frequent homogeneous nucleation in orographic cirrus and generates a higher number concentration of ice crystals observed during the SPARTICUS campaign. Moreover, our evaluation of the model simulations against satellite observations indicates that the McFarlane scheme generates high in‐cloud ice number concentrations (>200 L−1) in the upper troposphere over mountains and high plateaus at mid‐ and high‐latitudes of the winter hemisphere as shown in the observations. More ice crystals with smaller sizes absorb more infrared radiation (+0.523 ± 0.125 W m−2). The net radiative cloud forcing change at the top of the atmosphere is +0.330 W m−2 due to the orographic cirrus.

Importance of ice nucleation and precipitation on climate with the Parameterization of Unified Microphysics Across Scales version 1 (PUMASv1)

Abstract. Cloud microphysics is critical for weather and climate prediction. In this work, we document updates and corrections to the cloud microphysical scheme used in the Community Earth System Model (CESM) and other models. These updates include a new nomenclature for the scheme, now called Parameterization of Unified Microphysics Across Scales (PUMAS), and the ability to run the scheme on graphics processing units (GPUs). The main science changes include refactoring an ice number limiter and associated changes to ice nucleation, adding vapor deposition onto snow, and introducing an implicit numerical treatment for sedimentation. We also detail the improvements in computational performance that can be achieved with GPU acceleration. We then show the impact of these scheme changes on the (a) mean state climate, (b) cloud feedback response to warming, and (c) aerosol forcing. We find that corrections are needed to the immersion freezing parameterization and that ice nucleation has important impacts on climate. We also find that the revised scheme produces less cloud liquid and ice but that this can be adjusted by changing the loss process for cloud liquid (autoconversion). Furthermore, there are few discernible effects of the PUMAS changes on cloud feedbacks but some reductions in the magnitude of aerosol–cloud interactions (ACIs). Small cloud feedback changes appear to be related to the implicit sedimentation scheme, with a number of factors affecting ACIs.

Extremely High Concentrations of Ice Particles in East Asian Dust‐Infused Baroclinic Storm (DIBS) Cirrus Shield: Dominant Role of Dust Ice Nucleation Effect

AbstractActing as efficient ice‐nucleating particles, dust particles can affect cirrus cloud properties through heterogeneous nucleation. A dust‐infused baroclinic storm (DIBS) over East Asia in May 2017 was examined using satellite observations and the Weather Research and Forecasting model coupled with chemistry (WRF‐Chem) to study the dust effects on cirrus clouds. Satellite observations showed that the ice number concentration (Nice) was extremely high, exceeding 1–10 cm−3 at the cirrus shield of the DIBS, and the ice effective radius (rei) was about 10–30 μm. An aerosol‐aware heterogeneous ice nucleation parameterization (Ullrich et al., 2017, https://doi.org/10.1175/jas-d-16-0074.1) was implemented into WRF‐Chem to explore the formation mechanism of such a high Nice. The simulated Nice from Ullrich17 was 1–2 orders of magnitude higher than that from the default scheme (Cooper, 1986, https://doi.org/10.1007/978-1-935704-17-1_4), and the rei was 2–3 times smaller. The results from Ullrich17 were more consistent with satellite observations, which can be attributed to the enhanced deposition nucleation from the high dust concentrations. Our results demonstrated that the extremely high Nice could be fully explained by the heterogeneous nucleation of dust particles in DIBS. A conceptual framework is proposed to explain the dust‐cloud interactions in DIBS. This study confirms the significant dust effects on cirrus clouds in DIBS and highlights the importance of using aerosol‐aware ice nucleation parameterizations in modeling cirrus clouds when the dust is present.

Sensitivity of convectively driven tropical tropopause cirrus properties to ice habits in high-resolution simulations

Abstract. Cirrus clouds that form in the tropical tropopause layer (TTL) can play a key role in vertical transport through the upper troposphere and lower stratosphere, which can significantly impact the radiative energy budget and stratospheric chemistry. However, the lack of realistic representation of natural ice cloud habits in microphysical parameterizations can lead to uncertainties in cloud-related processes and cloud–climate feedbacks. The main goal of this study is to investigate the role of different cloud regimes and the associated ice habits in regulating the properties of the TTL. We compare aircraft measurements from the StratoClim field campaign to a set of numerical experiments at the scale of large-eddy simulations (LESs) for the same case study that employ different microphysics schemes. Aircraft measurements over the southern slopes of the Himalayas captured high ice water content (HIWC) up to 2400 ppmv and ice particle aggregates exceeding 700 µm in size with unusually long residence times. The observed ice particles were mainly of liquid origin, with a small amount formed in situ. The corresponding profile of ice water content (IWC) from the ERA5 reanalysis corroborates the presence of HIWC detrained from deep-convective plumes in the TTL but underestimates HIWC by an order of magnitude. In the TTL, only the scheme that predicts ice habits can reproduce the observed HIWC, ice number concentration, and bimodal ice particle size distribution. The lower range of particle sizes is mostly represented by planar and columnar habits, while the upper range is dominated by aggregates. Large aggregates with sizes between 600 and 800 µm have fall speeds of less than 20 cm s−1, which explains the long residence time of the aggregates in the TTL. Planar ice particles of liquid origin contribute substantially to HIWC. The columnar and aggregate habits are in the in situ range with lower IWC and number concentrations. For all habits, the ice number concentration increases with decreasing temperature. For the planar ice habit, relative humidity is inversely correlated with fall speed. This correlation is less evident for the other two ice habits. In the lower range of supersaturation with respect to ice, the columnar habit has the highest fall speed. The difference in ice number concentration across habits can be up to 4 orders of magnitude, with aggregates occurring in much smaller numbers. We demonstrate and quantify the linear relationship between the differential sedimentation of pristine ice crystals and the size of the aggregates that form when pristine crystals collide. The slope of this relationship depends on which pristine ice habit sediments faster. Each simulated ice habit is associated with distinct radiative and latent heating rates. This study suggests that a model configuration nested down to LES scales with a microphysical parameterization that predicts ice shape evolution is crucial to provide an accurate representation of the microphysical properties of TTL cirrus and thus the associated (de)hydration process.

Important Ice Processes Are Missed by the Community Earth System Model in Southern Ocean Mixed‐Phase Clouds: Bridging SOCRATES Observations to Model Developments

AbstractGlobal climate models (GCMs) are challenged by difficulties in simulating cloud phase and cloud radiative effect over the Southern Ocean (SO). Some of the new‐generation GCMs predict too much liquid and too little ice in mixed‐phase clouds. This misrepresentation of cloud phase in GCMs results in weaker negative cloud feedback over the SO and a higher climate sensitivity. Based on a model comparison with observational data obtained during the Southern Ocean Cloud Radiation and Aerosol Transport Experimental Study, this study addresses a key uncertainty in the Community Earth System Model version 2 (CESM2) related to cloud phase, namely ice formation in pristine remote SO clouds. It is found that sea spray organic aerosols (SSOAs) are the most important type of ice nucleating particles (INPs) over the SO with concentrations 1 order of magnitude higher than those of dust INPs based on measurements and CESM2 simulations. Secondary ice production (SIP) which includes riming splintering, rain droplet shattering, and ice‐ice collisional fragmentation as implemented in CESM2 is the dominant ice production process in moderately cold clouds with cloud temperatures greater than −20°C. SIP enhances the in‐cloud ice number concentrations (Ni) by 1–3 orders of magnitude and predicts more mixed‐phase (with percentage occurrence increased from 15% to 21%), in better agreement with the observations. This study highlights the importance of accurately representing the cloud phase over the pristine remote SO by considering the ice nucleation of SSOA and SIP processes, which are currently missing in most GCM cloud microphysics parameterizations.

Observations of microphysical properties and radiative effects of a contrail cirrus outbreak over the North Atlantic

Abstract. Contrail cirrus constitute the largest radiative forcing (RF) component to the total aviation effect on climate. However, the microphysical properties and radiative effects of contrail cirrus and natural cirrus clouds in the same meteorological conditions are still not completely resolved. Motivated by these uncertainties, we investigate an extended cirrus region perturbed by aviation in the North Atlantic region (NAR) on 26 March 2014 during the Midlatitude Cirrus (ML-CIRRUS) experiment. On that day, high air traffic density in the NAR combined with large scale cold and humid ambient conditions favored the formation of a contrail cirrus outbreak situation. In addition, low coverage by low-level water clouds and the homogeneous oceanic albedo increased the sensitivity for retrieving cirrus properties and their radiative effect from satellite remote sensing. This allowed us to extend the current knowledge on contrail cirrus by combining airborne in situ, lidar and satellite observations. In the synoptic context of a ridge cirrus, an extended thin ice cloud with many persistent contrails and contrail cirrus has been observed for many hours with the geostationary Meteosat Second Generation (MSG)/Spinning Enhanced Visible and InfraRed Imager (SEVIRI) from the early morning hours until dissipation after noon. Airborne lidar observations aboard the German High Altitude and LOng Range Research Aircraft (HALO) suggest that this cirrus has a significant anthropogenic contribution from aviation. A new method based on in situ measurements was used to distinguish between contrails, contrail cirrus and natural cirrus based on ice number and gas phase NO concentrations. Results show that contrail effective radii (Reff) reach at most 11 µm, while contrail cirrus Reff can be as large as 51 µm. Contrail and contrail cirrus mean Reff is 18 % smaller than that of natural cirrus. We find that a difference in Reff between contrail cirrus and natural cirrus survives in this contrail cirrus outbreak event. As for radiative effects, a new method to estimate top-of-atmosphere instantaneous RF in the solar and thermal range is developed based on radiative transfer model simulations exploiting in situ and lidar measurements, satellite observations and ERA5 reanalysis data for both cirrus and cirrus-free regions. Broadband irradiances estimated from our simulations compare well with satellite observations from MSG, indicating that our method provides a good representation of the real atmosphere and can thus be used to determine the RF of ice clouds. For a larger spatial area around the flight path, we find that the contrail cirrus outbreak is warming in the early morning and cooling during the day. The methods presented here and the results will be valuable for future research to constrain uncertainties in the assessment of radiative impacts of contrail cirrus and natural cirrus and for the formulation and evaluation of contrail mitigation options.

Implementation of Prognostic Cloud Ice Number Concentrations for the Weather Research and Forecasting (WRF) Double‐Moment 6‐Class (WDM6) Microphysics Scheme

AbstractThe ice microphysical processes in the Weather Research and Forecasting (WRF) Double‐Moment 6‐class (WDM6) microphysics scheme are treated as a single‐moment approach, in which the number concentration of cloud ice is diagnosed based on its mixing ratio. This study develops the revised WDM6 scheme through the implementation of prognostic cloud ice number concentrations. The effect of the prognostic number concentration on the simulated precipitation is verified through simulations of short‐term winter snowfall cases during International Collaborative Experiments for the Pyeongchang 2018 Olympics and Paralympics (ICE‐POP 2018) winter games and a 1‐month regional climate case during the summer season, July 2009. For all cases, the revised WDM6 simulates higher cloud ice number concentrations and lower cloud ice mixing ratios than the original WDM6. The microphysics budget analysis for the snowfall cases shows that the inefficient deposition and vapor freezing nucleation processes of cloud ice reduce the available cloud ice mixing ratio. Consequently, the accretion processes with cloud ice decrease and the deposition into snow increases due to the surplus water vapor. The revised WDM6 alleviates the positive bias of surface precipitation consisting of snow over the region where the original WDM6 simulates excessive precipitation, compared to the observed data. For the regional climate case, the reduced cloud ice amount strengthens the Western North Pacific high‐pressure system by allowing more solar radiation to reach the surface, leading to simulated precipitation bands and synoptic environments that are more comparable with the observed data.

The Role of Mineral Dust Aerosol Particles in Aviation Soot‐Cirrus Interactions

AbstractPredicting cirrus cloud properties with confidence requires a sound understanding of the relative roles of homogeneous and heterogeneous ice formation. This study explores the effect of mineral dust and contrail‐processed aviation soot particles as ice‐nucleating particles (INPs) competing with liquid solution droplets in cirrus formation. We study aerosol‐cirrus interactions by accounting for atmospheric variability in updraft speeds and INP number concentrations. Our results confirm the dominant role of mineral dust in ice nucleation events in cirrus clouds. In addition, we show that pre‐existing thin cirrus may suppress ice formation when updraft speeds are low. We find that homogeneous freezing of liquid solution droplets dominates clear‐sky aerosol‐cirrus interaction above a threshold updraft speed determined by total number concentrations and ice nucleation abilities of INPs. When mineral dust particles exceed number concentrations of 10 L−1, they reduce homogeneously nucleated ice crystal numbers significantly and even prevent homogeneous freezing for frequently observed local updraft speeds between 10 and 20 cm s−1. When both mineral dust and aviation soot particles coexist with solution droplets, dust typically prevents ice nucleation by aviation soot. Aviation soot exerts a notable impact on cirrus ice numbers only if updrafts are weak, large soot particles are present in number concentrations that are considerably higher than typically observed in emission measurements, and/or number concentrations of mineral dust and other INPs are low. Overall, our results elucidate the role of aviation soot‐cirrus interactions in the presence of other INP types.

Improvements of the Double-Moment Bulk Cloud Microphysics Scheme in the Nonhydrostatic Icosahedral Atmospheric Model (NICAM)

Abstract This study revises the collisional growth, heterogeneous ice nucleation, and homogeneous ice nucleation processes in a double-moment bulk cloud microphysics scheme implemented in the Nonhydrostatic Icosahedral Atmospheric Model (NICAM). The revised cloud microphysical processes are tested by 10-day global simulations with a horizontal resolution of 14 km. It is found that both the aggregation of cloud ice with smaller diameters and the graupel production by riming are overestimated in the current schemes. A new method that numerically integrates the collection kernel solves this issue, and consequently, the lifetime of cloud ice is reasonably extended in reference to satellite observations. In addition, the results indicate that a reduction in graupel modulates the convective intensity, particularly in intense rainfall systems. The revision of both heterogeneous and homogeneous ice nucleation significantly increases the production rate of cloud ice number concentration. With these revisions, the new version of the cloud microphysics scheme successfully improves outgoing longwave radiation, particularly over the intertropical convergence zone, in reference to satellite observations. Therefore, the revisions are beneficial for both long-term climate simulations and representing the structure of severe storms. Significance Statement Very high-resolution global atmospheric models have been developed to simultaneously address global climate and regional weather. In general, cloud microphysics schemes used in such global models are introduced from regional weather forecasting models to realistically represent mesoscale cloud systems. However, a cloud microphysics scheme that was originally developed with the aim of weather forecasting can cause unexpected errors in global climate simulations because such a cloud microphysics scheme is not designed for interdisciplinary usage across spatiotemporal scales. This study focuses on systematic model biases in evaluating the terminal velocity of ice cloud particles and proposes a method to accurately calculate the growth rate of ice cloud particles. Improvements in ice cloud modeling successfully reduce model biases in the global energy budget. In addition, the internal structure of intense rainfall systems is modified using the new cloud model. Therefore, improvements in ice cloud modeling could further increase the reliability of weather forecasting, seasonal prediction, and climate projection.

Code for JAMES publication - Implementation of prognostic cloud ice number concentrations for the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme.

Code for JAMES publication - Implementation of prognostic cloud ice number concentrations for the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme.

The influence of multiple groups of biological ice nucleating particles on microphysical properties of mixed-phase clouds observed during MC3E

Abstract. A new empirical parameterization (EP) for multiple groups of primary biological aerosol particles (PBAPs) is implemented in the aerosol–cloud model (AC) to investigate their roles as ice nucleating particles (INPs). The EP describes the heterogeneous ice nucleation by (1) fungal spores, (2) bacteria, (3) pollen, (4) detritus of plants, animals, and viruses, and (5) algae. Each group includes fragments from the originally emitted particles. A high-resolution simulation of a midlatitude mesoscale squall line by AC is validated against airborne and ground observations. Sensitivity tests are carried out by varying the initial vertical profiles of the loadings of individual PBAP groups. The resulting changes in warm and ice cloud microphysical parameters are investigated. The changes in warm microphysical parameters, including liquid water content and cloud droplet number concentration, are minimal (<10 %). Overall, PBAPs have little effect on the ice number concentration (<6 %) in the convective region. In the stratiform region, increasing the initial PBAP loadings by a factor of 1000 resulted in less than 40 % change in ice number concentrations. The total ice concentration is mostly controlled by various mechanisms of secondary ice production (SIP). However, when SIP is intentionally shut down in sensitivity tests, increasing the PBAP loading by a factor of 100 has an effect of less than 3 % on the ice phase. Further sensitivity tests revealed that PBAPs have little effect on surface precipitation and on the shortwave and longwave flux (<4 %) for a 100-fold perturbation in PBAPs.

Evidence for Secondary Ice Production in Southern Ocean Maritime Boundary Layer Clouds

AbstractMaritime boundary‐layer clouds over the Southern Ocean (SO) have a large shortwave radiative effect. Yet, climate models have difficulties in representing these clouds and, especially, their phase in this observationally sparse region. This study aims to increase the knowledge of SO cloud phase by presenting in‐situ cloud microphysical observations from the Southern Ocean Clouds, Radiation, Aerosol, Transport Experimental Study (SOCRATES). We investigate the occurrence of ice in summertime marine stratocumulus and cumulus clouds in the temperature range between 6 and −25°C. Our observations show that in ice‐containing clouds, maximum ice number concentrations of up to several hundreds per liter were found. The observed ice crystal concentrations were on average one to two orders of magnitude higher than the simultaneously measured ice nucleating particle (INP) concentrations in the temperature range below −10°C and up to five orders of magnitude higher than estimated INP concentrations in the temperature range above −10°C. These results highlight the importance of secondary ice production (SIP) in SO summertime marine boundary‐layer clouds. Evidence for rime splintering was found in the Hallett‐Mossop (HM) temperature range but the exact SIP mechanism active at lower temperatures remains unclear. Finally, instrument simulators were used to assess simulated co‐located cloud ice concentrations and the role of modeled HM rime‐splintering. We found that CAM6 is deficient in simulating number concentrations across the HM temperature range with little sensitivity to the model HM process, which is inconsistent with the aforementioned observational evidence of highly active SIP processes in SO low‐level clouds.

Code for JAMES publication - Implementation of prognostic cloud ice number concentrations for the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme.

Code for JAMES publication - Implementation of prognostic cloud ice number concentrations for the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme.