網友們如有到 RAMMB 看熱帶氣旋衛星資料，可能都早發現有熱帶氣旋動能的一項. 今日看看其計算方法.
關於該 ACE 指數的定義和分級可看這文章 MWR，指數定義為熱帶氣旋 200 公里內的動能．純粹從影響力來說是非常合理(個人認為 200-250 都是不錯的），因這會忽略了外圍不構成災害的風場，但範圍又足夠大地考慮了一些因大風場而帶來的風暴潮，包括 Katrina:
數據都是從飛機量度得來，因此亦較可靠:
強度和 ACE 的平均比例為 ACE ~ wnd^1.872，較接近(但仍小於)平方．同一強度的熱帶氣旋其 ACE 差異可以很大，這亦正正是這指數具有意義的原因．唯一盲點是一些非常細小但強的熱帶氣旋會得分很低(上圖偏向右下角的那些數據點)．
昨天給出的三個例子，如果也以 200 公里計算 ACE 的話, 結果是 ACE ~ wnd^1.57 和 wnd^1.19, 明顯小於以上的觀測值．這是由於昨天的風場沒有考慮熱帶氣旋移速．如加入 20 km/h 移速的話，則會變成 wnd^1.80 和 wnd^1.30, 前者和觀測平均值接近．
前兩回提及時間和強度可以結合, 構成更能代表熱帶氣旋影響力的指標, 例如 hrs*wnd^2 和 hrs*wnd^3. 這概念不是由本人提出, 事實上國內外一些研究都會考慮到 Accumulated Cyclone Energy(ACE) 和 Power Dissipation Index(PDI), 在計算的時候正正就採用了以上的算法. 今日討論這算法的問題.
前兩回提到與 nos 或 hrs 或 hrs*wnd 相比, hrs*wnd^2 和 hrs*wnd^3 的預報潛在可能都不高. 原因可以包括: 1) 強度不是大尺度如 ENSO 可以定奪, 因此加入強度因素反而減低了相關系數; 2) 該兩指數的算法可能有問題, 違反了一些基本定律(如物理, 數學), 因此其與 ENSO 的相關較低. 以下只討論後者.
以下集中討論 ACE. 顧名思義, ACE 即熱帶氣旋的動能, 但具體計算 ACE 的時候還要有其他考慮, 包括怎樣定義, 和有沒有資料用作計算驗證. 由於歷史風場資料缺乏, 不少氣象學家們就把 ACE 用 hrs*wnd^2 來表示, 這亦間接地只考慮了風眼附近的範圍. 從影響力來說, 只考慮風眼附近範圍問題不大, 因熱帶氣旋帶來的災害還是在風眼附近. 但這會否忽略了某些物理因素? 以下定義 ACE 為整個熱帶氣旋的動能, 看看有何分別.
圖一共考慮了三個風場結構:
實線是一個最高風速 120 km/h (颱風下限), 最大風半徑 30 km, 熱帶氣旋環流 800 km. 可看到烈風圈約 200 km, 應為頗正常的熱帶氣旋. 另外考慮兩個較強的熱帶氣旋 (150 km/h), 最大風半徑分別為 30 km (長虛線) 和 25 km (短虛線), 400 km 外的風速和實線的一樣. 考慮 25 km 是因為較強熱帶氣旋其風眼一般較小.
由於動能為風速的二次方, 而較遠的距離包圍了更大的面積, 因此可知 ACE(R) ~ v(R)^2*R. 總 ACE 就是 ACE = SumOverAllR ACE(R). 因為這關係, 從下圖可看到對 ACE 最大貢獻的半徑不是在風速最大的 25-30 km, 而是外移到約 100 km, 約莫是暴風半徑的位置. 這樣我們可以定斷 ACE 對稍為外圍 (如 100-200 km) 的風速也很敏感, 而強度不足以完全決定 100-200 km 的風速亦即 ACE.
根據以上公式, 若假設整個熱帶氣旋的風速 v 都以某個比例上升, 則 ACE(R) 和 ACE 都應該以平方比例上升. 但這假設並不合理. 舉例, 一個超颱和一個颱風強度比例可能是 1.5, 但在 200 km 處, 颱風的 60 km/h 風速未必就對應為超颱的 60*1.5 = 90 km/h 的風速. 圖一考慮的風速分佈中, 較強的 150 km/h 熱帶氣旋外圍風速增長較少可能才較貼近現實.
換言之, 如以 hrs*wnd^2 來表示 ACE, 較強兩個熱帶氣旋的 ACE 會是原先的 (150/120)^2 = 1.5625 倍, 但這樣會高估. 如以整個風場來計算 ACE, 則兩者和原先的比例只是 1.26 和 1.18, 其量級和 hrs*wnd 比較接近.
結論: 如以整個熱帶氣旋的風場來定義 ACE, 則 hrs*wnd 可能比 hrs*wnd^2 較合理.
以下是 1986-2012 年天文台發出一號風球時各熱帶氣旋的位置:
位置和強度資料大大部分都是從每隔六小時的天文台最佳路徑插值而成, 因此可能有一點誤差. 並不包括新制下離岸仍吹強風時三號轉回一號的情況.
簡簡單單的一張圖也花了半天時間才繪製完成, 雖則無甚可以解讀, 但還是和各位分享. 風迷題:
* 為何從 1986 年起繪劃?
* 哪三個颱風在 800 公里外掛起一號風球?
* 最接近香港的藍點是哪個熱帶氣旋?
* 發出的距離和強度有多大關係? 有沒有臨界(如 600, 700)距離?
任何帶有時空隨機性的過程都有一共通點, 就是要準確預報某一細小區域很困難. 例如要預測香港的降雨必然比預測整個廣東省的困難. 同樣, 預測影響香港的熱帶氣旋也比預測整個西北太平洋的困難. 所以預測可能是一個較大的範圍(如 4-7 個), 而西北太平洋預測的幅度則會相對較小(如 28-32 個).
眾所周知, 影響香港的熱帶氣旋也受厄爾尼諾現象影響, 因此今回看看各"香港四百(HK400)"指數與 Nino 3.4 區水溫的關係. 另外亦參考"香港八百(HK800)", 看看預報可能性有否提高.
以下是水溫和香港四百的線性相關系數. X-軸為月份, 從之前一年的五月至本年的十二月. 各指數均為負值, 並於之前一年年底至本年年初最明顯, 即"香港四百"受年初 ENSO 活動影響. 上回提到 Hrs*Wnd^2(3) 原理上比 Nos 和 Hrs*Wnd 更有效地表示熱帶氣旋帶來的影響, 但其與水溫的相關最低(藍色和綠色), 因此未必是容易預報的項目.
上回圖表顯示自制指數 WMOHK 和 Hrs 接近, 相關達 0.969. 與水溫相關方面, WMOHK(紫色)稍優於 Hrs(橙色). Nos(紅色)與水溫的關係和 WMOHK/Hrs 相約, 於年中時尤其較佳. 但總體而言, 相關系數都並不高, 只略過了 -0.4, 要用作預報可能還要作更深入的分析或採用其他方法.
香港八百方面(下圖), 除 Nos 外, 其餘指數都沒有預期的進步. 與 HK400 時比較, WMOHK 比 Hrs 的優勢更明顯.
西北太平洋的研究顯示年初的 ENSO 現象對風季前期有較大影響, 因此以下只看一至七月的結果, 即水溫和 HK400(Jan-Jul), HK800(Jan-Jul)的線性相關系數. 結果顯示 WMOHK 比起其他指數有較大關系, HK400(Jan-Jul)超過了 -0.5. 而有趣的是 HK800(Jan-Jul)的相關反而沒那麼明顯.
八至十二月的熱帶氣旋活動可能受年中 ENSO 影響, 但只得 Nos(紅色)有較理想結果:
總結:
* Jan-Jul 和 Aug-Dec 結果的分別, 可能顯示 ENSO 的"遙時"(去年年底至年初)和即時(年中)效應對熱帶氣旋有不一樣的影響, 主要展現為遙時效應與自制指數 WMOHK (即熱帶氣旋有多強大)有較大相關, 而即時效應只對熱帶氣旋數目 Nos 有較大影響. 其原因有待進一步探索.
* 由於資料不足, 自制指數只考慮了再分析資料距離熱帶氣旋 300 公里外的風速. 以後如取得更多和更準確的風場數據, 有望進一步提高其代表性和可預測性.
香港天文台每年都會在世界氣象日前後發佈一年的熱帶氣旋預測, 今年預測影響香港的數目為 4-7 個. 心水清或做過相關工作的朋友可能有一些疑問.
影響香港是否就是指發出熱帶氣旋警告? 但(一號)風球的準則數年前也曾作修改, 加上近十年預報誤差大幅減少, 預報員處理風球時可能(理應)和以往不同. 預報這 4-7 時, 如要參考過往發出風球的數據, 會否出現困難?
如不以過往發出風球的數目而改以其他的客觀指標, 有沒有可能? 天文台網頁的統計資料包括進入香港五百公里範圍熱帶氣旋的數目. 五百這數字是隨意, 還是暗藏玄機?
以個數來表示對香港的影響雖則是國際社會多年來的做法, 也廣為普羅大眾所理解和接納, 但其缺點也不少. 第一是沒有考慮熱帶氣旋強度. 4-7 個本身是個頗濶的範圍, 4-7 個熱低和 4-7 個超颱的分別更是可想而知. 第二是沒有考慮氣旋影響的時間. 第三沒有考慮熱帶氣旋的風圈大小.
去年五月時曾探討相關問題(http://hkmet.blogspot.hk/2012/05/blog-post.html), 今回再較詳細地剖析. 資料還是 1979-2010 年 CMA 的最佳颱風路徑, 以下是"香港四百", 即每年進入香港四百公里熱帶氣旋的統計:
圖中 Nos 是個數, Hrs 是在 400 公里內出現的時間, Hrs*Wnd 是時間乘以風速, Hrs*Wnd^2(3) 是時間乘以風速的 2(3) 次方. WMOHK 是自制指數, 大致上是影響時間乘以距離熱帶氣旋中心三百公里的風速.
個數(Nos)和時間(Hrs)上面已作出陳述. Hrs*Wnd 則除了時間外考慮了熱帶氣旋的強度, 越強的熱帶氣旋得到更大權重. Hrs*Wnd^2(3) 則會給強熱帶氣旋更大的權重.
Hrs*Wnd^2 甚至是 Hrs*Wnd^3, 原理上都應該比 Nos, Hrs 和 Hrs*Wnd 更為合理地描述熱帶氣旋帶來的影響. 一來在物理學上, 我們有動能 K.E. ~ v^2. 二來從直觀上, 60 海浬風速帶來的破壞和 30 海浬風速帶來的破壞也不會只是一倍, 即 Hrs*Wnd 給予強度的權重還是不夠. 從上表所看, 個數和 Hrs*Wnd^2(3) 的相關系數只為 0.711 和 0.629, 很明顯單以個數 Nos 來表示熱帶氣旋帶來的影響會有很大的缺憾.
再舉一個例子:
以上是 Nos 和 Hrs*Wnd^2 的線性關係. 1981 和 1986 年同樣有 8 個熱帶氣旋進入香港 400 公里範圍, 但 Hrs*Wnd^2 1986 是 1981 的 2.58 倍. 而時數 Hrs 為 204 和 354, 即只為 1.74 倍. 因此可解讀為 1986 年熱帶氣旋逗留在 400 公里範內的時間較長(可以是移動較慢, 也可以是更為接近香港), 而且強度較大.
香港十一月通常是所謂秋高氣爽的好日子, 但不少人都注意到今年天氣反常, 潮濕多雲有雨, 又是時候淺談當中可能的原因.
先看究竟有多異常. 判別秋冬季是否好天, 香港上空是否吹偏南風為一重要因素, 因此先抓了 NCEP 再分析資料 1948-2012 共 65 個 11 月的 850 百帕南風數據查看. 65 年當中只有 9 年(14%) 為偏南風. 今年 11 月果真是 65 年來最大偏南風的一年, 緊接其後的是 1997 (厄爾尼諾年), 但 97 強度還是差了一大截. 如以標準差表示, 則今年為 +2.6, 應是頗為罕有的個案.
天氣/氣候異常定必聯想到 ENSO. 從 CPC 網站下載 1950 年至今的海水溫度數據(Nino3.4),並看看今年和以往哪些年比較相似. 發現最接近的為 2006 年:
2006: -0.9,-0.7,-0.5,-0.3,0.0,0.1,0.2,0.3,0.5,0.8
2012: -0.9,-0.7,-0.5,-0.3,-0.1,0.0,0.1,0.3,0.4,0.6
再從黃家香港氣象局網頁劃圖, 比較 2012 年和其他各年的差異. 肉眼發現和 2012 年北半球的異常形態最接近的也是 2006 年(共查看了 1971-2011). 注意左上方北半球的異常形態:
2006 年 11 月:
2012 年 11 月:
結論: 雖然今年最終並沒有發展為明顯的厄爾尼諾年, 但以上資料顯示今年可能和以往厄爾尼諾年(包括最像的 2006 微弱偏北風年)偏暖偏濕相似, ENSO 或赤道水溫或可解釋到一部分原因. 還未看各中心之前的月天氣/氣候預測, 估計今年 11 月的情況應該都在模式的掌握之中.
西北太平洋一片寧靜之際(預測下周有野橫過南海中部), 美國東北岸登陸的 Sandy 可算為風迷帶來一點生氣, 也是時候更新網誌了.
看了一些分析和文件後, 暫未看到 Sandy 的結構有何特殊之處, 反而溫帶氣旋轉回偏西登陸東岸可能比較罕有. 這或許和東面的阻塞現象有關: Rex Block, ECMWF
從歐中網站下載了一些資料, 看看 Sandy 的結構和轉變:
高度場顯示在登陸前的十二小時, 各層上的氣旋都加深了, 在中上層(200 至 500 毫巴)似乎和西北面的槽混合. 登陸後和熱帶氣旋一樣, 中心填塞消散(圖略).
風場亦顯示氣旋加深及和槽的混合. 例如在 500 毫巴, 三十日零時的環流中心明顯較十二小時前清晰, 最大風速亦較高. 在低層 850 百帕也隱約看到環流整合, 最大風速雖然下跌, 但風圈擴大了, 整體上還是稍為加深.
值得一提是 warm seclusion 現象, 即當溫帶氣旋發展到最後階段時, 低層中心附近被冷空氣包圍, 相對來說是暖的. 翻了一些文件, 才發現自己對溫帶氣旋的認識仍停留在上世紀初挪威學派 Bjerknes 的模型. Bjerknes 的模型雖然在任何氣象學教材上還會提及, 但原來早於 1990 年已有更準確的 Shapiro-Keyser 模型. Bjerknes 模型雖也有 warm seclusion 這現象, 但和 Shapiro-Keyser 模型還是有較大分別.
溫度場可看到低層冷空氣約在登陸前十二小時入侵氣旋中心. 在 925 和 850 百帕上, 二十九日十二時仍可看到氣旋上空的暖中心, 三十日零時則不見了. 再早一點(從二十五日開始, 圖略)可看到低層暖中心逐步減弱, 因此歐中分析圖上似乎看不到明顯的 warm seclusion.
最後送上位渦圖(這圖是等位溫面, 高度約在 300 毫巴), 可很清楚地看到氣旋和槽的混合, 並開始出現在登陸前十二小時:
