Abstract:There are a wide variety of bulk microphysical parameterizations (BMPs) available in the Weather Research and Forecasting (WRF) model, ranging from the simple Kessler BMP to the latest Thompson scheme (Thompson et al. 2008). One of the critical parameters in a BMP is the treatment of snow and graupel. Most BMPs treat snow and graupel by assuming spherical particles with constant density. However, in reality, ice particles have a variety of habits and complex structures, while riming changes and modifies the properties of ice particles. Both ice habit and riming should be considered in BMPs to better represent the ice processes and precipitation fallout. In general, there have been two approaches to include habits in BMPs. One is to increase the numbers of prognosed variables to simulate the different habits explicitly (Straka and Mansell 2005), while the other is to estimate the preferred habit at each grid point based on the local conditions (Meyers et al. 1997), which can be tracked by partitioning the snow category into different habit types (Woods et al. 2007). Snow and graupel (or hail) are kept often as separate categories in traditional BMPs. Due to the large fall velocity difference between snow and graupel (or hail), the residence time and fallout of precipitation depends highly on the partition of snow and graupel, which has a large sensitivity on the surface precipitation distribution and precipitation efficiency (Lin and Colle 2008). By separately predicting the ice mixing ratios acquired through deposition versus riming, the riming effect on the mass and fall speed of snow can be considered (Stoelinga et al. 2007). Using a similar approach, Morrison and Grabowski (2008) considered the variation of particle properties with particle size and rimed mass fraction. Due to the size dependent properties, the method can be easily applied in a bin scheme. However, the approach has its shortcomings: only one habit type is used for nonspherical particles; it is difficult to derive a continuous mass-diameter (m-D) and area-diamater (A-D) relationships over the whole size spectrum; graupel size is smaller than partially-rimed particles; and analytical integration of total mass is impossible due to the size dependent m-D and velocity-diameter (V-D) relationships. In this paper we propose a new bulk microphysical scheme for WRF (SUNY scheme), which relates the crystal habits and degree of riming directly to the ice properties (m-D, A-D, and V-D). By diagnosing a varying __________________________________________ *Corresponding author address: Yanluan Lin, SoMAS, Stony Brook Univ. NY, 11794 email: yanlin@ic.sunysb.edu degree of riming, the continuous spectrum from pristine ice to slightly rimed ice to heavily rimed snow and graupel is represented.

Steady-state suspension of snow

Unified Model of Sediment Transport Threshold and Rate Across Weak and Intense Subaqueous Bedload, Windblown Sand, and Windblown Snow

On the relative importance of saltation and suspension for aeolian snow mass transport

Saltation of snow

Unified Model of Threshold and Rate of Equilibrium Nonsuspended Sediment Transport Across Most Aeolian and Fluvial Environments

Numerical method used for the simulation of fluid–solid transitions based on the VOF and SAM models and prediction of snow distributions on large-span roofs

A variational method for estimating surface fluxes with mass conservation constraint

Vertical Profiles of Mass Flux for Windblown Sand Movement at Steady State

The effect of turbulence on drifting snow sublimation

Settling and Clustering of Snow Particles in Atmospheric Turbulence

A Drifting and Blowing Snow Scheme in the Weather Research and Forecasting Model

Large eddy simulation of near-surface boundary layer dynamics over patchy snow

The Seasonal Snow Cover Dynamics: Review on Wind-Driven Coupling Processes

Field investigation of 3-D snow settling dynamics under weak atmospheric turbulence

Sublimation of Snow

Development of a new bulk microphysical scheme for WRF with varying snow characteristics and riming intensity

SnowPappus v1.0, a blowing-snow model for large-scale applications of the Crocus snow scheme

Field investigation of 3D snow settling dynamics under weak atmospheric turbulence

Simulation and Statistical Analysis of Ground Snow Loads Based on a Multi-Layer Snow Accumulation and Melt Model

Quantifying the Aerodynamic Roughness Length of Snow Surfaces with Time‐Lapse Structure‐From‐Motion

Numerical study on slipstream-induced snow drifting and accumulation in the bogie region of a high-speed train passing the snowy ballast bed