High Ice Water Content at Low Radar Reflectivity near Deep Convection: Part II. Evaluation of Microphysical Pathways in Updraft Parcel Simulations

The aeronautics industry has established that a threat to aircraft is posed by atmospheric conditions of substantial ice water content (IWC) where equivalent radar reflectivity (Ze) does not exceed 20-30 dBZ and supercooled water is not present; these conditions are encountered almost exclusively in the vicinity of deep convection. Part 1 (Fridlind et al., 2015) of this two-part study presents in situ measurements of such conditions sampled by Airbus in three tropical regions, commonly near 11 km and -43 C, and concludes that the measured ice particle size distributions are broadly consistent with past literature with profiling radar measurements of Z(sub e) and mean Doppler velocity obtained within monsoonal deep convection in one of the regions sampled. In all three regions, the Airbus measurements generally indicate variable IWC that often exceeds 2 gm (exp -3) with relatively uniform mass median area-equivalent diameter (MMD(sub eq) of 200-300 micrometers. Here we use a parcel model with size-resolved microphysics to investigate microphysical pathways that could lead to such conditions. Our simulations indicate that homogeneous freezing of water drops produces a much smaller ice MMD(sub eq) than observed, and occurs only in the absence of hydrometeor gravitational collection for the conditions considered. Development of a mass mode of ice aloft that overlaps with the measurements requires a substantial source of small ice particles at temperatures of about -10 C or warmer, which subsequently grow from water vapor. One conceivable source in our simulation framework is Hallett-Mossop ice production; another is abundant concentrations of heterogeneous ice freezing nuclei acting together with copious shattering of water drops upon freezing. Regardless of the production mechanism, the dominant mass modal diameter of vapor-grown ice is reduced as the ice-multiplication source strength increases and as competition for water vapor increases. Both mass and modal diameter are reduced by entrainment and by increasing aerosol concentrations. Weaker updrafts lead to greater mass and larger modal diameters of vapor-grown ice, the opposite of expectations regarding lofting of larger ice particles in stronger updrafts. While stronger updrafts do loft more dense ice particles produced primarily by raindrop freezing, we find that weaker updrafts allow the warm rain process to reduce competition for diffusional growth of the less dense ice expected to persist in convective outflow

Remote sensing of ice crystal asymmetry parameter using multi-directional polarization measurements – Part 1: Methodology and evaluation with simulated measurements

We present a new remote sensing technique to infer the average asymmetry parameter of ice crystals near cloud top from multi-directional polarization measurements. The method is based on previous findings that (a) complex aggregates of hexagonal crystals generally have scattering phase matrices resembling those of their components; and (b) scattering phase matrices systematically vary with aspect ratios of crystals and their degree of microscale surface roughness. Ice cloud asymmetry parameters are inferred from multi-directional polarized reflectance measurements by searching for the closest fit in a look-up table of simulated polarized reflectances computed for cloud layers that contain individual, randomly oriented hexagonal columns and plates with varying aspect ratios and roughness values. The asymmetry parameter of the hexagonal particle that leads to the best fit with the measurements is considered the retrieved value. For clouds with optical thickness less than 5, the cloud optical thickness must be retrieved simultaneously with the asymmetry parameter, while for optically thicker clouds the asymmetry parameter retrieval is independent of cloud optical thickness. Evaluation of the technique using simulated measurements based on the optical properties of a number of complex particles and their mixtures shows that the ice crystal asymmetry parameters are generally retrieved to within 5%, or about 0.04 in absolute terms. The retrieval scheme is largely independent of calibration errors, range and sampling density of scattering angles and random noise in the measurements. The approach can be applied to measurements of past, current and future airborne and satellite instruments that measure multi-directional polarized reflectances of ice-topped clouds

How the Assumed Size Distribution of Dust Minerals Affects the Predicted Ice Forming Nuclei

The formation of ice in clouds depends on the availability of ice forming nuclei (IFN). Dust aerosol particles are considered the most important source of IFN at a global scale. Recent laboratory studies have demonstrated that the mineral feldspar provides the most efficient dust IFN for immersion freezing and together with kaolinite for deposition ice nucleation, and that the phyllosilicates illite and montmorillonite (a member of the smectite group) are of secondary importance.A few studies have applied global models that simulate mineral specific dust to predict the number and geographical distribution of IFN. These studies have been based on the simple assumption that the mineral composition of soil as provided in data sets from the literature translates directly into the mineral composition of the dust aerosols. However, these tables are based on measurements of wet-sieved soil where dust aggregates are destroyed to a large degree. In consequence, the size distribution of dust is shifted to smaller sizes, and phyllosilicates like illite, kaolinite, and smectite are only found in the size range 2 m. In contrast, in measurements of the mineral composition of dust aerosols, the largest mass fraction of these phyllosilicates is found in the size range 2 m as part of dust aggregates. Conversely, the mass fraction of feldspar is smaller in this size range, varying with the geographical location. This may have a significant effect on the predicted IFN number and its geographical distribution.An improved mineral specific dust aerosol module has been recently implemented in the NASA GISS Earth System ModelE2. The dust module takes into consideration the disaggregated state of wet-sieved soil, on which the tables of soil mineral fractions are based. To simulate the atmospheric cycle of the minerals, the mass size distribution of each mineral in aggregates that are emitted from undispersed parent soil is reconstructed. In the current study, we test the null-hypothesis that simulating the presence of a large mass fraction of phyllosilicates in dust aerosols in the size range 2 m, in comparison to a simple model assumption where this is neglected, does not yield a significant effect on the magnitude and geographical distribution of the predicted IFN number. Results from sensitivity experiments are presented as well

High Ice Water Content at Low Radar Reflectivity near Deep Convection

Occurrences of jet engine power loss and damage have been associated with flight through fully glaciated deep convection at -10 to -50 degrees Centigrade. Power loss events commonly occur during flight through radar reflectivity (Zeta (sub e)) less than 20-30 decibels relative to Zeta (dBZ - radar returns) and no more than moderate turbulence, often overlying moderate to heavy rain near the surface. During 2010-2012, Airbus carried out flight tests seeking to characterize the highest ice water content (IWC) in such low-radar-reflectivity regions of large, cold-topped storm systems in the vicinity of Cayenne, Darwin, and Santiago. Within the highest IWC regions encountered, at typical sampling elevations (circa 11 kilometers), the measured ice size distributions exhibit a notably narrow concentration of mass over area-equivalent diameters of 100-500 micrometers. Given substantial and poorly quantified measurement uncertainties, here we evaluate the consistency of the Airbus in situ measurements with ground-based profiling radar observations obtained under quasi-steady, heavy stratiform rain conditions in one of the Airbus-sampled locations. We find that profiler-observed radar reflectivities and mean Doppler velocities at Airbus sampling temperatures are generally consistent with those calculated from in situ size-distribution measurements. We also find that column simulations using the in situ size distributions as an upper boundary condition are generally consistent with observed profiles of radar reflectivity (Ze), mean Doppler velocity (MDV), and retrieved rain rate. The results of these consistency checks motivate an examination of the microphysical pathways that could be responsible for the observed size-distribution features in Ackerman et al. (2015)

How the Emitted Size Distribution and Mixing State of Feldspar Affect Ice Nucleating Particles in a Global Model

The effect of aerosol particles on ice nucleation and, in turn, the formation of ice and mixed phase clouds is recognized as one of the largest sources of uncertainty in climate prediction. We apply an improved dust mineral specific aerosol module in the NASA GISS Earth System ModelE, which takes into account soil aggregates and their fragmentation at emission as well as the emission of large particles. We calculate ice nucleating particle concentrations from K-feldspar abundance for an active site parameterization for a range of activation temperatures and external and internal mixing assumption. We find that the globally averaged INP concentration is reduced by a factor of two to three, compared to a simple assumption on the size distribution of emitted dust minerals. The decrease can amount to a factor of five in some geographical regions. The results vary little between external and internal mixing and different activation temperatures, except for the coldest temperatures. In the sectional size distribution, the size range 24 micrometer contributes the largest INP number

Fast Physics Testbed for the FASTER Project

This poster describes the Fast Physics Testbed for the new FAst-physics System Testbed and Research (FASTER) project. The overall objective is to provide a convenient and comprehensive platform for fast turn-around model evaluation against ARM observations and to facilitate development of parameterizations for cloud-related fast processes represented in global climate models. The testbed features three major components: a single column model (SCM) testbed, an NWP-Testbed, and high-resolution modeling (HRM). The web-based SCM-Testbed features multiple SCMs from major climate modeling centers and aims to maximize the potential of SCM approach to enhance and accelerate the evaluation and improvement of fast physics parameterizations through continuous evaluation of existing and evolving models against historical as well as new/improved ARM and other complementary measurements. The NWP-Testbed aims to capitalize on the large pool of operational numerical weather prediction products. Continuous evaluations of NWP forecasts against observations at ARM sites are carried out to systematically identify the biases and skills of physical parameterizations under all weather conditions. The highresolution modeling (HRM) activities aim to simulate the fast processes at high resolution to aid in the understanding of the fast processes and their parameterizations. A four-tier HRM framework is established to augment the SCM- and NWP-Testbeds towards eventual improvement of the parameterizations

Use of Cloud Radar Doppler Spectra to Evaluate Stratocumulus Drizzle Size Distributions in Large-Eddy Simulations with Size-Resolved Microphysics

A case study of persistent stratocumulus over the Azores is simulated using two independent large-eddy simulation (LES) models with bin microphysics, and forward-simulated cloud radar Doppler moments and spectra are compared with observations. Neither model is able to reproduce the monotonic increase of downward mean Doppler velocity with increasing reflectivity that is observed under a variety of conditions, but for differing reasons. To a varying degree, both models also exhibit a tendency to produce too many of the largest droplets, leading to excessive skewness in Doppler velocity distributions, especially below cloud base. Excessive skewness appears to be associated with an insufficiently sharp reduction in droplet number concentration at diameters larger than ~200 μm, where a pronounced shoulder is found for in situ observations and a sharp reduction in reflectivity size distribution is associated with relatively narrow observed Doppler spectra. Effectively using LES with bin microphysics to study drizzle formation and evolution in cloud Doppler radar data evidently requires reducing numerical diffusivity in the treatment of the stochastic collection equation; if that is accomplished sufficiently to reproduce typical spectra, progress toward understanding drizzle processes is likely

Contemplating Synergistic Algorithms for the NASA ACE Mission

ACE is a proposed Tier 2 NASA Decadal Survey mission that will focus on clouds, aerosols, and precipitation as well as ocean ecosystems. The primary objective of the clouds component of this mission is to advance our ability to predict changes to the Earth's hydrological cycle and energy balance in response to climate forcings by generating observational constraints on future science questions, especially those associated with the effects of aerosol on clouds and precipitation. ACE will continue and extend the measurement heritage that began with the A-Train and that will continue through Earthcare. ACE planning efforts have identified several data streams that can contribute significantly to characterizing the properties of clouds and precipitation and the physical processes that force these properties. These include dual frequency Doppler radar, high spectral resolution lidar, polarimetric visible imagers, passive microwave and submillimeter wave radiometry. While all these data streams are technologically feasible, their total cost is substantial and likely prohibitive. It is, therefore, necessary to critically evaluate their contributions to the ACE science goals. We have begun developing algorithms to explore this trade space. Specifically, we will describe our early exploratory algorithms that take as input the set of potential ACE-like data streams and evaluate critically to what extent each data stream influences the error in a specific cloud quantity retrieval

Confronting the Challenge of Modeling Cloud and Precipitation Microphysics

In the atmosphere, microphysics refers to the microscale processes that affect cloud and precipitation particles and is a key linkage among the various components of Earth\u27s atmospheric water and energy cycles. The representation of microphysical processes in models continues to pose a major challenge leading to uncertainty in numerical weather forecasts and climate simulations. In this paper, the problem of treating microphysics in models is divided into two parts: (i) how to represent the population of cloud and precipitation particles, given the impossibility of simulating all particles individually within a cloud, and (ii) uncertainties in the microphysical process rates owing to fundamental gaps in knowledge of cloud physics. The recently developed Lagrangian particle‐based method is advocated as a way to address several conceptual and practical challenges of representing particle populations using traditional bulk and bin microphysics parameterization schemes. For addressing critical gaps in cloud physics knowledge, sustained investment for observational advances from laboratory experiments, new probe development, and next‐generation instruments in space is needed. Greater emphasis on laboratory work, which has apparently declined over the past several decades relative to other areas of cloud physics research, is argued to be an essential ingredient for improving process‐level understanding. More systematic use of natural cloud and precipitation observations to constrain microphysics schemes is also advocated. Because it is generally difficult to quantify individual microphysical process rates from these observations directly, this presents an inverse problem that can be viewed from the standpoint of Bayesian statistics. Following this idea, a probabilistic framework is proposed that combines elements from statistical and physical modeling. Besides providing rigorous constraint of schemes, there is an added benefit of quantifying uncertainty systematically. Finally, a broader hierarchical approach is proposed to accelerate improvements in microphysics schemes, leveraging the advances described in this paper related to process modeling (using Lagrangian particle‐based schemes), laboratory experimentation, cloud and precipitation observations, and statistical methods