Ice Supersaturation Research Articles

Cloud Phase and Relative Humidity Distributions over the Southern Ocean in Austral Summer Based on In Situ Observations and CAM5 Simulations

AbstractCloud phase and relative humidity (RH) distributions at −67° to 0°C over the Southern Ocean during austral summer are compared between in situ airborne observations and global climate simulations. A scale-aware comparison is conducted using horizontally averaged observations from 0.1 to 50 km. Cloud phase frequencies, RH distributions, and liquid mass fraction are found to be less affected by horizontal resolutions than liquid and ice water content (LWC and IWC, respectively), liquid and ice number concentrations (Ncliq and Ncice, respectively), and ice supersaturation (ISS) frequency. At −10° to 0°C, observations show 27%–34% and 17%–37% of liquid and mixed phases, while simulations show 60%–70% and 3%–4%, respectively. Simulations overestimate (underestimate) LWC and Ncliq in liquid (mixed) phase, overestimate Ncice in mixed phase, underestimate IWC in ice and mixed phases, and underestimate (overestimate) liquid mass fraction below (above) −5°C, indicating that observational constraints are needed for different cloud phases. RH frequently occurs at liquid saturation in liquid and mixed phases for all datasets, yet the observed RH in ice phase can deviate from liquid saturation by up to 20%–40% at −20° to 0°C, indicating that the model assumption of liquid saturation for coexisting ice and liquid is inaccurate for low liquid mass fractions (<0.1). Simulations lack RH variability for partial cloud fractions (0.1–0.9) and underestimate (overestimate) ISS frequency for cloud fraction <0.1 (≥0.6), implying that improving RH subgrid-scale parameterizations may be a viable path to account for small-scale processes that affect RH and cloud phase heterogeneities. Two sets of simulations (nudged and free-running) show very similar results (except for ISS frequency) regardless of sample sizes, corroborating the statistical robustness of the model–observation comparisons.

High Depolarization Ratios of Naturally Occurring Cirrus Clouds Near Air Traffic Regions Over Europe

AbstractCirrus clouds have a large influence on the Earth's climate and anthropogenic activities such as aviation can alter their properties. Besides the formation of contrails, indirect effects on naturally occurring cirrus like increased heterogeneous freezing due to exhaust soot particles are discussed in the literature. However, hardly any observational study exists. In this work we present cirrus optical properties measured by an airborne lidar over Europe during the Midlatitude Cirrus experiment (ML‐CIRRUS). One half of the cloud cases showed elevated depolarization ratios with a mode difference of 10 percentage points indicating differences in the clouds microphysical properties. Their origin can be traced back to highly frequented air traffic regions, and they show lower in‐cloud ice supersaturations. Our analysis reveals no influence of embedded contrails and temperature. These results could be explained by an indirect aerosol effect where heterogeneous freezing is caused by aviation exhaust particles.

Heterogeneous Ice Nucleation in the Tropical Tropopause Layer

AbstractRecent laboratoryexperiments have advanced our understanding of heterogeneous ice nucleation at low temperatures. We use these laboratory results, along with field measurements of aerosol composition, to constrain a parameterization of heterogeneous nucleation rate dependence on temperature and supersaturation. We apply this nucleation parameterization in detailed simulations of transport and cloud formation in the Tropical Tropopause Layer (TTL, ≃14–18 km) constrained by high‐altitude aircraft measurements of TTL cirrus microphysical properties and supersaturation. The model results indicate that if the abundant TTL organic‐containing aerosols were effective glassy ice nuclei (IN), as indicated by laboratory experiments with simple sugars/acids, then heterogeneous nucleation would prevent the occurrences of large ice concentrations and large ice supersaturations, both of which are clearly indicated by the in situ observations. If glassy organic aerosols are relatively poor IN, as indicated by laboratory experiments with more atmospherically relevant secondary organic aerosol‐like composition, then the simulated ice concentrations and supersaturations are in reasonable agreement with the observations. Concentrations of effective mineral dust or ammonium sulfate IN larger than about 50/L can drive significant changes in cirrus microphysical properties and occurrence frequency. The cloud occurrence frequency can either increase or decrease, depending on the efficacy and abundance of IN added to the TTL. Relatively low TTL concentrations of mineral dust particles (≤10 L−1 indicated by the limited available field measurements) have negligible impact on cirrus frequencies or microphysical properties. Given the assortment of uncertainties, meaningful estimates of TTL aerosol‐cirrus indirect effects on radiative forcing are likely not possible at this time.

Multiyear observations of deposition-mode ice nucleating particles at two high-altitude stations in India

Ice nucleating particle (INP) measurements were made at two high-altitude stations in India. Aerosols collected on filter paper at Girawali Observatory, Inter University Center for Astronomy & Astrophysics (IGO), and at the Radio Astronomy Center, Ooty (RAC), were activated in deposition mode using a thermal gradient diffusion chamber to determine the INP concentrations. The measurement campaigns at IGO were conducted during 2011, 2013 and 2014, and at RAC during 2013 and 2014. When the aerosol samples were exposed to an ice supersaturation of between 5% and 23% in the temperature range −17.6°C to −22°C, the maximum INP number concentration at IGO and RAC was 1.0 L−1 and 1.6 L−1, respectively. A maximum correlation coefficient of 0.76 was observed between the INP number concentration and ice supersaturation. The airmass trajectories analyzed for the measurement campaigns showed that the Arabian Desert and arid regions were the main INP contributors. Elemental analysis of particles showed the presence of Na, Cl, Si, Al, Fe, Cu, Co, Cd, S, Mn and K, as well as some rare-Earth elements like Mo, Ru, La, Ce, V and Zr. When aerosols in the size range 0.5–20 μm were considered, the fraction that acted as INPs was 1: 104 to 1: 106 at IGO, and 1: 103 to 1: 104 at RAC. The higher ratio of INPs to aerosols at RAC than IGO may be attributable to the presence of rare-Earth elements observed in the aerosol samples at RAC, which were absent at IGO.

Ice Nucleation Parameterization and Relative Humidity Distribution in Idealized Squall-Line Simulations

AbstractOutput from idealized simulations of a squall line are compared with in situ aircraft-based observations from the Deep Convective Clouds and Chemistry campaign. Relative humidity distributions around convection are compared between 1-Hz aircraft observations (≈250-m horizontal scale) and simulations using a double-moment bulk microphysics scheme at three horizontal grid spacings: Δx = 0.25, 1, and 4 km. The comparisons focus on the horizontal extent of ice supersaturated regions (ISSRs), the maximum and average relative humidity with respect to ice (RHi) in ISSRs, and the ice microphysical properties during cirrus cloud evolution, with simulations at 0.25 and 1 km providing better results than the 4-km simulation. Within the ISSRs, all the simulations represent the dominant contributions of water vapor horizontal heterogeneities to ISSR formation on average, but with larger variabilities in such contributions than the observations. The best results are produced by a Δx = 0.25-km simulation with the RHi threshold for initiating ice nucleation increased to 130%, which improves almost all the ISSR characteristics and allows for larger magnitude and frequency of ice supersaturation (ISS) > 8%. This simulation also allows more occurrences of clear-sky ISSRs and a higher spatial fraction of ISS for in-cloud conditions, which are consistent with the observations. These improvements are not reproduced by modifying other ice microphysical processes, such as a factor-of-2 reduction in the ice nuclei concentration; a factor-of-10 reduction in the vapor deposition rate; turning off heterogeneous contact and immersion freezing; or turning off homogeneous freezing of liquid water.

Susceptibility of contrail ice crystal numbers to aircraft soot particle emissions

AbstractWe develop an idealized, physically based model describing combined effects of ice nucleation and sublimation on ice crystal number during persistent contrail formation. Our study represents the first effort to predict ice numbers at the point where contrails transition into contrail cirrus—several minutes past formation—by connecting them to aircraft soot particle emissions and atmospheric supersaturation with respect to ice. Results averaged over an observed exponential distribution of ice supersaturation (mean value 15%) indicate that large reductions in soot particle numbers are needed to lower contrail ice crystal numbers significantly for soot emission indices around 1015 (kg fuel)−1, because reductions in nucleated ice number are partially compensated by sublimation losses. Variations in soot particle (−50%) and water vapor (+10%) emission indices at threefold lower soot emissions resulting from biofuel blending cause ice crystal numbers to change by −35% and <5%, respectively. The efficiency of reduction depends on ice supersaturation and the size distribution of nucleated ice crystals in jet exhaust plumes and on atmospheric ice supersaturation, making the latter another key factor in contrail mitigation. We expect our study to have important repercussions for planning airborne measurements targeting contrail formation, designing parameterization schemes for use in large‐scale models, reducing uncertainties in predicting contrail cirrus, and mitigating the climate impact of aviation.

Modeling Condensation in Deep Convection

Abstract Cloud-scale models apply two drastically different methods to represent condensation of water vapor to form and grow cloud droplets. Maintenance of water saturation inside liquid clouds is assumed in the computationally efficient saturation adjustment approach used in most bulk microphysics schemes. When super- or subsaturations are allowed, condensation/evaporation can be calculated using the predicted saturation ratio and (either predicted or prescribed) mean droplet radius and concentration. The study investigates differences between simulations of deep unorganized convection applying a saturation adjustment condensation scheme (SADJ) and a scheme with supersaturation prediction (SPRE). A double-moment microphysics scheme with CCN activation parameterized as a function of the local vertical velocity is applied to compare cloud fields simulated applying SPRE and SADJ. Clean CCN conditions are assumed to demonstrate upper limits of the SPRE and SADJ difference. Microphysical piggybacking is used to extract the impacts with confidence. Results show a significant impact on deep convection dynamics, with SADJ featuring more cloud buoyancy and thus stronger updrafts. This leads to around a 3% increase of the surface rain accumulation in SADJ. Upper-tropospheric anvil cloud fractions are much larger in SPRE than in SADJ because of the higher ice concentrations and thus longer residence times of anvil particles in SPRE, as demonstrated by sensitivity tests. Higher ice concentrations in SPRE come from significantly larger ice supersaturations in strong convective updrafts that feature water supersaturations of several percent.

Influence of Bulk Microphysics Schemes upon Weather Research and Forecasting (WRF) Version 3.6.1 Nor'easter Simulations.

This study evaluated the impact of five, single- or double- moment bulk microphysics schemes (BMPSs) on Weather Research and Forecasting model (WRF) simulations of seven, intense winter time cyclones impacting the Mid-Atlantic United States. Five-day long WRF simulations were initialized roughly 24 hours prior to the onset of coastal cyclogenesis off the North Carolina coastline. In all, 35 model simulations (5 BMPSs and seven cases) were run and their associated microphysics-related storm properties (hydrometer mixing ratios, precipitation, and radar reflectivity) were evaluated against model analysis and available gridded radar and ground-based precipitation products. Inter-BMPS comparisons of column-integrated mixing ratios and mixing ratio profiles reveal little variability in non-frozen hydrometeor species due to their shared programming heritage, yet their assumptions concerning snow and graupel intercepts, ice supersaturation, snow and graupel density maps, and terminal velocities lead to considerable variability in both simulated frozen hydrometeor species and radar reflectivity. WRF-simulated precipitation fields exhibit minor spatio-temporal variability amongst BMPSs, yet their spatial extent is largely conserved. Compared to ground-based precipitation data, WRF-simulations demonstrate low-to-moderate (0.217-0.414) threat scores and a rainfall distribution shifted toward higher values. Finally, an analysis of WRF and gridded radar reflectivity data via contoured frequency with altitude (CFAD) diagrams reveals notable variability amongst BMPSs, where better performing schemes favored lower graupel mixing ratios and better underlying aggregation assumptions.

Technical note: On the intercalibration of HIRS channel 12 brightness temperatures following the transition from HIRS 2 to HIRS 3/4 for ice saturation studies

Abstract. In the present study we explore the capability of the intercalibrated HIRS brightness temperature data at channel 12 (the HIRS water vapour channel; T12) to reproduce ice supersaturation in the upper troposphere during the period 1979–2014. Focus is given on the transition from the HIRS 2 to the HIRS 3 instrument in the year 1999, which involved a shift of the central wavelength in channel 12 from 6.7 to 6.5 µm. It is shown that this shift produced a discontinuity in the time series of low T12 values ( < 235 K) and associated cases of high upper-tropospheric humidity with respect to ice (UTHi > 70 %) in the year 1999 which prevented us from maintaining a continuous, long-term time series of ice saturation throughout the whole record (1979–2014). We show that additional corrections are required to the low T12 values in order to bring HIRS 3 levels down to HIRS 2 levels. The new corrections are based on the cumulative distribution functions of T12 from NOAA 14 and 15 satellites (that is, when the transition from HIRS 2 to HIRS 3 occurred). By applying these corrections to the low T12 values we show that the discontinuity in the time series caused by the transition of HIRS 2 to HIRS 3 is not apparent anymore when it comes to calculating extreme UTHi cases. We come up with a new time series for values found at the low tail of the T12 distribution, which can be further exploited for analyses of ice saturation and supersaturation cases. The validity of the new method with respect to typical intercalibration methods such as regression-based methods is presented and discussed.

Dynamical conditions of ice supersaturation and ice nucleation in convective systems: A comparative analysis between in situ aircraft observations and WRF simulations

AbstractOccurrence frequency and dynamical conditions of ice supersaturation (ISS, where relative humidity with respect to ice (RHi) > 100%) are examined in the upper troposphere around convective activity. Comparisons are conducted between in situ airborne observations and the Weather Research and Forecasting model simulations using four double‐moment microphysical schemes at temperatures ≤ −40°C. All four schemes capture both clear‐sky and in‐cloud ISS conditions. However, the clear‐sky (in‐cloud) ISS conditions are completely (significantly) limited to the RHi thresholds of the Cooper parameterization. In all of the simulations, ISS occurrence frequencies are higher by ~3–4 orders of magnitude at higher updraft speeds (>1 m s−1) than those at the lower updraft speeds when ice water content (IWC) > 0.01 g m−3, while observations show smaller differences up to ~1–2 orders of magnitude. The simulated ISS also occurs less frequently at weaker updrafts and downdrafts than observed. These results indicate that the simulations have a greater dependence on stronger updrafts to maintain/generate ISS at higher IWC. At lower IWC (≤0.01 g m−3), simulations unexpectedly show lower ISS frequencies at stronger updrafts. Overall, the Thompson aerosol‐aware scheme has the closest magnitudes and frequencies of ISS >20% to the observations, and the modified Morrison has the closest correlations between ISS frequencies and vertical velocity at higher IWC and number density. The Cooper parameterization often generates excessive ice crystals and therefore suppresses the frequency and magnitude of ISS, indicating that it should be initiated at higher ISS (e.g., ≥25%).

An observational study of atmospheric ice nuclei number concentration during three fog-haze weather periods in Shenyang, northeastern China

Characteristics of ice nuclei (IN) number concentrations during three fog-haze weather periods from November 2010 to January 2012 in Shenyang were presented in this paper. A static diffusion chamber was used and sampling of IN aerosols was conducted using a membrane filter method. Sampling membrane filter processing conditions were unified in the activation temperature at −15°C under conditions of 20% ice supersaturation and 3% water supersaturation. The variations of natural IN number concentrations in different weather conditions were investigated. The relations between the meteorological factors and the IN number concentrations were analyzed, and relationships between pollutants and IN number concentrations were also studied. The results showed that mean IN number concentration were 38.68L−1 at −20°C in Shenyang, for all measurements. Mean IN number concentrations are higher during haze days (55.92L−1 at −20°C) and lower after rain. Of all meteorological factors, wind speed, boundary stability, and airflow direction appeared to influence IN number concentrations. IN number concentrations were positively correlated with particulate matters PM1, PM2.5, and PM10 during haze weather.

On the Life Cycle of Individual Contrails and Contrail Cirrus

Abstract The life cycle of individual (initially line shaped) contrails behind aircraft and of contrail cirrus (aged contrails mixed with other ice clouds) is described. The full contrail life cycle is covered, from ice formation for given water, heat, and particulate emissions; to changes in the jet, wake, and dispersion phases; through final sublimation or sedimentation. Contrail properties are deduced from various in situ, remote sensing, and model studies. Aerodynamically induced contrails and distrails are explained briefly. Contrails form both in clear air and inside cirrus. Young contrails consume most of the ambient ice supersaturation. Optical properties of contrails are age and humidity dependent. Contrail occurrence and radiative forcing depends on the ambient Earth–atmosphere conditions. Contrail cirrus seems to be optically thicker than assessed previously and may not only increase cirrus coverage but also thicken existing cirrus. Some observational constraints for contrail cirrus occurrence and radiative forcing are derived. Key parameters controlling contrail properties—besides aircraft and fuel properties, ambient pressure, temperature, and humidity—are the number of ice particles per flight distance surviving the wake vortex phase, the contrail depth, and particle sedimentation, wind shear, turbulence, and vertical motions controlling contrail dispersion. The climate impact of contrails depends among other things on the ratio of shortwave to longwave radiative forcing (RF) and on the efficacy with which contrail RF contributes to surface warming. Several open issues are identified, including renucleation from residuals of sublimated contrail ice particles.

Remote sensing ice supersaturation inside and near cirrus clouds: a case study in the subtropics

AbstractCombining vertically resolved lidar retrievals of water vapor and cloud detection, we document a 2‐day subtropical cirrus case study over La Réunion (20.9°S–55.5°E) in March 2005, focusing on the conditions of ice supersaturation inside and near the observed cloud. Using satellite observations, we describe the synoptic conditions leading to cloud formation. Supersaturation occurs 25% of the time within the cirrus, up to 35% in its middle segment, where relative humidity goes beyond 150%. In clear‐sky areas, relative humidity stays consistently low, especially in profiles without clouds. High‐troposphere atmospheric waves could initiate the formation of supersaturation conditions, especially on 16 March.

An assessment of the radiative effects of ice supersaturation based on in situ observations

AbstractWe use aircraft observations combined with the reanalysis data to investigate the radiative effects of ice supersaturation (ISS). Our results show that although the excess water vapor over ice saturation itself has relatively small radiative effects, mistaking it as ice crystals in climate models would lead to considerable impacts: on average, +2.49 W/m2 change in the top of the atmosphere (TOA) radiation, −2.7 W/m2 change in surface radiation, and 1.47 K/d change in heating rates. The radiative effects of ISS generally increase with the magnitudes of supersaturation. However, there is a strong dependence on the preexisting ice water path, which can even change the sign of the TOA radiative effect. It is therefore important to consider coexistence between ISS and ice clouds and to validate their relationship in the parameterizations of ISS in climate models.

Development and characterization of an ice-selecting pumped counterflow virtual impactor (IS-PCVI) to study ice crystal residuals

Abstract. Separation of particles that play a role in cloud activation and ice nucleation from interstitial aerosols has become necessary to further understand aerosol-cloud interactions. The pumped counterflow virtual impactor (PCVI), which uses a vacuum pump to accelerate the particles and increase their momentum, provides an accessible option for dynamic and inertial separation of cloud elements. However, the use of a traditional PCVI to extract large cloud hydrometeors is difficult mainly due to its small cut-size diameters (< 5 µm). Here, for the first time we describe a development of an ice-selecting PCVI (IS-PCVI) to separate ice in controlled mixed-phase cloud system based on the particle inertia with the cut-off diameter ≥ 10 µm. We also present its laboratory application demonstrating the use of the impactor under a wide range of temperature and humidity conditions. The computational fluid dynamics simulations were initially carried out to guide the design of the IS-PCVI. After fabrication, a series of validation laboratory experiments were performed coupled with the Aerosol Interaction and Dynamics in the Atmosphere (AIDA) expansion cloud simulation chamber. In the AIDA chamber, test aerosol particles were exposed to the ice supersaturation conditions (i.e., RHice > 100 %), where a mixture of droplets and ice crystals was formed during the expansion experiment. In parallel, the flow conditions of the IS-PCVI were actively controlled, such that it separated ice crystals from a mixture of ice crystals and cloud droplets, which were of diameter ≥ 10 µm. These large ice crystals were passed through the heated evaporation section to remove the water content. Afterwards, the residuals were characterized with a suite of online and offline instruments downstream of the IS-PCVI. These results were used to assess the optimized operating parameters of the device in terms of (1) the critical cut-size diameter, (2) the transmission efficiency and (3) the counterflow-to-input flow ratio. Particle losses were characterized by comparing the residual number concentration to the rejected interstitial particle number concentration. Overall results suggest that the IS-PCVI enables inertial separation of particles with a volume-equivalent particle size in the range of ~ 10–30 µm in diameter with small inadvertent intrusion (~ 5 %) of unwanted particles.

Levitation Diffusion Chamber Measurements of the Mass Growth of Small Ice Crystals from Vapor

Abstract A levitation diffusion chamber designed to examine the mass growth from the vapor of small ice particles (diameter < 100 μm) at ambient pressure (≃970 hPa) and low temperature (T < −30°C) is presented. The diffusion chamber is unique in that charged ice particles are levitated by an opposing voltage on the lower copper plate with lateral stability provided by button quadrupole electrodes attached to the upper copper plate. The button electrodes are far from the ice particle growth region, allowing ice particles to grow free of substrate influences. Experiments have been conducted for temperatures from −30° to −35.7°C, ice supersaturations from 2.5% to 28.6%, and over growth times ranging from 5 to 15 min. The experiments indicate that mass varies nonlinearly in time and exhibits a dependence on initial particle radius and ice supersaturation in accord with expectations from theory. In contrast to expectations from spherical capacitance theory, the derived mass growth rates do not scale linearly with radius, and derived effective shape factors (capacitance normalized with radius) are approximately 0.5. Fitting the growth data with a theoretical model indicates that growth is limited by surface kinetics with deposition coefficients ranging from 0.003 to 0.02.

Sulfate Aerosols from Non-Explosive Volcanoes: Chemical-Radiative Effects in the Troposphere and Lower Stratosphere

SO2 and H2S are the two most important gas-phase sulfur species emitted by volcanoes, with a global amount from non-explosive emissions of the order 10 Tg-S/yr. These gases are readily oxidized forming SO42− aerosols, which effectively scatter the incoming solar radiation and cool the surface. They also perturb atmospheric chemistry by enhancing the NOx to HNO3 heterogeneous conversion via hydrolysis on the aerosol surface of N2O5 and Br-Cl nitrates. This reduces formation of tropospheric O3 and the OH to HO2 ratio, thus limiting the oxidation of CH4 and increasing its lifetime. In addition to this tropospheric chemistry perturbation, there is also an impact on the NOx heterogeneous chemistry in the lower stratosphere, due to vertical transport of volcanic SO2 up to the tropical tropopause layer. Furthermore, the stratospheric O3 formation and loss, as well as the NOx budget, may be slightly affected by the additional amount of upward diffused solar radiation and consequent increase of photolysis rates. Two multi-decadal time-slice runs of a climate-chemistry-aerosol model have been designed for studying these chemical-radiative effects. A tropopause mean global net radiative flux change (RF) of −0.23 W·m−2 is calculated (including direct and indirect aerosol effects) with a 14% increase of the global mean sulfate aerosol optical depth. A 5–15 ppt NOx decrease is found in the mid-troposphere subtropics and mid-latitudes and also from pole to pole in the lower stratosphere. The tropospheric NOx perturbation triggers a column O3 decrease of 0.5–1.5 DU and a 1.1% increase of the CH4 lifetime. The surface cooling induced by solar radiation scattering by the volcanic aerosols induces a tropospheric stabilization with reduced updraft velocities that produce ice supersaturation conditions in the upper troposphere. A global mean 0.9% decrease of the cirrus ice optical depth is calculated with an indirect RF of −0.08 W·m−2.

The temporal evolution of a long‐lived contrail cirrus cluster: Simulations with a global climate model

AbstractThe representation of contrail cirrus in climate models has advanced in the last years tremendously. Nevertheless, uncertainties in particular regarding the representation of contrail microphysics still remain. Properties of young contrail cirrus differ from those of natural cirrus due to the large ice crystal number concentration common in contrails. Consequently, microphysical process rates in contrail cirrus, which control its lifetime, can be very different to those in natural cirrus. We extend a contrail cirrus scheme within a climate model by implementing a microphysical two‐moment scheme and study the life cycle of a contrail cirrus cluster. In an idealized experiment we study the properties and microphysical process rates of a contrail cirrus cluster in a large and long‐lived ice supersaturated region. We find that at flight level contrail cirrus display their typical high ice crystal number concentration (of about 10–100 cm−3) for a few hours with far lower densities in lower levels caused by sedimentation. After about 7 h contrail cirrus have spread considerably so that even at flight level associated ice crystal number concentrations have dropped to values that prohibit fast relaxation of ice supersaturation. The reduced ice crystal number and the resulting limited water uptake in the contrail cirrus limit the lifetime of the contrail cirrus cluster to about 10 h even though surrounding conditions would be still favorable for contrail cirrus persistence. In our case studies, contrail cirrus resembles natural cirrus regarding their ice crystal number concentration and size after 5–7 h.

Dynamics of Cloud-Top Generating Cells in Winter Cyclones. Part I: Idealized Simulations in the Context of Field Observations

Abstract This paper assesses the influence of radiative forcing and latent heating on the development and maintenance of cloud-top generating cells (GCs) in high-resolution idealized Weather Research and Forecasting Model simulations with initial conditions representative of the vertical structure of a cyclone observed during the Profiling of Winter Storms campaign. Simulated GC kinematics, structure, and ice mass are shown to compare well quantitatively with Wyoming Cloud Radar, cloud probe, and other observations. Sensitivity to radiative forcing was assessed in simulations with longwave-only (nighttime), longwave-and-shortwave (daytime), and no-radiation parameterizations. The domain-averaged longwave cooling rate exceeded 0.50 K h−1 near cloud top, with maxima greater than 2.00 K h−1 atop GCs. Shortwave warming was weaker by comparison, with domain-averaged values of 0.10–0.20 K h−1 and maxima of 0.50 K h−1 atop GCs. The stabilizing influence of cloud-top shortwave warming was evident in the daytime simulation’s vertical velocity spectrum, with 1% of the updrafts in the 6.0–8.0-km layer exceeding 1.20 m s−1, compared to 1.80 m s−1 for the nighttime simulation. GCs regenerate in simulations with radiative forcing after the initial instability is released but do not persist when radiation is not parameterized, demonstrating that radiative forcing is critical to GC maintenance under the thermodynamic and vertical wind shear conditions in this cyclone. GCs are characterized by high ice supersaturation (RHice > 150%) and latent heating rates frequently in excess of 2.00 K h−1 collocated with vertical velocity maxima. Ice precipitation mixing ratio maxima of greater than 0.15 g kg−1 were common within GCs in the daytime and nighttime simulations.

High-resolution NU-WRF simulations of a deep convective-precipitation system during MC3E: Further improvements and comparisons between Goddard microphysics schemes and observations.

The Goddard microphysics was recently improved by adding a fourth ice class (frozen drops/hail). This new 4ICE scheme was developed and tested in the Goddard Cumulus Ensemble (GCE) model for an intense continental squall line and a moderate, less organized continental case. Simulated peak radar reflectivity profiles were improved in intensity and shape for both cases, as were the overall reflectivity probability distributions versus observations. In this study, the new Goddard 4ICE scheme is implemented into the regional-scale NASA Unified-Weather Research and Forecasting (NU-WRF) model, modified and evaluated for the same intense squall line, which occurred during the Midlatitude Continental Convective Clouds Experiment (MC3E). NU-WRF simulated radar reflectivities, total rainfall, propagation, and convective system structures using the 4ICE scheme modified herein agree as well as or significantly better with observations than the original 4ICE and two previous 3ICE (graupel or hail) versions of the Goddard microphysics. With the modified 4ICE, the bin microphysics-based rain evaporation correction improves propagation and in conjunction with eliminating the unrealistic dry collection of ice/snow by hail can replicate the erect, narrow, and intense convective cores. Revisions to the ice supersaturation, ice number concentration formula, and snow size mapping, including a new snow breakup effect, allow the modified 4ICE to produce a stronger, better organized system, more snow, and mimic the strong aggregation signature in the radar distributions. NU-WRF original 4ICE simulated radar reflectivity distributions are consistent with and generally superior to those using the GCE due to the less restrictive domain and lateral boundaries.