A nonlocal ob ser va tion op er a tor has been de vel oped to as sim i late GPS ra dio occultation (RO) refractivity with WRF 3DVAR. For sim plic ity, in the past GPS RO refractivity was of ten as sim i lated us ing a lo cal ob ser va tion op er a tor with the as sump tion that the GPS RO ob ser va tion was rep re sen ta tive of a model lo cal point. Such an op er a tor did not take into ac count the ef fects of hor i zon tal inhomogeneity on the de rived GPS RO refractivity. In or der to more ac cu rately model the observables, Sokolovskiy et al. (2005a) de vel oped a nonlocal ob ser va tion op er a tor, which would take into ac count the ef fects of hor i zon tal inhomogeneity on GPS RO mea sure ments. This nonlocal ob ser va tion op er a tor cal cu lates the in te grated amount of the model refractivity along the ray paths cen tered at the peri gee points. For com par a tive pur poses, the nonlocal ob ser va tion op er a tor can be sim pli fied by lim it ing the length of in te gra tion near the RO point. This is called the "lo cal op er a tor vari ant", which is equiv a lent to the orig i nal lo cal op er a tor ex cept that the orig i nal one is per formed with fixed tan gent points at ob ser va tion lev els. For com pu ta tional ef fi ciency, as sim i la tion us ing both the nonlocal op er a tor and lo cal op er a tor vari ant now is per formed with smear tan gent points at the mean height of each model ver ti cal level. In this study, the sta tis tics of ob ser va tion er rors us ing both lo cal and nonlocal op er a tors were es ti mated based on WRF sim u la tions. The ob ser va tion er rors produced by the nonlocal operator are about two times smaller than those generated by the local operator and in agreement with Sokolovskiy et al. (2005b).Each of the three op er a tors is used to as sim i late GPS RO refractivity soundings from the FORMOSAT-3/COS MIC mis sion us ing the WRF 3DVAR sys tem. The WRF model then is ap plied to sim u late Ty phoons Kaemi (July 2006) which struck Tai wan with sig nif i cant tor ren tial rain fall. The anal y sis in cre ments pro duced by the nonlocal op er a tor and lo cal op er a tor vari ant are quite sim i lar in hor i zon tal and ver ti cal dis tri bu tions; whereas, the for mer is slightly stretched along the ray's di rec tion, as a re sult of the lon ger in te gra tion length. The sim u lated ty phoon tracks prior to land fall are quite sim i lar for the three op er a tors. Both the nonlocal op er a tor and lo cal op er a tor vari ant im prove the de toured track af ter land fall as pre dicted by the lo cal op er a tor. The nonlocal op er a tor out per forms the two lo cal op er a tors in rain fall pre dic tion at later times. The per for mances of the nonlocal op er a tor in gen eral are promising and can replace the local operator at no marked cost of computational efficiency.Key words: FORMOSAT-3/COS MIC, Data as sim i la tion, Typhoon Ci ta tion: Chen, S. Y., C. Y. Huang, Y. H. Kuo, Y. R. Guo, and S. Sokolovskiy, 2009: As sim i la tion of GPS refractivity from FORM...