Comparison of simulated with co-located satellite millimeter-wave observations
Contributed by D. H. Staelin, MIT, Cambridge, MA, USA

Accurate retrievals of precipitation rates and water paths require consistency between reality and the many physical assumptions underlying the retrieval method.  One test of this consistency involves comparisons of predicted and observed brightness temperatures over the same ensemble of storms.  Reasonable agreement has been reported between: (a) AMSU observations of brightness temperature histograms 50-190 GHz over 122 globally and seasonally distributed storms, and (b) 15-km predictions made using NCEP-initialized 4-6 hour MM5 forecasts.  The particular cloud-resolving form of MM5 used and the radiative transfer model are described in the attached preprint entitled "Comparison of AMSU millimeter-wave satellite observations, MM5/TBSCAT predicted radiances, and electromagnetic models for hydrometeors" by C. Surussavadee and D. H. Staelin (IEEE Trans. Geosci. and Remote Sensing, in press, 2006).

The Goddard explicit cloud physics model was used, which provides profiles for snow, graupel, cloud ice, and rain water.  Although Mie scattering was assumed, the density of the icy spheres was assumed to have a wavelength-dependent "ice factor" value F that was determined separately for snow and graupel using electromagnetic scattering computations (DDSCAT) for hexagonal plates and 6-pointed rosettes, respectively.  Reasonable agreement between observed and predicted brightness temperature histograms was also separately found for each of five classes of precipitating AMSU pixels (convective, stratiform, snowy, non-snowy, or non-glaciated), although the brightness temperature histograms for window channels over snow were degraded by the unmodeled surface emissivity spectrum of snow (an omission that can be remedied).  The validity of the DDSCAT computations was tested by varying F for each wavelength so as to minimize the discrepancy between the AMSU observations and the simulated radiances; the resulting two sets of F values agreed within ~0.1 or less, where 0 < F < 1.  We have separately noted that our F values derived using DDSCAT are a weak function of the size distribution of the icy hydrometeors for frequencies below ~190 GHz, but that size matters more at higher frequencies.

We have also documented additional tests of the robustness of these MM5 and radiative transfer comparisons with AMSU that involve varying several of the physical assumptions such as the altitudes and abundances of the snow and graupel predicted by MM5, the size distributions of the various ice habits, and the hydrometeor loss tangents, F values, back-scattering ratios, etc.

This work (Surussavadee, C., and D. H. Staelin, 2007: Millimeter-Wave Precipitation Retrievals and Observed-versus-Simulated Radiance Distributions: Sensitivity to Assumptions. J. Atmos. Sci., 64, 3808-3826) is discussed further here under the theme Retrieval Techniques.