Ice Crystal Number Research Articles

Investigating secondary ice production in a deep convective cloud with a 3D bin microphysics model: Part I - Sensitivity study of microphysical processes representations

Secondary ice production (SIP) is a crucial phenomenon for explaining the formation of ice crystal clouds, especially when addressing the discrepancies between observed ice crystal number concentrations and ice nucleating particles (INPs). In this study, we investigate parameterizations of three SIP processes (Hallett-Mossop, fragmentation of freezing drops, and fragmentation due to ice–ice collision) by simulating a deep convective cloud observed during the HAIC/HIWC campaign with the 3D bin microphysics scheme DESCAM (DEtailed SCAvening and Microphysics model). The simulated mean cloud properties, including particle size distributions and ice crystal number concentration are compared with in situ probe observations obtained during the campaign. Simulation excluding SIP show a large underestimation of small ice crystals (< 1 mm diameter) for temperatures warmer than ‐30∘C. In our results, incorporating Hallett-Mossop and fragmentation due to ice–ice collision processes leads to ice crystal number concentrations close to observed values, thereby reducing discrepancies by two orders of magnitude. Our simulations also indicates that fragmentation of freezing drops affect minimally the properties of the cloud at its mature stage. Furthermore, we investigate the impact of fragments sizes resulting from SIP processes and show that the size of fragments generated from fragmentation due to ice–ice collision significantly influences the shape of ice particle size distribution. Employing various parameterizations of the ice crystal sticking efficiency reveals a notable impact on cloud properties. This study shows that SIP mechanisms are important and have to be considered for cold and mixed-phase clouds. However their parameterization lack reliability, highlighting the need for better quantifying these mechanisms. The companion paper, investigates the effects of SIP processes on the formation and the evolution of the deep convective system.

Impact of Asian Dust on Cirrus Formation Over the Central Pacific: CALIOP‐ and CloudSat‐Observation‐Based Case Studies

AbstractCirrus clouds are of great importance to the global climate, with their net radiative forcing strongly dependent on the microphysical properties that are related to the ice‐nucleating regime. However, the influence of long‐range transport of dust on primary ice formation in cirrus clouds is limitedly understood, specifically over the clean remote ocean regions. Here, two case studies show that transpacific Asian dust can impact the ice formation of cirrus clouds over the central Pacific based on Cloud‐Aerosol Lidar with Orthogonal Polarization and Cloud Profiling Radar (CPR, CloudSat) observations. One case shows a well‐developed horizontally extended cirrus embedded in a pure dust layer, with an average dust‐related ice‐nucleating particle concentration (INPC) of 7 L−1 and 96 L−1 for an ice saturation ratio Si of 1.15 and 1.25, respectively; ice crystal number concentration (ICNC) with diameters &gt;25 and 100 μm (denoted as nice,25 μm and nice,100 μm) are 64 L−1 and 7 L−1, respectively. Another case shows that cirrus clouds with a much smaller horizontal extent appeared in the vicinity of polluted dust, with an average INPC of 42–310 L−1 for the typical higher Si of 1.25–1.35 by considering a tenfold reduction of the ice nucleation efficiency of ice crystals; nice,25 μm and nice,100 μm are 168 L−1 and 20 L−1, respectively. The estimated INPC and ICNC values suggest the dominance of ice formation by dust‐induced heterogeneous nucleation, proving that the long‐range transport of dust toward the upper troposphere and the potential influence on cirrus formation over the central Pacific should be well considered in atmospheric models.

Repurposing weather modification for cloud research showcased by ice crystal growth.

The representation of cloud processes in models is one of the largest sources of uncertainty in weather forecast and climate projections. While laboratory settings offer controlled conditions for studying cloud processes, they cannot reproduce the full range of conditions and interactions present in natural cloud systems. To bridge this gap, here we leverage weather modification, specifically glaciogenic cloud seeding, to investigate ice growth rates within natural clouds. Seeding experiments were conducted in supercooled stratus clouds (at to C) using an uncrewed aerial vehicle, and the created ice crystals were measured 4-10 min downwind by in situ and ground-based remote sensing instrumentation. We observed substantial variability in ice crystal growth rates within natural clouds, attributed to variations in ice crystal number concentrations and in the supersaturation, which is difficult to reproduce in the laboratory and which implies faster precipitation initiation than previously thought. We found that for the experiments conducted at C, the ice crystal populations grew nearly linearly during the time interval from 6 to 10 min. Our results demonstrate that the targeted use of weather modification techniques can be employed for fundamental cloud research (e.g. ice growth processes, aerosol-cloud interactions), helping to advance cloud microphysics parameterizations and to improve weather forecasts and climate projections.

In-situ aircraft observations of aerosol and cloud microphysical characteristics of mixed-phase clouds over the North China Plain

Aerosol-cloud interactions play a vital role in climate change. This study leverages observations from the King-350 aircraft over the North China Plain on November 29, 2019, to examine aerosol and cloud microphysical characteristics of mixed-phase clouds. Through detailed vertical and spectral distributions, we investigate aerosol, cloud droplet, and ice crystal distributions in seeded clouds (SC) and natural precipitating clouds (NPC) within the same cloud system. From the vertical profile, SC and NPC have similar vertical distributions of aerosol and cloud droplets, with over 95 % of aerosols concentrated below 1600 m near the ground. Cloud droplets are more evenly distributed within the two clouds, cloud droplet number concentrations (Nc) in SC were three orders of magnitude higher than in NPC. Ice water content (IWC) and ice crystal number concentration (Ni) show distinct layer preferences—accumulating predominantly in SC's top layer and NPC's middle layer. From spectral distribution, a smaller proportion of cloud droplets (40–50 μm in diameter, the same hereafter) in SC compared to NPC. Rimed ice crystals and globular graupel (1325–1550 μm in diameter) were in SC, while plate and irregular ice crystals (300–450 μm) were in NPC with an order of magnitude higher than in SC. These microphysical differences highlight the complexity of cloud seeding efficacy, which varies based on cloud conditions and microphysical properties. In the first seeding, Ni increased by 1–2 orders of magnitude (125–300 μm) in the high Nc (Nc > 1.11 × 105 L−1) region. Seeding in low Nc (Nc < 1.11 × 105 L−1) regions was hard to be effective, especially in low Nc and low liquid water content (LWC) (LWC < 0.122 g/m3) regions. In the second seeding, ice crystals (125–250 μm) produced by the first seeding enhance the seeding efficiency. The responded regions were more sensitive to subsequent seeding, resulting in stronger reactions or longer duration.

Simulations of primary and secondary ice production during an Arctic mixed-phase cloud case from the Ny-Ålesund Aerosol Cloud Experiment (NASCENT) campaign

Abstract. The representation of Arctic clouds and their phase distributions, i.e., the amount of ice and supercooled water, influences predictions of future Arctic warming. Therefore, it is essential that cloud phase is correctly captured by models in order to accurately predict the future Arctic climate. Ice crystal formation in clouds happens through ice nucleation (primary ice production) and ice multiplication (secondary ice production). In common weather and climate models, rime splintering is the only secondary ice production process included. In addition, prescribed number concentrations of cloud condensation nuclei or cloud droplets and ice-nucleating particles are often overestimated in Arctic environments by standard model configurations. This can lead to a misrepresentation of the phase distribution and precipitation formation in Arctic mixed-phase clouds, with important implications for the Arctic surface energy budget. During the Ny-Ålesund Aerosol Cloud Experiment (NASCENT), a holographic probe mounted on a tethered balloon took in situ measurements of number and mass concentrations of ice crystals and cloud droplets in Svalbard, Norway, during fall 2019 and spring 2020. In this study, we choose one case study from this campaign that shows evidence of strong secondary ice production and use the Weather Research and Forecasting (WRF) model to simulate it at a high vertical and spatial resolution. We test the performance of different microphysical parametrizations and apply a new state-of-the-art secondary ice parametrization. We find that agreement with observations highly depends on the prescribed cloud condensation nuclei/cloud droplet and ice-nucleating particle concentrations and requires an enhancement of secondary ice production processes. Lowering mass mixing ratio thresholds for rime splintering inside the Morrison microphysics scheme is crucial to enable secondary ice production and thereby match observations for the right reasons. In our case, rime splintering is required to initiate collisional breakup. The simulated contribution from collisional breakup is larger than that from droplet shattering. Simulating ice production correctly for the right reasons is a prerequisite for reliable simulations of Arctic mixed-phase cloud responses to future temperature or aerosol perturbations.

The Preliminary Application of Spectral Microphysics in Numerical Study of the Effects of Aerosol Particles on Thunderstorm Development

Progress in numerical models and improved computational capabilities have significantly advanced our comprehension of how aerosol particles impact thunderstorm clouds. Yet, much of this research has focused on employing bulk microphysics models to explain the impacts of aerosol particles acting as cloud condensation nuclei (CCN) on electrical activities in thunderstorm clouds. The bulk thunderstorm models use mean sizes of particles and terminal-fall velocities. This causes calculation deviation in the electrification simulation, which in turn leads to deviations in the simulation of lightning processes. Developing this further, we established a three-dimensional high-resolution cloud–aerosol bin thunderstorm model with electrification and lightning to provide more accurate microphysics and dynamic fields for studying electrical activities. For evaluating the impacts of aerosol particles, specifically CCN, on the properties of continental thunderclouds, aerosols from both clean and polluted continental environments were selected. Cloud simulations indicate that droplets develop a narrower spectrum in polluted continental conditions, and weakened ice crystal growth increases the number of small ice crystals compared to clean conditions. Smaller droplets and ice crystals result in less effective riming and decreased graupel concentration and mass. Consequently, a significant decrease in large ice particles leads to a weakened process of charge separation under conditions of pollution. As a direct result, there is about a 43% reduction in lightning frequency and a delay of approximately 5 min in the lightning process under polluted conditions.

RaFSIP: Parameterizing Ice Multiplication in Models Using a Machine Learning Approach

AbstractAccurately representing mixed‐phase clouds (MPCs) in global climate models (GCMs) is critical for capturing climate sensitivity and Arctic amplification. Secondary ice production (SIP), can significantly increase ice crystal number concentration (ICNC) in MPCs, affecting cloud properties and processes. Here, we introduce a machine‐learning (ML) approach, called Random Forest SIP (RaFSIP), to parameterize SIP in stratiform MPCs. RaFSIP is trained on 16 grid points with 10‐km horizontal spacing derived from a 2‐year simulation with the Weather Research and Forecasting (WRF) model, including explicit SIP microphysics. Designed for a temperature range of 0 to −25°C, RaFSIP simplifies the description of rime splintering, ice‐ice collisional break‐up, and droplet‐shattering using only a limited set of inputs. RaFSIP was evaluated offline before being integrated into WRF, demonstrating its stable online performance in a 1‐year simulation keeping the same model setup as during training. Even when coupled with the 50‐km grid spacing domain of WRF, RaFSIP reproduces ICNC predictions within a factor of 3 when compared to simulations with explicit SIP microphysics. The coupled WRF‐RaFSIP scheme replicates regions of enhanced SIP and accurately maps ICNCs and liquid water content, particularly at temperatures above −10°C. Uncertainties in RaFSIP minimally impact surface cloud radiative forcing in the Arctic, resulting in radiative biases under 3 Wm−2 compared to simulations with detailed microphysics. Although the performance of RaFSIP in convective clouds remains untested, its adaptable nature allows for data set augmentation to address this aspect. This framework opens possibilities for GCM simplification and process description through physics‐guided ML algorithms.

Simulating the seeder–feeder impacts on cloud ice and precipitation over the Alps

Abstract. The ice phase impacts many cloud properties as well as cloud lifetime. Ice particles that sediment into a lower cloud from an upper cloud (external seeder–feeder process) or into the mixed-phase region of a deep cloud from cirrus levels (internal seeder–feeder process) can influence the ice phase of the lower cloud, amplify cloud glaciation and enhance surface precipitation. Recently, numerical weather prediction modeling studies have aimed at representing the ice crystal number concentration in mixed-phase clouds more accurately by including secondary ice formation processes. The increase in the ice crystal number concentration can impact the number of ice particles that sediment into the lower cloud and alter its composition and precipitation formation. In the Swiss Alps, the orography permits the formation of orographic clouds, making it ideal for studying the occurrence of multi-layered clouds and the seeder–feeder process. We present results from a case study on 18 May 2016, showing the occurrence frequency of multi-layered clouds and the seeder–feeder process. About half of all observed clouds were categorized as multi-layered, and the external seeder–feeder process occurred in 10 % of these clouds. Between cloud layers, ≈60 % of the ice particle mass was lost due to sublimation or melting. The external seeder–feeder process was found to be more important than the internal seeder–feeder process with regard to the impact on precipitation. In the case where the external seeder–feeder process was inhibited, the average surface precipitation and riming rate over the domain were both reduced by 8.5 % and 3.9 %, respectively. When ice–graupel collisions were allowed, further large reductions were seen in the liquid water fraction and riming rate. Inhibiting the internal seeder–feeder process enhanced the liquid water fraction by 6 % compared to a reduction of 5.8 % in the cloud condensate, therefore pointing towards the de-amplification in cloud glaciation and a reduction in surface precipitation. Adding to the observational evidence of frequent seeder–feeder situations, at least over Switzerland, our study highlights the extensive influence of sedimenting ice particles on the properties of feeder clouds as well as on precipitation formation.

Secondary Ice Production in Simulated Deep Convective Clouds: A Sensitivity Study

Abstract Multiple mechanisms have been proposed to explain secondary ice production (SIP), and SIP has been recognized to play a vital role in forming cloud ice crystals. However, most weather and climate models do not consider SIP in their cloud microphysical schemes. In this study, in addition to the default rime splintering (RS) process, two SIP processes, namely, shattering/fragmentation during freezing of supercooled rain/drizzle drops (DS) and breakup upon ice–ice collisions (BR), were implemented into a two-moment cloud microphysics scheme. Besides, two different parameterization schemes for BR were introduced. A series of sensitivity experiments were performed to investigate how SIP impacts cloud microphysics and cloud phase distributions in warm-based deep convective clouds developed in the central part of Europe. Simulation results revealed that cloud microphysical properties were significantly influenced by the SIP processes. Ice crystal number concentrations (ICNCs) increased up to more than 20 times and surface precipitation was reduced by up to 20% with the consideration of SIP processes. Interestingly, BR was found to dominate SIP, and the BR process rate was larger than the RS and DS process rates by four and three orders of magnitude, respectively. Liquid pixel number fractions inside clouds and at the cloud top decreased when implementing all three SIP processes, but the decrease depended on the BR scheme. Peak values of ice enhancement factors (IEFs) in the simulated deep convective clouds were 102–104 and located at −24°C with the consideration of all three SIP processes, while the temperature dependency of IEF was sensitive to the BR scheme. However, if only RS or RS and DS processes were included, the IEFs were comparable, with peak values of about 6, located at −7°C. Moreover, switching off the cascade effect led to a remarkable reduction in ICNCs and ice crystal mass mixing ratios. Significance Statement The cloud phase is found to have a significant impact on cloud evolution, radiative properties, and precipitation formation. However, the simulation of the cloud phase is a big challenge for cloud research because multiple processes are not well described or missing in numerical models. In this study, we implemented two secondary ice production (SIP) processes, namely, shattering/fragmentation during the freezing of supercooled rain/drizzle drops and breakup upon ice–ice collisions, which are missing in most numerical models. Sensitivity experiments were conducted to investigate how SIP impacts cloud microphysics and cloud phase in deep convective clouds. We found that SIP significantly impacts in-cloud and cloud-top phase distribution. We also identified that the collisional breakup of ice particles is the dominant SIP process in the simulated deep convective clouds.

Powering aircraft with 100 % sustainable aviation fuel reduces ice crystals in contrails

Abstract. Powering aircraft by sustainable aviation fuels (SAFs) is a pathway to reduce the climate impact of aviation by lowering aviation lifecycle CO2 emissions and by reducing ice crystal numbers and radiative forcing from contrails. While the effect of SAF blends on contrails has been measured previously, here we present novel measurements on particle emission and contrails from 100 % SAF combustion. During the ECLIF3 (Emission and CLimate Impact of alternative Fuels) campaign, a collaboration between the Deutsches Zentrum für Luft- und Raumfahrt (DLR), Airbus, Rolls-Royce, and Neste, the DLR Falcon 20 research aircraft performed in situ measurements following an Airbus A350-941 source aircraft powered by Rolls-Royce Trent XWB-84 engines in 1 to 2 min old contrails at cruise altitudes. Apparent ice emission indices of 100 % HEFA-SPK (hydro-processed esters and fatty acids–synthetic paraffinic kerosene) were measured and compared to Jet A-1 fuel contrails at similar engine and ambient ice-supersaturated conditions within a single flight. A 56 % reduction in ice particle numbers per mass of burned fuel was measured for 100 % HEFA-SPK compared to Jet A-1 under engine cruise conditions. The measured 35 % reduction in soot particle numbers suggests reduced ice activation by the low-sulfur HEFA fuel. Contrail properties are consistently modeled with a contrail plume model. Global climate model simulations for the 2018 fleet conservatively estimate a 26 % decrease in contrail radiative forcing and stronger decreases for larger particle reductions. Our results indicate that higher hydrogen content fuels as well as clean engines with low particle emissions may lead to reduced climate forcing from contrails.

Contrail formation on ambient aerosol particles for aircraft with hydrogen combustion: a box model trajectory study

Abstract. Future air traffic using (green) hydrogen (H2) promises zero carbon emissions, but the effects of contrails from this new technology have hardly been investigated. We study contrail formation behind aircraft with H2 combustion by means of the particle-based Lagrangian Cloud Module (LCM) box model. Assuming the absence of soot and ultrafine volatile particle formation, contrail ice crystals form solely on atmospheric background particles mixed into the plume. While a recent study extended the original LCM with regard to the contrail formation on soot particles, we further advance the LCM to cover the contrail formation on ambient particles. For each simulation, we perform an ensemble of box model runs using the dilution along 1000 different plume trajectories. The formation threshold temperature of H2 contrails is around 10 K higher than for conventional contrails (which form behind aircraft with kerosene combustion). Then, contrail formation becomes primarily limited by the homogeneous freezing temperature of the water droplets such that contrails can form at temperatures down to around 234 K. The number of ice crystals formed varies strongly with ambient temperature even far away from the contrail formation threshold. The contrail ice crystal number clearly increases with ambient aerosol number concentration and decreases significantly for ambient particles with mean dry radii ⪅ 10 nm due to the Kelvin effect. Besides simulations with one aerosol particle ensemble, we analyze contrail formation scenarios with two co-existing aerosol particle ensembles with different mean dry sizes or hygroscopicity parameters. We compare them to scenarios with a single ensemble that is the average of the two aerosol ensembles. We find that the total ice crystal number can differ significantly between the two cases, in particular if nucleation-mode particles are involved. Due to the absence of soot particle emissions, the ice crystal number in H2 contrails is typically reduced by more than 80 %–90 % compared to conventional contrails. The contrail optical thickness is significantly reduced, and H2 contrails either become visible later than kerosene contrails or are not visible at all for low ambient particle number concentrations. On the other hand, H2 contrails can form at lower flight altitudes where conventional contrails would not form.

Technical note: Bimodal parameterizations of in situ ice cloud particle size distributions

Abstract. The cloud particle size distribution (PSD) is a key parameter for the retrieval of microphysical and optical properties from remote-sensing instruments, which in turn are necessary for determining the radiative effect of clouds. Current representations of PSDs for ice clouds rely on parameterizations that were largely based on aircraft in situ measurements where the distribution of small ice crystals were uncertain. This makes current parameterizations deficient to simulate remote-sensing observations sensitive to small ice, such as from lidar and thermal infrared instruments. In this study we fit the in situ PSDs of ice crystals from the JULIA (JÜLich In situ Aircraft data set) database, which consists of 11 campaigns covering the tropics, midlatitudes and the Arctic, consistently processed and considered more robust in their measurements of small ice. For the fitting, we implement an established approach to PSD parameterizations, which consists of finding an adequate set of parameters for a modified gamma function after normalization of both PSD axes. These parameters are constrained to match in situ measurements when predicting microphysical properties from the PSDs, via a cost function minimization method. We selected the ice water content and the ice crystal number concentration, which are currently key parameters for modern satellite retrievals and model microphysics schemes. We found that a bimodal parameterization yields better results than a monomodal one. The bimodal parameterization has a lower spread for almost all ice crystal sizes over the entire range of analyzed temperatures and fits better the observations, especially for particles between 20 and about 110 µm at temperatures between −60 and −20 ∘C. For this temperature range, the root mean square error for the retrieved Nice is reduced from 0.36 to 0.20. This demonstrates a clear advantage to considering the bimodality of PSDs, e.g., for satellite retrievals.

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Abstract Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column. Significance Statement Secondary ice production (SIP) remains a key challenge in our understanding of boundary layer mixed-phase clouds. Here, we use sensitivity experiments performed with the ICON model at the cloud-resolving scale to explore potential interactions between primary nucleation, SIP, and mesoscale cloud organization. We simulate an Arctic single-layer cold air outbreak stratocumulus deck that was sampled during the M-PACE campaign. We find that once established, SIP alone is sufficient to maintain the mixed-phase cloud state until the end of the simulation. Our sensitivity analysis also shows that numerically more efficient treatments of immersion freezing are statistically invariant from simulations with a full prognostic INP budget.

Two new multirotor uncrewed aerial vehicles (UAVs) for glaciogenic cloud seeding and aerosol measurements within the CLOUDLAB project

Abstract. Uncrewed aerial vehicles (UAVs) have become widely used in a range of atmospheric science research applications. Because of their small size, flexible range of motion, adaptability, and low cost, multirotor UAVs are especially well-suited for probing the lower atmosphere. However, their use so far has been limited to conditions outside of clouds, first because of the difficulty of flying beyond visual line of sight and second because of the challenge of flying in icing conditions in supercooled clouds. Here, we present two UAVs for cloud microphysical research: one UAV (the measurement UAV) equipped with a Portable Optical Particle Spectrometer (POPS) and meteorological sensors to probe the aerosol and meteorological properties in the boundary layer and one UAV (the seeding UAV) equipped with seeding flares to produce a plume of particles that can nucleate ice in supercooled clouds. A propeller heating mechanism on both UAVs allows for operating in supercooled clouds with icing conditions. These UAVs are an integral part of the CLOUDLAB project in which glaciogenic cloud seeding of supercooled low stratus clouds is utilized for studying aerosol–cloud interactions and ice crystal formation and growth. In this paper, we first show validations of the POPS on board the measurement UAV, demonstrating that the rotor turbulence has a small effect on measured particle number concentrations. We then exemplify the applicability for profiling the planetary boundary layer, as well as for sampling and characterizing aerosol plumes, in this case, the seeding plume. We also present a new method for filtering out high-concentration data to ensure good data quality of POPS. We explain the different flight patterns that are possible for both UAVs, namely horizontal or vertical leg patterns or hovering, with an extensive and flexible parameter space for designing the flight patterns according to our scientific goals. Finally, we show two examples of seeding experiments: first characterizing an out-of-cloud seeding plume with the measurement UAV flying horizontal transects through the plume and, second, characterizing an in-cloud seeding plume with downstream measurements from a POPS and a holographic imager mounted on a tethered balloon. Particle number concentrations and particle number size distributions of the seeding plume from the experiments reveal that we can successfully produce and measure the seeding plume, both in-cloud (with accompanying elevated ice crystal number concentrations) and out-of-cloud. The methods presented here will be useful for probing the lower atmosphere, for characterizing aerosol plumes, and for deepening our cloud microphysical understanding through cloud seeding experiments, all of which have the potential to benefit the atmospheric science community.

Tropical Cirrus Are Highly Sensitive to Ice Microphysics Within a Nudged Global Storm‐Resolving Model

AbstractCirrus dominate the longwave radiative budget of the tropics. For the first time, the variability in cirrus properties and longwave cloud radiative effects (CREs) that arises from using different microphysical schemes within nudged global storm‐resolving simulations from a single model, is quantified. Nudging allows us to compute radiative biases precisely using coincident satellite measurements and to fix the large‐scale dynamics across our set of simulations to isolate the influence of microphysics. We run 5‐day simulations with four commonly‐used microphysics schemes of varying complexity (SAM1MOM, Thompson, M2005 and P3) and find that the tropical average longwave CRE varies over 20 W m−2 between schemes. P3 best reproduces observed longwave CRE. M2005 and P3 simulate cirrus with realistic frozen water path but unrealistically high ice crystal number concentrations which commonly hit limiters and lack the variability and dependence on frozen water content seen in aircraft observations. Thompson and SAM1MOM have too little cirrus.

Interactions between Gravity Waves and Cirrus Clouds: Asymptotic Modeling of Wave-Induced Ice Nucleation

Abstract We present an asymptotic approach for the systematic investigation of the effect of gravity waves (GWs) on ice clouds formed through homogeneous nucleation. In particular, we consider high- and midfrequency GWs in the tropopause region driving the formation of ice clouds, modeled with a double-moment bulk ice microphysics scheme. The asymptotic approach allows for identifying reduced equations for self-consistent description of the ice dynamics forced by GWs including the effects of diffusional growth and nucleation of ice crystals. Further, corresponding analytical solutions for a monochromatic GW are derived under a single-parcel approximation. The results provide a simple expression for the nucleated number of ice crystals in a nucleation event. It is demonstrated that the asymptotic solutions capture the dynamics of the full ice model and accurately predict the nucleated ice crystal number. The present approach is extended to allow for superposition of GWs, as well as for variable ice crystal mean mass in the deposition. Implications of the results for an improved representation of GW variability in cirrus parameterizations are discussed.

Wildfire smoke triggers cirrus formation: lidar observations over the eastern Mediterranean

Abstract. The number of intense wildfires may increase further in upcoming years as a consequence of climate change. It is therefore necessary to improve our knowledge about the role of smoke in the climate system, with emphasis on the impact of smoke particles on the evolution of clouds, precipitation, and cloud radiative properties. Presently, one key aspect of research is whether or not wildfire smoke particles can initiate cirrus formation. In this study, we present lidar observations over Limassol, Cyprus, from 27 October to 3 November 2020, when extended wildfire smoke fields crossed the Mediterranean Basin from Portugal to Cyprus. We found strong evidence that aged smoke (organic aerosol particles) originating from wildfires in North America triggered significant ice nucleation at temperatures from −47 to −53 ∘C and caused the formation of extended cirrus layers. The observations suggest that the ice crystals were nucleated just below the tropopause in the presence of smoke particles serving as ice-nucleating particles (INPs). The main part of the 2–3 km thick smoke layer was, however, in the lower stratosphere just above the tropopause. With actual radiosonde observations of temperature and relative humidity and lidar-derived smoke particle surface area concentrations used as starting values, gravity wave simulations show that the lofting of air by 100–200 m is sufficient to initiate significant ice nucleation on the smoke particles, leading to ice crystal number concentrations of 1–100 L−1.

Distinct secondary ice production processes observed in radar Doppler spectra: insights from a case study

Abstract. Secondary ice production (SIP) has an essential role in cloud and precipitation microphysics. In recent years, substantial insights were gained into SIP by combining experimental, modeling, and observational approaches. Remote sensing instruments, among them meteorological radars, offer the possibility of studying clouds and precipitation in extended areas over long time periods and are highly valuable to understand the spatiotemporal structure of microphysical processes. Multi-modal Doppler spectra measured by vertically pointing radars reveal the coexistence, within a radar resolution volume, of hydrometeor populations with distinct properties; as such, they can provide decisive insight into precipitation microphysics. This paper leverages polarimetric radar Doppler spectra as a tool to study the microphysical processes that took place during a snowfall event on 27 January 2021 in the Swiss Jura Mountains during the ICE GENESIS campaign. A multi-layered cloud system was present, with ice particles sedimenting through a supercooled liquid water (SLW) layer in a seeder–feeder configuration. Building on a Doppler peak detection algorithm, we implement a peak labeling procedure to identify the particle type(s) that may be present within a radar resolution volume. With this approach, we can visualize spatiotemporal features in the radar time series that point to the occurrence of distinct mechanisms during different stages of the event. By focusing on three 30 min phases of the case study and by using the detailed information contained in the Doppler spectra, together with dual-frequency radar measurements, aircraft in situ images, and simulated profiles of atmospheric variables, we narrow down the possible processes that could be responsible for the observed signatures. Depending on the availability of SLW and the droplet sizes, on the temperature range, and on the interaction between the liquid and ice particles, various SIP processes are identified as plausible, with distinct fingerprints in the radar Doppler spectra. A simple modeling approach suggests that the ice crystal number concentrations likely exceed typical concentrations of ice-nucleating particles by 1 to 4 orders of magnitude. While a robust proof of occurrence of a given SIP mechanism cannot be easily established, the multi-sensor data provide various independent elements each supporting the proposed interpretations.

Dust effects on mixed-phase clouds and precipitation during a super dust storm over northern China

Dust aerosols have been generally regarded as efficient ice nuclei (IN), but their influences on mixed-phase clouds and precipitation are poorly quantified. In this study, combining satellite observations, reanalysis data and the Weather Research and Forecasting (WRF) model, we investigate the impacts of dust aerosols on mixed-phase clouds and precipitation during a super dust storm over northern China. The Suomi National Polar-Orbiting Partnership (S-NPP) satellite observed a super dust storm from March 15–19, 2021 in northern China, followed by heavy precipitation between March 17 and 19. The results showed that the occurrence of super dust storms may play an important role in cloud formation and precipitation processes. Furthermore, numerical modeling revealed that dust aerosols can increase the ice crystal number concentration (QNi) and decrease the cloud droplet number concentration (QNc) in mixed-phase clouds through the Wegener–Bergeron–Findeisen (WBF) process, in which the maximum increase in QNi can reach 21%. Simultaneously, the mass mixing ratios of rain, graupel and snow increased due to dust aerosols. Consequently, precipitation could increase by up to 9.8% in northern China. This study could provide evidence for understanding the mechanisms of dust effects on mixed-phase clouds over northern China.

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.