Attendee Response to Guiding Questions of the September 2014 DOE Workshop on North Slope of Alaska Priorities

We advocate for several modeling approaches suitable to prove understanding of and capability to model aerosol behavior in nucleating droplets and ice and to model regional cloud systems at the mesoscale. The full complement of data that is ideal for particular study types is increasingly well understood; the ASR ISDAC Science Plan provides an excellent example of listing the data required for each study type targeted by that campaign. Whereas discussion here focused on the microscale and mesoscale, other studies should obviously focus on the regional, seasonal and global scale. When considering ice microphysics alone, laser focus on unconstrained model components could be usefully intensified (e.g., collision-coalescence kernels). Similarly, a focus on boundary layer structure would be deeply valuable to a wide range of model types. The expense of observations motivates efficient planning for observational data sets to support specific study approaches; the concept of closure is often useful (e.g., the surface energy budget)

An Evaluation of Size-Resolved Cloud Microphysics Scheme Numerics for Use with Radar Observations. Part I: Collision-Coalescence

This study evaluates some available schemes designed to solve the stochastic collection equation (SCE) for collision-coalescence of hydrometeors using a size-resolved (bin) microphysics approach, and documents their numerical properties within the framework of a box model. Comparing three widely used SCE schemes, we find that all converge to almost identical solutions at sufficiently fine mass grids. However, one scheme converges far slower than the other two and shows pronounced numerical diffusion at the large-drop tail of the size distribution. One of the remaining two schemes is recommended on the basis that it is well-converged on a relatively coarse mass grid, stable for large time steps, strictly mass-conservative, and computationally efficient. To examine the effects of SCE scheme choice on simulating clouds and precipitation, two of the three schemes are compared in large-eddy simulations of a drizzling stratocumulus field. A forward simulator that produces Doppler spectra from the large-eddy simulation results is used to compare the model output directly with radar observations. The scheme with pronounced numerical diffusion predicts excessively large mean Doppler velocities and overly broad and negatively skewed spectra compared with observations, consistent with numerical diffusion demonstrated in the box model. Statistics obtained using the recommended scheme are closer to observations, but notable differences remain, indicating that factors other than SCE scheme accuracy are limiting simulation fidelity

A Flexible Parameterization for Shortwave Optical Properties of Ice Crystals

A parameterization is presented that provides extinction cross section sigma (sub e), single-scattering albedo omega, and asymmetry parameter (g) of ice crystals for any combination of volume, projected area, aspect ratio, and crystal distortion at any wavelength in the shortwave. Similar to previous parameterizations, the scheme makes use of geometric optics approximations and the observation that optical properties of complex, aggregated ice crystals can be well approximated by those of single hexagonal crystals with varying size, aspect ratio, and distortion levels. In the standard geometric optics implementation used here, sigma (sub e) is always twice the particle projected area. It is shown that omega is largely determined by the newly defined absorption size parameter and the particle aspect ratio. These dependences are parameterized using a combination of exponential, lognormal, and polynomial functions. The variation of (g) with aspect ratio and crystal distortion is parameterized for one reference wavelength using a combination of several polynomials. The dependences of g on refractive index and omega are investigated and factors are determined to scale the parameterized (g) to provide values appropriate for other wavelengths. The parameterization scheme consists of only 88 coefficients. The scheme is tested for a large variety of hexagonal crystals in several wavelength bands from 0.2 to 4 micron, revealing absolute differences with reference calculations of omega and (g) that are both generally below 0.015. Over a large variety of cloud conditions, the resulting root-mean-squared differences with reference calculations of cloud reflectance, transmittance, and absorptance are 1.4%, 1.1%, and 3.4%, respectively. Some practical applications of the parameterization in atmospheric models are highlighted

Systematic Relationships Between Lidar Observables and Sizes And Mineral Composition Of Dust Aerosols

The physical and chemical properties of soil dust aerosol particles fundamentally affect their interaction with climate, including shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, heterogeneous formation of sulfates and nitrates on the surface of dust particles, and atmospheric processing of iron into bioavailable forms that increase the productivity of marine phytoplankton. Lidar measurements, such as extinction-to-backscatter, color and depolarization ratios, are frequently used to distinguish between aerosol types with different physical and chemical properties. The chemical composition of aerosol particles determines their complex refractive index, hence affecting their backscattering properties. Here we present a study on how dust aerosol backscattering and depolarization properties at wavelengths of 355, 532 and 1064 nm are related to size and complex refractive index, which varies with the mineral composition of the dust. Dust aerosols are represented by collections of spheroids with a range of prolate and oblate aspect ratios and their optical properties are obtained using T-matrix calculations. We find simple, systematic relationships between lidar observables and the dust size and complex refractive index that may aid the use of space-based or airborne lidars for direct retrieval of dust properties or for the evaluation of chemical transport models using forward simulated lidar variables. In addition, we present first results on the spatial variation of forward-simulated lidar variables based on a dust model that accounts for the atmospheric cycle of eight different mineral types plus internal mixtures of seven mineral types with iron oxides, which was recently implemented in the NASA GISS Earth System ModelE2

Highly supercooled riming and unusual triple-frequency radar signatures over Antarctica

Riming of ice crystals by supercooled water droplets is an efficient ice growth process, but its basic properties are still poorly known. While it has been shown to contribute significantly to surface precipitation at mid-latitudes, little is known about its occurrence at high latitudes. In Antarctica, two competing effects can influence the occurrence of riming: the scarcity of supercooled liquid water clouds due to the extremely low tropospheric temperatures and the low aerosol concentration, which may lead to the formation of fewer and larger supercooled drops potentially resulting in an enhanced riming efficiency. In this work, by exploiting the deployment of an unprecedented number of multi-wavelength active and passive remote sensing systems (including triple-frequency radar measurements) in West Antarctica, during the Atmospheric Radiation Measurements West Antarctic Radiation Experiment (AWARE) field campaign, we evaluate the importance of riming incidence in Antarctica and find that riming occurs at much lower temperatures compared to the mid-latitudes. We then focus on a case study featuring a persistent layer of unexpectedly pronounced triple-frequency radar signatures but only a relatively modest amount of supercooled liquid water. In-depth analysis of the radar observations suggests that such signatures can only be explained by the combined effects of moderately rimed aggregates or similarly shaped florid polycrystals and a narrow particle size distribution (PSD). Simulations of this case study performed with a 1D bin model %by introducing an additional class corresponding to rimed ice indicate that similar triple frequency radar observations can be reproduced when narrow PSDs are simulated. Such narrow PSDs can in turn be explained by two key factors: (i) the presence of a shallow homogeneous droplet or humidified aerosol freezing layer aloft seeding an underlying supercooled liquid layer, and (ii) the absence of turbulent mixing throughout a stable polar atmosphere that sustains narrow PSDs, as hydrometeors grow from the nucleation region aloft to several millimeter ice particles, by vapor deposition and then riming

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

ARM/GCSS/SPARC TWP-ICE CRM Intercomparison Study

Specifications are provided for running a cloud-resolving model (CRM) and submitting results in a standardized format for inclusion in a n intercomparison study and archiving for public access. The simulated case study is based on measurements obtained during the 2006 Tropical Warm Pool - International Cloud Experiment (TWP-ICE) led by the U. S. department of Energy Atmospheric Radiation Measurement (ARM) program. The modeling intercomparison study is based on objectives developed in concert with the Stratospheric Processes And their Role in Climate (SPARC) program and the GEWEX cloud system study (GCSS) program. The Global Energy and Water Cycle Experiment (GEWEX) is a core project of the World Climate Research PRogramme (WCRP)

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

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