Final Report on the Development of an Improved Cloud Microphysical Product for Model and Remote Sensing Evaluation using RACORO Observations

We proposed to analyze data collected during the Routine Aerial Facilities (AAF) Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations (RACORO) in order to develop an integrated product of cloud microphysical properties (number concentration of drops in different size bins, total liquid drop concentration integrated over all bin sizes, liquid water content LWC, extinction of liquid clouds bw, effective radius of water drops re, and radar reflectivity factor) that could be used to evaluate large-eddy simulations (LES), general circulation models (GCMs) and ground-based remote sensing retrievals, and to develop cloud parameterizations with the end goal of improving the modeling of cloud processes and properties and their impact on atmospheric radiation. We have completed the development of this microphysical database and have submitted it to ARM for consideration of its inclusion on the ARM database as a PI product. This report describes the development of this database, and also describes research that has been conducted on cloud-aerosol interactions using the data obtained during RACORO. A list of conference proceedings and publications is also included

Multi-Layer Arctic Mixed-Phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity Experiments

A cloud-resolving model (CRM) is used to simulate the multiple-layer mixed-phase stratiform (MPS) clouds that occurred during a three-and-a-half day subperiod of the Department of Energy-Atmospheric Radiation Measurement Program s Mixed-Phase Arctic Cloud Experiment (M-PACE). The CRM is implemented with an advanced two-moment microphysics scheme, a state-of-the-art radiative transfer scheme, and a complicated third-order turbulence closure. Concurrent meteorological, aerosol, and ice nucleus measurements are used to initialize the CRM. The CRM is prescribed by time-varying large-scale advective tendencies of temperature and moisture and surface turbulent fluxes of sensible and latent heat. The CRM reproduces the occurrences of the single- and double-layer MPS clouds as revealed by the M-PACE observations. However, the simulated first cloud layer is lower and the second cloud layer thicker compared to observations. The magnitude of the simulated liquid water path agrees with that observed, but its temporal variation is more pronounced than that observed. As in an earlier study of single-layer cloud, the CRM also captures the major characteristics in the vertical distributions and temporal variations of liquid water content (LWC), total ice water content (IWC), droplet number concentration and ice crystal number concentration (nis) as suggested by the aircraft observations. However, the simulated mean values differ significantly from the observed. The magnitude of nis is especially underestimated by one order of magnitude. Sensitivity experiments suggest that the lower cloud layer is closely related to the surface fluxes of sensible and latent heat; the upper cloud layer is probably initialized by the large-scale advective cooling/moistening and maintained through the strong longwave (LW) radiative cooling near the cloud top which enhances the dynamical circulation; artificially turning off all ice-phase microphysical processes results in an increase in LWP by a factor of 3 due to interactions between the excessive LW radiative cooling and extra cloud water; heating caused by phase change of hydrometeors could affect the LWC and cloud top height by partially canceling out the LW radiative cooling. It is further shown that the resolved dynamical circulation appears to contribute more greatly to the evolution of the MPS cloud layers than the parameterized subgrid-scale circulation

Ice and Supercooled Liquid Water Distributions Over the Southern Ocean Based on In Situ Observations and Climate Model Simulations

Three climate models are evaluated using in situ airborne observations from the Southern Ocean Clouds, Radiation, Aerosol Transport Experimental Study (SOCRATES) campaign. The evaluation targets cloud phases, microphysical properties, thermodynamic conditions, and aerosol indirect effects from −40°C to 0°C. Compared with 580-s averaged observations (i.e., 100 km horizontal scale), the Community Atmosphere Model version 6 (CAM6) shows the most similar result for cloud phase frequency distribution and allows more liquid-containing clouds below −10°C compared with its predecessor—CAM5. The Energy Exascale Earth System Model (E3SM) underestimates (overestimates) ice phase frequencies below (above) −20°C. CAM6 and E3SM show liquid and ice water contents (i.e., LWC and IWC) similar to observations from −25°C to 0°C, but higher LWC and lower IWC than observations at lower temperatures. Simulated in-cloud RH shows higher minimum values than observations, possibly restricting ice growth during sedimentation. As number concentrations of aerosols larger than 500 nm (Na500) increase, observations show increases of LWC, IWC, liquid, and ice number concentrations (Nliq, Nice). Number concentrations of aerosols larger than 100 nm (Na100) only show positive correlations with LWC and Nliq. From −20°C to 0°C, higher aerosol number concentrations are correlated with lower glaciation ratio and higher cloud fraction. From −40°C to −20°C, large aerosols show positive correlations with glaciation ratio. CAM6 shows small increases of LWC and Nliq with Na500 and Na100. E3SM shows small increases of Nice with Na500. Overall, CAM6 and E3SM underestimate aerosol indirect effects on ice crystals and supercooled liquid droplets over the Southern Ocean

A New Parameterization of Single Scattering Solar Radiative Properties for Tropical Anvils Using Observed Ice Crystal Size and Shape Distributions

Parameterizations of single scattering properties currently used in cloud resolving and general circulation models are somewhat limited in that they typically assume the presence of single particle habits, do not adequately account for the numbers of ice crystals with diameters smaller than 100 mm, and contain no information about the variance of parameterization coefficients. Here, new parameterizations of mean single scattering properties (e.g., single scatter albedo, asymmetry parameter, and extinction efficiency) for distributions of ice crystals in tropical anvils are developed. Using information about the size and shape of ice crystals acquired by a two-dimensional cloud probe during the Central Equatorial Pacific Experiment (CEPEX), a self-organized neural network defines shape based on simulations of how the particle maximum dimension and area ratio (ratio of projected area to that of circumscribed circle with maximum dimension) vary for random orientations of different idealized shapes (i.e., columns, bullet rosettes, rough aggregates, and particles represented by Chebyshev poly-nomials). The size distributions for ice crystals smaller than 100 mm are based on parameterizations developed using representative samples of 11 633 crystals imaged by a video ice particle sampler (VIPS). The mean-scattering properties for distributions of ice crystals are then determined by weighting the single scattering properties of individual ice crystals, determined using an improved geometric ray-tracing method, according t

ARM/GCSS/SPARC TWP-ICE CRM Intercomparison Study

Specifications are provided for running a cloud-resolving model (CRM) and submitting results in a standardized format for inclusion in a n intercomparison study and archiving for public access. The simulated case study is based on measurements obtained during the 2006 Tropical Warm Pool - International Cloud Experiment (TWP-ICE) led by the U. S. department of Energy Atmospheric Radiation Measurement (ARM) program. The modeling intercomparison study is based on objectives developed in concert with the Stratospheric Processes And their Role in Climate (SPARC) program and the GEWEX cloud system study (GCSS) program. The Global Energy and Water Cycle Experiment (GEWEX) is a core project of the World Climate Research PRogramme (WCRP)

Dynamics of Cloud-Top Generating Cells in Winter Cyclones. Part III: Shear and Convective Organization

Cloud-top generating cells (GCs) are a common feature atop stratiform clouds within the comma head of winter cyclones. The dynamics of cloud-top GCs are investigated using very high-resolution idealized WRF Model simulations to examine the role of shear in modulating the structure and intensity of GCs. Simulations were run for the same combinations of radiative forcing and instability as in Part II of this series, but with six different shear profiles ranging from 0 to 10ms21 km21 within the layer encompassing the GCs. The primary role of shear was to modulate the organization of GCs, which organized as closed convective cells in simulations with radiative forcing and no shear. In simulations with shear and radiative forcing, GCs organized in linear streets parallel to the wind. No GCs developed in the initially stable simulations with no radiative forcing. In the initially unstable and neutral simulations with no radiative forcing or shear, GCs were exceptionally weak, with no clear organization. In moderate-shear (Du/Dz 5 2, 4ms21 km21) simulations with no radiative forcing, linear organization of the weak cells was apparent, but this organization was less coherent in simulations with high shear (Du/Dz 5 6, 8, 10ms21 km21). The intensity of the updrafts was primarily related to the mode of radiative forcing but was modulated by shear. The more intense GCs in nighttime simulations were either associated with no shear (closed convective cells) or strong shear (linear streets). Updrafts within GCs under conditions with radiative forcing were typically ;1–2 ms21 with maximum values , 4ms21