本研究建置一套都卜勒雷達資料同化系統，結合局地系集轉換卡爾曼濾波器與天氣研究預報模式，並針對在台灣中、南部造成降雨紀錄和嚴重災害的莫拉克颱風（2009）個案，進行觀測系統模擬實驗與真實觀測實驗，以評估此雷達資料同化系統對於定量降水即時預報的助益。而研究的最大意義，為在台灣地形複雜且四周海域缺乏觀測等不利預報的條件下，提供以雷達資料同化改善颱風降雨即時預報的有效方案。

在觀測系統模擬實驗中，三維風場和雨水混合比為同化徑向風和回波之後改善最多的預報變數，可歸因於觀測算符內的直接關係。在降雨預報方面，雷達資料的正面影響可長達6小時。同化回波的改善主要在預報初期，同化徑向風的改善則較為延遲，而兩者皆同化的降雨預報表現最好。增加上游對流區域的觀測覆蓋量，亦可大幅提升降雨預報表現。另外，本研究針對颱風環流下所發展的對流雨帶，提出使用混合局地法進行雷達資料同化，以進一步改善降雨預報。

在真實觀測實驗中，此雷達資料同化系統仍能有效改善定量降水即時預報。同化回波時須使用變數局地化法，只用來更新雨水混合比。使用觀測空間的統計方法，能診斷預報偏差和理想系集離散度。混合局地化法在真實觀測實驗的效益更加明顯，尤其能提升觀測資料稀疏或破碎處的風場準確度，進而改善降雨預報。

This study develops a Doppler radar data assimilation system, which couples the local ensemble transform Kalman filter with the Weather Research and Forecasting model. Its benefits to quantitative precipitation nowcasting (QPN) are evaluated with observing system simulation experiments (OSSEs) and real observation experiments on Typhoon Morakot (2009), which brought record-breaking rainfall and extensive damage to central and southern Taiwan. The purpose is to provide a useful plan of radar data assimilation for improving typhoon rainfall nowcasts in Taiwan, which are challenges due to complex terrain and the lack of in-situ observations over the surrounding sea.

In the OSSEs, the assimilation of radial velocity and reflectivity improves the three-dimensional winds and rain-mixing ratio most significantly because of the direct relations in the observation operator. For QPN, the positive impact of radar data lasts for 6 hours; the performance responds to reflectivity assimilation more quickly than radial velocity assimilation while assimilating both is most recommended. Increasing the observation coverage over upstream convection areas also largely enhances the QPN performance. For multi-scale interactions, we propose a mixed localization method, which yields further improvement.

Our system also improves QPN effectively with real observations. When real reflectivity data are assimilated, the variable localization method must be used to update only the rain mixing ratio. With observation-space statistics, the model bias and ideal ensemble spread can be diagnosed. The mixed localization method, which is more beneficial in the real case, enhances the accuracy of the wind field especially for the areas with sparse or discontinuous radar observations and also improves QPN.

Aksoy, A., D. C. Dowell, and C. Snyder, 2009: A multicase comparative assessment of the ensemble Kalman filter for assimilation of radar observations. Part I: Storm-scale analyses. Mon. Wea. Rev., 137, 1805-1824.
Aksoy, A., D. C. Dowell, and C. Snyder, 2010: A multicase comparative assessment of the ensemble Kalman filter for assimilation of radar observations. Part II: Short-range ensemble forecasts. Mon. Wea. Rev., 138, 1273-1292.
Anderson, J. L., 2001: An ensemble adjustment Kalman filter for data assimilation. Mon. Wea. Rev., 129, 2884-2903.
Anderson, J. L., 2010: A non-Gaussian ensemble filter update for data assimilation. Mon. Wea. Rev., 138, 4186-4198.
Arnold Jr., C. P., and C. H. Dey, 1986: Observing-systems simulation experiments: Past, present, and future. B. Am. Meteor. Soc., 67, 687-695.
Berenguer, M., and I. Zawadzki, 2008: A study of the error covariance matrix of radar rainfall estimates in stratiform rain. Wea. Forecasting, 23, 1085-1101.
Bishop, C. H., B. J. Etherton, and S. J. Majumdar, 2001: Adaptive sampling with the ensemble transform Kalman filter. Part I: Theoretical aspects. Mon. Wea. Rev., 129, 420-436.
Browning, K. A., C. G. Collier, P. R. Larke, P. Menmuir, G. A. Monk, and R. G. Owens, 1982: On the forecasting of frontal rain using a weather radar network. Mon. Wea. Rev., 110, 534-552.
Caine, N., 1980: The rainfall intensity-duration control of shallow landslides and debris flows. Geogr. Ann., 62A, 23-27.
Caya, A., J. Sun, and C. Snyder, 2005: A comparison between the 4DVAR and the ensemble Kalman filter techniques for radar data assimilation. Mon. Wea. Rev., 133, 3081-3094.
Chang, P.-L., P.-F. Lin, B. J.-D. Jou, and J. Zhang, 2009: An application of reflectivity climatology in constructing radar hybrid scans over complex terrain. J. Atmos. Oceanic Technol., 26, 1315-1327.
Chen, B., M. Mu, and X. Qin, 2013: The impact of assimilating dropwindsonde data deployed at different sites on typhoon track forecasts. Mon. Wea. Rev., 141, 2669-2682.
Chen, C.-S., and Y.-L. Chen, 2003: The rainfall characteristics of Taiwan. Mon. Wea. Rev., 131, 1323-1341.
Chen, F., and J. Dudhia, 2001: Coupling an advanced land surface-hydrology model with the Penn State-NCAR MM5 modeling system. Part I: Model implementation and sensitivity. Mon. Wea. Rev., 129, 569-585.
Chien, F.‐C., and H.‐C. Kuo, 2011: On the extreme rainfall of Typhoon Morakot (2009). J. Geophys. Res., 116, D05104, doi: 10.1029/2010JD015092.
Chung, K.-S., I. Zawadzki, M. K. Yau, and L. Fillion, 2009: Short-term forecasting of a midlatitude convective storm by the assimilation of single-Doppler radar observations. Mon. Wea. Rev., 137, 4115-4135.
Dowell, D. C., F. Zhang, L. J. Wicker, C. Snyder, and N. A. Crook, 2004: Wind and temperature retrievals in the 17 May 1981 Arcadia, Oklahoma, supercell: Ensemble Kalman filter experiments. Mon. Wea. Rev., 132, 1982-2005.
Evensen, G., 1994: Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics. J. Geophys. Res., 99, 10143-10162.
Fang, X., and Y.-H. Kuo, 2013: Improving ensemble-based quantitative precipitation forecasts for topography-enhanced typhoon heavy rainfall over Taiwan with a modified probability-matching technique. Mon. Wea. Rev., 141, 3908-3932.
Fang, X., Y.-H. Kuo, and A. Wang, 2011: The impacts of Taiwan topography on the predictability of Typhoon Morakot’s record-breaking rainfall: A high-resolution ensemble simulation. Wea. Forecasting, 26, 613-633.
Gao, J., and M. Xue, 2008: An efficient dual-resolution approach for ensemble data assimilation and tests with simulated Doppler radar data. Mon. Wea. Rev., 136, 945-963.
Germann, U., and I. Zawadzki, 2002: Scale-dependence of the predictability of precipitation from continental radar images. Part I: Description of the methodology. Mon. Wea. Rev., 130, 2859-2873.
Gourley, J. J., R. A. Maddox, K. W. Howard, and D. W. Burgess, 2002: An exploratory multisensor technique for quantitative estimation of stratiform rainfall. J. Hydrometeor., 3, 166-180.
Greybush, S. J., E. Kalnay, T. Miyoshi, K. Ide, and B. R. Hunt, 2011: Balance and ensemble Kalman filter localization techniques. Mon. Wea. Rev., 139, 511-522.
Guzzetti, F., S. Peruccacci, M. Rossi, and C. P. Stark, 2008: The rainfall intensity-duration control of shallow landslides and debris flows: An update. Landslides, 5, 3-17.
Hamill, T. M., and C. Snyder, 2000: A hybrid ensemble Kalman filter-3D variational analysis scheme. Mon. Wea. Rev., 128, 2905-2919.
Hong, S.-Y., Y. Noh, and J. Dudhia, 2006: A new vertical diffusion package with an explicit treatment of entrainment processes. Mon. Wea. Rev., 134, 2318-2341.
Houtekamer, P. L., and H. L. Mitchell, 2005: Ensemble Kalman filtering. Quart. J. Roy. Meteor. Soc., 131, 3269-3289.
Houze, R. A., S. S. Chen, W.-C. Lee, R. F. Rogers, J. A. Moore, G. J. Stossmeister, M. M. Bell, J. Cetrone, W. Zhao, and S. R. Brodzik, 2006: The hurricane rainband and intensity change experiment: Observations and modeling of Hurricanes Katrina, Ophelia, and Rita. Bull. Amer. Meteor. Soc., 87, 1503-1521.
Hunt, B. R., E. J. Kostelich, and I. Szunyogh, 2007: Efficient data assimilation for spatiotemporal chaos: a local ensemble transform Kalman filter. Physica D, 230, 112-126.
Jung, Y., G. Zhang, and M. Xue, 2008a: Assimilation of simulated polarimetric radar data for a convective storm using the ensemble Kalman filter. Part I: Observation operators for reflectivity and polarimetric variables. Mon. Wea. Rev., 136, 2228-2245.
Jung, Y., M. Xue, G. Zhang, and J. M. Straka, 2008b: Assimilation of simulated polarimetric radar data for a convective storm using the ensemble Kalman filter. Part II: Impact of polarimetric data on storm analysis. Mon. Wea. Rev., 136, 2246-2260.
Kain, J. S., 2004: The Kain-Fritsch convective parameterization: An update. J. Appl. Meteor., 43, 170-181.
Kalnay, E., 2003: Atmospheric Modeling, Data Assimilation and Predictability. Cambridge University Press, New York, NY, USA, 340 pp.
Kalnay, E., H. Li, T. Miyoshi, S.-C. Yang, and J. Ballabrera-Poy, 2007: 4-D-Var or ensemble Kalman filter? Tellus A, 59, 758-773.
Kang, J.-S., E. Kalnay, J. Liu, I. Fung, T. Miyoshi, and K. Ide, 2011: “Variable localization” in an ensemble Kalman filter: Application to the carbon cycle data assimilation. J. Geophys. Res., 116, D09110.
Kawabata, T., T. Kuroda, H. Seko, and K. Saito, 2011: A cloud-resolving 4DVAR assimilation experiment for a local heavy rainfall event in the Tokyo metropolitan area. Mon. Wea. Rev., 139, 1911-1931.
Krzysztofowicz, R., W. J. Drzal, T. R. Drake, J. C. Weyman, and L. A. Giordano, 1993: Probabilistic quantitative precipitation forecasts for river basins. Wea. Forecasting, 8, 424-439.
Le Dimet, F.-X., and O. Talagrand, 1986: Variational algorithms for analysis and assimilation of meteorological observations: Theoretical aspects. Tellus A, 38, 97-110.
Liang, J., L. Wu, X. Ge, and C.-C. Wu, 2011: Monsoonal influence on Typhoon Morakot (2009). Part II: Numerical study. J. Atmos. Sci., 68, 2222-2235.
Lin, Y.-L., R. D. Farley, and H. D. Orville, 1983: Bulk parameterization of the snow field in a cloud model. J. Climate Appl. Meteor., 22, 1065-1092.
Lindskog, M., K. Salonen, H. Järvinen, and D. B. Michelson, 2004: Doppler radar wind data assimilation with HIRLAM 3DVAR. Mon. Wea. Rev., 132, 1081-1092.
Liou, Y.-C., T.-C. Chen Wang, Y.-C. Tsai, Y.-S. Tang, P.-L. Lin, and Y.-A. Lee, 2013: Structure of precipitating systems over Taiwan’s complex terrain during Typhoon Morakot (2009) as revealed by weather radar and rain gauge observations. J. Hydrology, 506, 14-25.
Mandapaka, P. V., U. Germann, L. Panziera, and A. Hering, 2012: Can Lagrangian extrapolation of radar fields be used for precipitation nowcasting over complex Alpine orography?. Wea. Forecasting, 27, 28-49.
Marshall, J. S., and W. McK. Palmer, 1948: The distribution of raindrops with size. J. Meteor., 5, 165-166.
Miyoshi, T., 2011: The Gaussian approach to adaptive covariance inflation and its implementation with the local ensemble transform Kalman filter. Mon. Wea. Rev., 139, 1519-1535.
Miyoshi, T., Y. Sato, and T. Kadowaki, 2010: Ensemble Kalman filter and 4D-Var intercomparison with the Japanese operational global analysis and prediction system. Mon. Wea. Rev., 138, 2846-2866.
Montmerle, T., and C. Faccani, 2009: Mesoscale assimilation of radial velocities from Doppler radars in a preoperational framework. Mon. Wea. Rev., 137, 1939-1953.
Murphy, J. M., 1988: The impact of ensemble forecasts on predictability. Quart. J. Roy. Meteor. Soc., 114, 463-493.
NCDR, 2010: Disaster Survey and Analysis of Morakot Typhoon (in Chinese). National Science and Technology Center for Disaster Reduction, New Taipei City, Taiwan, 109 pp.
Nguyen, H. V., and Y.-L. Chen, 2011: High-resolution initialization and simulations of Typhoon Morakot (2009). Mon. Wea. Rev., 139, 1463-1491.
Ott, E., B. R. Hunt, I. Szunyogh, A. V. Zimin, E. J. Kostelich, M. Corazza, E. Kalnay, D. J. Patil, and J. A. Yorke, 2004: A local ensemble Kalman filter for atmospheric data assimilation. Tellus A, 56, 415-428.
Ruiz, J., C. Saulo, and E. Kalnay, 2009: Comparison of methods used to generate probabilistic quantitative precipitation forecasts over South America. Wea. Forecasting, 24, 319-336.
Sasaki, Y., 1958: An objective analysis based on the variational method. J. Meteor. Soc. Japan, 36, 77-88.
Schwartz, C. S., Z. Liu, Y. Chen, and X.-Y. Huang, 2012: Impact of assimilating microwave radiances with a limited-area ensemble data assimilation system on forecasts of Typhoon Morakot. Wea. Forecasting, 27, 424-437.
Skamarock, W. C., J. B. Klemp, J. Dudhia, D. O. Gill, D. M. Barker, M. G. Duda, X.-Y. Huang, W. Wang, and J. G. Powers, 2008: A description of the Advanced Research WRF version 3. National Center for Atmospheric Research, Boulder, CO, USA, 113 pp.
Smith, P. L., Jr., C. G. Myers, and H. D. Orville, 1975: Radar reflectivity factor calculations in numerical cloud models using bulk parameterization of precipitation processes. J. Appl. Meteor., 14, 1156-1165.
Snyder, C., and F. Zhang, 2003: Assimilation of simulated Doppler radar observations with an ensemble Kalman filter. Mon. Wea. Rev., 131, 1663-1677.
Sugimoto, S., N. A. Crook, J. Sun, Q. Xiao, and D. M. Barker, 2009: An examination of WRF 3DVAR radar data assimilation on its capability in retrieving unobserved variables and forecasting precipitation through Observing System Simulation Experiments. Mon. Wea. Rev., 137, 4011-4029.
Sun, J., 2005: Initialization and numerical forecasting of a supercell storm observed during STEPS. Mon. Wea. Rev., 133, 793-813.
Sun, J., and N. A. Crook, 1997: Dynamical and microphysical retrieval from Doppler radar observations using a cloud model and its adjoint. Part I: Model development and simulated data experiments. J. Atmos. Sci., 54, 1642-1661.
Sun, J., M. Xue, J. W. Wilson, I. Zawadzki, S. P. Ballard, J. Onvlee-Hooimeyer, P. Joe, D. Barker, P.-W. Li, B. Golding, M. Xu, and J. Pinto, 2013: Use of NWP for nowcasting convective precipitation: Recent progress and challenges. B. Am. Meteor. Soc., doi: 10.1175/BAMS-D-11-00263.1.
Sun, J., and Y. Zhang, 2008: Analysis and prediction of a squall line observed during IHOP using multiple WSR-88D observations. Mon. Wea. Rev., 136, 2364-2388.
Tai, S.-L., Y.-C. Liou, J. Sun, S.-F. Chang, and M.-C. Kuo, 2011: Precipitation forecasting using Doppler radar data, a cloud model with adjoint, and the Weather Research and Forecasting model: Real case studies during SoWMEX in Taiwan. Wea. Forecasting, 26, 975-992.
Tong, M., and M. Xue, 2005: Ensemble Kalman filter assimilation of Doppler radar data with a compressible nonhydrostatic model: OSS experiments. Mon. Wea. Rev., 133, 1789-1807.
Tsai, C.-C., S.-C. Yang, and Y.-C. Liou, 2014: Improving quantitative precipitation nowcasting with a local ensemble transform Kalman filter radar data assimilation system: Observing system simulation experiments. Tellus A, 66, 21804, doi: 10.3402/tellusa.v66.21804.
Wang, C.-C., H.-C. Kuo, Y.-H. Chen, H.-L. Huang, C.-H. Chung, and K. Tsuboki, 2012: Effects of asymmetric latent heating on typhoon movement crossing Taiwan: The Case of Morakot (2009) with extreme rainfall. J. Atmos. Sci., 69, 3172-3196.
Wang, X., C. Snyder, and T. M. Hamill, 2007: On the theoretical equivalence of differently proposed ensemble-3DVAR hybrid analysis schemes. Mon. Wea. Rev., 135, 222-227.
Weng, Y., and F. Zhang, 2012: Assimilating airborne Doppler radar observations with an ensemble Kalman filter for convection-permitting hurricane initialization and prediction: Katrina (2005). Mon. Wea. Rev., 140, 841-859.
Whitaker, J. S., and T. M. Hamill, 2002: Ensemble data assimilation without perturbed observations. Mon. Wea. Rev., 130, 1913-1924.
Wilks, D. S., 1990: Probabilistic quantitative precipitation forecasts derived from PoPs and conditional precipitation amount climatologies. Mon. Wea. Rev., 118, 874-882.
Willoughby, H. E., F. D. Marks, and R. J. Feinberg, 1984: Stationary and moving convective bands in hurricanes. J. Atmos. Sci., 22, 3189-3211.
Wu, C.-C., T.-H. Yen, Y.-H. Kuo, and W. Wang, 2002: Rainfall simulation associated with Typhoon Herb (1996) near Taiwan. Part I: The topographic effect. Wea. Forecasting, 17, 1001-1015.
Wu, L., J. Liang, and C.-C. Wu, 2011: Monsoonal influence on Typhoon Morakot (2009). Part I: Observational analysis. J. Atmos. Sci., 68, 2208-2221.
Xiao, Q., Y.-H. Kuo, J. Sun, W.-C. Lee, E. Lim, Y.-R. Guo, and D. M. Barker, 2005: Assimilation of Doppler radar observations with a regional 3DVAR system: Impact of Doppler velocities on forecasts of a heavy rainfall case. J. Appl. Meteor., 44, 768-788.
Xiao, Q., and J. Sun, 2007: Multiple-radar data assimilation and short-range quantitative precipitation forecasting of a squall line observed during IHOP_2002. Mon. Wea. Rev., 135, 3381-3404.
Xie, B., and F. Zhang, 2012: Impacts of typhoon track and island topography on the heavy rainfalls in Taiwan associated with Morakot (2009). Mon. Wea. Rev., 140, 3379-3394.
Xie, B., F. Zhang, Q. Zhang, J. Poterjoy, and Y. Weng, 2013: Observing strategy and observation targeting for tropical cyclones using ensemble-based sensitivity analysis and data assimilation. Mon. Wea. Rev., 141, 1437-1453.
Xue, M., M. Tong, and K. K. Droegemeier, 2006: An OSSE framework based on the ensemble square root Kalman filter for evaluating the impact of data from radar networks on thunderstorm analysis and forecasting. J. Atmos. Oceanic Technol., 23, 46-66.
Yang, S.-C., M. Corazza, A. Carrassi, E. Kalnay, and T. Miyoshi, 2009: Comparison of local ensemble transform Kalman filter, 3DVAR, and 4DVAR in a quasigeostrophic model. Mon. Wea. Rev., 137, 693-709.
Yang, S.-C., E. Kalnay, and T. Miyoshi, 2012: Accelerating the EnKF spinup for typhoon assimilation and prediction. Wea. Forecasting, 27, 878-897.
Yang, S.-C., K.-J. Lin, T. Miyoshi, and E. Kalnay, 2013: Improving the spin-up of regional EnKF for typhoon assimilation and forecasting with Typhoon Sinlaku (2008). Tellus A, 65, 20804, doi: 10.3402/tellusa.v65i0.20804.
Zhang, F., C. Snyder, and J. Sun, 2004: Impacts of initial estimate and observation availability on convective-scale data assimilation with an ensemble Kalman filter. Mon. Wea. Rev., 132, 1238-1253.
Zhang, F., Y. Weng, Y.-H. Kuo, J. S. Whitaker, and B. Xie, 2010: Predicting Typhoon Morakot’s catastrophic rainfall with a convection-permitting mesoscale ensemble system. Wea. Forecasting, 25, 1816-1825.
Zhang, F., Y. Weng, J. A. Sippel, Z. Meng, and C. H. Bishop, 2009: Cloud-resolving hurricane initialization and prediction through assimilation of Doppler radar observations with an ensemble Kalman filter. Mon. Wea. Rev., 137, 2105-2125.
Zhang, J., K. Howard, C. Langston, S. Vasiloff, B. Kaney, A. Arthur, S. Van Cooten, K. Kelleher, D. Kitzmiller, F. Ding, D.-J. Seo, E. Wells, and C. Dempsey, 2011: National Mosaic and Multi-Sensor QPE (NMQ) system: Description, results, and future plans. B. Am. Meteor. Soc., 92, 1321-1338.