Engineering topics

Fridlind, Ann M.

Publications and source records attributed to Fridlind, Ann M..

At least 37 records · Page 2

DOE OSTI2020

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’s 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.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2020

A New Look at the Environmental Conditions Favorable to Secondary Ice Production

This study attempts a new identification of mechanisms of secondary ice production (SIP) based on the observation of small faceted ice crystals (hexagonal plates or columns) with typical sizes smaller than 100µm. Due to their young age, such small ice crystals can be used as tracers for identifying the conditions for SIP. Observations reported here were conducted in oceanic tropical mesoscale convective systems (MCSs) and midlatitude frontal clouds in the temperature range from 0 to −15◦C and heavily seeded by aged ice particles. It was found that in both MCSs and frontal clouds, SIP was observed right above the melting layer and extended to higher altitudes with colder temperatures. The roles of six possible mechanisms to generate the SIP particles are assessed using additional observations. In most observed SIP cases, small secondary ice particles spatially correlated with liquid-phase, vertical updrafts and aged rimed ice particles. However, in many cases, neither graupel nor liquid drops were observed in the SIP regions, and therefore, the conditions for an active Hallett–Mossop process were not met. In many cases, large concentrations of small pristine ice particles were observed right above the melting layer, starting at temperatures as warm as −0.5◦C. It is proposed that the initiation of SIP above the melting layer is stimulated by the recirculation of large liquid drops through the melting layer with convective turbulent updrafts. After re-entering a super cooled environment above the melting layer, they impact with aged ice, freeze, and shatter. The size of the splinters generated during SIP was estimated as 10µm or less. A principal conclusion of this work is that only the freezing drop-shattering mechanism could be clearly supported by the airborne in situ observations.

Korolev, Alexei↗

NASA NTRS2019

Use of Polarimetric Radar Measurements to Constrain Simulated Convective Cell Evolution: A Pilot Study with Lagrangian Tracking

To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (K(sub DP)) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (Z(sub DR)) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived K(sub DP) and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that ∼ 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign.

Fridlind, Ann M.↗

NASA NTRS2019

Evaluating Models' Response of Tropical Low Clouds to SST Forcings Using CALIPSO Observations

Recent studies have shown that, in response to a surface warming, the marine tropical low-cloud cover (LCC) as observed by passive-sensor satellites substantially decreases, therefore generating a smaller negative value of the top-of-the-atmosphere (TOA) cloud radiative effect (CRE). Here we study the LCC and CRE interannual changes in response to sea surface temperature (SST) forcings in the GISS model E2 climate model, a developmental version of the GISS model E3 climate model, and in 12 other climate models, as a function of their ability to represent the vertical structure of the cloud response to SST change against 10 years of CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) observations. The more realistic models (those that satisfy the observational constraint) capture the observed interannual LCC change quite well ([delta]LCC]/[delta]SST= -3.49Â±1.01%K [negative 1 superscript] vs. [delta]LCC/[delta]SST[subscript obs]= -3.59Â±0.28%K[negative 1 superscript]) while the others largely underestimate it ([delta]LCC/[delta]SST= -1.32Â±1.28%K[negative 1 superscript]). Consequently, the more realistic models simulate more positive shortwave (SW) feedback ([delta]CRE/[delta]SST=2.60Â±1.13Wm[negative 2 superscript] K[negative 1 superscript]) than the less realistic models (CRE/SST=0.87Â±2.63Wm2K1), in better agreement with the observations ([delta]CRE/[delta]SST[subscript obs]=3Â±0.26Wm[negative 2 superscript] K[negative 1 superscript] ), although slightly underestimated. The ability of the models to represent moist processes within the planetary boundary layer (PBL) and produce persistent stratocumulus (Sc) decks appears crucial to replicating the observed relationship between clouds, radiation and surface temperature. This relationship is different depending on the type of low clouds in the observations. Over stratocumulus regions, cloud-top height increases slightly with SST, accompanied by a large decrease in cloud fraction, whereas over trade cumulus (Cu) regions, cloud fraction decreases everywhere, to a smaller extent.

low clouds↗

NASA NTRS2018

(GO)2-SIM: a GCM-Oriented Ground-Observation Forward-Simulator Framework for Objective Evaluation of Cloud and Precipitation Phase

General circulation model (GCM) evaluation using ground-based observations is complicated by inconsistencies in hydrometeor and phase definitions. Here we describe (GO)2-SIM, a forward simulator designed for objective hydrometeor-phase evaluation, and assess its performance over the North Slope of Alaska using a 1-year GCM simulation. For uncertainty assessment, 18 empirical relationships are used to convert model grid-average hydrometeor (liquid and ice, cloud, and precipitation) water contents to zenith polarimetric micropulse lidar and Ka-band Doppler radar measurements, producing an ensemble of 576 forward-simulation realizations. Sensor limitations are represented in forward space to objectively remove from consideration model grid cells with undetectable hydrometeor mixing ratios, some of which may correspond to numerical noise.Phase classification in forward space is complicated by the inability of sensors to measure ice and liquid signals distinctly. However, signatures exist in lidarâ€“radar space such that thresholds on observables can be objectively estimated and related to hydrometeor phase. The proposed phase-classification technique leads to misclassification in fewer than 8% of hydrometeor-containing grid cells. Such misclassifications arise because, while the radar is capable of detecting mixed-phase conditions, it can mistake water- for ice-dominated layers. However, applying the same classification algorithm to forward-simulated and observed fields should generate hydrometeor-phase statistics with similar uncertainty. Alternatively, choosing to disregard how sensors define hydrometeor phase leads to frequency of occurrence discrepancies of up to 40%. So, while hydrometeor-phase maps determined in forward space are very different from model "reality" they capture the information sensors can provide and thereby enable objective model evaluation.

Lamer, K.↗

NASA NTRS2018

Simulations of Arctic Mixed-Phase Boundary Layer Clouds: Advances in Understanding and Outstanding Questions

An overview of mixed-phase boundary layer cloud simulations emphasizes what detailed studies reveal about the level of understanding of microphysical processes, using analogous liquid-phase boundary layer clouds as a reference for the dynamical conditions. We focus primarily on modeling results from three major field campaign case studies, spanning a range of liquid water path, aerosol loading, cloud temperatures, and active processes (including drizzle, aggregation, and riming). Simulations with detailed microphysics are able to reproduce basic aspects of all three case studies, including continuous ice formation within the context of a well-mixed liquid-cloud-topped layer. However, simulations constrained with measurements of ice nuclei generally underpredict the observed amount of ice present, indicating a possible role for ice multiplication that is poorly established. Assumed ice properties may also bear a significant influence on the water vapor budget and reflectivity properties, but remain imprecisely quantified on a case study basis.

Arctic clouds; Mixed-phase clouds; Mixed-phase mic↗

NASA NTRS2017

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.

Perlwitz, Jan P.↗

NASA NTRS2017

Impacts of Solar-Absorbing Aerosol Layers on the Transition of Stratocumulus to Trade Cumulus Clouds

The effects of an initially overlying layer of solar-absorbing aerosol on the transition of stratocumulus to trade cumulus clouds are examined using large-eddy simulations. For lightly drizzling cloud the transition is generally hastened, resulting mainly from increased cloud droplet number concentration (Nc) induced by entrained aerosol. The increased Nc slows sedimentation of cloud droplets and shortens their relaxation time for diffusional growth, both of which accelerate entrainment of overlying air and thereby stratocumulus breakup. However, the decrease in albedo from cloud breakup is more than offset by redistributing cloud water over a greater number of droplets, such that the diurnal-average shortwave forcing at the top of the atmosphere is negative. The negative radiative forcing is enhanced by sizable longwave contributions, which result from the greater cloud breakup and a reduced boundary layer height associated with aerosol heating. A perturbation of moisture instead of aerosol aloft leads to a greater liquid water path and a more gradual transition. Adding absorbing aerosol to that atmosphere results in substantial reductions in liquid water path (LWP) and cloud cover that lead to positive short-wave and negative longwave forcings on average canceling each other. Only for heavily drizzling clouds is the breakup delayed, as inhibition of precipitation overcomes cloud water loss from enhanced entrainment. Considering these simulations as an imperfect proxy for biomass burning plumes influencing Namibian stratocumulus, we expect regional indirect plus semi-direct forcings to be substantially negative to negligible at the top of the atmosphere, with its magnitude sensitive to background and perturbation properties.

Zhou, Xiaoli↗

NASA NTRS2017

Derivation of Aerosol Profiles for MC3E Convection Studies and Use in Simulations of the 20 May Squall Line Case

Advancing understanding of deep convection microphysics via mesoscale modeling studies of well-observed case studies requires observation-based aerosol inputs. Here, we derive hygroscopic aerosol size distribution input profiles from ground-based and airborne measurements for six convection case studies observed during the Midlatitude Continental Convective Cloud Experiment (MC3E) over Oklahoma. We demonstrate use of an input profile in simulations of the only well-observed case study that produced extensive stratiform outflow on 20 May 2011. At well-sampled elevations between −11 and −23 C over widespread stratiform rain, ice crystal number concentrations are consistently dominated by a single mode near approx. 400 microm in randomly oriented maximum dimension (Dmax). The ice mass at −23 C is primarily in a closely collocated mode, whereas a mass mode near Dmax approx. 1000 microns becomes dominant with decreasing elevation to the −11 C level, consistent with possible aggregation during sedimentation. However, simulations with and without observation-based aerosol inputs systematically overpredict mass peak Dmax by a factor of 3-5 and underpredict ice number concentration by a factor of 4-10. Previously reported simulations with both two-moment and size-resolved microphysics have shown biases of a similar nature. The observed ice properties are notably similar to those reported from recent tropical measurements. Based on several lines of evidence, we speculate that updraft microphysical pathways determining outflow properties in the 20 May case are similar to a tropical regime, likely associated with warm-temperature ice multiplication that is not well understood or well represented in models.

Fridlind, Ann M.↗

NASA NTRS2017

Derivation of Aerosol Profiles for MC3E Convection Studies and Use in Simulations of the 20 May Squall Line Case

Advancing understanding of deep convection microphysics via mesoscale modeling studies of well-observed case studies requires observation-based aerosol inputs. Here, we derive hygroscopic aerosol size distribution input profiles from ground-based and airborne measurements for six convection case studies observed during the Midlatitude Continental Convective Cloud Experiment (MC3E) over Oklahoma. We demonstrate use of an input profile in simulations of the only well-observed case study that produced extensive stratiform outflow on 20 May 2011. At well-sampled elevations between -11 and -23 degree C over widespread stratiform rain, ice crystal number concentrations are consistently dominated by a single mode near 400 micrometer in randomly oriented maximum dimension (D[superscript max] ). The ice mass at -23 degree C is primarily in a closely collocated mode, whereas a mass mode near D[superscript max] -1000 micrometer becomes dominant with decreasing elevation to the -11 degree C level, consistent with possible aggregation during sedimentation. However, simulations with and without observation-based aerosol inputs systematically overpredict mass peak D[superscript max] by a factor of 3-5 and underpredict ice number concentration by a factor of 4-10. Previously reported simulations with both two-moment and size-resolved microphysics have shown biases of a similar nature. The observed ice properties are notably similar to those reported from recent tropical measurements. Based on several lines of evidence, we speculate that updraft microphysical pathways determining outflow properties in the 20 May case are similar to a tropical regime, likely associated with warm-temperature ice multiplication that is not well understood or well represented in models.

Fridlind, Ann M.↗

NASA NTRS2016

Derivation of Physical and Optical Properties of Midlatitude Cirrus Ice Crystals for a Size-Resolved Cloud Microphysics Model

Single-crystal images collected in mid-latitude cirrus are analyzed to provide internally consistent ice physical and optical properties for a size-resolved cloud microphysics model, including single-particle mass, projected area, fall speed, capacitance, single-scattering albedo, and asymmetry parameter. Using measurements gathered during two flights through a widespread synoptic cirrus shield, bullet rosettes are found to be the dominant identifiable habit among ice crystals with maximum dimension (Dmax) greater than 100μm. Properties are therefore first derived for bullet rosettes based on measurements of arm lengths and widths, then for aggregates of bullet rosettes and for unclassified (irregular) crystals. Derived bullet rosette masses are substantially greater than reported in existing literature, whereas measured projected areas are similar or lesser, resulting in factors of 1.5-2 greater fall speeds, and, in the limit of large Dmax, near-infrared single-scattering albedo and asymmetry parameter (g) greater by approx. 0.2 and 0.05, respectively. A model that includes commonly imaged side plane growth on bullet rosettes exhibits relatively little difference in microphysical and optical properties aside from approx. 0:05 increase in mid-visible g primarily attributable to plate aspect ratio. In parcel simulations, ice size distribution, and g are sensitive to assumed ice properties.

asymmetry↗

NASA NTRS2016

Vertical Variation of Ice Particle Size in Convective Cloud Tops

A novel technique is used to estimate derivatives of ice effective radius with respect to height near convective cloud tops (dr(sub e)/dz) from airborne shortwave reflectance measurements and lidar. Values of dr(sub e)/dz are about −6 micrometer/km for cloud tops below the homogeneous freezing level, increasing to near 0 micrometer/km above the estimated level of neutral buoyancy. Retrieved dr(sub e)/dz compares well with previously documented remote sensing and in situ estimates. Effective radii decrease with increasing cloud top height, while cloud top extinction increases. This is consistent with weaker size sorting in high, dense cloud tops above the level of neutral buoyancy where fewer large particles are present and with stronger size sorting in lower cloud tops that are less dense. The results also confirm that cloud top trends of effective radius can generally be used as surrogates for trends with height within convective cloud tops. These results provide valuable observational targets for model evaluation.

cloud height indicators↗

NASA NTRS2015

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.

Dust↗

NASA NTRS2015

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.

Dust↗

NASA NTRS2015

The Dependence of Cirrus Gamma Size Distributions Expressed as Volumes in N(sub 0)-Lambda-Mu Phase Space and Bulk Cloud Properties on Environmental Conditions: Results from the Small Ice Particles in Cirrus Experiment (SPARTICUS)

The variability of cirrus ice microphysical properties is investigated using observations obtained during the Small Particles in Cirrus (SPARTICUS) campaign. An existing approach that represents a size distribution (SD) as a single gamma function using an ellipsoid of equally realizable solutions in (N(sub 0), lambda, mu) phase space is modified to automatically identify multiple modes in SDs and characterize each mode by such an ellipsoid. The modified approach is applied to ice crystals with maximum dimension D greater than15 micrometers collected by the 2-D stereo and 2-D precipitation probes on the Stratton Park Engineering Company Learjet. The dependencies of N(sub 0), mu, and lambda from each mode, total number concentration, bulk extinction, ice water content (IWC), and mass median maximum dimension D(sub mm) as a function of temperature T and cirrus type are then analyzed. The changes in the observed codependencies between N(sub 0), mu, and lambda, bulk extinction, IWC, and D(sub mm) with environmental conditions indicate that particles were larger at higher T during SPARTICUS. At most two modes were observed in any SD during SPARTICUS, with the average boundary between them at 115 micrometers, similar to past studies not using probes with shatter mitigating tips and artifact removal algorithms. The bimodality of the SDs increased with T. This and the differences in N(sub 0), mu, and lambda between the modes suggest that particles with smaller D nucleated more recently than particles with larger D, which grew via vapor deposition and aggregation. Because smaller crystals, whose concentrations are uncertain, make marginal contributions to higher order moments, the use of higher moments for evaluating model fields is suggested.

Ellipsoids↗

NASA NTRS2015

RACORO Continental Boundary Layer Cloud Investigations: 1. Case Study Development and Ensemble Large-Scale Forcings

Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60 h case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in situ measurements from the Routine AAF (Atmospheric Radiation Measurement (ARM) Aerial Facility) CLOWD (Clouds with Low Optical Water Depth) Optical Radiative Observations (RACORO) field campaign and remote sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, kappa, are derived from observations to be approximately 0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing data sets are derived from the ARM variational analysis, European Centre for Medium-Range Weather Forecasts, and a multiscale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in "trial" large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.

Boundary layers↗

NASA NTRS2014

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.

parameterization↗

NASA NTRS2014

Intercomparison of Large-Eddy Simulations of Arctic Mixed-Phase Clouds: Importance of Ice Size Distribution Assumptions

Large-eddy simulations of mixed-phase Arctic clouds by 11 different models are analyzed with the goal of improving understanding and model representation of processes controlling the evolution of these clouds. In a case based on observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC), it is found that ice number concentration, Ni, exerts significant influence on the cloud structure. Increasing Ni leads to a substantial reduction in liquid water path (LWP), in agreement with earlier studies. In contrast to previous intercomparison studies, all models here use the same ice particle properties (i.e., mass-size, mass-fall speed, and mass-capacitance relationships) and a common radiation parameterization. The constrained setup exposes the importance of ice particle size distributions (PSDs) in influencing cloud evolution. A clear separation in LWP and IWP predicted by models with bin and bulk microphysical treatments is documented and attributed primarily to the assumed shape of ice PSD used in bulk schemes. Compared to the bin schemes that explicitly predict the PSD, schemes assuming exponential ice PSD underestimate ice growth by vapor deposition and overestimate mass-weighted fall speed leading to an underprediction of IWP by a factor of two in the considered case. Sensitivity tests indicate LWP and IWP are much closer to the bin model simulations when a modified shape factor which is similar to that predicted by bin model simulation is used in bulk scheme. These results demonstrate the importance of representation of ice PSD in determining the partitioning of liquid and ice and the longevity of mixed-phase clouds.

clouds↗