Abstract-Falling snow is an important component of global precipitation in extratropical regions. This paper describes the methodology and results of physically based retrievals of snow falling over land surfaces. Because microwave brightness temperatures emitted by snow-covered surfaces are highly variable, precipitating snow above such surfaces is difficult to observe using window channels that occur at low frequencies ( 100 GHz). Furthermore, at frequencies 37 GHz, sensitivity to liquid hydrometeors is dominant. These problems are mitigated at high frequencies ( 100 GHz) where water vapor screens the surface emission, and sensitivity to frozen hydrometeors is significant. However, the scattering effect of snowfall in the atmosphere at those higher frequencies is also impacted by water vapor in the upper atmosphere. The theory of scattering by randomly oriented dry snow particles at high microwave frequencies appears to be better described by regarding snow as a concatenation of "equivalent" ice spheres rather than as a sphere with the effective dielectric constant of an air-ice mixture. An equivalent sphere snow scattering model was validated against high-frequency attenuation measurements. Satellite-based high-frequency observations from an Advanced Microwave Sounding Unit (AMSU-B) instrument during the March 5-6, 2001 New England blizzard were used to retrieve snowfall over land. Vertical distributions of snow, temperature, and relative humidity profiles were derived from the Mesoscale Model (MM5) cloud model. Those data were applied and modified in a radiative transfer model that derived brightness temperatures consistent with the AMSU-B observations. The retrieved snowfall distribution was validated with radar reflectivity measurements obtained from a ground-based radar network.

This study seeks to evaluate the impact of several newly available sources of meteorological data on mesoscale model forecasts of the extratropical cyclone that struck Florida on 2 February 1998. Intermittent measurements of precipitation and integrated water vapor (IWV) distributions were obtained from Special Sensor Microwave/ Imager (SSM/I) and Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) observations. The TMI also provided sea surface temperatures (SSTs) with structural detail of the Loop Current and Gulf Stream. Continuous lightning distributions were measured with a network of very low frequency radio receivers. Lightning data were tuned with intermittent spaceborne microwave radiometer data through a probability matching technique to continuously estimate convective rainfall rates. A series of experiments were undertaken to evaluate the effect of those data on mesoscale model forecasts produced after assimilating processed rainfall and IWV for 6 h. Assimilating processed rainfall, IWV, and SSTs from TMI measurements in the model yielded improved forecasts of precipitation distributions and vertical motion fields. Assimilating those data also produced an improved 9-h forecast of the radar reflectivity cross section that was validated with a coincident observation from the TRMM spaceborne precipitation radar. Sensitivity experiments showed that processed rainfall information had greater impact on the rainfall forecast than IWV and SST information. Assimilating latent heating in the correct location of the forecast model was found to be more important than an accurate determination of the rainfall intensity.

In this study, we present an improved physical model to retrieve snowfall rate over land using brightness temperature observations from the National Oceanic and Atmospheric Administration's (NOAA) Advanced Microwave Sounder Unit-B (AMSU-B) at 89 GHz, 150 GHz, 183.3±1 GHz, 183.3±3 GHz, and 183.3±7 GHz. The retrieval model is applied to the New England blizzard of March 5, 2001 which deposited about 75 cm of snow over much of Vermont, New Hampshire, and northern New York.In this improved physical model, prior retrieval assumptions about snowflake shape, particle size distributions, environmental conditions, and optimization methodology have been updated. Here, single scattering parameters for snow particles are calculated with the Discrete-Dipole Approximation (DDA) method instead of assuming spherical shapes.Five different snow particle models (hexagonal columns, hexagonal plates, and three different kinds of aggregates) are considered. Snow particle size distributions are assumed to vary with air temperature and to follow aircraft measurements described by previous studies.Brightness temperatures at AMSU-B frequencies for the New England blizzard are calculated using these DDA calculated single scattering parameters and particle size distributions. The vertical profiles of pressure, temperature, relative humidity and hydrometeors are provided by MM5 model simulations. These profiles are treated as the a priori data base in the Bayesian retrieval algorithm. In algorithm applications to the blizzard data, calculated brightness temperatures associated with selected database profiles agree with AMSU-B observations to within about ±5 K at all five frequencies.Retrieved snowfall rates compare favorably with the near-concurrent National Weather Service (NWS) radar reflectivity measurements. The relationships between the NWS radar measured reflectivities Z e and retrieved snowfall rate R for a given snow particle model are derived by a histogram matching technique. All of these Z e -R relationships fall in the range of previously established Z e -R relationships for snowfall. This suggests that the current physical model developed in this study can reliably estimate the snowfall rate over land using the AMSU-B measured brightness temperatures.1

Abstract. In this study, optimal parameter estimations are performed for both physical and computational parameters in a mesoscale meteorological model, and their impacts on the quantitative precipitation forecasting (QPF) are assessed for a heavy rainfall case occurred at the Korean Peninsula in June 2005. Experiments are carried out using the PSU/NCAR MM5 model and the genetic algorithm (GA) for two parameters: the reduction rate of the convective available potential energy in the Kain-Fritsch (KF) scheme for cumulus parameterization, and the Asselin filter parameter for numerical stability. The fitness function is defined based on a QPF skill score. It turns out that each optimized parameter significantly improves the QPF skill. Such improvement is maximized when the two optimized parameters are used simultaneously. Our results indicate that optimizations of computational parameters as well as physical parameters and their adequate applications are essential in improving model performance.