Search

The Tropical Warm Pool.International Cloud Experiment (TWP ]ICE) provided extensive observational data sets designed to initialize, force, and constrain atmospheric model simulations. In this first of a two ]part study, precipitation and cloud structures within nine cloud ]resolving model simulations are compared with scanning radar reflectivity and satellite infrared brightness temperature observations during an active monsoon period from 19 to 25 January 2006. Seven of nine simulations overestimate convective area by 20% or more leading to general overestimation of convective rainfall. This is balanced by underestimation of stratiform rainfall by 5% to 50% despite overestimation of stratiform area by up to 65% because of a preponderance of very low stratiform rain rates in all simulations. All simulations fail to reproduce observed radar reflectivity distributions above the melting level in convective regions and throughout the troposphere in stratiform regions. Observed precipitation ]sized ice reaches higher altitudes than simulated precipitation ]sized ice despite some simulations that predict lower than observed top ]of ]atmosphere infrared brightness temperatures. For the simulations that overestimate radar reflectivity aloft, graupel is the cause with one ]moment microphysics schemes whereas snow is the cause with two ]moment microphysics schemes. Differences in simulated radar reflectivity are more highly correlated with differences in mass mean melted diameter (Dm) than differences in ice water content. Dm is largely dependent on the mass ]dimension relationship and gamma size distribution parameters such as size intercept (N0) and shape parameter (m). Having variable density, variable N0, or m greater than zero produces radar reflectivities closest to those observed.

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).

Convective systems are an important mechanism in the transport of boundary layer air into the upper troposphere. The Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) campaign, in July 2002, was developed as a comprehensive atmospheric mission to improve knowledge of subtropical cirrus systems and their roles in regional and global climate. In situ measurements of carbon monoxide (CO), water vapor (H2Ov), and total water (H2Ot) aboard NASA's WB-57F aircraft and CO aboard the U.S. Navy's Twin Otter aircraft were obtained to study the role of convective transport. Three flights sampled convective outflow on 11, 16 and 29 July found varying degrees of CO enhancement relative to the free troposphere. A cloud-resolving model used the in situ observations and meteorological fields to study these three systems. Several methods of filtering the observations were devised here using ice water content, relative humidity with respect to ice, and particle number concentration as a means to statistically sample the model results to represent the flight tracks. A weighted histogram based on ice water content observations was then used to sample the simulations for the three flights. In addition, because the observations occurred in the convective outflow cirrus and not in the storm cores, the model was used to estimate the maximum CO within the convective systems. In general, anvil-level air parcels contained an estimated 20-40% boundary layer air in the analyzed storms.

The transport of H2O and HDO within deep convection is investigated with 3-D large eddy simulations (LES) using bin microphysics. The lofting and sublimation of HDO-rich ice invalidate the Rayleigh fractionation model of isotopologue distribution within deep convection. Bootstrapping the correlation of the ratio of HDO to H2O (deltaD) to water vapor mixing ratio (q(sub v)) through a sequence of convective events produced non-Rayleigh correlations resembling observations. These results support two mechanisms for stratospheric entry. Deep convection can inject air with water vapor of stratospheric character directly into the tropical transition layer (TTL). Alternatively, moister air detraining from convection may be dehydrated via cirrus formation n the TTL to produce stratospheric water vapor. Significant production of subsaturated air in the TTL via convective dehydration is not observed in these simulations, nor is it necessary to resolve the stratospheric isotope paradox.

The 2004 Multi-Phase Arctic Cloud Experiment (M-PACE) IOP at the ARM NSA site focused on measuring the properties of autumn transition-season arctic stratus and the environmental conditions controlling them, including concentrations of heterogeneous ice nuclei. Our work aims to use a large-eddy simulation (LES) code with embedded size-resolved aerosol and cloud microphysics to identify factors controlling multi-phase arctic stratus. Our preliminary simulations of autumn transition-season clouds observed during the 1994 Beaufort and Arctic Seas Experiment (BASE) indicated that low concentrations of ice nuclei, which were not measured, may have significantly lowered liquid water content and thereby stabilized cloud evolution. However, cloud drop concentrations appeared to be virtually immune to changes in liquid water content, indicating an active Bergeron process with little effect of collection on drop number concentration. We will compare these results with preliminary simulations from October 8-13 during MPACE. The sensitivity of cloud properties to uncertainty in other factors, such as large-scale forcings and aerosol profiles, will also be investigated. Based on the LES simulations with M-PACE data, preliminary results from the NASA GlSS single-column model (SCM) will be used to examine the sensitivity of predicted cloud properties to changing cloud drop number concentrations for multi-phase arctic clouds. Present parametrizations assumed fixed cloud droplet number concentrations and these will be modified using M-PACE data.

Activation of some fraction of increased concentrations of sub-micron soluble aerosol particles lead to enhanced cloud droplet concentrations and hence smaller droplets, increasing their total cross sectional area and thus reflecting solar radiation more efficiently (the Twomey, or first indirect, effect). However, because of competition during condensational growth, droplet distributions tend to broaden as numbers increase, reducing the sensitivity of cloud albedo to droplet concentration on the order of 10%. Also, smaller droplets less effectively produce drizzle through collisions and coalescence, and it is widely expected (and found in large-scale models) that decreased precipitation leads to clouds with more cloud water on average (the so-called cloud lifetime, or second indirect, effect). Much of the uncertainty regarding the overall indirect aerosol effect stems from inadequate understanding of such changes in cloud water. Detailed simulations based on FIRE-I, ASTEX, and DYCOMS-II conditions show that suppression of precipitation from increased droplet concentrations leads to increased cloud water only when sufficient precipitation reaches the surface, a condition favored when the overlying air is-humid or droplet concentrations are very low. Otherwise, aerosol induced suppression of precipitation enhances entrainment of overlying dry air, thereby reducing cloud water and diminishing the indirect climate forcing.

The intense heat from forest fires can generate explosive deep convective cloud systems that inject pollutants to high altitudes. Both satellite and high-altitude aircraft measurements have documented cases in which these pyrocumulonimbus clouds inject large amounts of smoke well into the stratosphere (Fromm and Servranckx 2003; Jost et al. 2004). This smoke can remain in the stratosphere, be transported large distances, and affect lower stratospheric chemistry. In addition recent in situ measurements in pyrocumulus updrafts have shown that the high concentrations of smoke particles have significant impacts on cloud microphysical properties. Very high droplet number densities result in delayed precipitation and may enhance lightning (Andrew et al. 2004). Presumably, the smoke particles will also lead to changes in the properties of anvil cirrus produces by the deep convection, with resulting influences on cloud radiative forcing. In situ sampling near the tops of mature pyrocumulonimbus is difficult due to the high altitude and violence of the storms. In this study, we use large eddy simulations (LES) with size-resolved microphysics to elucidate physical processes in pyrocumulonimbus clouds.

Numerical models that resolve cloud particles into discrete mass size distributions on an Eulerian grid provide a uniquely powerful means of studying the closely coupled interaction of aerosols, cloud microphysics, and transport that determine cloud properties and evolution. However, such models require many experimentally derived paramaterizations in order to properly represent the complex interactions of droplets within turbulent flow. Many of these parameterizations remain poorly quantified, and the numerical methods of solving the equations for temporal evolution of the mass size distribution can also vary considerably in terms of efficiency and accuracy. In this work, we compare results from two size-resolved microphysics models that employ various widely-used parameterizations and numerical solution methods for several aspects of stochastic collection.

Increases in cloud cover and condensed water contribute more than half of the indirect aerosol effect in an ensemble of general circulation model (GCM) simulations estimating the global radiative forcing of anthropogenic aerosols. We use detailed simulations of marine stratocumulus clouds and airborne observations of ship tracks to show that increases in cloud cover and condensed water in reality are far less than represented by the GCM ensemble. Our results offer an explanation for recent simplified inverse climate calculations indicating that indirect aerosol effects are greatly exaggerated in GCMs.

The indirect radiative forcing of aerosols is poorly constrained by the observational data underlying the simple cloud parameterizations in GCMs. signal of cloud response to increased aerosol concentrations from meteorological noise. Recent satellite observations indicate a significant decrease of cloud water in ship tracks, in contrast to an ensemble of in situ measurements showing no average change in cloud water relative to the surrounding clouds. Both results contradict the expectation of cloud water increasing in polluted clouds. We find through large-eddy simulations of stratocumulus that the trend in the satellite data is likely an artifact of sampling only overcast clouds. The simulations instead show cloud cover increasing with droplet concentrations. The simulations also show that increases in cloud water from suppressing drizzle by increased droplet concentrations are favored at night or at extremely low droplet concentrations. At typical droplet concentrations we find that the Twomey effect on cloud albedo is amplified very little by the secondary indirect effect of drizzle suppression, largely because the absorption of solar radiation by cloud water reduces boundary-layer mixing in the daytime and thereby restricts any possible increase in cloud water from drizzle suppression. The cloud and boundary layer respond to radiative heating variations on a time scale of hours, and on longer time scales respond to imbalances between large-scale horizontal advection and the entrainment of inversion air. We analyze the co-varying response of cloud water, cloud thickness, width of droplet size distributions, and dispersion of the optical depth, as well as the overall response of cloud albedo, to changes in droplet concentrations. We also dissect the underlying physical mechanisms through sensitivity studies. Ship tracks represent an ideal natural laboratory to extricate the

Recent satellite observations of ship tracks surprisingly indicate that cloud water decreases with increasing droplet concentrations. However, we find by analyzing detailed simulations of stratocumulus that the reported trend is likely an artifact of sampling, only overcast clouds. The simulations instead show cloud coverage increasing with droplet concentrations, accounting for 25% of cloud albedo increase at moderate droplet concentrations. Our simulations also show that increases in cloud water from drizzle suppression (by increasing droplet concentrations) are favored only at night or at extremely low droplet concentrations, suggesting that the indirect aerosol forcing is overestimated in climate change projections by many general circulation models.

Recent satellite observations indicate a significant decrease of cloud water in ship tracks, in contrast to an ensemble of in situ ship-track measurements that show no average change in cloud water relative to the surrounding clouds. We find through large-eddy simulations of stratocumulus that the trend in the satellite data is likely an artifact of sampling only overcast clouds. The simulations instead show cloud cover increasing with droplet concentrations. Our simulations also show that increases in cloud water from drizzle suppression (by increasing droplet concentrations) are favored at night or at extremely low droplet concentrations.

The processes controlling production of ice crystals in deep, rapidly ascending convective columns are poorly understood due to the difficulties involved with either modeling or in situ sampling of these violent clouds. A large number of ice crystals are no doubt generated when droplets freeze at about -40 C. However, at higher levels, these crystals are likely depleted due to precipitation and detrainment. As the ice surface area decreases, the relative humidity can increase well above ice saturation, resulting in bursts of ice nucleation. We will present simulations of these processes using a large-eddy simulation model with detailed microphysics. Size bins are included for aerosols, liquid droplets, ice crystals, and mixed-phase (ice/liquid) hydrometers. Microphysical processes simulated include droplet activation, freezing, melting, homogeneous freezing of sulfate aerosols, and heterogeneous ice nucleation. We are focusing on the importance of ice nucleation events in the upper part of the cloud at temperatures below -40 C. We will show that the ultimate evolution of the cloud in this region (and the anvil produced by the convection) is sensitive to these ice nucleation events, and hence to the composition of upper tropospheric aerosols that get entrained into the convective column.

We present a method to calculate vertical profiles of particle size distributions in condensation clouds of giant planets and brown dwarfs. The method assumes a balance between turbulent diffusion and precipitation in horizontally uniform cloud decks. Calculations for the Jovian ammonia cloud are compared with previous methods. An adjustable parameter describing the efficiency of precipitation allows the new model to span the range of predictions from previous models. Calculations for the Jovian ammonia cloud are found to be consistent with observational constraints. Example calculations are provided for water, silicate, and iron clouds on brown dwarfs and on a cool extrasolar giant planet.

A detailed 1D model of the stratocumulus-topped marine boundary layer is described. The model has three coupled components: a microphysics module that resolves the size distributions of aerosols and cloud droplets, a turbulence module that treats vertical mixing between layers, and a multiple wavelength radiative transfer module that calculates radiative heating rates and cloud optical properties. The results of a 12-h model simulation reproduce reasonably well the bulk thermodynamics, microphysical properties, and radiative fluxes measured in an approx. 500-m thick, summertime marine stratocumulus cloud layer by Nicholls. However, in this case, the model predictions of turbulent fluxes between the cloud and subcloud layers exceed the measurements. Results of model simulations are also compared to measurements of a marine stratus layer made under gate conditions and with measurements of a high, thin marine stratocumulus layer. The variations in cloud properties are generally reproduced by the model, although it underpredicts the entrainment of overlying air at cloud top under gale conditions. Sensitivities of the model results are explored. The vertical profile of cloud droplet concentration is sensitive to the lower size cutoff of the droplet size distribution due to the presence of unactivated haze particles in the lower region of the modeled cloud. Increases in total droplet concentrations do not always produce less drizzle and more cloud water in the model. The radius of the mean droplet volume does not correlate consistently with drizzle, but the effective droplet radius does. The greatest impacts on cloud properties predicted by the model are produced by halving the width of the size distribution of input condensation nuclei and by omitting the effect of cloud-top radiative cooling on the condensational growth of cloud droplets. The omission of infrared scattering produces noticeable changes in cloud properties. The collection efficiencies for droplets less than 30-micron radius, and the value of the accommodation coefficient for condensational droplet growth, have noticeable effects on cloud properties. The divergence of the horizontal wind also has a significant effect on a 12-h model simulation of cloud structure. Conclusions drawn from the model are tentative because of the limitations of the 1D model framework. A principal simplification is that the model assumes horizontal homogeneity, and, therefore, does not resolve updrafts and downdrafts. Likely consequences of this simplification include overprediction of the growth of droplets by condensation in the upper region of the cloud, underprediction of droplet condensational growth in the lower region of the cloud, and under-prediction of peak supersaturations.