Assimilation of TRMM Multisatellite Precipitation Analysis with a Low-Resolution NCEP Global Forecast System

Current methods of assimilation of precipitation into numerical weather prediction models are able to make the model precipitation become similar to the observed precipitation during the assimilation, but the model forecasts tend to return to their original solution after a few hours. To facilitate the precipitation assimilation, a logarithm transformation has been used in several past studies. Lien et al. proposed instead to assimilate precipitation using the local ensemble transform Kalman filter (LETKF) with a Gaussian transformation technique and succeeded in improving the model forecasts in perfect-model observing system simulation experiments (OSSEs).In this study, the method of Lien et al. is tested within a more realistic configuration: the TRMM Multisatellite Precipitation Analysis (TMPA) data are assimilated into a low-resolution version of the NCEP Global Forecast System (GFS). With guidance from a statistical study comparing the GFS model background precipitation and the TMPA data, some modifications of the assimilation methods proposed in Lien et al. are made, including 1) applying separate Gaussian transformations to model and to observational precipitation based on their own cumulative distribution functions; 2) adopting a quality control criterion based on the correlation between the long-term model and observed precipitation data at the observation location; and 3) proposing a new method to define the transformation of zero precipitation that takes into account the zero precipitation probability in the background ensemble rather than the climatology. With these modifications, the assimilation of the TMPA precipitation data improves both the analysis and 5-day model forecasts when compared with a control experiment assimilating only rawinsonde data.