A new bulk microphysical parameterization (BMP) has been developed for use with the Weather Research and Forecasting (WRF) Model or other mesoscale models. As compared with earlier single-moment BMPs, the new scheme incorporates a large number of improvements to both physical processes and computer coding, and it employs many techniques found in far more sophisticated spectral/bin schemes using lookup tables. Unlike any other BMP, the assumed snow size distribution depends on both ice water content and temperature and is represented as a sum of exponential and gamma distributions. Furthermore, snow assumes a nonspherical shape with a bulk density that varies inversely with diameter as found in observations and in contrast to nearly all other BMPs that assume spherical snow with constant density. The new scheme’s snow category was readily modified to match previous research in sensitivity experiments designed to test the sphericity and distribution shape characteristics. From analysis of four idealized sensitivity experiments, it was determined that the sphericity and constant density assumptions play a major role in producing supercooled liquid water whereas the assumed distribution shape plays a lesser, but nonnegligible, role. Further testing using numerous case studies and comparing model results with in situ and other observations confirmed the results of the idealized experiments and are briefly mentioned herein, but more detailed, microphysical comparisons with observations are found in a companion paper in this series (Part III, forthcoming).

A new two-moment cloud microphysics scheme predicting the mixing ratios and number concentrations of five species (i.e., cloud droplets, cloud ice, snow, rain, and graupel) has been implemented into the Weather Research and Forecasting model (WRF). This scheme is used to investigate the formation and evolution of trailing stratiform precipitation in an idealized two-dimensional squall line. Results are compared to those using a one-moment version of the scheme that predicts only the mixing ratios of the species, and diagnoses the number concentrations from the specified size distribution intercept parameter and predicted mixing ratio. The overall structure of the storm is similar using either the one-or two-moment schemes, although there are notable differences. The two-moment (2-M) scheme produces a widespread region of trailing stratiform precipitation within several hours of the storm formation. In contrast, there is negligible trailing stratiform precipitation using the one-moment (1-M) scheme. The primary reason for this difference are reduced rain evaporation rates in 2-M compared to 1-M in the trailing stratiform region, leading directly to greater rain mixing ratios and surface rainfall rates. Second, increased rain evaporation rates in 2-M compared to 1-M in the convective region at midlevels result in weaker convective updraft cells and increased midlevel detrainment and flux of positively buoyant air from the convective into the stratiform region. This flux is in turn associated with a stronger mesoscale updraft in the stratiform region and enhanced ice growth rates. The reduced (increased) rates of rain evaporation in the stratiform (convective) regions in 2-M are associated with differences in the predicted rain size distribution intercept parameter (which was specified as a constant in 1-M) between the two regions. This variability is consistent with surface disdrometer measurements in previous studies that show a rapid decrease of the rain intercept parameter during the transition from convective to stratiform rainfall.

This study evaluates the sensitivity of winter precipitation to numerous aspects of a bulk, mixed-phase microphysical parameterization found in three widely used mesoscale models [the fifth-generation Pennsylvania State University-National Center for Atmospheric Research Mesoscale Model (MM5), the Rapid Update Cycle (RUC), and the Weather Research and Forecast (WRF) model]. Sensitivities of the microphysics to primary ice initiation, autoconversion, cloud condensation nuclei (CCN) spectra, treatment of graupel, and parameters controlling the snow and rain size distributions are tested. The sensitivity tests are performed by simulating various cloud depths (with different cloud-top temperatures) using flow over an idealized two-dimensional mountain. The height and width of the two-dimensional barrier are designed to reproduce an updraft pattern with extent and magnitude consistent with documented freezing-drizzle cases. By increasing the moisture profile to saturation at low temperatures, a deep, precipitating snow cloud is also simulated. Upon testing the primary sensitivities of the microphysics scheme in two dimensions as reported in the present study, the MM5 with the modified scheme will be tested in multiple case studies and the results will be compared to observations in a forthcoming companion paper, Part II. The key results of this study are 1) the choice of ice initiation schemes is relatively unimportant for deep precipitating snow clouds but more important for shallow warm clouds having cloud-top temperature greater than Ϫ13ЊC, 2) the assumed snow size distribution and associated snow diffusional growth along with the assumed graupel size distribution and method of transforming rimed snow into graupel have major impacts on the mass of cloud water and formation of freezing drizzle, and 3) a proper simulation of drizzle using a singlemoment scheme and exponential size distribution requires an increase in the rain intercept parameter, thereby reducing rain terminal velocities to values more characteristic of drizzle.

No abstract

Aerosols influence cloud and precipitation development in complex ways due to myriad feedbacks at a variety of scales from individual clouds through entire storm systems. This paper describes the implementation, testing, and results of a newly modified bulk microphysical parameterization with explicit cloud droplet nucleation and ice activation by aerosols. Idealized tests and a high-resolution, convection-permitting, continental-scale, 72-h simulation with five sensitivity experiments showed that increased aerosol number concentration results in more numerous cloud droplets of overall smaller size and delays precipitation development. Furthermore, the smaller droplet sizes cause the expected increased cloud albedo effect and more subtle longwave radiation effects. Although increased aerosols generally hindered the warm-rain processes, regions of mixed-phase clouds were impacted in slightly unexpected ways with more precipitation falling north of a synoptic-scale warm front. Aerosol impacts to regions of light precipitation, less than approximately 2.5 mm h−1, were far greater than impacts to regions with higher precipitation rates. Comparisons of model forecasts with five different aerosol states versus surface precipitation measurements revealed that even a large-scale storm system with nearly a thousand observing locations did not indicate which experiment produced a more correct final forecast, indicating a need for far longer-duration simulations due to the magnitude of both model forecast error and observational uncertainty. Last, since aerosols affect cloud and precipitation phase and amount, there are resulting implications to a variety of end-user applications such as surface sensible weather and aircraft icing.