Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cloud microphysics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗

Stratus and Stratocumulus Cloud Microphysics and Drizzle Relationships With CCN Modality

High resolution extended-range cloud condensation nuclei (CCN) spectral comparisons with cloud microphysics and drizzle of the Physics of Stratocumulus Tops (POST) field experiment confirmed results in the Marine Stratus/Stratocumulus Experiment (MASE). Both of these stratus cloud projects demonstrated that bimodal CCN spectra typically caused by cloud processing were associated with clouds that exhibited higher concentrations of smaller droplets with narrower distributions and less drizzle than clouds associated with unimodal CCN spectra. Resulting brighter clouds and increased cloudiness could enhance both indirect aerosol effects (IAE). These stratus findings are opposite of analogous measurements in two cumulus cloud projects, which showed bimodal CCN associated with fewer larger droplets more broadly distributed and with more drizzle than clouds associated with unimodal CCN. Resulting reduced cumulus brightness and cloudiness could reduce both IAE. Physics of Stratocumulus Tops (POST) flights in air masses with higher CCN concentrations, N CCN , showed more extremes of the stratus characteristics. However, POST flights with lower N CCN showed opposite droplet characteristics similar to the cumulus clouds, yet still showed less drizzle in clouds associated with bimodal CCN, but not as much less as the flights with higher N CCN . Since all MASE clouds were in polluted air masses, while the two cumulus projects were in clean air masses we deduce from these four projects that both the dynamic stratus/cumulus differences (vertical wind) and N CCN are responsible for the microphysics and drizzle differences among these projects. This is because the clean POST characteristics are a hybrid between MASE/POST high N CCN and the two cumulus projects.

bimodality↗

Karhunen-Loeve expansion analysis of uncertainties in cloud microphysical property retrievals

Here, this study proposes a methodology of quantifying uncertainties for cloud retrievals on model resolution to facilitate the comparison with model outputs. Primary component analysis is applied to reduce the dimension of random variables (up to a factor of 50) and reveal the cross correlations in the input data, making large sampling computationally feasible and uncertainty quantification accurate and reliable. Our approach has the capability of parameterizing input uncertainties and attributing the uncertainties in the retrieval output to each individual source, which allows sensitivity analysis of cloud retrieval algorithms and provides directions for improving observation instruments as well as strategies. We applied the method to characterize uncertainties in cloud ice water content (IWC) retrieved from the DOE Atmospheric Radiation Measurement (ARM) programs baseline cloud microphysical retrieval algorithm (MICROBASE). We test it with a selected ice cloud case observed on 9 March 2000 at the ARM Southern Great Plains site during its 2000 cloud intensive observing period. The test results indicate that (1) uncertainties in the output retrieved by MICROBASE are comparable amongst different retrievals; (2) The mean values obtained by our UQ method are closer to the aircraft data with less errors compared to the direct ensemble average; (3) Ice water path (IWP) generally incurred larger uncertainty in optically thin ice clouds and there was more variability in vertical in the retrieved IWC; and (4) Uncertainties in the output are mainly due to the interactions among different modes of ARM radar pro files.

54 ENVIRONMENTAL SCIENCES↗

Cloud microphysical response to entrainment and mixing is locally inhomogeneous and globally homogeneous: Evidence from the lab

Entrainment of dry air into clouds strongly influences cloud optical and precipitation properties and the response of clouds to aerosol perturbations. The response of cloud droplet size distributions to entrainment–mixing is examined in the Pi convection-cloud chamber that creates a turbulent, steady-state cloud. The experiments are conducted by injecting dry air with temperature (T e ) and flow rate (Q e ) through a flange in the top boundary, into the otherwise well-mixed cloud, to mimic the entrainment–mixing process. Due to the large-scale circulation, the downwind region is directly affected by entrained dry air, whereas the upwind region is representative of the background conditions. Droplet concentration (Cn) and liquid water content (L) decrease in the downwind region, but the difference in the mean diameter of droplets (D m ) is small. The shape of cloud droplet size distributions relative to the injection point is unchanged, to within statistical uncertainty, resulting in a signature of inhomogeneous mixing, as expected for droplet evaporation times small compared to mixing time scales. As T e and Q e of entrained air increase, however, Cn, L, and D m of the whole cloud system decrease, resulting in a signature of homogeneous mixing. The apparent contradiction is understood as the cloud microphysical responses to entrainment and mixing differing on local and global scales: locally inhomogeneous and globally homogeneous. This implies that global versus local sampling of clouds can lead to seemingly contradictory results for mixing, which informs the long-standing debate about the microphysical response to entrainment and the parameterization of this process for coarse-resolution models.

54 ENVIRONMENTAL SCIENCES↗

Version 4 CALIPSO Imaging Infrared Radiometer ice and liquid water cloud microphysical properties – Part I: The retrieval algorithms

Following the release of the version 4 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a new version (version 4; V4) of the CALIPSO Imaging Infrared Radiometer (IIR) Level 2 data products has been developed. The IIR Level 2 data products include cloud effective emissivities and cloud microphysical properties such as effective diameter and ice or liquid water path estimates. Dedicated retrievals for water clouds were added in V4, taking advantage of the high sensitivity of the IIR retrieval technique to small particle sizes. This paper (Part I) describes the improvements in the V4 algorithms compared to those used in the version 3 (V3) release, while results will be presented in a companion (Part II) paper. The IIR Level 2 algorithm has been modified in the V4 data release to improve the accuracy of the retrievals in clouds of very small (close to 0) and very large (close to 1) effective emissivities. To reduce biases at very small emissivities that were made evident in V3, the radiative transfer model used to compute clear-sky brightness temperatures over oceans has been updated and tuned for the simulations using Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data to match IIR observations in clear-sky conditions. Furthermore, the clear-sky mask has been refined compared to V3 by taking advantage of additional information now available in the V4 CALIOP 5 km layer products used as an input to the IIR algorithm. After sea surface emissivity adjustments, observed and computed brightness temperatures differ by less than ±0.2 K at night for the three IIR channels centered at 08.65, 10.6, and 12.05 µm, and inter-channel biases are reduced from several tens of Kelvin in V3 to less than 0.1 K in V4. We have also improved retrievals in ice clouds having large emissivity by refining the determination of the radiative temperature needed for emissivity computation. The initial V3 estimate, namely the cloud centroid temperature derived from CALIOP, is corrected using a parameterized function of temperature difference between cloud base and top altitudes, cloud absorption optical depth, and CALIOP multiple scattering correction factor. As shown in Part II, this improvement reduces the low biases at large optical depths that were seen in V3 and increases the number of retrievals. As in V3, the IIR microphysical retrievals use the concept of microphysical indices applied to the pairs of IIR channels at 12.05 and 10.6 µm and at 12.05 and 08.65 µm. The V4 algorithm uses ice look-up tables (LUTs) built using two ice habit models from the recent “TAMUice2016” database, namely the single-hexagonal-column model and the eight-element column aggregate model, from which bulk properties are synthesized using a gamma size distribution. Four sets of effective diameters derived from a second approach are also reported in V4. Here, the LUTs are analytical functions relating microphysical index applied to IIR channels 12.05 and 10.6 µm and effective diameter as derived from in situ measurements at tropical and midlatitudes during the Tropical Composition, Cloud, and Climate Coupling (TC4) and Small Particles in Cirrus Science and Operations Plan (SPARTICUS) field experiments.

54 ENVIRONMENTAL SCIENCES↗

Signatures of aerosol-cloud interactions in GiOcean: a coupled global reanalysis with two-moment cloud microphysics

Aerosols influence the Earth's radiative balance through direct interactions with radiation and by affecting cloud properties. Anthropogenic aerosols have led to cooling during the industrial era through aerosol–cloud interactions (ACI), including aerosol effects on cloud microphysical properties and the subsequent adjustments. However, large uncertainties remain in Earth system models (ESMs) regarding the magnitude of this cooling. In part, ESMs substantially disagree on cloud properties, thermodynamics, the hydrological cycle, and general circulation. Reanalysis provides a useful avenue for exploring the impact of ACI on clouds and radiation because its atmosphere is forced to match realistic conditions through the assimilation of observations. Here, we explore the impact of ACI on clouds in the GiOcean reanalysis – the first to incorporate aerosol-cloud interactions. We contrast variables important for ACI between GiOcean and satellite observations and develop 2-dimensional lookup tables of ACI for both using a source-sink budget perspective to attribute the changes in cloud droplet number (Nd) and liquid water path (LWP) to aerosol and meteorology. A compositing analysis using lookup tables shows that GiOcean captures key aspects of aerosol–cloud–precipitation interactions, including (1) activation of aerosol into cloud droplets, (2) effective precipitation scavenging of Nd, (3) suppression of precipitation by high Nd in regions with heavy aerosol emissions. In contrast, satellite observations do not exhibit clear patterns for processes (2) and (3). Random Forest analysis shows that interannual variability in Nd and LWP over the Northern Hemisphere ocean in GiOcean is primarily driven by precipitation, consistent with satellite observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Sensitivity Study of Radiative Fluxes at the Top of Atmosphere to Cloud-Microphysics and Aerosol Parameters in the Community Atmosphere Model CAM5

In this study, we investigated the sensitivity of net radiative fluxes (FNET) at the top of atmosphere (TOA) to 16 selected uncertain parameters mainly related to the cloud microphysics and aerosol schemes in the Community Atmosphere Model version 5 (CAM5). We adopted a quasi-Monte Carlo (QMC) sampling approach to effectively explore the high dimensional parameter space. The output response variables (e.g., FNET) were simulated using CAM5 for each parameter set, and then evaluated using generalized linear model analysis. In response to the perturbations of these 16 parameters, the CAM5-simulated global annual mean FNET ranges from -9.8 to 3.5 W m-2 compared to the CAM5-simulated FNET of 1.9 W m-2 with the default parameter values. Variance-based sensitivity analysis was conducted to show the relative contributions of individual parameter perturbation to the global FNET variance. The results indicate that the changes in the global mean FNET are dominated by those of cloud forcing (CF) within the parameter ranges being investigated. The size threshold parameter related to auto-conversion of cloud ice to snow is confirmed as one of the most influential parameters for FNET in the CAM5 simulation. The strong heterogeneous geographic distribution of FNET variation shows parameters have a clear localized effect over regions where they are acting. However, some parameters also have non-local impacts on FNET variance. Although external factors, such as perturbations of anthropogenic and natural emissions, largely affect FNET variations at the regional scale, their impact is weaker than that of model internal parameters in terms of simulating global mean FNET in this study. The interactions among the 16 selected parameters contribute a relatively small portion of the total FNET variations over most regions of the globe. This study helps us better understand the CAM5 model behavior associated with parameter uncertainties, which will aid the next step of reducing model uncertainty via calibration of uncertain model parameters with the largest sensitivity.

aerosol↗

An evaluation of size-resolved cloud microphysics scheme numerics for use with radar observations Part II: Condensation and evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part I of this study addressed collision–coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision–coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Here, large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

54 ENVIRONMENTAL SCIENCES↗

Investigation of Liquid Cloud Microphysical Properties of Deep Convective Systems: 2. Parameterization of Raindrop Size Distribution and its Application for Convective Rain Estimation

The liquid cloud microphysical properties for stratiform rain (SR) have been investigated in Part 1 of this series. Since convective rain (CR) has characteristics in raindrop size distribution (DSD) and precipitation that are distinct from SR, we investigate the CR properties in this study by using 20 hr of CR samples collected by 17 Automatic Parsivel Units disdrometers during the Midlatitude Continental Convective Clouds Experiment over the Atmospheric Radiation Measurement Southern Great Plains site. A full spectrum of DSD is constructed based on a total of 23 size channels (0.321 to 9.785 mm), and both Gamma and Exponential fitted functions are applied to extract the DSD shape parameters. Compared to SR properties, CR has distinct features including broader size range and narrower exponential slope parameter (λ E ). These results indicate that the assumption of constant exponential intercept parameter ( N 0E ) is inappropriate for CR rain rate estimates. Additionally, the subsetting scheme of the CR DSD spectra also has a strong impact on the Gamma/Exponential functions. Therefore, a new CR DSD parameterization scheme is developed by choosing the appropriate spectra subset with the constraint of rain rate and using the constant λ E instead of constant N 0E . With the input of the radar reflectivity at the lowest observed height, the newly calculated CR rain rates match well with the collocated surface rain gauge measurements (127 Mesonet stations and 17 Automatic Parsivel Units), while the rain rates calculated using traditional Z-R relationship are 3–4 times larger, indicating that constant λ E is a better assumption for CR DSD.

58 GEOSCIENCES↗

Bulk cloud microphysical properties as seen from numerical simulation and remote sensing products: case study of a hailstorm event over the La Plata Basin

Hailstorms develop over the La Plata Basin, in south-eastern South America, more often during later winter and early austral spring, between September and October. These systems have significant socioeconomic impacts over the region. Thus, a better understanding of how atmospheric drivers modulate the formation of hailstorms is important to improve the forecast of such phenomena. In this study, we selected a hailstorm event observed over the eastern La Plata Basin during 14–15 July 2016 to evaluate the performance of the Brazilian developments on the Regional Atmospheric Modelling System (BRAMS) model. The ability of the model in simulating cloud microphysical properties was evaluated by comparing simulations driven by different global forcings against in situ and remote sensing observations. The model results showed good skill in capturing the basic characteristics of the thunderstorm, particularly in terms of the spatial distribution of hydrometeors. The simulated spatial distribution of hail covers locations where hail fall was reported. The BRAMS simulations suggest that, despite relatively low values of the convective available potential energy (CAPE) (700–1000 J kg -1 ), environments with strong 0–8-km bulk shear (60–70 kt, ~30.9–36.0 m s –1 ) can promote the formation of ice clouds and hail fall over the eastern La Plata Basin. To be more conclusive, however, further research is needed to understand how different combinations of CAPE and shear affect hail formation over the region.

54 ENVIRONMENTAL SCIENCES↗

DOE Early Career Award # DE-SC0012152: Understanding the Roles of Cloud Microphysics and Land Surface Coupling Feedbacks in Multi-Scale Predictions of Central US Summer Hydroclimate (Final Report)

This summarizes technical highlights from the DOE Early Career Project entitled “Understanding the Roles of Cloud Microphysics and Land Surface Coupling Feedbacks in Multi-Scale Predictions of Central US Summer Hydroclimate” led by PI Mike Pritchard at the University of California. The project’s original aims were to investigate trade-offs of cloud superparameterization such as the effects of microphysical assumptions on mesoscale convective systems (Section I) and the representation of land-atmosphere coupling under explicit convection (Section II). Under the same broader theme of water cycle dynamics, superparameterization, and land interaction, the project evolved to further investigate the climate dynamics of irrigation, vegetation and streamflow (Section III), the dynamics of the Madden- Julian Oscillation (Section IV), as well as the use of modern machine learning to replace traditional sub grid approaches for climate simulation (Section V). Twenty peer-reviewed journal publications were produced as a result of this project (Section VI).

54 ENVIRONMENTAL SCIENCES↗

Two-moment Bulk Stratiform Cloud Microphysics in the Grid-point Atmospheric Model of IAP LASG (GAMIL)

A two-moment bulk stratiform microphysics scheme, including recently developed physically-based droplet activation/ice nucleation parameterizations has been implemented into the Grid-point Atmospheric Model of IAP LASG (GAMIL) as an effort to enhance the model capability for studying aerosol indirect effects. Unlike the previous one-moment cloud microphysics scheme, the new scheme produces reasonable representation of cloud particle size and number concentration. This scheme captures the observed spatial variations in cloud droplet number concentrations. Simulated ice crystal number concentrations in cirrus clouds qualitatively agree with in-situ observations. The longwave and shortwave cloud forcing are in better agreement with observations. Sensitivity tests show that the column cloud droplet number concentrations calculated from two different droplet activation parameterizations are similar. However, ice crystal number concentration in mixed-phased clouds is sensitive to different heterogeneous freezing formulations. The simulation with high ice crystal number concentration in mixed-phase clouds has less liquid water path and weaker cloud forcing. Furthermore, ice crystal number concentration in cirrus clouds is sensitive to different ice nucleation parameterizations. Sensitivity tests also suggest that impact of pre-existing ice crystals on homogeneous freezing in old clouds should be taken into account.

(GAMIL)↗

Asymmetries in cloud microphysical properties ascribed to sea ice leads via water vapour transport in the central Arctic

To investigate the influence of sea ice openings like leads on wintertime Arctic clouds, the air mass transport is exploited as a heat and humidity feeding mechanism which can modify Arctic cloud properties. Cloud microphysical properties in the central Arctic are analysed as a function of sea ice conditions during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in 2019–2020. The Cloudnet classification algorithm is used to characterize the clouds based on remote sensing observations and the atmospheric thermodynamic state from the observatory on board the research vessel (RV) Polarstern. To link the sea ice conditions around the observational site with the cloud observations, the water vapour transport (WVT) being conveyed towards RV Polarstern has been utilized as a mechanism to associate upwind sea ice conditions with the measured cloud properties. This novel methodology is used to classify the observed clouds as coupled or decoupled to the WVT based on the location of the maximum vertical gradient of WVT height relative to the cloud-driven mixing layer. Only a conical sub-sector of sea ice concentration (SIC) and the lead fraction (LF) centred on the RV Polarstern location and extending up to 50 km in radius and with an azimuth angle governed by the time-dependent wind direction measured at the maximum WVT is related to the observed clouds. We found significant asymmetries for cases when the clouds are coupled or decoupled to the WVT and selected by LF regimes. Liquid water path of low-level clouds is found to increase as a function of LF, while the ice water path does so only for deep precipitating systems. Clouds coupled to WVT are found to generally have a lower cloud base and larger thickness than decoupled clouds. Thermodynamically, for coupled cases the cloud-top temperature is warmer and accompanied by a temperature inversion at the cloud top, whereas the decoupled cases are found to be closely compliant with the moist adiabatic temperature lapse rate. The ice water fraction within the cloud layer has been found to present a noticeable asymmetry when comparing coupled versus decoupled cases. This novel approach of coupling sea ice to cloud properties via the WVT mechanism unfolds a new tool to study Arctic surface–atmosphere processes. With this formulation, long-term observations can be analysed to enforce the statistical significance of the asymmetries. Furthermore, our results serve as an opportunity to better understand the dynamic linkage between clouds and sea ice and to evaluate its representation in numerical climate models for the Arctic system.

54 ENVIRONMENTAL SCIENCES↗

Derivation of physical and optical properties of mid-latitude cirrus ice crystals for a size-resolved cloud microphysics model

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

54 ENVIRONMENTAL SCIENCES↗

Investigation of ice cloud microphysical properties of DCSs using aircraft in situ measurements during MC3E over the ARM SGP site

In this study, six deep convective systems (DCSs) with a total of 5589 five-second samples and a range of temperatures from -41°C to 0°C during the Midlatitude Continental Convective Clouds Experiment (MC3E) were selected to investigate the ice cloud microphysical properties of DCSs over the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site. The ice cloud measurements of the DCS cases were made by the University of North Dakota Citation II research aircraft, and the ice cloud properties were derived through the following processes. First, the instances of supercooled liquid water in the ice-dominated cloud layers of DCSs have been eliminated using multisensor detection, including the Rosemount Icing Detector, King and Cloud Droplet Probes, as well as 2DC and Cloud Imaging Probe images. Then the Nevzorov-measured ice water contents (IWCs) at maximum diameter D max < 4000 µm are used as the best estimation to determine a new mass-dimensional relationship. Finally, the newly derived mass-dimensional relationship (a = 0.00365, b = 2.1) has been applied to a full spectrum of particle size distributions (PSDs, 120–30,000 µm) constructed from both 2DC and High-Volume Precipitation Spectrometer measurements to calculate the best-estimated IWCs of DCSs during MC3E. The averages of the total number concentrations (N t ), median mass diameter (D m ), maximum diameter (D max ), and IWC from six selected cases are 0.035 cm -3 , 1666 µm, 8841 µm, and 0.45 g m -3 , respectively. The gamma-type-size distributions are then generated matching the observed PSDs (120–30,000 µm), and the fitted gamma parameters are compared with the observed PSDs through multimoment assessments including first moment (D m ), third moment (IWC), and sixth moment (equivalent radar reflectivity, Z e ). Lastly, for application of observed PSDs to the remote sensing community, a series of empirical relationships between fitted parameters and Z e values has been derived, and the bullet rosette ice crystal backscattering relationship has been suggested for ground-based remote sensing.

58 GEOSCIENCES↗

Parameterizations of Cloud Microphysics and Indirect Aerosol Effects

Aerosols and especially their effect on clouds are one of the key components of the climate system and the hydrological cycle [Ramanathan et al., 2001]. Yet, the aerosol effect on clouds remains largely unknown and the processes involved not well understood. A recent report published by the National Academy of Science states "The greatest uncertainty about the aerosol climate forcing - indeed, the largest of all the uncertainties about global climate forcing - is probably the indirect effect of aerosols on clouds [NRC, 2001]." The aerosol effect on clouds is often categorized into the traditional "first indirect (i.e., Twomey)" effect on the cloud droplet sizes for a constant liquid water path [Twomey, 1977] and the "semi-direct" effect on cloud coverage [e.g., Ackerman et al., 2000]. Enhanced aerosol concentrations can also suppress warm rain processes by producing a narrow droplet spectrum that inhibits collision and coalescence processes [e.g., Squires and Twomey, 1961; Warner and Twomey, 1967; Warner, 1968; Rosenfeld, 1999]. The aerosol effect on precipitation processes, also known as the second type of aerosol indirect effect [Albrecht, 1989], is even more complex, especially for mixed-phase convective clouds. Table 1 summarizes the key observational studies identifying the microphysical properties, cloud characteristics, thermodynamics and dynamics associated with cloud systems from high-aerosol continental environments. For example, atmospheric aerosol concentrations can influence cloud droplet size distributions, warm-rain process, cold-rain process, cloud-top height, the depth of the mixed phase region, and occurrence of lightning. In addition, high aerosol concentrations in urban environments could affect precipitation variability by providing an enhanced source of cloud condensation nuclei (CCN). Hypotheses have been developed to explain the effect of urban regions on convection and precipitation [van den Heever and Cotton, 2007 and Shepherd, 2005]. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated.

54 ENVIRONMENTAL SCIENCES↗

Deriving Arctic Cloud Microphysics at Barrow, Alaska. Algorithms, Results, and Radiative Closure

Cloud phase and microphysical properties control the radiative effects of clouds in the climate system and are therefore crucial to characterize in a variety of conditions and locations. An Arctic-specific, ground-based, multi-sensor cloud retrieval system is described here and applied to two years of observations from Barrow, Alaska. Over these two years, clouds occurred 75% of the time, with cloud ice and liquid each occurring nearly 60% of the time. Liquid water occurred at least 25% of the time even in the winter, and existed up to heights of 8 km. The vertically integrated mass of liquid was typically larger than that of ice. While it is generally difficult to evaluate the overall uncertainty of a comprehensive cloud retrieval system of this type, radiative flux closure analyses were performed where flux calculations using the derived microphysical properties were compared to measurements at the surface and top-of-atmosphere. Radiative closure biases were generally smaller for cloudy scenes relative to clear skies, while the variability of flux closure results was only moderately larger than under clear skies. The best closure at the surface was obtained for liquid-containing clouds. Radiative closure results were compared to those based on a similar, yet simpler, cloud retrieval system. These comparisons demonstrated the importance of accurate cloud phase classification, and specifically the identification of liquid water, for determining radiative fluxes. Enhanced retrievals of liquid water path for thin clouds were also shown to improve radiative flux calculations.

97 MATHEMATICS AND COMPUTING↗