台风双眼墙形成机制之探讨.doc

球面二維亂流探討 蔡禹明 郭鴻基 國立台灣大學大氣科學研究所 摘要 太陽系中的類木行星，由可見光照片可看出其表面有明顯緯向帶狀的結構，以及大渦旋存在；地球中緯度的高層區域存在著西風噴流，顯示帶狀噴流是行星風系的重要特徵之一。此動力特徵和球面上二維亂流合併的過程，以及背景的行星渦度梯度有關，系統經向的尺度受到萊茵斯尺度（Rhines scale）限制，緯向上可無限制成長，最終形成帶狀的結構，甚至出現如大紅斑的大渦旋。 地球上存在著渦旋與地形，二者皆會造成局部的位渦擾動，不同的地方在於：渦旋會移動，地形則靜止。本文延續前人的研究，使用全球波譜模式地形規則一、前言 航海家2號（Voyager 2）以及位於南半球反時鐘旋轉的大紅斑（great red spot）球面上二維亂流，在物理空間及波譜中，均看出渦度及能量場有呈緯向分佈的現象。圖2（a）年月日0Z的紅外線衛星雲圖僅有些零星對流擾動年月日Z的紅外線衛星雲圖颱風因此我們想上，正壓強旋轉的二維亂流，在考慮黏滯作用的情況下，其基本方程式如下： （1） （2） （1）式為渦度預報方程，（）式則為渦度場和流函數關係式。很小的情況下，動能將近似保守，而渦度擬能將往小尺度串跌（cascade）。這種現象稱為「二維亂流之選擇性衰減」（selective decay of 2-D turbulence）。Kuo et al.（2004）由模式驗證了這個理論。 二維亂流隨時間會合併成較大尺度的渦旋（McWilliams, 1984）。當系統尺度不斷增加，將逐漸感受到背景的行星渦度梯度，而激發羅士比波 Rossby wave）將能量由系統中帶走，則系統尺度無法再繼續成長Rhines（1975）定義這個最大的尺度稱做Rhines scale），並經由尺度分析，得到： （3） 其中是波數，是整個區域中速度場的均方根（root mean square）球體來說，由於只有經向方向上存在行星渦度梯度，因此僅存在於經向方向上，是不等方向性的（aisotropic）。 二維亂流使用一維能量波譜法分析其特性。為單位質量動能的串跌率，由因次分析，可推導出： （4） 接著，定義為單位質量渦度擬能的串跌率，同樣由因次分析，可推導出： （5） 其中與為常數。由（）式及（）式，在雙對數座標（log-log coordinate）上，及的斜率呈現。分析系統的一維能量波譜可幫助我們瞭解系統是否已達成穩定平衡。 Vallis and Maltrud（1993），在淺水模式中植入無經向變化的餘弦函數地形，實驗結果顯示在地形西側出現南風，東側出現北風，此現象可由位渦守恆解釋。Grose and Hoskins（1979），在線性化後的淺水模式中植入理想的緯向風場及地形，可看出在地形下游處出現正渦度，以及一連串向東南方傳遞的波列。 三、模式介紹 本研究使用的模式，球面諧函數全球波譜模式（spherical harmonics global spectral model）。選用全球模式來作模擬，是可以不需要處理邊界條件。 （）是緯向波數，是經向波數，是勒壤得多項式（Legendre polynomials）。球面諧函數是二維拉普拉斯運算子（Laplacian operator）的特徵函數（eigenfunctions），如下所示： （）法（finite difference method）。數值模擬，考慮模擬正確性及計算效率，採用T64，三角形截斷於第64個波。東西方向微分運算使用快速傅立葉轉換（fast Fourier transform, FFT），以32個正弦波和32個餘弦波當作基本函數（basis functions）南北方向微分運算使用，以64個當作基本函數。同時，南北方向的積分方法使用高斯求積法（Gaussian Quadrature），所以南北方向 如圖3所示，T64全球共有個網格點網格點間距（grid spacing）為1.875度，在赤道上約為208公里。時間積分格式使用跳蛙法時間積分間距為00秒圖之間。圖（b）為積分80天後的結果，可清楚看出能量聚集在一個類似漏斗形狀的區域中，顯示動能反串跌到緯向結構中。實驗結果與HR98的實驗結果相同，顯示我們的模式有能力模擬這種能量聚集在緯向帶狀結構的現象。如此我們可進一步探討渦旋及地形對於噴流結構的影響。 四、植入渦旋之實驗 氣候統計資料顯示，每年全球平均有80個颱風，假設每個颱風的生命期大約是8天，則平均每天全球有1.75個颱風。因此，我們想要探討渦旋這種高渦度聚集的系統對於緯向噴流結構的影響。 模式中所使用的正壓方程式為： （8） 等號右手邊第二項為線性曳力項，加入這個項的原因是為了避免動能往大尺度反串跌的過程中，在萊茵斯尺度的地方累積。給定，其對於渦度