喬曉燕, 馬舒慶, 陶士偉, 趙培濤, 官福順

(中國氣象科學(xué)研究院,北京100081)

1 引言

目前,商用飛機(jī)氣象觀測用AMDAR(Aircraft Meteorological Data Relay)來代指。AMDAR即航空器氣象數(shù)據(jù)中繼是飛機(jī)對地通訊系統(tǒng)的泛指,用來代表對現(xiàn)代商用飛機(jī)上進(jìn)行的自動氣象觀測資料的收集系統(tǒng)。AMDAR是國際民航組織和世界氣象組織積極推動的全球性的項(xiàng)目。我國于2001年底在中國氣象局和中國民航總局的協(xié)商下啟動了AMDAR項(xiàng)目,此后AMDAR報告的數(shù)量逐漸增加。并且開始逐步應(yīng)用于氣象預(yù)報以及機(jī)場臨近預(yù)報業(yè)務(wù)。

AMDAR資料的質(zhì)量好壞關(guān)系著其在天氣預(yù)報業(yè)務(wù)中的應(yīng)用效果,就這一點(diǎn)而言,研究AMDAR資料的質(zhì)量是一項(xiàng)非常有意義的工作。目前,美國NOAA在這方面開展了大量的工作。例如1995年,通過與機(jī)場附近的探空資料對比分析了飛機(jī)氣象資料的誤差特性。1999年,選取滿足一定時間間隔(不超過10min)和空間間隔(水平距離不超過10km、垂直距離不超過30m)的飛機(jī)觀測要素數(shù)據(jù)對(數(shù)據(jù)對取自不同飛機(jī))進(jìn)行了對比研究,利用統(tǒng)計學(xué)方法計算出了飛機(jī)觀測隨機(jī)誤差均方根。國內(nèi)在這方面開展的工作較少,2006年國家氣象中心利用統(tǒng)計學(xué)方法對北京54511站2005年7月早00：00時探空資料與同時間北京首都國際機(jī)場飛機(jī)起飛時的AMDAR數(shù)據(jù)分5個層次進(jìn)行對比,各要素取月平均值。結(jié)果表明AMDAR資料和探空資料的差別很小。但是沒有計算AMDAR資料和探空資料的差值,以及均方根誤差。

AMDAR資料與探空資料對比時,兩者要素差值的均方根包括來自AMDAR觀測的誤差和來自探空觀測的誤差,但無法從標(biāo)準(zhǔn)差中分離出AMDAR觀測誤差。借鑒1999年美國NOAA所用的分析方法,對北京首都國際機(jī)場2007年1-8月份飛機(jī)起降時的氣象數(shù)據(jù)進(jìn)行了研究,計算了不同飛機(jī)的AMDAR數(shù)據(jù)之間的差值以及該差值的均方根,并利用統(tǒng)計學(xué)方法估計了飛機(jī)觀測的誤差。

2 AMDAR氣象觀測

目前國內(nèi)商業(yè)飛機(jī)上氣象觀測要素包括時間、高度(有氣壓換算得到)、經(jīng)緯度、飛機(jī)飛行狀態(tài)、氣溫、風(fēng)向、風(fēng)速。飛機(jī)在各個階段的觀測頻次見表1。

飛機(jī)在降落階段垂直分辨率大約在600m左右,在飛機(jī)起飛階段垂直分辨率要好一些,在起飛階段0-1400m數(shù)據(jù)空間密度最大。

表1 飛機(jī)各階段觀測頻率

AMDAR氣象要素觀測原理為：

(1)氣溫觀測：溫度傳感器探頭觀測值是大氣總溫,空氣靜溫(自由氣流的溫度)與大氣總溫不同,由馬赫數(shù)對大氣總溫修正得到空氣靜溫,其表達(dá)式如下：

其中T0是大氣靜溫、T1是大氣總溫、M 是馬赫數(shù)。馬赫數(shù)由空氣總壓和大氣靜壓的計算得到,總壓和靜壓均有皮托管測得。由此可知,氣溫觀測誤差源包括來自總溫觀測、總壓觀測、靜壓觀測以及計算和訂正誤差。

(2)風(fēng)速風(fēng)向觀測：風(fēng)矢量的測量是非常復(fù)雜的。由飛機(jī)相對于地面的速度和飛機(jī)相對空氣速度相減獲得。利用從飛機(jī)導(dǎo)航系統(tǒng)(一般是慣性導(dǎo)航系統(tǒng))和空速系統(tǒng)(一般是皮托管)獲得的數(shù)據(jù),加上從溫度傳感器獲得的數(shù)據(jù),可以計算出具有很高準(zhǔn)確度的飛機(jī)相對于地面的速度 Vg和空氣相對于飛機(jī)的速度Va,從而可計算出風(fēng)速 V,表達(dá)式如下：

對于大多數(shù)應(yīng)用,只測量風(fēng)的水平分量。這時要求輸入的數(shù)據(jù)縮減為只需空速、航向和地速。航向和地速取自導(dǎo)航系統(tǒng),真空速需根據(jù)空速指示器的校正空速計算出來。水平風(fēng)的分量(u,v)為：

其中,|va|是真空速的絕對值,φ為相對正北的航向,ug和vg是地速的分量。

3 分析方法及資料獲取

3.1 誤差分析方法

首先利用統(tǒng)計學(xué)知識計算出數(shù)據(jù)對間要素差值的均方根 σ。假設(shè)：不同飛機(jī)的系統(tǒng)誤差是一樣的或沒有系統(tǒng)誤差;不同飛機(jī)的隨機(jī)誤差是相互獨(dú)立的。則此均方根滿足下面的公式

其中σ1、σ2分別代表兩飛機(jī)的觀測隨機(jī)誤差,σ0代表兩飛機(jī)所處位置的中尺度過程引起的均方根誤差。如果兩飛機(jī)的時空間隔足夠小,中尺度過程引起的均方根誤差可以忽略。在大樣本統(tǒng)計的情況下,可以認(rèn)為兩飛機(jī)的觀測誤差相同,則飛機(jī)隨機(jī)誤差為

3.2 AMDAR資料獲取

研究數(shù)據(jù)是北京首都國際機(jī)場2007年1-8月飛機(jī)在航站起飛和降落的AMDAR數(shù)據(jù)。所選數(shù)據(jù)通過質(zhì)量控制去除了含有粗大誤差的數(shù)據(jù)點(diǎn),所用數(shù)據(jù)僅包括飛機(jī)系統(tǒng)誤差和觀測隨機(jī)誤差。表2給出了溫度觀測隨機(jī)誤差隨數(shù)據(jù)對時空間隔的變化情況。由表2可看出隨著數(shù)據(jù)對時間間隔和空間間隔的變小,溫度觀測隨機(jī)誤差明顯變小,也就是說隨著時空間隔的減小,由于中尺度變化造成的隨機(jī)誤差減小。由此可以用時空間隔足夠小的數(shù)據(jù)對來估計飛機(jī)氣象觀測隨機(jī)誤差。

表2 AMDAR溫度資料觀測誤差均方根隨樣本數(shù)據(jù)對(來自不同飛機(jī))時間間隔和水平間隔的變化

因此,選取數(shù)據(jù)對滿足的條件設(shè)定為：時間間隔小于10min、水平間隔小于10km、垂直高度差不超過30m。該條件足以反映兩飛機(jī)觀測的隨機(jī)誤差,由于中尺度變化造成的誤差已經(jīng)很小。

4 結(jié)果分析

4.1 溫度誤差特性分析

表3給出溫度觀測隨機(jī)誤差在各層的情況?？梢钥闯鲭S著氣壓的降低(高度的增加)飛機(jī)觀測隨機(jī)誤差總體趨勢是下降的。在大氣邊界層溫度觀測誤差均方根最大,可達(dá)0.91℃。在300hPa-500hPa層溫度觀測誤差均方根最小,主要是由于對流層上層湍流運(yùn)動減弱,由于湍流運(yùn)動引起的氣象要素的變化減小。

表3 各氣壓層的溫度觀測隨機(jī)誤差

數(shù)據(jù)對溫度差值的分布也與溫度本身有關(guān)(如圖1)。由圖1可見,在溫度較高時溫度差的分布比較離散,大差值較多。較低溫度時大差值較少且越來越集中在較小值。溫度小于-20℃時,其差不超過 2℃。溫度大于 20℃時,其差可達(dá)6℃。

4.2 風(fēng)矢量誤差特性分析

風(fēng)矢量誤差特性分析中與溫度相同,風(fēng)速觀測隨機(jī)誤差也隨著氣壓的降低而減小。如表4所示,風(fēng)速差觀測隨機(jī)誤差在300hPa-500hPa氣壓層為最小值(1.03m/s)。因?yàn)樵诟邔哟髿鈱α鲗?湍流運(yùn)動引起的風(fēng)速空間變化減小。在300hPa-500hPa氣壓層可以認(rèn)為中尺度變化引起的誤差達(dá)到最小值。

圖1 溫度差與平均溫度關(guān)系

表4 各氣壓層風(fēng)速觀測隨機(jī)誤差

風(fēng)速差的分布也與風(fēng)速的大小有關(guān)。圖2反映了風(fēng)速差隨風(fēng)速的變化情況。所用的樣本數(shù)為1651個,由于風(fēng)速、風(fēng)速差都是整數(shù),很多點(diǎn)在圖中重合在一起?？傮w趨勢是：風(fēng)速較大時風(fēng)速差都較小,風(fēng)速較小時會出現(xiàn)較大的風(fēng)速差值。特別是風(fēng)速大于30m/s時,風(fēng)速差值不超過4m/s。

風(fēng)向是比較難測量的要素,影響因素很多,所以風(fēng)向差一般比較大,但是風(fēng)向的測量精度與風(fēng)速的大小也有關(guān)。圖3給出了樣本風(fēng)向差隨風(fēng)速的變化,風(fēng)速較小時風(fēng)向差較大,當(dāng)風(fēng)速達(dá)到一定值時風(fēng)向差較小。風(fēng)速大于10m/s時,不同飛機(jī)測得的風(fēng)向差值很小。

圖2 風(fēng)速差隨風(fēng)速的變化

圖3 風(fēng)向差隨風(fēng)速的變化

5 與美國NOAA統(tǒng)計結(jié)果的對比

1999年美國NOAA對美國中部地區(qū)1996年8月-1997年8月的AMDAR數(shù)據(jù)進(jìn)行了統(tǒng)計。所選數(shù)據(jù)對與本文相同：時間間隔小于10min、水平間隔小于10km、垂直高度差不超過30m。其結(jié)果如表5所示。對比表3和表5可知：我們的溫度觀測誤差均方根略大于美國NOAA的統(tǒng)計結(jié)果,但差別不大,兩溫度觀測隨機(jī)誤差隨氣壓的變化趨勢基本一致。對比表4和表5可以看出,我國的商用飛機(jī)風(fēng)速觀測隨機(jī)誤差小于美國NOAA統(tǒng)計的風(fēng)速觀測隨機(jī)誤差,兩者隨氣壓的變化趨勢一致。由此可見：我國2007年1-8月份的飛機(jī)觀測數(shù)據(jù)和美國1996年8月-1997年8月的飛機(jī)觀測數(shù)據(jù)相比,溫度觀測資料隨機(jī)誤差相近,我國飛機(jī)風(fēng)速觀測隨機(jī)誤差要小于美國的。

表5 美國NOAA統(tǒng)計結(jié)果

6 結(jié)論

通過選取滿足一定時空間隔不同飛機(jī)觀測數(shù)據(jù)并對每對觀測風(fēng)速差、風(fēng)向差、溫度差進(jìn)行統(tǒng)計分析,揭示了商用飛機(jī)氣象觀測誤差的一些特性,結(jié)論如下：

(1)溫度觀測隨機(jī)誤差、風(fēng)速觀測隨機(jī)誤差均隨著氣壓的減小(即高度的增加)而變小。在對流層底部(＞850hPa)時觀測隨機(jī)誤差最大,溫度觀測隨機(jī)誤差均方根為0.90℃,風(fēng)速觀測隨機(jī)誤差均方根為1.24m/s。在300hPa-500hPa層溫度觀測誤差均方根最小。與1999年美國NOAA統(tǒng)計結(jié)果相比,風(fēng)速觀測隨機(jī)誤差均方根在各氣壓層均比美國NOAA的統(tǒng)計結(jié)果要小一些,而溫度觀測隨機(jī)誤差均方根與美國NOAA的統(tǒng)計結(jié)果差別不大。

(2)溫度觀測誤差與溫度本身有關(guān),溫度較低時,溫度差較小,精度較高。同樣風(fēng)速觀測誤差、風(fēng)向觀測誤差也與風(fēng)速有關(guān),風(fēng)速較大時,風(fēng)向差、風(fēng)速差均較小。風(fēng)向差隨風(fēng)速的變化較明顯,風(fēng)速超過10m/s時,風(fēng)向差都很小。

致謝：感謝中國氣象局氣象探測中心高空室給予的幫助!

[1]Stanley G Benjamin,Barry E Schwartz.Accuracy of ACARS Wind and Temperature Observations Determined by Collocation[J].Weather and Forecasting,1999,14：1032-1038.

[2]Barry Schwartz,Stanley G Benjamin.A Comparsion of Temperature and Wind Measurements from ACARS-E-quipped Aircraft and Rawinsondes[J].Weather and Forecasting,1995,10：528-543.

[3]Geneva Switzerland.Sensors and Measurements[M].AMDAR Reference Manual,2003.

[4]Fleming,R J.The use of commercial aircraft as platforms for environmental measurements[J].American meteorological society,1996,77：2229-2242.

[5]劉小魏,曹之玉,蘭海波.AMDAR資料特征及質(zhì)量分析[J].氣象科技,2007,35(4)：480-483.

[6]黃卓,李延香,王慧.AMDAR資料在天氣預(yù)報中的應(yīng)用[J].氣象,2006,32(9)：42-48.

[7]拓瑞芳,金山,丁葉鳳.AMDAR資料在機(jī)場天氣預(yù)報中的應(yīng)用[J].氣象,2006,32(3)：44-48.

[8]賈朋群,胡英,王金星.民用航空氣象觀測綜述[J].氣象科技,2004,32(4)：213-218.