An Evaluation of Size-Resolved Cloud Microphysics Scheme Numerics for Use with Radar Observations. Part I: Collision-Coalescence

This study evaluates some available schemes designed to solve the stochastic collection equation (SCE) for collision-coalescence of hydrometeors using a size-resolved (bin) microphysics approach, and documents their numerical properties within the framework of a box model. Comparing three widely used SCE schemes, we find that all converge to almost identical solutions at sufficiently fine mass grids. However, one scheme converges far slower than the other two and shows pronounced numerical diffusion at the large-drop tail of the size distribution. One of the remaining two schemes is recommended on the basis that it is well-converged on a relatively coarse mass grid, stable for large time steps, strictly mass-conservative, and computationally efficient. To examine the effects of SCE scheme choice on simulating clouds and precipitation, two of the three schemes are compared in large-eddy simulations of a drizzling stratocumulus field. A forward simulator that produces Doppler spectra from the large-eddy simulation results is used to compare the model output directly with radar observations. The scheme with pronounced numerical diffusion predicts excessively large mean Doppler velocities and overly broad and negatively skewed spectra compared with observations, consistent with numerical diffusion demonstrated in the box model. Statistics obtained using the recommended scheme are closer to observations, but notable differences remain, indicating that factors other than SCE scheme accuracy are limiting simulation fidelity

Establishing cross-discipline consensus on contraception, pregnancy and breast feeding-related educational messages and clinical practices to support women with rheumatoid arthritis: An Australian Delphi study

Objective: Recognising the need for a best-practice and consistent approach in providing care to women with rheumatoid arthritis (RA) in relation to (1) general health, (2) contraception, (3) conception and pregnancy, (4) breast feeding and (5) early parenting, we sought to achieve cross-discipline, clinical consensus on key messages and clinical practice behaviours in these 5 areas. Design: 3-round eDelphi study. In round 1, panellists provided free-text responses to open-ended questions about care for women with RA across the 5 areas. Subsequently, panellists refined and scored the synthesised responses, presented as metathemes, themes and detailed elements. Where =5% of panellists did not support a theme in a given round, it was removed. Setting: Panel of practicing Australian rheumatologists (n=22), obstetricians/obstetric medicine physicians (n=9) and pharmacists (n=5). Results: 34 (94.4%) panellists participated in all 3 rounds. The panel supported 18 themes across the 5 areas (support/strongly support: 88.2-100%) underpinned by 5 metathemes. Metathemes focused on coordination in information delivery, the mode and timing of information delivery, evidence underpinning information, engagement of the right health professionals at the right time and a non-judgemental approach to infant feeding. Themes included practices for primary prevention of chronic disease and their sequelae, the importance of contraception and planning pregnancy and breast feeding, close monitoring of medications, supporting mental wellbeing, managing disease activity and providing practical support for early parenting. Conclusions: A cross-disciplinary clinical panel highly supported key information and clinical practices in the care for women with RA across the continuum of contraception to early parenting within a whole-person, chronic disease management approach

A Model for Particle Microphysics, Turbulent Mixing, and Radiative Transfer in the Stratocumulus-Topped Marine Boundary Layer and Comparisons with Measurements

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 gale 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-micrometers 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 underprediction of peak supersaturations

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

Spectral Dependence of MODIS Cloud Droplet Effective Radius Retrievals for Marine Boundary Layer Clouds

Low-level warm marine boundary layer (MBL) clouds cover large regions of Earth's surface. They have a significant role in Earth's radiative energy balance and hydrological cycle. Despite the fundamental role of low-level warm water clouds in climate, our understanding of these clouds is still limited. In particular, connections between their properties (e.g. cloud fraction, cloud water path, and cloud droplet size) and environmental factors such as aerosol loading and meteorological conditions continue to be uncertain or unknown. Modeling these clouds in climate models remains a challenging problem. As a result, the influence of aerosols on these clouds in the past and future, and the potential impacts of these clouds on global warming remain open questions leading to substantial uncertainty in climate projections. To improve our understanding of these clouds, we need continuous observations of cloud properties on both a global scale and over a long enough timescale for climate studies. At present, satellite-based remote sensing is the only means of providing such observations

Accuracy Assessments of Cloud Droplet Size Retrievals from Polarized Reflectance Measurements by the Research Scanning Polarimeter

We present an algorithm for the retrieval of cloud droplet size distribution parameters (effective radius and variance) from the Research Scanning Polarimeter (RSP) measurements. The RSP is an airborne prototype for the Aerosol Polarimetery Sensor (APS), which was on-board of the NASA Glory satellite. This instrument measures both polarized and total reflectance in 9 spectral channels with central wavelengths ranging from 410 to 2260 nm. The cloud droplet size retrievals use the polarized reflectance in the scattering angle range between 135deg and 165deg, where they exhibit the sharply defined structure known as the rain- or cloud-bow. The shape of the rainbow is determined mainly by the single scattering properties of cloud particles. This significantly simplifies both forward modeling and inversions, while also substantially reducing uncertainties caused by the aerosol loading and possible presence of undetected clouds nearby. In this study we present the accuracy evaluation of our algorithm based on the results of sensitivity tests performed using realistic simulated cloud radiation fields

Effects of Cloud Horizontal Inhomogeneity and Drizzle on Remote Sensing of Cloud Droplet Effective Radius: Case Studies Based on Large-eddy Simulations

This study investigates effects of drizzle and cloud horizontal inhomogeneity on cloud effective radius (re) retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS). In order to identify the relative importance of various factors, we developed a MODIS cloud property retrieval simulator based on the combination of large-eddy simulations (LES) and radiative transfer computations. The case studies based on synthetic LES cloud fields indicate that at high spatial resolution (100 m) 3-D radiative transfer effects, such as illumination and shadowing, can induce significant differences between retrievals ofre based on reflectance at 2.1 m (re,2.1) and 3.7 m (re,3.7). It is also found that 3-D effects tend to have stronger impact onre,2.1 than re,3.7, leading to positive difference between the two (re,3.72.1) from illumination and negative re,3.72.1from shadowing. The cancellation of opposing 3-D effects leads to overall reasonable agreement betweenre,2.1 and re,3.7 at high spatial resolution as far as domain averages are concerned. At resolutions similar to MODIS, however, re,2.1 is systematically larger than re,3.7when averaged over the LES domain, with the difference exhibiting a threshold-like dependence on bothre,2.1and an index of the sub-pixel variability in reflectance (H), consistent with MODIS observations. In the LES cases studied, drizzle does not strongly impact reretrievals at either wavelength. It is also found that opposing 3-D radiative transfer effects partly cancel each other when cloud reflectance is aggregated from high spatial resolution to MODIS resolution, resulting in a weaker net impact of 3-D radiative effects onre retrievals. The large difference at MODIS resolution between re,3.7 and re,2.1 for highly inhomogeneous pixels with H 0.4 can be largely attributed to what we refer to as the plane-parallelrebias, which is attributable to the impact of sub-pixel level horizontal variability of cloud optical thickness onre retrievals and is greater for re,2.1 than re,3.7. These results suggest that there are substantial uncertainties attributable to 3-D radiative effects and plane-parallelre bias in the MODIS re,2.1retrievals for pixels with strong sub-pixel scale variability, and theH index can be used to identify these uncertainties

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)Surface Heat Budget of the Arctic Ocean (SHEBA) campaign provide a unique opportunity 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 N(sub IN) measured above cloud top, but details of IN activation behavior are unknown. If activated rapidly (in deposition, condensation, or immersion modes), as commonly assumed, IN are depleted from the well-mixed boundary layer within minutes. Quasi-equilibrium ice number concentration N(sub i) is then limited to a small fraction of overlying N(sub IN) that is determined by the cloud-top entrainment rate w(sub e) divided by the number-weighted ice fall speed at the surface v(sub f). Because w(sub c) 10 cm/s, N(sub i)/N(sub IN)<< 1. Such conditions may be common for this cloud type, which has implications for modeling IN diagnostically, interpreting measurements, and quantifying sensitivity to increasing N(sub IN) (when w(sub e)/v(sub f)< 1, entrainment rate limitations serve to buffer cloud system response). To reproduce observed ice crystal size distributions and cloud radar reflectivities with rapidly consumed IN in this case, the measured above-cloud N(sub IN) must be multiplied by approximately 30. However, results are sensitive to assumed ice crystal properties not constrained by measurements. In addition, simulations do not reproduce the pronounced mesoscale heterogeneity in radar reflectivity that is observed