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)

High-Resolution NU-WRF Simulations of a Deep Convective-Precipitation System During MC3E

The Goddard microphysics scheme was recently improved by adding a 4th ice class (frozen dropshail). This new 4ICE scheme was implemented and tested in the Goddard Cumulus Ensemble model (GCE) for an intense continental squall line and a moderate,less-organized continental case. Simulated peak radar reflectivity profiles were improved both 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 model (NU-WRF) and tested on an intense mesoscale convective system that occurred during the Midlatitude Continental Convective Clouds Experiment (MC3E). The NU42WRF simulated radar reflectivities, rainfall intensities, and vertical and horizontal structure using the new 4ICE scheme agree as well as or significantly better with observations than when using previous versions of the Goddard 3ICE (graupel or hail) schemes. In the 4ICE scheme, the bin microphysics-based rain evaporation correction produces more erect convective cores, while modification of the unrealistic collection of ice by dry hail produces narrow and intense cores, allowing more slow-falling snow to be transported rearward. Together with a revised snow size mapping, the 4ICE scheme produces a more horizontally stratified trailing stratiform region with a broad, more coherent light rain area. In addition, the NU-WRF 4ICE simulated radar reflectivity distributions are consistent with and generally superior to those using the GCE due to the less restrictive open lateral boundarie

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

Evaluation of Intercomparisons of Four Different Types of Model Simulating TWP-ICE

Four model intercomparisons were run and evaluated using the TWP-ICE field campaign, each involving different types of atmospheric model. Here we highlight what can be learnt from having single-column model (SCM), cloud-resolving model (CRM), global atmosphere model (GAM) and limited-area model (LAM) intercomparisons all based around the same field campaign. We also make recommendations for anyone planning further large multi-model intercomparisons to ensure they are of maximum value to the model development community. CRMs tended to match observations better than other model types, although there were exceptions such as outgoing long-wave radiation. All SCMs grew large temperature and moisture biases and performed worse than other model types for many diagnostics. The GAMs produced a delayed and significantly reduced peak in domain-average rain rate when compared to the observations. While it was shown that this was in part due to the analysis used to drive these models, the LAMs were also driven by this analysis and did not have the problem to the same extent. Based on differences between the models with parametrized convection (SCMs and GAMs) and those without (CRMs and LAMs), we speculate that that having explicit convection helps to constrain liquid water whereas the ice contents are controlled more by the representation of the microphysics