劉曉燕，陳莉，邢建勇，高元勇，崔利民

（1.國家海洋環(huán)境預(yù)報中心，北京 100081；2.國家衛(wèi)星海洋應(yīng)用中心，北京 100081）

1 引言

我國東部毗鄰西北太平洋，因此我國是受臺風(fēng)影響最為嚴(yán)重的國家之一[1]。西北太平洋平均每年可生成28.2個臺風(fēng)，在我國登陸的約有7個[2-3]，其伴隨的狂風(fēng)和暴雨等災(zāi)害現(xiàn)象，給所影響的省市和地區(qū)帶來了嚴(yán)重的損失。提高臺風(fēng)預(yù)報水平，及時做好防范措施，可有效降低災(zāi)害天氣帶來的損失，因此臺風(fēng)的準(zhǔn)確預(yù)報對國家和社會具有重要意義[4-5]。

近些年，隨著衛(wèi)星和雷達(dá)等遙感技術(shù)的發(fā)展，非常規(guī)觀測資料迅速增多，有效彌補(bǔ)了常規(guī)觀測資料的不足。這些非常規(guī)觀測資料，不僅可以對臺風(fēng)進(jìn)行及時準(zhǔn)確的監(jiān)測，而且可以用于數(shù)值模式中，通過同化技術(shù)改善模式初始場，提高臺風(fēng)預(yù)報效果。黃嘉宏等[6]利用四維變分（Four-Dimensional VARiational，4DVAR）同化技術(shù)，對南海臺風(fēng)“Vongfong”進(jìn)行了模擬試驗，證明了QuikSCAT 衛(wèi)星海面風(fēng)場資料可以使模擬的臺風(fēng)路徑與實況更為接近；王業(yè)桂等[5]探究了同化不同衛(wèi)星的同一種微波遙感資料對臺風(fēng)路徑模擬效果的影響，證明了同化衛(wèi)星微波遙感資料可以有效改善臺風(fēng)路徑模擬，但是同化不同資料的效果不一；希爽等[7]針對2008 年影響我國大陸的8個典型臺風(fēng)開展了批量衛(wèi)星資料同化和數(shù)值預(yù)報試驗，證明了衛(wèi)星資料對臺風(fēng)預(yù)報有改進(jìn)作用。

海洋2B（HY-2B）衛(wèi)星和中法海洋衛(wèi)星（China-France Oceanography SATellite，CFOSAT）均為2018年底發(fā)射，截至2020 年底，已經(jīng)有一年多的數(shù)據(jù)可供使用，也形成了一些相關(guān)研究。陳克海等[8]利用歐洲中期天氣預(yù)報中心（European Centre for Medium-Range Weather Forecasts，ECMWF）再分析風(fēng)場數(shù)據(jù)、熱帶大氣海洋觀測計劃和美國國家數(shù)據(jù)浮標(biāo)中心（National Data Buoy Center，NDBC）浮標(biāo)獲取的海面風(fēng)矢量實測數(shù)據(jù)對HY-2B 散射計海面風(fēng)場數(shù)據(jù)產(chǎn)品的質(zhì)量進(jìn)行了統(tǒng)計分析。結(jié)果表明，不論是在開闊海域還是近海海域，HY-2B 散射計海面風(fēng)資料都能較好地滿足業(yè)務(wù)應(yīng)用的精度要求，即風(fēng)速優(yōu)于2 m/s，風(fēng)向優(yōu)于20°。國家衛(wèi)星海洋應(yīng)用中心的中法衛(wèi)星在軌測試評審報告指出，CFOSAT對ECMWF 再分析風(fēng)場數(shù)據(jù)的風(fēng)速精度在1.26～1.48 m/s 之間，風(fēng)向精度在11°～15°之間，優(yōu)于風(fēng)速測量誤差2 m/s 和風(fēng)向誤差20°的指標(biāo)要求[9-10]。Xu等[11]利用CFOSAT同步觀測到了由臺風(fēng)引起的風(fēng)場和海浪場，實現(xiàn)了臺風(fēng)的有效監(jiān)測。到目前為止，這兩種衛(wèi)星散射計數(shù)據(jù)用于同化的工作還比較少。本文基于新一代高分辨率中尺度預(yù)報模式（Weather Research and Forecasting model，WRF）及其三維變分（Three-Dimensional Variational，3DVAR）同化系統(tǒng)，開展這兩種海洋衛(wèi)星散射計海面風(fēng)資料在同化中的應(yīng)用，并考察同化不同衛(wèi)星散射計海面風(fēng)資料對臺風(fēng)模擬的影響。

2 臺風(fēng)個例概述

2019年第18號熱帶氣旋“米娜”（Mitag）于2019年9月28日00時（世界時，下同）在臺灣省臺北市東南方約1540 km 的西北太平洋上（132.1°E，15.4°N）生成，中心附近最大風(fēng)力8 級（18 m/s），中心最低氣壓為998 hPa；隨后向西北方向移動，強(qiáng)度逐漸加強(qiáng)，28 日21 時加強(qiáng)為強(qiáng)熱帶風(fēng)暴級，29 日09 時加強(qiáng)為臺風(fēng)級；10 月1 日00 時開始減弱，并于1 日12 時30分在浙江舟山普陀沿海登陸，登陸時中心附近最大風(fēng)力為11級（30 m/s）），中心最低氣壓980 hPa；之后向東北方向移動，強(qiáng)度逐漸減弱，并于2 日12 時10分前后在韓國全羅南道沿海再次登陸，登陸時為熱帶風(fēng)暴級，最大風(fēng)力9 級（23 m/s），中心最低氣壓為988 hPa；之后繼續(xù)減弱，3 日06 時中央氣象臺對其停止編號[12-13]。

3 資料介紹、方法介紹、模式設(shè)置及試驗設(shè)計

3.1 資料介紹

本文所使用的數(shù)據(jù)主要有美國全球預(yù)報系統(tǒng)（Global Forecast System，GFS）的模式輸出產(chǎn)品、我國自主研制的HY-2B 散射計海面風(fēng)數(shù)據(jù)以及與法國聯(lián)合研制的CFOSAT散射計海面風(fēng)數(shù)據(jù)。GFS數(shù)據(jù)每6 h 輸出一次，取其第12～84 h 共計72 h 的預(yù)報結(jié)果為本文試驗提供初值（背景場）和邊界條件，該數(shù)據(jù)空間分辨率為0.5°×0.5°。HY-2B 是我國第二顆海洋動力環(huán)境衛(wèi)星，于2018年10月25日發(fā)射，其上搭載的散射計主要功能為測量海面風(fēng)矢量場，工作頻率為13.256 GHz，極化方式包含HH和VV兩種，其中H 極化的刈幅優(yōu)于1350 km，V 極化的刈幅優(yōu)于1700 km，分辨率為25 km×25 km，風(fēng)速測量范圍為2～24 m/s，風(fēng)速精度為±2 m/s 或10%，風(fēng)向精度 為±20°[14]。CFOSAT 于2018 年10 月29 日 發(fā)射，是一顆具有海風(fēng)海浪聯(lián)合探測能力的衛(wèi)星，其上搭載的散射計是國際上首次采用旋轉(zhuǎn)扇形波束掃描方式測量海洋風(fēng)場的微波散射計，刈幅大于1000 km，工作頻率為13.256 GHz，空間分辨率有25 km×25 km 和12.5 km×12.5 km 兩種，本文使用前者，風(fēng)速的有效測量范圍為4～24 m/s，精度為±2 m/s或10%，風(fēng)向精度為±20°[15]。此外，臺風(fēng)路徑和強(qiáng)度觀測數(shù)據(jù)來自中國氣象局熱帶氣旋最佳路徑數(shù)據(jù)集（網(wǎng)址：https：//tcdata.typhoon.org.cn）[16]。

3.2 方法介紹

本文采用WRFV3.7.1 及其3DVAR 同化系統(tǒng)。3DVAR 是變分法的一種方法，它是通過迭代函數(shù)，產(chǎn)生一個分析時刻的大氣真實狀態(tài)的最優(yōu)估計，通過目標(biāo)泛函表示為：

式中，x為大氣狀態(tài)變量；xb為背景場資料；y0為觀測資料；H為觀測算子，它將分析變量投影到觀測空間；B為背景場誤差協(xié)方差矩陣；O為觀測誤差協(xié)方差矩陣。三維變分問題就是求解上式極小值的問題。

文中使用的3DVAR同化方法依賴于WRF模式的變分同化系統(tǒng)，該同化系統(tǒng)根據(jù)設(shè)置的不同，可以識別little_r和bufr兩種數(shù)據(jù)格式。文中采用的是little_r 格式，即一種該同化系統(tǒng)所識別的有著固定格式的文本格式，將觀測資料輸出為此格式，進(jìn)入同化系統(tǒng)，該同化系統(tǒng)在進(jìn)行質(zhì)量控制后（剔除超過閾值的數(shù)據(jù)），使用上述公式計算，便可得到優(yōu)化后的初始場。

3.3 模式設(shè)置及試驗設(shè)計

本文模擬區(qū)域的模式水平分辨率為25 km，垂直為不均勻35 層，以（128°E，23°N）為中心，共400×300 個格點。試驗采用的模式物理過程方案設(shè)置見表1。模擬時間為2019 年9 月29 日12 時—10 月2 日12 時。同 化分析時刻為2019 年9 月29 日12時，同化時間窗口為2019年9月29日09—15時。

表1 物理過程參數(shù)化方案設(shè)置Tab.1 Parameterization scheme setting of physical processes

文中共設(shè)計了4 組試驗，試驗HY2B、試驗CFO和試驗HY2B+CFO 分別為同化HY-2B 衛(wèi)星散射計海面風(fēng)數(shù)據(jù)、CFO 衛(wèi)星散射計海面風(fēng)數(shù)據(jù)和同化以上兩種數(shù)據(jù)，試驗CTRL為控制試驗（不同化任何資料），用于評估其他3組同化試驗的效果。具體設(shè)計見表2。

表2 試驗設(shè)計Tab.2 Design of numerical experiments

4 試驗結(jié)果分析

4.1 同化效果分析

圖1 為3 組同化試驗中吸收的觀測資料的分布情況。圖中可以看到兩顆衛(wèi)星的散射計海面風(fēng)數(shù)據(jù)的部分區(qū)域有重合，且都完整地捕捉到了臺風(fēng)位置信息，不同的是兩顆衛(wèi)星散射計捕捉臺風(fēng)信息的時間存在差異，HY-2B 衛(wèi)星散射計捕捉到臺風(fēng)的時間在09—10時之間，而CFO衛(wèi)星散射計捕捉到臺風(fēng)的時間在10—11時之間，更接近于同化時間。

圖1 3組同化試驗的觀測資料分布情況Fig.1 Observation data distribution of three groups of assimilation experiments

圖2 是3 組同化試驗使用不同衛(wèi)星觀測資料進(jìn)行同化后，觀測資料的U/V 分量分別與背景場和分析場進(jìn)行比較的結(jié)果，表3 更直觀地展示了圖2 中同化前后的U 分量的分析結(jié)果。由圖2 和表3 均可以看出，不論是使用HY-2B 衛(wèi)星散射計資料、CFO衛(wèi)星散射計資料亦或二者和，其效果都是一致的，相對于背景場，分析場的平均誤差（mean）、均方根值（rms）和標(biāo)準(zhǔn)差（std）均更接近觀測場，即OMA

圖2 3組同化試驗的OMB和OMA（OBS為觀測資料，BAK為背景場資料，ANA為分析場資料，單位：m/s）Fig.2 OMB and OMA of three groups of assimilation experiments（OBS is observation data，BAK is background field data，ANA is analysis field data，unit：m/s）

表3 同化前后U分量的背景場與分析場對比結(jié)果Tab.3 Comparison of background field and analysis field of U component

此外，我們還可以看到HY2B 試驗中觀測資料與背景場的差異小于其他兩組試驗的差異，說明HY-2B 衛(wèi)星散射計海面風(fēng)資料與背景場資料接近。從初始增量場圖（略）也可以看出，HY2B 試驗的增量小于其他兩組試驗增量，這與圖2 反映出的特征吻合。

4.2 臺風(fēng)路徑72 h模擬結(jié)果分析

圖3 是4 組試驗的臺風(fēng)路徑72 h 模擬結(jié)果與臺風(fēng)實況的比較。由圖可以看出，4 組試驗結(jié)果整體趨勢比較接近，與實況相比，CTRL 和HY2B 試驗?zāi)M的路徑偏西，其他兩組試驗?zāi)M的路徑偏東，且4組試驗?zāi)M的臺風(fēng)移動速度都要慢于實況。臺風(fēng)移動既受臺風(fēng)內(nèi)部因子的作用又受環(huán)境條件的影響，環(huán)境條件主要是指高空引導(dǎo)氣流，但4組試驗的500 hPa 位勢高度場（圖略）未能將路徑不同的原因體現(xiàn)出來，說明同化過程對該高度場的影響不太明顯，同時還說明造成4 組試驗臺風(fēng)移動路徑差異的原因可能是臺風(fēng)內(nèi)部因子的作用。

圖3 臺風(fēng)“米娜”的模擬路徑和實況路徑Fig.3 Simulated tracks and observed track of typhoon"Mitag"

圖4 是4 組試驗在72 h 內(nèi)逐6 h 的路徑誤差變化。由圖可以看出，HY2B 試驗的路徑誤差趨勢跟CTRL 試驗結(jié)果比較相似，均為在模擬后期誤差明顯增大；HY2B 試驗的結(jié)果除個別時刻與CTRL 試驗結(jié)果比較一致外，其他時刻的效果是優(yōu)于CTRL 試驗的；CFO 試驗結(jié)果明顯優(yōu)于控制試驗CTRL，在整個模擬時段內(nèi)，誤差均小于CTRL；最后一組同化試驗HY2B+CFO 結(jié)果的模擬誤差隨時間變化不大，在預(yù)報前36 h 略微優(yōu)于CTRL，在42 h 時誤差大于CTRL，之后隨著CTRL 的誤差快速增大將二者之間的差距拉開，HY2B+CFO 試驗結(jié)果在后期明顯優(yōu)于CTRL。從圖中左上角的統(tǒng)計結(jié)果也可以得知，3 組同化試驗對于臺風(fēng)路徑的模擬均起到了較好的改善效果。

圖4 4組試驗的臺風(fēng)路徑模擬誤差逐6 h變化Fig.4 Temporal variations of errors in typhoon track simulations at 6-hour interval for 4 groups of experiments

從圖4 中我們可以明顯地看到，在模擬時段的前54 h 內(nèi)，路徑誤差均控制在150 km 以內(nèi)，54 h 之后，CTRL 和HY2B 兩組試驗的誤差開始明顯增大，72 h 時誤差達(dá)到300 km 左右，而另外兩組試驗的誤差隨模擬時效增大不明顯。

4.3 臺風(fēng)強(qiáng)度72 h模擬結(jié)果分析

4.3.1 臺風(fēng)中心海平面最低氣壓分析

圖5為臺風(fēng)中心海平面最低氣壓模擬結(jié)果和實況對比情況。由圖可以看出，4 組試驗結(jié)果走勢比較一致，均為氣壓先減小后增大，與實況的變化趨勢也基本一致，但是4 組試驗結(jié)果在模擬前60 h 均大于實況，即模擬的臺風(fēng)強(qiáng)度均偏弱，在模擬時段的后12 h 里4 組試驗出現(xiàn)了分支，CFO 和HY2B+CFO兩組試驗?zāi)M的臺風(fēng)強(qiáng)度偏強(qiáng)。

圖5 臺風(fēng)中心海平面最低氣壓模擬結(jié)果和實況對比Fig.5 Comparison of simulation results of lowest sea-level pressure in typhoon center and CMA best track data

4.3.2 臺風(fēng)中心附近最大風(fēng)速分析

圖6為臺風(fēng)中心附近最大風(fēng)速模擬結(jié)果和實況對比情況。由圖可以看出，4 組試驗結(jié)果走勢比較一致（起報時刻除外），均為風(fēng)速先增大再減小，與實況的變化趨勢也基本一致，只是4 組試驗的臺風(fēng)強(qiáng)度模擬結(jié)果均比實況弱。

圖6 臺風(fēng)中心附近最大風(fēng)速模擬結(jié)果和實況對比Fig.6 Comparison of simulation results of maximum wind speed in typhoon center and CMA best track data

不論是從臺風(fēng)中心最低氣壓還是中心附近最大風(fēng)速情況分析，4 組試驗?zāi)M的臺風(fēng)強(qiáng)度均偏弱。分析其原因，可能因為臺風(fēng)強(qiáng)度主要和熱力結(jié)構(gòu)直接相關(guān)，本文的觀測資料是衛(wèi)星散射計海面風(fēng)資料，更多的是直接改進(jìn)臺風(fēng)風(fēng)場等動力結(jié)構(gòu)，雖然利用背景誤差協(xié)方差，風(fēng)場的調(diào)整可以傳遞到其他與熱力結(jié)構(gòu)有關(guān)的控制變量，但影響有限，因此對臺風(fēng)的強(qiáng)度改進(jìn)不明顯。

5 總結(jié)

本文進(jìn)行了HY-2B 和CFO 兩顆衛(wèi)星散射計海面風(fēng)資料在同化中的應(yīng)用?；?019 年第18 號臺風(fēng)“米娜”，利用WRF 模式及其3DVAR 同化系統(tǒng)對不同的海洋衛(wèi)星散射計海面風(fēng)資料開展同化試驗，考察了同化不同散射計海面風(fēng)對臺風(fēng)模擬效果的影響。結(jié)論如下：

（1）對于臺風(fēng)路徑，4 組試驗的模擬和實況較為一致，且3 組同化試驗的結(jié)果均優(yōu)于CTRL 試驗，CFO 試驗效果最優(yōu)，其次是HY2B+CFO，HY2B 試驗效果相對弱一些，可能是因為CFOSAT 散射計海面風(fēng)資料的時間更接近模擬時刻。

（2）對于臺風(fēng)強(qiáng)度，4 組試驗的模擬情況與實況的變化趨勢較為一致，但是模擬結(jié)果比實況偏弱，3組同化試驗未在臺風(fēng)強(qiáng)度模擬中體現(xiàn)出同化的優(yōu)勢，只有在海平面最低氣壓模擬中，CFO 和HY2B+CFO兩組試驗的結(jié)果略微優(yōu)于CTRL。

基于上述結(jié)論可知，同化不同海洋衛(wèi)星散射計海面風(fēng)資料對于此次臺風(fēng)路徑的模擬都具有不同效果的改善，但是對于臺風(fēng)強(qiáng)度的模擬效果比較微弱。由于本文結(jié)論只是基于個例試驗的衛(wèi)星散射計海面風(fēng)場數(shù)據(jù)的一次應(yīng)用，具有一定的局限性，其結(jié)論并不能直接代表這兩種衛(wèi)星散射計海面風(fēng)場數(shù)據(jù)的真實效果，因而下一步將考慮進(jìn)行一年或更長時間內(nèi)的臺風(fēng)模擬試驗，得出更具代表性和說服力的結(jié)論。

致謝：所有CFOSAT 數(shù)據(jù)由中國國家航天局（CNSA）和法國國家空間研究中心（CNES）提供。CFOSAT 和HY-2B 數(shù)據(jù)來源于https://osdds.nsoas.org.cn.