使用系集卡曼濾波器同化都卜勒雷達資料之研究

Doppler Radar Data Assimilation Using Ensemble Kalman Filter

黃國禎(Guo-Jhen Huang)

指導教授 ： 廖宇慶

國立中央大學/地球科學學院/大氣物理研究所/碩士(2007年)

摘要

系集卡曼濾波器（EnKF）（Ensemble Kalman Filter）是一種資料同化方法。本研究是利用EnKF技術以及觀測系統模擬實驗（OSSE）（Observation System Simulation Experiments）來測試EnKF同化都卜勒雷達資料的表現。本研究使用的模式為ARPS（Advanced Regional Prediction System），此模式是ㄧ個包含多種微物理參數的非靜力可壓縮模式。模擬個案為1977年5月20日發生在美國中部的Del City Storm。本研究在無地形及有地形情況下，設計多組實驗，如「同化系集數目」、「同化時間間隔」、「同化次數」與「同化區域」等若干實驗，探討其對模式各氣象場預報誤差的影響，並且著重在定量降水預報度的改善。結果顯示，以40個系集樣本，每5分鐘同化一次，可得到足夠準確的結果。配合雷達觀測，盡量多次同化，可有助於誤差的降低，並避免誤差擴大而無法修正回來。在有地形時，只要在過山之前就同化初期暴風發展的資訊，便能快速地掌握住暴風的位置與型態，使後續的預報更加準確。同化一個小時雷達資料以後，此系集模式的降水預報能力約可維持至少一小時的準確度。使用EnKF同化雷達資料的確對降低預報誤差具有顯著的改進。

關鍵字

觀測系統模擬實驗 ； 都卜勒雷達 ； 資料同化 ； 系集卡曼濾波器

並列摘要

Ensemble Kalman filter (EnKF) is a method for data assimilation. In this study, we apply the Observation System Simulation Experiments (OSSE) type of experimental designs to explore the performance of assimilating Doppler radar data using EnKF. A general purpose non-hydrostatic compressible model, the ARPS (Advanced Regional Prediction System), with complex multi-class microphysics, is employed for conducting all the experiments. Artificial data sets are from a simulated classic storm case that occurred on 20 May 1977 in Del City, Oklahoma. With and without terrain, we investigate the impact of several factors on the model forecasts, with the emphasis on the issue of quantitative precipitation forecast (QPF). These factors consider the number of ensembles, the time interval and frequency of data injection, the area of data availability, and so on. The major results show that using 40 members, and assimilating the radar data once every 5 minutes, can effectively produce the forecasts with sufficient accuracy. Assimilating as many data sets as possible can help to reduce the errors, and prevent the errors from growing to an uncontrollable scale. When the terrain is present and becomes a potential blockage to the radar beams, and if one can assimilate into the model the information of the initial storm development before the storm reaches the lee side of the mountain, then it is still possible to catch the location and pattern of the storm. Such a measure makes the following model forecast maintain its accuracy, even after the storm passes the mountain, and arrives at a region where the radar beams are completely blocked. Finally, to obtain an accurate one-hour QPF also requires an one-hour of radar data assimilation. Overall speaking, the assimilation of Doppler radar data does reveal significant improvements on reducing the forecast errors.

並列關鍵字

ensemble Kalman filter ； Doppler radar ； observation system simulation experiments ； data assimilation

參考文獻

紀博庭，2005：利用中央大學雙偏極化雷達資料反求雨滴粒徑分佈及降雨率方法的研究，國立中央大學大氣物理碩士論文，70頁。

呂崇華，2006：雙偏極化雷達資料分析梅雨鋒面雨滴粒徑分佈的物理特性，國立中央大學大氣物理研究所碩士論文，100頁。

Crook, N.A., and J. Sun, 2002: Assimilating radar, surface, and profiler data for the Sydney 2000 Forecast Demonstration Project. J. Atmos. ceanic Technol., 19, 888–898.

——, ——, 2004: Analysis and forecasting of the low-level wind during the Sydney 2000 Forecast Demonstration Project. Wea. Forecasting, 19, 151–167.

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.