張釗揚(yáng) 鐘劍 周炯 王國云

摘要利用COSMIC大氣主要?dú)庀笠乩€資料，經(jīng)過數(shù)據(jù)標(biāo)準(zhǔn)化處理、質(zhì)量控制、數(shù)據(jù)可視化，使其成為預(yù)報(bào)員能直接使用的可視化產(chǎn)品，以期改善包括廣大海洋地區(qū)缺乏實(shí)況氣象觀測資料的現(xiàn)狀。通過與再分析資料的對比發(fā)現(xiàn)，該資料雖然存在隨機(jī)性和不確定性等缺點(diǎn)，但其覆蓋面廣、精度高，能夠較好地反映大范圍的環(huán)流形勢和一定區(qū)域內(nèi)的天氣系統(tǒng)發(fā)展演變，尤其在測站稀少的大洋地區(qū)，在可預(yù)見的將來沒有大范圍氣象站點(diǎn)布設(shè)的情況下，COSMIC資料能夠起到部分替代實(shí)況觀測的作用。

關(guān)鍵詞COSMIC；標(biāo)準(zhǔn)化；質(zhì)量控制；Cressman客觀分析；遠(yuǎn)洋氣象保障

中圖分類號U666.135文獻(xiàn)標(biāo)識碼A文章編號0517-6611（2017）35-0190-03

AbstractBased on the COSMIC profile data of the main meteorological elements，through standardization，quality control and data visualization，make it a visualized product which can be directly used by the forecasters，in order to improve the current situation of the lack of the actual meteorological observational data in the vast ocean area.Comparing with reanalysis data found that despite the shortcomings of randomness and uncertainty，its wide coverage and high precision can well reflect the largescale circulation situation and the development of the weather system of a certain area.Especially in the ocean，in the foreseeable future without a wide range of meteorological site layout，COSMIC data can play a part of the role of alternative reality observation.

Key wordsCOSMIC；Standardization；Quality control；Cressman objective analysis；Ocean meteorological support

隨著中國經(jīng)濟(jì)的發(fā)展和貿(mào)易全球化的持續(xù)推進(jìn)，遠(yuǎn)洋航運(yùn)業(yè)務(wù)也日趨繁榮。面對廣闊的大洋和陌生的海域，氣象保障的重要性也日趨凸顯。由于大洋地區(qū)觀測資料稀少，目前遠(yuǎn)洋氣象保障只能依靠數(shù)值預(yù)報(bào)產(chǎn)品。雖然數(shù)值預(yù)報(bào)的準(zhǔn)確性越來越高，但其也經(jīng)常存在錯(cuò)報(bào)漏報(bào)的情況，尤其是在缺乏實(shí)況檢驗(yàn)的大洋地區(qū)。

20世紀(jì)90年代開始發(fā)展起來的空基GPS氣象學(xué)就是利用GPS掩星資料反演大氣主要?dú)庀笠?。其原理主要是根?jù)低軌衛(wèi)星LEO（Low Earth Orbit）在掩星事件過程中接收到的GPS衛(wèi)星發(fā)出的無線電信號的相位和振幅，利用GPS衛(wèi)星和LEO衛(wèi)星的位置和速度數(shù)據(jù)，計(jì)算出GPS信號經(jīng)過地球大氣時(shí)發(fā)生的彎曲角，再對局部球?qū)ΨQ的假設(shè)下利用Abel積分求得折射率的廓線；最后利用折射率與大氣溫度、濕度和氣壓的關(guān)系公式，結(jié)合一維變分等方法反演出大氣的溫、壓、濕等要素[1-2]。目前，GPS掩星反演資料主要應(yīng)用于資料同化領(lǐng)域，用于改善數(shù)值模式初始場，直接應(yīng)用于業(yè)務(wù)的較少。筆者選取較為成熟的GPS掩星反演大氣主要?dú)庀笠財(cái)?shù)據(jù)——COSMIC數(shù)據(jù)，經(jīng)過標(biāo)準(zhǔn)化處理、質(zhì)量控制及可視化，使其成為預(yù)報(bào)員能直接使用的可視化產(chǎn)品，并考察其在氣象預(yù)報(bào)保障中的可用性。

1COSMIC資料介紹

中國臺灣地區(qū)和美國UCAR、JPL聯(lián)合的氣象、氣候和電離層星座觀測系統(tǒng)COSMIC（Constellation Observing System for Meteorology，Ionosphere and Climate）計(jì)劃由6顆LEO衛(wèi)星組成，其主要目的有：改善全球和局地的天氣預(yù)報(bào)；改善全球空間天氣的檢測和預(yù)報(bào)；為氣候和全球變化研究提供數(shù)據(jù)；提高對地球重力場的認(rèn)識[2-4]。

從2006年開始，COSMIC系統(tǒng)每天可提供全球2 500個(gè)掩星點(diǎn)的大氣和電離層的觀測資料和反演的氣象要素廓線，資料基本達(dá)到實(shí)時(shí)要求，最新的觀測資料2 h內(nèi)就可以通過網(wǎng)上獲得。

COSMIC系統(tǒng)發(fā)布的中性大氣要素的廓線資料最常用的有3種類型：atmPrf資料為干模式下反演的大氣要素廓線，不包括水汽信息；wetPrf資料是考慮了大氣中水汽含量的濕模式下，利用AVN等預(yù)報(bào)模式的一維同化反演的大氣主要要素廓線，主要包括不同高度上的溫度、氣壓、水汽壓和折射率；bfrPrf資料綜合了atmPrf和wetPrf的結(jié)果，同時(shí)降低垂直分辨率，主要用于模式的資料同化[2]。在此選取的是wetPrf資料作為研究對象，其探測高度范圍為0～40 km，垂直分辨率為100 m，分層數(shù)為40層，為表述方便，下文中簡稱為COSMIC資料。

2COSMIC資料的處理

COSMIC資料可從網(wǎng)頁直接查看或下載，也可在UNIX或LINUX環(huán)境下通過cshell語言編程下載獲得所需資料文件。需要注意的是，COSMIC反演的廓線數(shù)據(jù)采用的是NetCDF數(shù)據(jù)格式。因此，在數(shù)據(jù)使用之前需要對其進(jìn)行解碼。該研究主要使用FORTRAN程序讀取NetCDF數(shù)據(jù)并將其轉(zhuǎn)換為方便使用的文本數(shù)據(jù)格式。

對解碼后的數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理、質(zhì)量控制、圖形化顯示后，得到預(yù)報(bào)員可以直接使用的天氣圖，具體流程如圖1所示。

2.1數(shù)據(jù)的標(biāo)準(zhǔn)化

COSMIC每一條氣象廓線數(shù)據(jù)代表的是一個(gè)掩星點(diǎn)上空的氣象情況，由于掩星事件發(fā)生的隨機(jī)性和不確定性，因此COSMIC資料的時(shí)空分布也是沒有規(guī)律的。如圖2所示，每天2 000多個(gè)資料點(diǎn)代表互相獨(dú)立的掩星事件，每個(gè)資料垂直層達(dá)400層，對于一般業(yè)務(wù)使用來說，過高的垂直分辨率不僅沒有太大價(jià)值，而且還給計(jì)算和圖形化顯示帶來不便。

針對COSMIC資料的上述特征，按照常規(guī)探空資料的規(guī)范對其進(jìn)行了標(biāo)準(zhǔn)化處理。為了更加有利于預(yù)報(bào)員對天氣形勢轉(zhuǎn)變的掌握，標(biāo)準(zhǔn)化過程中觀測時(shí)次（6 h間隔）和垂直層（11 層）進(jìn)行了加密，即世界時(shí)00Z、06Z、12Z、18Z的 925、900、850、800、700、600、500、400、300、200、100 hPa 氣象要素情況。

垂直層的標(biāo)準(zhǔn)化相對簡單，就是將所需要的垂直層的資料從原始資料中直接提取出來。時(shí)間的標(biāo)準(zhǔn)化過程中，引入一個(gè)大小為±3 h的時(shí)間窗口，將這個(gè)時(shí)間窗口內(nèi)的所有COSMIC資料作為標(biāo)準(zhǔn)時(shí)刻的數(shù)據(jù)，即將T-3 h→T+3 h的數(shù)據(jù)都作為T時(shí)刻數(shù)據(jù)。以世界時(shí)12：00為標(biāo)準(zhǔn)時(shí)為例，09：00—15：00窗口期內(nèi)所有資料都作為標(biāo)準(zhǔn)時(shí)刻資料（圖3）。

2.2質(zhì)量控制

國內(nèi)外專家對COSMIC資料的可用性進(jìn)行了深入研究[5-8]，發(fā)現(xiàn)其總體精度較高，其中溫度的精度最高，整層均方差不超過2 K；500 hPa以下中低層位勢高度具有較高精度，均方差小于3 dagpm，500 hPa以上均方差迅速增大到5 dagpm以上。因此，COSMIC資料總體上是可以滿足一般業(yè)務(wù)需要。但COSMIC資料的精度受大氣條件、儀器誤差和資料分辨率等因素影響，資料中存在少量噪點(diǎn)。因此，該研究用設(shè)定閾值的方法對資料進(jìn)行了噪點(diǎn)剔除處理，判定為噪點(diǎn)的數(shù)據(jù)標(biāo)記為缺測，例如500 hPa高度場的閾值就設(shè)為492和592 dagpm，將不在這個(gè)范圍內(nèi)的數(shù)據(jù)點(diǎn)標(biāo)記為缺測。

2.3數(shù)據(jù)可視化準(zhǔn)備

可視化是氣象資料業(yè)務(wù)應(yīng)用的重要一環(huán)，COSMIC資料作為一種較新的氣象數(shù)據(jù)，目前直接將其圖形化顯示，應(yīng)用于業(yè)務(wù)還比較少。由于掩星事件的隨機(jī)性和不確定性，COSMIC資料的分布也具有相應(yīng)的特點(diǎn)。同時(shí)，目前COSMIC的6顆LEO衛(wèi)星都是極軌衛(wèi)星，因此中高緯度地區(qū)COSMIC掩星觀測點(diǎn)明顯多于低緯度地區(qū)。相對于陸上稠密的常規(guī)觀測網(wǎng)來說，COSMIC資料沒有固定的站點(diǎn)，呈稀疏的不均勻分布，這給其業(yè)務(wù)應(yīng)用帶來了困難。

針對以上問題，該研究參照常規(guī)站點(diǎn)資料格式，提取溫、壓、濕、風(fēng)等氣象要素并重新編碼為二進(jìn)制站點(diǎn)數(shù)據(jù)。為了克服資料稀疏和不均勻等問題，在此采用Cressman客觀分析方法將標(biāo)準(zhǔn)化后的資料插值為格點(diǎn)數(shù)據(jù)。

Cressman客觀分析通過逐步訂正將離散點(diǎn)內(nèi)插到規(guī)則格點(diǎn)，實(shí)現(xiàn)最優(yōu)化插值，被廣泛應(yīng)用于氣象領(lǐng)域。該方法基本原理是先給定初猜測場，然后用實(shí)際觀測場逐步修正初猜場，直到訂正后的值逼近觀測值為止。用公式表示為

a′=a0+Δaij

其中：

Δaij=Kk=1（W2ijkΔak）Kk=1W2ijk

式中，ak為任意氣象要素；a0是變量a在格點(diǎn)（i，j）上的初猜值；a′是變量ak在格點(diǎn)（i，j）上的訂正值；Δak是觀測點(diǎn)k上的觀測值和初猜值之差；Wijk是權(quán)重因子，為0～1.0；K是影響半徑R內(nèi)的臺站數(shù)。

Cressman客觀分析方法的關(guān)鍵是確定權(quán)重函數(shù)Wijk，它的形式為

Wijk=R2-d2ijkR2+d2ijkdijk0dijk≥R

影響半徑R一般取為1、2、4、7和10。dijk為格點(diǎn)（i，j）到觀測點(diǎn)k的距離。

3可用性檢驗(yàn)

以美國NECP （National Centers for Environmental Prediction）中心GFS （The Global Forecast System） 0.5°×0.5°分析場資料作參照，通過大氣環(huán)流形勢和區(qū)域天氣形勢對比，檢驗(yàn)COSMIC資料在業(yè)務(wù)應(yīng)用中的可用性。

3.1環(huán)流形勢檢驗(yàn)

大氣環(huán)流一般是大范圍大氣基本運(yùn)動狀態(tài)，是各種天氣系統(tǒng)發(fā)生、發(fā)展和移動的背景。了解大氣環(huán)流形勢及其變化，對于制作短、中期乃至長期天氣預(yù)報(bào)有重要指導(dǎo)作用。

安徽農(nóng)業(yè)科學(xué)2017年

從GFS資料北半球500 hPa環(huán)流形勢（圖4）可以看出，北半球秋季冷暖氣團(tuán)勢力相當(dāng)，天氣系統(tǒng)較為復(fù)雜，西北太平洋副熱帶高壓雖勢力減弱，但仍控制我國南方地區(qū)；亞歐大陸北部呈兩槽一脊型，槽線位置分別位于烏拉爾山附近和日本列島北部；菲律賓群島東部的低緯度地區(qū)有熱帶低壓系統(tǒng)發(fā)展。COSMIC資料和GFS資料的北半球環(huán)流形勢場基本一致，能較為準(zhǔn)確地體現(xiàn)天氣尺度以上系統(tǒng)，例如北半球極渦多中心分布和亞歐大陸兩槽一脊的環(huán)流形勢都有很好的刻畫。但是，相對于北半球稠密的觀測網(wǎng)來說，COSMIC資料的空間分辨率偏低，尺度偏小的系統(tǒng)往往會被濾掉，有些天氣尺度以上的系統(tǒng)在強(qiáng)度的體現(xiàn)上還有差異。 這在天氣圖上的表現(xiàn)就是COSMIC資料形勢場為等值線較為平滑的大槽大脊形態(tài)，而GFS資料形勢場呈不同尺度系統(tǒng)疊加的多中心分布。不同于COSMIC資料形勢場中高緯大槽大脊，副熱帶高壓呈不規(guī)則帶狀分布的環(huán)流形勢，GFS資料中高緯低槽中有發(fā)展比較成熟的氣旋系統(tǒng)，副熱帶高壓在大陸和海上有多個(gè)強(qiáng)度不等的中心。

3.2區(qū)域天氣形勢檢驗(yàn)

西南太平洋洋面開闊，氣象海況惡劣，溫帶氣旋年平均15.0個(gè)，主要集中在30°S附件及以南海域；熱帶氣旋年平均發(fā)生數(shù)8.8個(gè)，占全球熱帶氣旋總數(shù)的11%[9]。它們伴隨有大風(fēng)、強(qiáng)降水及低能見度等惡劣天氣，對船舶航行和經(jīng)濟(jì)生產(chǎn)活動常常造成嚴(yán)重威脅。選取西南太平洋為關(guān)注區(qū)域，通過對比COSMIC資料和GFS再分析資料500 hPa形勢場，檢驗(yàn)COSMIC資料在區(qū)域天氣形勢分析中的應(yīng)用效果。從南太平洋一次西風(fēng)帶槽脊系統(tǒng)東移過程GFS的形勢場（圖5a）可以看出，該區(qū)域中緯度基本是兩槽一脊形態(tài)，槽線分別位于155°E、165°W，脊線位于174°E附近。如果以540 dagpm為參考高度線，兩槽線與540 dagpm 交點(diǎn)在40°S，脊線與其交點(diǎn)在52°S。從COSMIC資料

500 hPa形勢場（圖5b）來看，槽脊線位置與GFS資料基本一致，同樣以540 dagpm為參考，新西蘭西側(cè)高空槽較弱，與540 dagpm高度線交點(diǎn)在45°S左右，新西蘭以東高空槽較深，與540 dagpm高度線相交于35°S左右，脊線與之相交在50°S附近。COSMIC資料形勢場圖的溫度場與高度場配合較好，與GFS資料500 hPa溫度場基本一致。

COSMIC資料可以較好地顯示西南太平洋區(qū)域的天氣形勢，但系統(tǒng)強(qiáng)度還有一定的誤差。掩星觀測點(diǎn)在高緯度地區(qū)多于低緯度地區(qū)，赤道附近部分區(qū)域由于原始數(shù)據(jù)過于稀疏，導(dǎo)致插值后仍然存在誤差較大甚至是數(shù)據(jù)空白。

4小結(jié)與討論

通過對COSMIC系統(tǒng)下發(fā)的GPS掩星反演大氣主要?dú)庀笠刭Y料進(jìn)行數(shù)據(jù)處理和可視化，使其轉(zhuǎn)化為能供預(yù)報(bào)員直接使用的天氣圖，并從環(huán)流形勢和區(qū)域天氣形勢2個(gè)方面探討COSMIC資料直接應(yīng)用于氣象預(yù)報(bào)保障效果。結(jié)論如下：

（1）COSMIC資料能較好地反映大氣環(huán)流形勢，北半球天氣尺度以上槽脊系統(tǒng)分布和形態(tài)均與再分析資料基本一致。在缺乏實(shí)況資料的情況下，COSMIC資料為中長期預(yù)報(bào)的制作提供了依據(jù)，同時(shí)也給短期天氣預(yù)報(bào)的制作提供環(huán)流背景。

（2）選取實(shí)況資料稀少的南太平洋作為參考區(qū)域，對比再分析資料發(fā)現(xiàn)，COSMIC資料天氣形勢場能準(zhǔn)確定位中尺度到天氣尺度槽脊系統(tǒng)的位置，同時(shí)溫度場分布也與再分析資料基本一致。

（3）雖然COSMIC資料在大尺度環(huán)流形勢和資料稀疏的大洋區(qū)域與再分析資料基本一致，但由于其空間分辨率偏低，客觀上對較小尺度系統(tǒng)起到了濾波作用，部分區(qū)域由于觀測點(diǎn)偏少，造成插值后誤差偏大。

COSMIC掩星觀測資料具有覆蓋面廣、精度高、全天候、穩(wěn)定性好等優(yōu)點(diǎn)，但其目前2 500個(gè)左右掩星點(diǎn)的密度分布還不夠，未來隨著低軌衛(wèi)星觀測星座的不斷加入，可以預(yù)見在不久的將來，掩星點(diǎn)的分布密度將大大增加，其應(yīng)用前景也更加廣闊。

參考文獻(xiàn)

[1]

丁金才，郭英華，郭永潤，等.利用 COSMIC 資料對 17 個(gè)臺風(fēng)熱力結(jié)構(gòu)的合成分析[J].熱帶氣象學(xué)報(bào)，2011，27（1）：31-43.

[2] 丁金才.GPS氣象學(xué)及其應(yīng)用[M].北京：氣象出版社，2009：198-211.

[3] ROCKEN C，KUO Y H，SCHREINER W，et al.COSMIC system Description[J].Special issue of terrestrial，atmospheric and oceanic science，2000，11（1）：21-52.

[4] HAJJ G A，LEE L C，PI X，et al.COSMIC GPS ionospheric sensing and space weather[J].Terrestrial，atmospheric and oceanic sciences，2000，11（1）：235-272.

[5] KURSINSKI E R，HAJJ G A，LEROY S S，et al.The GPS radio occultation technique[J].Terrestrial，atmospheric and oceanic sciences，2000，11（1）：53-114.

[6] POLI P，HEALY S B，DEE D P.Assimilation of global positioning system radio occultation data in the ECMWF ERA：Interim reanalysis[J].Quarterly journal of the royal meteorological society，2010，136：1972-1990.

[7] GUO P，KUO Y H，SOKOLOVSKIY S，et al.Estimating atmospheric boundary layer depth using COSMIC radio occultation data[J].Journal of the atmospheric sciences，2011，68（8）：1703-1713.

[8] 王也英，杜曉勇，袁勇.我國天基GPS掩星探測現(xiàn)狀及發(fā)展趨勢[J].氣象科技，2011，39（2）：202-206.

[9] 陳信雄.航天測量船氣象保障技術(shù)[M].北京：國防工業(yè)出版社，2009：86-99.