劉 樂，李耀東

(1.中國人民解放軍理工大學氣象海洋學院，江蘇南京 211101；2.北京航空氣象研究所，北京 100085)

?

RTTOV快速輻射傳輸模式應(yīng)用研究進展

劉 樂1，2，李耀東2

(1.中國人民解放軍理工大學氣象海洋學院，江蘇南京 211101；2.北京航空氣象研究所，北京 100085)

摘要RTTOV(Radiative Transfer for TOVS)是近期出現(xiàn)的一個比較優(yōu)秀的快速輻射傳輸模式。該研究介紹了RTTOV工作的原理和方法，從提高數(shù)值預(yù)報模式準確度和模擬衛(wèi)星輻射亮溫2個方面出發(fā)，綜述了RTTOV模式的應(yīng)用研究進展。結(jié)果顯示，利用RTTOV將衛(wèi)星資料直接同化到數(shù)值預(yù)報模式中，能顯著提高數(shù)值預(yù)報的準確度；同時，RTTOV能夠較好地模擬晴空條件下衛(wèi)星接收到的輻射亮溫，有云條件下的模擬值與實測值之間的誤差大于晴空條件。利用RTTOV模擬得到的亮溫進行仿云圖制作，在一定程度上可以判斷云體和云系的移動及發(fā)展變化，進而改善氣象預(yù)報和保障業(yè)務(wù)。

關(guān)鍵詞RTTOV；快速輻射傳輸模式；同化；正演；亮溫

The Application and Research Advances of RTTOV Fast Radiative Transfer Mode

LIU Le1,2, LI Yao-dong2(1.Meteorological Oceanographic Institute, PLA University of Science and Technology, Nanjing, Jiangsu 211101;2. Beijing Aviation Meteorological Institute, Beijing 100085)

AbstractRTTOV (Radiative Transfer for TOVS) is a good fast radiative transfer model. The working principle and method of RTTOV was introduced, the application and research advances were reviewed from two aspects of improving the numerical forecast’s accuracy and simulating satellite radiation TBB. The results showed that using RTTOV to assimilate satellite data directly into the numerical prediction model can obviously improve the accuracy of numerical prediction; meanwhile, RTTOV can better simulate satellite accepted TBB under clean sky condition, the error between simulated value and measured value under the cloud condition is greater that in clear sky condition. Using TBB can make satellite cloud images, determine the movement and change of clouds, cloud system in a certain extent, further improve weather forecasting and security business.

Key wordsRTTOV; Fast radiative transfer mode;Assimilation; Forward modeling; TBB

數(shù)值天氣預(yù)報(Numerical Weather Prediction)的準確度主要由數(shù)值模式本身的完善程度和數(shù)值模式初始場的準確度來決定。目前，數(shù)值預(yù)報模式本身已經(jīng)不斷趨于完善，而模式初始場的準確度還有很大的提升空間。數(shù)值模式的初始場通常由觀測資料提供，觀測資料包括常規(guī)觀測資料和非常規(guī)觀測資料，其中，衛(wèi)星資料就屬于非常規(guī)觀測資料中很重要的一項。目前世界上各個主要天氣預(yù)報中心的全球資料同化系統(tǒng)所采用的觀測資料中，衛(wèi)星資料的數(shù)量超過了90%[1]。如果數(shù)值預(yù)報中不使用衛(wèi)星資料，預(yù)報準確率將下降23.4%[2]。在世界先進的中期數(shù)值天氣預(yù)報中，衛(wèi)星資料同化已起到十分關(guān)鍵的作用[3]。

與常規(guī)觀測不同，衛(wèi)星探測量一般為輻射率，它不是數(shù)值模式的變量，無法為模式所直接使用。衛(wèi)星資料同化有間接同化和直接同化2種方式，間接同化即先由衛(wèi)星探測資料反演出大氣溫度和濕度等大氣參數(shù)的垂直廓線，然后同化反演結(jié)果，其反演過程與數(shù)據(jù)同化過程是相互獨立、分別進行的；直接同化即利用快速輻射傳輸模式正演得到的輻射亮溫值與直接引入數(shù)值預(yù)報模式變分同化系統(tǒng)中的實際衛(wèi)星輻射亮溫值進行對比，要使正演的亮溫與實際值趨于接近，就必須不斷調(diào)整初始場的大氣溫度、濕度等大氣廓線，經(jīng)過調(diào)整的大氣廓線更為準確，這樣不僅能從觀測中獲取更多的信息，且避免了復雜的衛(wèi)星資料反演計算及其帶來的誤差。所有基于衛(wèi)星探測輻射信息的同化系統(tǒng)均利用預(yù)報場或無線電探空廓線作為初始場，正演模擬輻射亮溫的方案。作為同化過程的一部分，利用輻射傳輸模式進行正演計算是必須的。因此，運用RTTOV(Radiative Transfer for TOVS)快速輻射傳輸模式和數(shù)值天氣預(yù)報模式產(chǎn)品資料，可以模擬衛(wèi)星輻射亮溫，進行仿云圖制作，從衛(wèi)星云圖上得到云的有關(guān)信息，在天氣預(yù)報業(yè)務(wù)上有重要作用?；谝陨侠碚撘罁?jù)，RTTOV快速輻射傳輸模式的應(yīng)用主要分為2個部分：①在數(shù)值預(yù)報初始場中直接同化衛(wèi)星輻射亮溫，提高預(yù)報準確率；②將快速輻射傳輸模式單獨使用，用于模擬衛(wèi)星輻射亮溫，進行仿云圖的制作。筆者介紹了RTTOV工作的原理和方法，從提高數(shù)值預(yù)報模式準確度和模擬衛(wèi)星輻射亮溫2個方面出發(fā)，綜述了RTTOV模式的應(yīng)用研究進展。

1RTTOV模型介紹

RTTOV系列是20世紀90年代初由歐洲中期天氣預(yù)報中心(European Centre for Medium-Range Weather Forecasts，ECMWF)開發(fā)的快速輻射傳輸模式，最初用于模擬TOVS(NOAA’s TIROS Operational Vertical Sounder，TIROS 業(yè)務(wù)垂直探測儀)衛(wèi)星，經(jīng)過多年的發(fā)展，現(xiàn)在已經(jīng)能夠模擬全世界幾十種氣象衛(wèi)星探測的地球環(huán)境紅外和微波輻射。RTTOV7以后的版本均可以模擬大氣紅外探測器(Atmospheric InfraRed Sounder,AIRS)等新一代紅外高光譜分辨率衛(wèi)星探測器資料[4]。目前已經(jīng)更新至RTTOV11，該版本包含了云散射模塊(RTTOV_SCATT)，能夠更好地模擬有云條件下的微波和紅外輻射[5]。

假設(shè)大氣為局地熱力、動力平衡的平行平面大氣，在不考慮散射的條件下，衛(wèi)星接收到的輻射表示為：

L(υ，θ)=(1-N)LClr(υ，θ)+NLCld(υ，θ)

(1)

式中，υ為觀測通道的中心波數(shù)，θ為衛(wèi)星觀測角度，N為云量，LClr(υ，θ)為晴空時的大氣層頂向上輻射，LCld(υ，θ)為有云時的大氣層頂向上輻射。LClr(υ，θ)包含了表面輻射和大氣輻射兩部分，定義為：

(2)

LCld(υ，θ)的定義為：

(3)

式中，τCld(υ，θ)為云頂?shù)酵鈱涌臻g的透射率，TCld為云頂溫度。對于厚而密的云層，其紅外出射率假設(shè)為均一值(ε=1)。

由(2)式和(3)式得出，要計算有云條件下衛(wèi)星接收到的輻射，關(guān)鍵是要知道透射率。

模式相應(yīng)氣壓層j到大氣層頂?shù)目偟耐ǖ劳干渎蕿?種氣體透射率的乘積：

(4)

式中，上標mix、wv、oz分別表示均勻混合氣體、水汽、臭氧3種氣體種類，而透射率和光學厚度之間可以相互轉(zhuǎn)化：

τv,j=exp(-|σv,j|)

(5)

對于給定的大氣狀態(tài)(大氣溫度、水汽、臭氧廓線)和表面參量(表面發(fā)射率、氣壓、溫度)，沿著儀器的觀測方向，可以計算出3種氣體的光學厚度：

(6)

式中，Xj，k為依賴于廓線的預(yù)報因子，m為預(yù)報因子的個數(shù)。av，j，k是分別按3組氣體預(yù)先計算好的快速透過率系數(shù)，RTTOV中有專門的系數(shù)文件來儲存這些系數(shù)。

利用快速輻射傳輸模式，依據(jù)輸入的預(yù)報因子、回歸系數(shù)和式(6)計算各模式層到大氣層頂?shù)墓鈱W厚度，再根據(jù)式(4)和式(5)計算大氣透射率，最后根據(jù)式(1)正演計算出衛(wèi)星觀測到的輻射率。若云層很厚而稠密，則在紅外譜段可以將云視作黑體，這時衛(wèi)星在紅外大氣窗測量的輻射來自云頂表面的溫度，根據(jù)普朗克黑體輻射公式就能推得云頂表面的溫度，再根據(jù)大氣廓線就可推得云頂高度。

RTTOV模式要求的初始參量包含了從地面到0.1hPa的大氣溫度和水汽廓線、表面溫度、水汽和氣壓以及一些云參數(shù)。為了滿足O3探測通道的需求，RTTOV將O3廓線也作為模式的初始輸入變量。RTTOV8以后的版本中包含了云散射模塊(RTTOV_SCATT)，相比RTTOV的晴空輻射模塊，輸入變量中增加了對大氣水凝物含量的要求，這些變量包括云水混合比、云冰混合比、雨水流量和云量。

在確定表面出射率時，針對微波波段，RTTOV提供FASTEM1～5共5種參數(shù)化方案。對于紅外輻射，RTTOV有默認的系統(tǒng)值，陸面出射率為0.98，海-冰表面出射率為0.99，而海面出射率由ISEM參數(shù)化方案計算得出。

2RTTOV的應(yīng)用研究進展

2.1RTTOV在數(shù)值天氣預(yù)報同化中的應(yīng)用直接同化衛(wèi)星資料一般采用變分法來實現(xiàn)。直接同化衛(wèi)星資料時，衛(wèi)星探測的輻射率資料將被直接引入變分同化系統(tǒng)，通過快速輻射傳輸模式模擬衛(wèi)星探測到的大氣層頂向上的輻射率，經(jīng)變分系統(tǒng)中定義目標函數(shù)表征該模擬值與相應(yīng)衛(wèi)星觀測輻射率之間的差，將衛(wèi)星探測信息轉(zhuǎn)換為大氣狀態(tài)變量，再利用極小化算法迭代計算出最接近實際大氣狀態(tài)大氣溫度、濕度等變量[2]。

利用RTTOV模式將衛(wèi)星資料同化到數(shù)值預(yù)報模式中，能顯著提高預(yù)報的準確度。Derber等[6]將TOVS晴空輻射資料直接同化到美國國家環(huán)境預(yù)報中心(NationalCentersforEnvironmentalPrediction，NCEP)同化系統(tǒng)之中，改進了模式的預(yù)報技巧。Okamoto等[7]在NCEP的全球模式中同化了特別微波成像輻射計(SpecialSensorMicrowaveImager，SSM/I)輻射資料，提高了預(yù)報結(jié)果的準確度。暴雨、龍卷和雷暴等強對流天氣系統(tǒng)，由于其空間尺度小、時間尺度短的特點，用常規(guī)資料很難觀測到，但衛(wèi)星卻可以觀測到中小尺度云系的發(fā)生、發(fā)展、成熟和消亡的全過程，因此，同化衛(wèi)星資料對中小尺度系統(tǒng)的預(yù)報具有重要意義。丁偉鈺等[8]基于全球區(qū)域同化預(yù)報系統(tǒng)(GlobalandRegionalAssimilationandPredictionEnhancedSystem，GRAPES)，根據(jù)RTTOV的特點，增加云水含量、云冰水含量和云量作為同化系統(tǒng)控制變量，同化了中分辨率成像光譜儀(Moderate-resolutionImagingSpectroradiometer，MODIS)資料，模擬了一次暴雨過程，結(jié)果表明，同化衛(wèi)星資料后可以改進初始場水汽和溫度分布，對短時降水預(yù)報有正面影響。張華等[9]基于GRAPES三維變分同化系統(tǒng)，利用輻射傳輸方程作為觀測算子直接同化了改進的微波探測裝置(AdvancedMicrowaveSoundingUnit，AMSU)資料，使預(yù)報的臺風路徑更加準確。

2.2RTTOV在衛(wèi)星亮溫模擬中的應(yīng)用除了在同化中的應(yīng)用，還可將RTTOV獨立于同化系統(tǒng)，用來模擬衛(wèi)星輻射亮溫。而云的存在對模式的模擬結(jié)果有很大的影響。隨著數(shù)值模式技術(shù)的發(fā)展，現(xiàn)在已經(jīng)能夠預(yù)報云物理量的三維結(jié)構(gòu)。不少學者針對云的影響展開了大量研究，如Thomas等[10]利用輻射傳輸模式、氣體消光模式和云模式構(gòu)造了全天氣條件下的輻射資料同化的觀測算子，由云模式來提供輻射傳輸模式需要的云物理量信息。Chevallier等[11]利用數(shù)值模式輸出的云物理量信息及RTTOV快速輻射傳輸模式成功模擬了有云條件下靜止衛(wèi)星11-μm通道亮溫。馬剛等[12]利用RTTOV5模擬1998年8月NOAA14極軌氣象衛(wèi)星上紅外探測器(High-resolutionInfraredRadiationSounder，HIRS)的輻射亮溫，與同時段實測TOVS探測資料對比，結(jié)果表明晴空模擬的水汽探測通道和長波窗區(qū)探測通道的亮溫誤差小于有云時的誤差。丁偉鈺等[13]結(jié)合天氣研究和預(yù)報模式(WeatherResearchandForecastingModel，WRF)輸出的氣象場資料，利用RTTOV快速輻射傳輸模式針對臺風“珍珠”進行了研究，模擬了HIRS/3紅外輻射亮溫，結(jié)果表明，晴空條件下，模擬的亮溫顯著偏高，且無法得到臺風的特征信息，考慮了云的影響之后，對應(yīng)低層通道(能量峰值在云頂以下)的模擬結(jié)果與觀測比較一致。馬剛等[14]綜合逐線模式GENLN2計算的精確大氣透過率，以及FY-2B紅外和水汽通道的光譜響應(yīng)函數(shù)，重新計算了RTTOV透過率系數(shù)，并將其用于FY-2B紅外和水汽通道亮溫的模擬；采用溫度、水汽和臭氧廓線進行敏感性試驗，得出2個通道對溫度均是敏感的，而對水汽的敏感度比溫度小一個量級，臭氧的變化對通道的模擬結(jié)果影響不大；對比模擬的亮溫值與實測值發(fā)現(xiàn)，除沙漠地區(qū)外，大部分地區(qū)的偏差較小。自RTTOV8以來，微波的散射模塊被開發(fā)出來并逐漸完善。蘇捷等[15]分析指出在微波亮溫的資料同化中，改用RTTOV-SCATT將大氣中水凝物的信息納入同化的計算中，可以有效地解決使用以往的方法所導致的對大氣中水汽含量的高估現(xiàn)象。

若云層很厚而稠密，則在紅外譜段可以將云視作黑體，這時衛(wèi)星在紅外大氣窗測量的輻射來自云頂表面的溫度，即模擬得到云頂亮溫?；谶@一原理，丁偉鈺[16]利用Grapes模式的模擬結(jié)果作為輸入量進行云頂亮溫模擬，計算RTTOV模式模擬得到的總云量與FY-2D反演得到的實測值之間的逐時相關(guān)系數(shù)和均方根誤差，結(jié)果表明，Grapes中尺度模式具備一定的云參數(shù)模擬能力；RTTOV在模擬的前12h模擬能力相對較高，之后逐漸減弱，云參數(shù)對FY-2衛(wèi)星云頂亮溫的模擬有很大影響。

對于超過幾百米厚而密實的云層，可以看成是黑體，模擬出來的輻射值主要來自云頂表面，由此可以直接估算云頂溫度或高度。但對于薄云或未能充滿視場的云單體，輻射值是云體與地表面發(fā)射輻射的總和。在這種情況下模擬得到的衛(wèi)星云頂亮溫存在較大誤差，由此輻射推算的云頂溫度比實際的暖，估算的云頂高度比實際的低。由于常規(guī)觀測資料經(jīng)數(shù)值預(yù)報模式得到的預(yù)報場中包含的云水信息有限，對云參數(shù)的預(yù)報能力還有待提高，導致在利用輻射傳輸模式模擬衛(wèi)星亮溫的初始階段，即初始參數(shù)的輸入上有誤差，使得模擬結(jié)果與實測值之間存在較大差異。

3結(jié)語

該研究綜述了多年來國內(nèi)外針對RTTOV的應(yīng)用研究進展，結(jié)果顯示，利用RTTOV將衛(wèi)星資料直接同化到數(shù)值預(yù)報模式中，能顯著提高數(shù)值預(yù)報的準確度，尤其是對于中小尺度的強對流天氣系統(tǒng)，無法用常規(guī)手段進行觀測，將衛(wèi)星輻射亮溫同化到數(shù)值預(yù)報模式初始場中顯得尤為重要；同時，RTTOV可以獨立于同化系統(tǒng)，能夠較好地模擬晴空條件下衛(wèi)星接收到的輻射亮溫，有云條件下的模擬值與實測值之間的誤差大于晴空條件。云的存在對輻射亮溫模擬效果有很大影響，模擬時要著重考慮云的三維物理結(jié)構(gòu)和云中水成物含量分布，從而使模擬結(jié)果更加準確。

研究RTTOV快速輻射傳輸模式主要有3個方面的意義：①因為衛(wèi)星資料直接同化需要利用快速輻射傳輸模式作為觀測算子，RTTOV模擬的結(jié)果是否準確影響了同化的效果。只有提高RTTOV模擬的準確度，才能確保提高數(shù)值預(yù)報模式的準確度；②在利用RTTOV進行衛(wèi)星輻射亮溫模擬時，要求輸入包括云和大氣溫度、濕度等由數(shù)值天氣預(yù)報得到的氣象要素，RTTOV的模擬效果在一定程度上也體現(xiàn)了數(shù)值模式的預(yù)報能力；③在實際的業(yè)務(wù)工作中，衛(wèi)星云圖反映的是過去或?qū)崟r的天氣過程，并不能預(yù)報未來的天氣形勢，而將RTTOV模擬得到的輻射亮溫進行仿云圖制作，則可以得到未來時刻的衛(wèi)星云圖，從衛(wèi)星云圖中可以直觀地判斷云體和云系的移動及發(fā)展變化，對暴雨、臺風、龍卷等強對流天氣系統(tǒng)進行預(yù)報，進而改善氣象預(yù)報和保障業(yè)務(wù)，減少災(zāi)害性天氣造成的損失。

參考文獻

[1] 王恬.氣象衛(wèi)星遙感資料同化及其在數(shù)值預(yù)報中的應(yīng)用[C]//創(chuàng)新驅(qū)動發(fā)展 提高氣象災(zāi)害防御能力：S4衛(wèi)星資料在中小尺度天氣分析中的應(yīng)用.中國氣象學會，2013.

[2] 董佩明，薛繼善，黃兵，等.數(shù)值天氣預(yù)報中衛(wèi)星資料同化應(yīng)用現(xiàn)狀和發(fā)展[J].氣象科技，2008，36(1)：1-7.

[3] 李俊，方宗義.氣象衛(wèi)星的發(fā)展——機遇與挑戰(zhàn)[J].氣象，2012，38(2)：129-146.

[4] SAUNDERS R，BRUNEL P.RTTOV-7 technical report[R].NWP-SAF document，2002.

[5] HOCKING J，SAUNDERS R，RAYE P，et al.RTTOV development status[C]//18th international TOVS study conference.Toulouse，F(xiàn)rance，2012.

[6] DERBER J C，WU W.The use of TOVS cloud-cleared radiances in the NCEP SSI analysis system[J].Monthly weather review，1998，126：2287-2299.

[7] OKAMOTO K，DERBER J C.Assimilation of SSM/I radiances in the NCEP global data assimilation system[J].Monthly weather review，2006，134：2612-2631.

[8] 丁偉鈺，萬齊林，黃燕燕，等.有云環(huán)境下 MODIS 亮溫資料的變分同化Ⅱ——對暴雨預(yù)報的影響[J].熱帶氣象學報，2010，26(1)：22-30.

[9] 張華，丑紀范，邱崇踐.西北太平洋威馬遜臺風結(jié)構(gòu)的衛(wèi)星觀測同化分析[J].科學通報，2004，49(5)：493-498.

[10] GREENWALD T J，HERTENSTEIN R，VUKICEVIC T.An all-weather observational operator for radiance data assimilation with mesoscale forecast models[J].Monthly weather review，2002，130：1882-1897.

[11] CHEVALLIER F，KELLY G.Model clouds as seen from space：Comparison with geostationary imagery in the 11-μm window channel[J].Monthly weather review，2002，130：712-722.

[12] 馬剛，方宗義，張鳳英.云參數(shù)對RTTOV5模式模擬誤差的影響分析[J].應(yīng)用氣象學報，2001，12(4)：385-392.

[13] 丁偉鈺，萬齊林.“珍珠”臺風衛(wèi)星紅外通道亮溫的數(shù)值模擬[J].大氣科學，2008，32(3)：572-580.

[14] 馬剛，邱崇踐，黎光清，等.利用RTTOV7快速輻射傳輸模式模擬風云二號紅外和水汽成像通道輻射率的研究[J].紅外與毫米波學報，2005，24(5)：37-40.

[15] 蘇捷，朱克云，張杰.RTTOV散射模塊在FY-3 MWHS資料同化中的應(yīng)用[J].氣象，2013，39(11)：1461-1472.

[16] 丁偉鈺.基于 Grapes 中尺度模式的 FY-2 衛(wèi)星 TBB 亮溫數(shù)值模擬[C]//第28屆中國氣象學會年會——S2風云衛(wèi)星定量應(yīng)用與數(shù)值.中國氣象學會，2011.

收稿日期2015-12-28

作者簡介劉樂(1990- )，女，廣西桂林人，碩士研究生，研究方向：衛(wèi)星資料正演模擬。

中圖分類號S 164

文獻標識碼A

文章編號0517-6611(2016)03-230-03