以作者查詢圖書館館藏 、以作者查詢臺灣博碩士 、以作者查詢全國書目 、勘誤回報 、線上人數:225 、訪客IP:52.14.17.240
姓名 陳盈文(Ying-wen Chen) 查詢紙本館藏 畢業系所 大氣物理研究所 論文名稱 使用WRF 3DVAR 及 4DVAR 同化虛擬位渦渦旋對颱風數值模擬之影響
(The Impact of Potential Vorticity Bogus Vortex Data Assimilation of Typhoon Using WRF 3DVAR and 4DVAR)相關論文 檔案 [Endnote RIS 格式] [Bibtex 格式] [相關文章] [文章引用] [完整記錄] [館藏目錄] [檢視] [下載]
- 本電子論文使用權限為同意立即開放。
- 已達開放權限電子全文僅授權使用者為學術研究之目的,進行個人非營利性質之檢索、閱讀、列印。
- 請遵守中華民國著作權法之相關規定,切勿任意重製、散佈、改作、轉貼、播送,以免觸法。
摘要(中) 在進行颱風預報時,初始場的誤差會導致颱風隨時間模擬而降低路徑的準確性,為了能夠更有效地提高颱風模擬的準確性,提高初始場的準確性為首當其衝的必要條件。位渦反演法可以有效地將颱風的位渦及切向風場模擬而出,且虛擬位渦擾動可依設定而反演出颱風的擾動位渦和切向風擾動場,因此,可將位渦的擾動場以虛擬位渦擾動來取代再進行位渦反演,則可獲得新的虛擬位渦渦旋(Potential Vorticity Bogus Vortex),並植入於所要探討的颱風初始場當中,來進行颱風初始場的修正。使用不同的背景場誤差和水平影響尺度,對於不同的颱風模擬也會有不同的結果產生。
在WRF 3DVAR的部分,對於輕度颱風凡那比而言,其初始場的颱風及駛流場強度較弱,因此在進行同化時必須同化較大的虛擬位渦渦旋,將環境流場也同化至初始場中,才能改變環境流場使颱風受駛流場的牽引讓颱風路徑產生偏移,使模擬路徑更接近於實際觀測,由於背景誤差cv3為使用NCEP GFS資料所得到的全球模式的誤差分析,會使得修正範圍較大使得環流較向外發散,能改變的環境流場也相對較多。因此,對於強度較弱的颱風而言,使用較大的同化半徑、較大的水平影響尺度和背景誤差cv3,則可以較有效地模擬出輕度颱風的行徑路徑。相對於初始場強度較強的強烈颱風梅姬而言,進行同化時虛擬位渦渦旋半徑不可同化太大,約只要同化RMW的大小,才不會讓颱風本身受到過多過強的環境流場影響而有多餘的偏移現象發生,而背景誤差cv5和cv3剛好相反,為針對實際模擬範圍網格點進行長達一個月的區域模式誤差統計結果,修正範圍較集中於植入渦旋的位置,使颱風流場較集中於颱風中心,讓颱風不受到過多的外圍流場來干涉路徑移動。因此,對於強度較強的颱風而言,使用較小的同化半徑、較小的水平影響尺度和cv5則可以較有效的模擬出強烈颱風的行徑路徑。
將凡那比颱風與梅姬颱風在WRF 3DVAR中,路徑模擬誤差較小的初始場設定改以WRF 4DVAR模擬其初始場及路徑。在輕度颱風凡那比的部分,初始場結果對於Ctrl Run而言其等壓線沒有明顯的變化,風場在同化半徑大小內為輻合,特別是眼牆的位置上最為強烈,當背景誤差為cv3時,則颱風的東北側會有明顯的增溫現象發生。和3DVAR相比,風場為輻散狀態,且最強處位於眼牆區域,等壓線向外擴散使得初始場模擬結果較為鬆散且範圍較為大,其路徑模擬結果均比3DVAR來的差、路徑誤差都較大,特別是在背景誤差為cv3時,雖然有將初期路徑特徵模擬而出,但其颱風行徑較緩慢且無登陸台灣,和實際觀測及3DVAR模擬有明顯的誤差。而強烈颱風梅姬的部分,和Ctrl Run相比,無論背景誤差為cv3或cv5,其初始場均有較強的氣旋式流場呈現,但氣壓場和溫度場均無明顯的變化產生。和3DVAR相比,在同化半徑大小內風場有些許的向外輻散,且當背景誤差為cv3時,颱風中心處有0.5度的降溫,使得颱風強度相較於3DVAR而言較微弱。而路徑模擬結果均無法勝過於3DVAR的初始場,整體路徑誤差均大於100公里,特別是在登陸呂宋島之後,其誤差更為明顯。
在此研究個案中,對於同化虛擬位渦渦旋於初始場中,使用WRF 3DVAR所模擬出的路徑不但可以修正Ctrl Run的誤差,且相對於WRF 4DVAR而言,在有限的電腦資源及有效的時間內可得到較好的模擬結果,有助於颱風預報的即時應用。摘要(英) On the typhoon prediction, the error of initial field will affect the accuracy of tracks along with the time. If we want to increase the accuracy of prediction, the first necessary factor is to improve the initial field accuracy of typhoon. By using different settings of potential vorticity Inversion could retrieved the potential vorticity perturbation and tangential wind perturbation field of typhoon. Therefore, we could get the new potential vorticity bogus vortex (PV Bogus Vortex) by using the potential vorticity perturbation field which is substituted for the original perturbation field. Then insert it to correct the initial field. It will have different simulate result by using different background error and horizontal scale.
In the WRF-3DVAR, the Fanapi typhoon’s intensity and steering flow are weak, so assimilate the bigger radius of PV Bogus Vortex and environmental field can correct the typhoon track so that the steering flow could be changed and the track would be more realistic. The outward divergence and larger correct range are using background error-cv3 to correct broad environment field. For the intensity of weaker typhoon, we can simulate the track efficiently by using the bigger PV Bogus Vortex and horizontal scale on the background error-cv3. It will over correct the strong typhoon’s track, because more environmental wind field will be alter when we assimilate the PV Bogus Vortex size of RMW and smaller horizontal scale. The background error-cv5 are opposed to the background error-cv3, it’s because cv5 is a regional model statistics results of a month on real simulate domain. The cv5 is focus on the typhoon center to correct the flow, which will make environment field not to be interfered with the insert position. For the intensity of stronger typhoon, we can simulate the track efficiently by using the smaller PV Bogus Vortex and horizontal scale on the background error-cv5.
The minimum track error of WRF-3DVAR of typhoon Fanapi and Megi will be replaced by WRF-4DVAR. In the part of Fanapi typhoon, to compare with the Ctrl Run on the initial field, the isobar don’t have obvious difference, but they will have stronger convergence wind field on the eyewall, and temperature of the north-east side of typhoon will increase. To compare with WRF-3DVAR on the initial field, the simulation result is looser and wider sphere of influence that isobar and wind field are divergence to effect the out region of typhoon, so it makes the track have larger error and doesn’t make landfall in Taiwan. In the part of Megi typhoon, to compare with the initial field of WRF-3DVAR, the simulation result shows the stronger cyclonic stream flow, but don’t have temperature change no matter background error-cv3 or cv5. To compare with WRF-3DVAR on the initial field, the simulation results is weaker than WRF-3DVAR that it have divergence wind pattern on the radius of typhoon, and decrease 0.5 Celsius degree on the typhoon center. The track simulation result can’t be better than WRF-3DVAR which has over 100 km track error when it land on Philippines.
On the part of assimilate PV Bogus Vortex, It can correct the track error of the Ctrl Run by using the WRF-3DVAR which makes the result better than WRF-4DVAR.
關鍵字(中) ★ 虛擬位渦渦旋
★ 背景場誤差關鍵字(英) ★ PV Bogus
★ Background error論文目次 中文摘要 …………………………………………………………………………… I
英文摘要 ………………………………………………………………………… III
誌謝 ………………………………………………………………………… IV
目錄 ………………………………………………………………………… VI
圖表目錄 …………………………………………………………………………VIII
符號說明 …………………………………………………………………………XVII
第一章 緒論 …………………………………………………………………………… 1
1-1 前言 …………………………………………………………………………… 1
1-2 前人研究 ……………………………………………………………………… 1
1-3 研究動機及目的 ……………………………………………………………… 4
第二章 個案與研究方法 ………………………………………………………………… 6
2-1 個案介紹 ……………………………………………………………………… 6
(1) 凡那比颱風 ……………………………………………………………… 6
(2) 梅姬颱風 ………………………………………………………………… 6
2-2 位渦反演方程式 ……………………………………………………………… 7
2-3 背景誤差 ……………………………………………………………………… 9
第三章 模式介紹與實驗設計 ………………………………………………………… 11
3-1 模式介紹 …………………………………………………………………… 11
3-2 模式設計 …………………………………………………………………… 12
3-3 實驗設計 …………………………………………………………………… 13
第四章 虛擬位渦反演與同化對颱風初始場之影響 ………………………………… 14
4-1 虛擬位渦擾動設定 ………………………………………………………… 14
4-2 虛擬位渦反演結果 ………………………………………………………… 14
4-3 虛擬位渦反演渦漩同化 …………………………………………………… 17
第五章 凡那比颱風模擬實驗 ………………………………………………………… 19
5-1 同化虛擬位渦渦旋 - 初始場修正與颱風模擬分析 ……………………… 19
(1) 三維變分資料同化 – cv3水平影響尺度實驗 ……………………… 20
(2) 三維變分資料同化 – cv5水平影響尺度實驗 ……………………… 22
5-2 四維變分資料同化和背景誤差cv3及cv5實驗 …………………………… 25
第六章 梅姬颱風模擬實驗 …………………………………………………………… 30
6-1 同化虛擬位渦渦旋 - 初始場修正與颱風模擬分析 …………………… 30
(1) 三維變分資料同化 – cv3水平影響尺度實驗 ……………………… 31
(2) 三維變分資料同化 – cv5水平影響尺度實驗 ……………………… 34
6-2 四維變分資料同化和背景誤差cv3及cv5實驗 …………………………… 39
第七章 總結與未來展望 …………………………………………………………… 44
參考文獻 ………………………………………………………………………………… 46
附錄一 …………………………………………………………………………………… 50
附錄二 …………………………………………………………………………………… 52
附錄三 …………………………………………………………………………………… 54
附錄四 …………………………………………………………………………………… 55
附表與附圖 ……………………………………………………………………………… 56參考文獻 林妗庭,2012:WRF 4DVAR 虛擬渦旋及雷達資料同化對於颱風凡那比(2010)數值模擬之影響。國立中央大學大氣物理研究所碩士論文,1-80。
黃子茂,2012:同化虛擬位渦反演渦旋對凡那比颱風初始場模擬及路徑遇到之影響。國立中央大學大氣物理研究所碩士論文,1-59。
洪于珺,2010:颱風辛樂克(2008) WRF 模擬及位渦反演之研究。國立中央大學大氣物
理研究所碩士論文,1-83。
劉豫臻,2009:聖帕颱風模擬的位渦反演之診斷分析。國立中央大學大氣物理研究所
碩士論文,1-85。
Davis, C. A., 1992: Piecewise potential vorticity inversion. J. Atmos. Sci., 49, 1397–1411.
Davis, C. A., and K. A. Emanuel, 1991: Potential vorticity diagnostics of cyclogenesis. Mon.
Wea. Rev., 119, 1929-1953.
Davis, C. A., E. D. Grell, and M. A. Shapiro, 1996: The balanced dynamical nature of a
rapidly intensifying oceanic cyclone. Mon. Wea. Rev., 124, 3-26.
Eliassen, A., 1952: Slow thermally or frictionally controlled meridional circulation in a
circular vortex. Astrophys. Norv., 5, 19-60.
Franklin, J. L., S. J. Lord, S. E. Feuer, and F. D. Marks, 1993: The kinematic structure of
Hurricane Gloria (1985) determined from nested analyses of dropwindsonde and
Doppler radar data. Mon. Wea. Rev., 121, 2433–2451.
Hoskins, B. J., and F. P. Bretherton, 1972: Atmospheric frontogenesis models: Mathematical
formulation and solution. J. Atmos. Sci., 27, 11-37.
Hoskins, B. J., M. E. McIntyre, and A. W. Robertson, 1985: On the use and significance of
isentropic potential vorticity maps. Quart. J. Roy. Meteor. Soc., 111, 877-946.
Huang, X.-Y., Q. Xiao, W. Huang, D. M Barker, J. Michalakes, J. Bray, Z. Ma, Y.-R., Guo,
H.-c.Lin, and Y.-H.Kuo,2006.Preliminary results of WRF 4D-Var. WRF user’s workshop, Boulder, Colorado, 19-22 June 2006.
Huang, X.-Y., Q. Xiao, D.M. Berker, X. Zhang, J. Michalakes, W. Huang, T. Henderson, J.
Bray, Y. Chen, Z. Ma, J. Dudhia, Y. Guo, X. Zhang, D.J. Won, H.C. Lin, and Y .H. Kuo,
2009: Four-Dimensional Variational Data Assimilation for WRF: Formulation and Preliminary Results. Mon. Wea. Rev., 137, 299-314.
Huo, Z.-H., D.-L. Zhang, and J. R. Gyakum, 1998: An Application of Potential Vorticity
Inversion to Improving the Numerical Prediction of the March 1993 Superstorm. Mon.
Wea. Rev., 126, 424-436.
Kieu, C. Q. 2005: Piecewise Potential Vorticity Inversion. Master thesis, pp. 81, Department
of Atmospheric and Oceanic Science, University of Maryland, USA.
Kieu, Chanh Q., Da-Lin Zhang, 2010: A Piecewise Potential Vorticity Inversion Algorithm
and Its Application to Hurricane Inner-Core Anomalies. J. Atmos. Sci., 67, 2616–2631.
Kurihara, Y., N. A. Bender, and R. J. Ross, 1993: An initialization scheme of hurricane models
by vortex specification. Mon. Wea. Rev., 121, 2030-2045.
Kurihara, Y., M. A. Bender, R. E. Tuleya, and R. J. Ross, 1995: Improvements in the GFDL
hurricane prediction system. Mon. Wea. Rev., 123, 2791-2801.
Liu, Y. D., D.-L. Zhang, and M. K. Yau, 1999: A multiscale numerical study of Hurricane
Andrew (1992). Part II: Kinematics and inner-core structures. Mon. Wea. Rev., 127,
2597–2616.
Marks, F. D., R. A. Houze, and J. F. Gamache, 1992: Dual-aircraft investigation of the inner
core of Hurricane Norbert. Part I: Kinematic structure. J. Atmos. Sci., 49, 919–942.
Marks, F. D., P. G. Black, M. T. Montgomery, and R. W. Burpee, 2008: Structure of the Eye
and Eyewall of Hurricane Hugo (1989). Mon. Wea. Rev., 136, 1237–1259.
Möller, J. D., and S. C. Jones, 1998: Potential vorticity inversion for tropical cyclones using
the asymmetric balance theory. J. Atmos. Sci., 55, 259–282.
Möller, J. D., and L. J. Shapiro, 2002: Balanced contributions to the intensification of
Hurricane Opal as diagnosed from a GFDL model forecast. Mon. Wea. Rev., 130,
1866–1881.
Olsson, P. Q., and W. R. Cotton, 1997: Balanced and unbalanced circulations in a primitive
equation simulation of a midlatitude MCC. Part II: Analysis of balance. J. Atmos. Sci.,
54, 479–497.
Reasor, P. D., M. T. Montgomery, F. D. Marks Jr., and J. F. Gamache, 2000:
Low-wavenumber structure and evolution of the hurricane inner core observed by
airborne dual-Doppler radar. Mon. Wea. Rev., 128, 1653–1680.
Schubert, W. H., S. A. Hausman, M. Garcia, K. V. Ooyama, and H.-C. Kuo, 2001: Potential
vorticity in a moist atmosphere. J. Atmos. Sci., 58, 3148–3157.
Shapiro, L. J., and J. L. Franklin, 1995: Potential vorticity in Hurricane Gloria. Mon. Wea.
Rev., 123, 1465–1475.
Shapiro, L. J., and M. T. Montgomery, 1993: A three-dimensional balance theory for rapidly
rotating vortices. J. Atmos. Sci., 50, 3322–3335.
Shapiro, L. J., 1996: The motion of Hurricane Gloria: A potential vorticity diagnosis. Mon.
Wea. Rev., 124, 2497–2508.
Wang, X., and D.-L. Zhang, 2003: Potential vorticity diagnosis of a simulation hurricane.
Part I: formulation and quasi-balanced flow. J. Atmos. Sci., 60, 1593–1607.
Willoughby, H. E., 1990: Gradient balance in tropical cyclones. J. Atmos. Sci., 47, 265-274.
Wu, C.-C., and K. A. Emanuel, 1995a: Potential vorticity diagnostics of hurricane
movement. Part I: A case study of Hurricane Bob (1991). Mon. Wea. Rev., 123,
69–92.
Wu, C.-C., and K. A. Emanuel, 1995b: Potential vorticity diagnostics of hurricane
movement. Part II: Tropical storm Ana (1991) and Hurricane Andrew (1992). Mon.
Wea. Rev., 123, 93–109.
Xiao, Q., X. Zou, and B. Wang, 2000: Initialization and simulation of a landfalling hurricane
using a variational bogus data assimilation scheme. Mon. Wea .Rev., 128, 2252-2269.
Zhang, D.-L., and C. Q. Kieu, 2006: Potential vorticity diagnosis of a simulation hurricane.
Part II: Quasi-balanced contributions to forced secondary circulations. J. Atmos. Sci.,
63, 2898-2914.
Zhang, D.-L., Y. Liu, and M.-K. Yau, 2000: A multiscale numerical study of hurricane
Andrew (1992). Part III: Dynamically induced vertical motion. Mon. Wea. Rev., 128,
3772-3788.
Zhang, D.-L., Y. Liu, and M. K. Yau, 2001: A multiscale numerical study of Hurricane
Andrew (1992). Part IV: Unbalanced flows. Mon. Wea. Rev., 129, 92-107.
Zhang, D.-L., Y. Liu, and M. K. Yau, 2002: A multiscale numerical study of Hurricane
Andrew (1992). Part V: Inner-core thermodynamics. Mon. Wea. Rev., 130, 2745–27指導教授 黃清勇(Ching-yuang Huang) 審核日期 2013-7-15 推文 facebook plurk twitter funp google live udn HD myshare reddit netvibes friend youpush delicious baidu 網路書籤 Google bookmarks del.icio.us hemidemi myshare