Ann M. Fridlind

Homogeneous Ice Nucleation in Subtropical and Tropical Convection and Its Influence on Cirrus Anvil Microphysics

Abstract This study uses a unique set of microphysical measurements obtained in a vigorous, convective updraft core at temperatures between −33° and −36°C, together with a microphysical model, to investigate the role of homogeneous ice nucleation in deep tropical convection and how it influences the microphysical properties of the associated cirrus anvils. The core and anvil formed along a sea-breeze front during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers–Florida Area Cirrus Experiment (CRYSTAL–FACE). The updraft core contained two distinct regions as traversed horizontally: the upwind portion of the core contained droplets of diameter 10–20 μm in concentrations of around 100 cm−3 with updraft speeds of 5–10 m s−1; the downwind portion of the core was glaciated with high concentrations of small ice particles and stronger updrafts of 10–20 m s−1. Throughout the core, rimed particles up to 0.6-cm diameter were observed. The anvil contained high concentrations of both small particles and ...

Intercomparison of cloud model simulations of Arctic mixed‐phase boundary layer clouds observed during SHEBA/FIRE‐ACE

An intercomparison of six cloud-resolving and large-eddy simulation models is presented. This case study is based on observations of a persistent mixed-phase boundary layer cloud gathered on 7 May, 1998 from the Surface Heat Budget of Arctic Ocean (SHEBA) and First ISCCP Regional Experiment - Arctic Cloud Experiment (FIRE-ACE). Ice nucleation is constrained in the simulations in a way that holds the ice crystal concentration approximately fixed, with two sets of sensitivity runs in addition to the baseline simulations utilizing different specified ice nucleus (IN) concentrations. All of the baseline and sensitivity simulations group into two distinct quasi-steady states associated with either persistent mixed-phase clouds or all-ice clouds after the first few hours of integration, implying the existence of multiple states for this case. These two states are associated with distinctly different microphysical, thermodynamic, and radiative characteristics. Most but not all of the models produce a persistent mixed-phase cloud qualitatively similar to observations using the baseline IN/crystal concentration, while small increases in the IN/crystal concentration generally lead to rapid glaciation and conversion to the all-ice state. Budget analysis indicates that larger ice deposition rates associated with increased IN/crystal concentrations have a limited direct impact on dissipation of liquid in these simulations. However, the impact of increased ice deposition is greatly enhanced by several interaction pathways that lead to an increased surface precipitation flux, weaker cloud top radiative cooling and cloud dynamics, and reduced vertical mixing, promoting rapid glaciation of the mixed-phase cloud for deposition rates in the cloud layer greater than about 122610 –5 gk g –1 s –1 for this case. These results indicate the critical importance of precipitation-radiative-dynamical interactions in simulating cloud phase, which have been neglected in previous fixed-dynamical parcel studies of the cloud phase parameter space. Large sensitivity to the IN/crystal concentration also suggests the need for improved understanding of ice nucleation and its parameterization in models.

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.

A FIRE-ACE/SHEBA Case Study of Mixed-Phase Arctic Boundary Layer Clouds: Entrainment Rate Limitations on Rapid Primary Ice Nucleation Processes

Observations of long-lived mixed-phase Arctic boundary layer clouds on 7 May 1998 during the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)‐Arctic Cloud Experiment(ACE)/SurfaceHeatBudgetoftheArcticOcean(SHEBA)campaignprovideauniqueopportunity to test understanding of cloud ice formation. Under the microphysically simple conditions observed (apparently negligible ice aggregation, sublimation, and multiplication), the only expected source of new ice crystals is activation of heterogeneous ice nuclei (IN) and the only sink is sedimentation. Large-eddy simulations with size-resolved microphysics are initialized with IN number concentration NIN measured above cloudtop,butdetailsofINactivationbehaviorareunknown.Ifactivatedrapidly(indeposition,condensation,

Evaluation of cloud‐resolving and limited area model intercomparison simulations using TWP‐ICE observations: 1. Deep convective updraft properties

Ten 3-D cloud-resolving model simulations and four 3-D limited area model simulations of an intense mesoscale convective system observed on 23–24 January 2006 during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) are compared with each other and with observed radar reflectivity fields and dual-Doppler retrievals of vertical wind speeds in an attempt to explain published results showing a high bias in simulated convective radar reflectivity aloft. This high-bias results from ice water content being large, which is a product of large, strong convective updrafts, although hydrometeor size distribution assumptions modulate the size of this bias. Making snow mass more realistically proportional to D2 rather than D3 eliminates unrealistically large snow reflectivities over 40 dBZ in some simulations. Graupel, unlike snow, produces high biased reflectivity in all simulations, which is partly a result of parameterized microphysics but also partly a result of overly intense simulated updrafts. Peak vertical velocities in deep convective updrafts are greater than dual-Doppler-retrieved values, especially in the upper troposphere. Freezing of liquid condensate, often rain, lofted above the freezing level in simulated updraft cores greatly contributes to these excessive upper tropospheric vertical velocities. The strongest simulated updraft cores are nearly undiluted, with some of the strongest showing supercell characteristics during the multicellular (presquall) stage of the event. Decreasing horizontal grid spacing from 900 to 100 m slightly weakens deep updraft vertical velocity and moderately decreases the amount of condensate aloft but not enough to match observational retrievals. Therefore, overly intense simulated updrafts may additionally be a product of unrealistic interactions between convective dynamics, parameterized microphysics, and large-scale model forcing that promote different convective strengths than observed.

Control of deep convection by sub-cloud lifting processes: the ALP closure in the LMDZ5B general circulation model

Recently, a new conceptual framework for deep convection scheme triggering and closure has been developed and implemented in the LMDZ5B general circulation model, based on the idea that deep convection is controlled by sub-cloud lifting processes. Such processes include boundary-layer thermals and evaporatively-driven cold pools (wakes), which provide an available lifting energy that is compared to the convective inhibition to trigger deep convection, and an available lifting power (ALP) at cloud base, which is used to compute the convective mass flux assuming the updraft vertical velocity at the level of free convection. While the ALP closure was shown to delay the local hour of maximum precipitation over land in better agreement with observations, it results in an underestimation of the convection intensity over the tropical ocean both in the 1D and 3D configurations of the model. The specification of the updraft vertical velocity at the level of free convection appears to be a key aspect of the closure formulation, as it is weaker over tropical ocean than over land and weaker in moist mid-latitudes than semi-arid regions. We propose a formulation making this velocity increase with the level of free convection, so that the ALP closure is adapted to various environments. Cloud-resolving model simulations of observed oceanic and continental case studies are used to evaluate the representation of lifting processes and test the assumptions at the basis of the ALP closure formulation. Results favor closures based on the lifting power of sub-grid sub-cloud processes rather than those involving quasi-equilibrium with the large-scale environment. The new version of the model including boundary-layer thermals and cold pools coupled together with the deep convection scheme via the ALP closure significantly improves the representation of various observed case studies in 1D mode. It also substantially modifies precipitation patterns in the full 3D version of the model, including seasonal means, diurnal cycle and intraseasonal variability.

Evaluation of Hydrometeor Phase and Ice Properties in Cloud-Resolving Model Simulations of Tropical Deep Convection Using Radiance and Polarization Measurements

Satellite measurements are used to evaluate the glaciation, particle shape, and effective radius in cloudresolvingmodel simulationsof tropicaldeepconvection.Multidirectional polarized reflectancesconstrainthe ice crystal geometry and the thermodynamic phase of the cloud tops, which in turn are used to calculate nearinfrared reflectances so as to constrain the simulated ice effective radius, thereby avoiding inconsistencies between retrieval algorithms and model simulations. Liquid index values derived from Polarization and Directionality of the Earth’s Reflectances (POLDER) measurements indicate only ice-topped clouds at brightness temperatures (BTs) lower than 2408C, only liquid clouds at BT .2 208C, and both phases occurringattemperaturesinbetween.Liquidindexvaluescalculatedfrommodelsimulationsgenerallyreveal too many ice-topped clouds at BT .2 208C. The model assumption of platelike ice crystals with an aspect ratioof0.7isfoundconsistentwithPOLDERmeasurementsforBT,2408Cwhenveryroughicecrystalsare assumed, leading to an asymmetry parameter of 0.74, whereas measurements indicate more extreme aspect ratios of ;0.15 at higher temperatures, yielding an asymmetry parameter of 0.84. MODIS-retrieved ice effective radii are found to be 18‐28 m ma t BT ,2 408C, but biased low by about 5 mm owing primarily to the assumption of pristine crystals in the retrieval. Simulated 2.13-mm reflectances at BT ,2 408C are found to be about 0.05‐0.1 too large comparedto measurements, suggesting that model-simulated effective radii are 7‐15 mm too small. Two simulations with contrasting ice nucleation schemes showed little difference in simulated effective radii at BT ,2 408C, indicating that homogeneous nucleation is dominating in the simulations. Changes around 2408C in satellite observations suggest a change in cloud-top ice shape and/or size in natural deep convection possibly related to a change in the freezing mechanism.

A Flexible Parameterization for Shortwave Optical Properties of Ice Crystals

A parameterization is presented that provides extinction cross section se, single-scattering albedo v, and asymmetryparametergof ice crystals for any combination of volume, projected area, aspect ratio, and crystal distortionatanywavelengthintheshortwave.Similartopreviousparameterizations,theschememakesuseof 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, se is always twice the particle projectedarea.Itisshownthatvislargelydeterminedbythenewlydefinedabsorptionsizeparameterandthe 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 v 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 4mm, revealing absolute differences with reference calculations of v 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,andabsorptanceare1.4%,1.1%,and3.4%,respectively.Somepracticalapplications of the parameterization in atmospheric models are highlighted.

Variation of ice crystal size, shape, and asymmetry parameter in tops of tropical deep convective clouds

The variation of ice crystal properties in the tops of deep convective clouds off the north coast of Australia is analyzed. Cloud optical thickness, ice effective radius, aspect ratio of ice crystal components, crystal distortion parameter and asymmetry parameter are simultaneously retrieved from combined measurements of the Moderate Resolution Imaging Spectroradiometer (MODIS) and Polarization and Directionality of the Earths Reflectances (POLDER) satellite instruments. The data are divided into periods with alternating weak and strong convection. Mostly plate-like particle components with aspect ratios closer to unity and lower asymmetry parameters characterize strongly convective periods, while weakly convective periods generally show lower aspect ratios, relatively more column-like shapes and somewhat greater asymmetry parameters. Results for strongly convective periods show that, with increasing cloud top temperature, the distortion parameter generally decreases, while the asymmetry parameter and effective radius increase. For one of the strongly convective periods, the rate at which effective radii increase with cloud top temperature is more than double that of the other periods, while the temperature dependence of the other microphysical quantities for this period is substantially weaker. Atmospheric state analysis indicates that these differences are concurrent with differences in middle-to-upper tropospheric zonal wind shear. The observed variation of microphysical properties may have significant effects on the shortwave radiative fluxes and cloud absorption associated with deep convection. Additionally, MODIS collection 5 effective radii are estimated to be biased small with an artificially narrow range. Collection 6 products are expected to have less severe biases that depend on cloud top temperature and atmospheric conditions.

Evaluation of cloud-resolving and limited area model intercomparison simulations using TWP-ICE observations: 2. Precipitation microphysics

Ten 3-D cloud-resolving model (CRM) simulations and four 3-D limited area model (LAM) simulations of an intense mesoscale convective system observed on 23-24 January 2006 during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) are compared with each other and with observations and retrievals from a scanning polarimetric radar, colocated UHF and VHF vertical profilers, and a Joss-Waldvogel disdrometer in an attempt to explain a low bias in simulated stratiform rainfall. Despite different forcing methodologies, similar precipitation microphysics errors appear in CRMs and LAMs with differences that depend on the details of the bulk microphysics scheme used. One-moment schemes produce too many small raindrops, which biases Doppler velocities low, but produces rainwater contents (RWCs) that are similar to observed. Two-moment rain schemes with a gamma shape parameter (mu) of 0 produce excessive size sorting, which leads to larger Doppler velocities than those produced in one-moment schemes but lower RWCs. Two-moment schemes also produce a convective median volume diameter distribution that is too broad relative to observations and, thus, may have issues balancing raindrop formation, collision-coalescence, and raindrop breakup. Assuming a mu of 2.5 rather than 0 for the raindrop size distribution improves one-moment scheme biases, and allowing mu to have values greater than 0 may improve excessive size sorting in two-moment schemes. Underpredicted stratiform rain rates are associated with underpredicted ice water contents at the melting level rather than excessive rain evaporation, in turn likely associated with convective detrainment that is too high in the troposphere and mesoscale circulations that are too weak. A limited domain size also prevents a large, well-developed stratiform region like the one observed from developing in CRMs, although LAMs also fail to produce such a region.